BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC SƯ PHẠM KỸ THUẬT THÀNH PHỐ HỒ CHÍ MINH LUẬN VĂN THẠC SĨ PHAN THÁI TỒN TỐI ƯU HÓA VỊ TRÍ CỦA THIẾT BỊ FACTS ĐỂ NÂNG CAO KHẢ NĂNG LÀM VIỆC CỦA HỆ THỐNG ĐIỆN NGÀNH: KỸ THUẬT ĐIỆN - 60520202 S K C0 3 Tp Hồ Chí Minh, tháng 10/2014 BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC SƯ PHẠM KỸ THUẬT THÀNH PHỐ HỒ CHÍ MINH LUẬN VĂN THẠC SĨ PHAN THÁI TỒN TỐI ƯU HÓA VỊ TRÍ CỦA THIẾT BỊ FACTS ĐỂ NÂNG CAO KHẢ NĂNG LÀM VIỆC CỦA HỆ THỐNG ĐIỆN NGÀNH: KỸ THUẬT ĐIỆN - 60520202 Tp Hồ Chí Minh, tháng 10/2014 BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC SƯ PHẠM KỸ THUẬT THÀNH PHỐ HỒ CHÍ MINH LUẬN VĂN THẠC SĨ PHAN THÁI TỒN TỐI ƯU HÓA VỊ TRÍ CỦA THIẾT BỊ FACTS ĐỂ NÂNG CAO KHẢ NĂNG LÀM VIỆC CỦA HỆ THỐNG ĐIỆN NGÀNH: KỸ THUẬT ĐIỆN - 60520202 Hướng dẫn khoa học: PGS.TS NGUYỄN HOÀNG VIỆT Tp Hồ Chí Minh, tháng 10/2014 LÝ LỊCH KHOA HỌC I LÝ LỊCH SƠ LƯỢC: Họ & tên: Phan Thái Tồn Giới tính: Nam Ngày, tháng, năm sinh: 03/03/1988 Nơi sinh: Tiền Giang Quê quán: Tiền giang Dân tộc: Kinh Chỗ riêng địa liên lạc: Q7, TPHCM E-mail: phanthaitonntt@yahoo.com II QUÁ TRÌNH ĐÀO TẠO: Cao Đẳng: Hệ đào tạo: Chính Quy Thời gian đào tạo từ 10/2006 đến 10/2009 Nơi học (trường, thành phố): Q4, Thành Phố Hồ Chí Minh Ngành học: Công nghệ Kỹ Thuật Điện-Điện Tử Đại học: Hệ đào tạo: Chính Quy Thời gian đào tạo từ 10/2009 đến 10/2011 Nơi học (trường, thành phố): Quận Bình Thạnh, Thành Phố Hồ Chí Minh Ngành học: Công nghệ Kỹ Thuật Điện-Điện Tử III QUÁ TRÌNH CÔNG TÁC CHUYÊN MÔN KỂ TỪ KHI TỐT NGHIỆP ĐẠI HỌC: Thời gian 11/200903/2011 03/2011- Nay Nơi công tác Công việc đảm nhiệm Công ty cổ phần khoa học công Kỹ thuật viên nghệ PE Công ty TNHH Intel Products Kỹ thuật viên Việt Nam LỜI CAM ĐOAN Tôi cam đoan công trình nghiên cứu Các số liệu, kết nêu luận văn trung thực chưa công bố công trình khác Tp Hồ Chí Minh, ngày … tháng … năm 201… (Ký tên ghi rõ họ tên) LỜI CẢM TẠ Trước tiên em xin bày tỏ lòng biết ơn sâu sắc đến cố PGS.TS Nguyễn Hoàng Việt dù điều kiện sức khỏe có phần sa sút thầy cố gắng tận tình hướng dẫn em hoàn thành luận văn tốt nghiệp Xin chân thành cảm ơn đến tất quý thầy cô Trường Đại Học Sư Phạm Kỹ Thuật trường Đại Học Bách Khoa TPHCM truyền đạt cho chúng em kiến thức vô quý báo làm hành trang cho em vững vàng công việc Xin cảm ơn tập thể lớp 2012b giúp đỡ, chia sẻ với khó khăn hoc tập, tạo cho em niềm tin nỗ lực phấn đấu hoàn thành khóa hoc Chân Thành Cảm Ơn! TÓM TẮT Luận văn trình bày phương pháp tìm vị trí thiết bị FACTS tối ưu hệ thống điện sử dụng giải thuật di truyền GA (Genetic Alglorithms), để từ tiến hành viết giao diện nhỏ ngôn ngữ lập trình Matlab tính toán tổn thất công suất phản kháng chi phí đầu tư hệ thống điện trước sau lắp thiết bị FACTS lưới điện IEEE 30 bus Thiết bị FACTS cụ thể luận văn SVC TCSC Kết tính toán chương trình xác ứng dụng cho nhiều lưới điện khác Trong chương giới thiệu dạng ổn định hệ thống điện, nguyên tắt điều khiển hệ thống điện Chương giới thiệu tổng quan thiết bị FACTS điều khiển thiết bị FACTS Chương giới thiệu giải thuật di truyền phân bố trào lưu công suất truyền thống làm sở cho lập trình chương trình tính toán Chương trình bày kết chương trình tính toán thay đổi trọng số toán tối ưu Chương đánh giá kết tính toán đề suất hướng phát triển để xây dựng nên công trình nghiên cứu tìm vị trí tối ưu thiết bị FACTS hệ thống điện hoàn chỉnh ABSTRACT This thesis presents methods to find the optimal location FACTS devices in power systems using genetic algorithm GA (Genetic Alglorithms), then proceed to write a small interface on the programming language Matlab computational costs reactive power losses and investment costs of the electricity system before and after installation of FACTS devices IEEE 30-bus grid FACTS devices specifically in this thesis is the SVC and TCSC Calculation results of this program is very accurate and can be applied to many different grids In Chapter introduces the types of power system stability, the material off of the power system control Chapters and present an overview of device controllers FACTS FACTS devices Chapter introduces genetic algorithms and distribution power flow as the basis for traditional programming calculation program Chapter presents the results of the calculation program to change the weights of the optimization problem Chapter reviews the results of calculations and proposal development continued to build up a study to find the optimal location of FACTS devices in power system complete MỤC LỤC TRANG Trang tựa i Quyết định giao đề tài ii Lý lịch cá nhân iii Lời cam đoan iv Cảm tạ v Tóm tắt vi Mục lục vii Danh sách chữ viết tắt viii Danh sách hình ix Danh sách bảng x PHẦN A: GIỚI THIỆU LUẬN VĂN Đặt vấn đề Mục tiêu nhiệm vụ luận văn Phạm vi nghiên cứu Phương pháp giải toán Điểm luận văn Giá trị thực tiễn Bố cục luận văn PHẦN B: NỘI DUNG CHƯƠNG 1: ĐIỀU KHIỂN VÀ CÁC GIỚI HẠN ỔN ĐỊNH CỦA HỆ THỐNG ĐIỆN 1.1 Đặc điểm lưới điện truyền tải 1.2 Cơ sở kiến thức 1.3 Điều khiển phân bố công suất 1.4 Giới hạn ổn định hệ thống 1.4.1Giới hạn nhiệt 1.4.2 Giới hạn điện áp 1.4.3 Giới hạn ổn định 1.4.3.1Ổn định độ 10 1.4.3.2 Ổn định dao động bé 10 1.4.4 Ổn định điện áp 11 1.5 Nguyên lý bù đường dây truyền tải 11 1.5.1 Bù song song 12 1.5.2 Bù nối tiếp 13 1.6 Tầm quan trọng việc giảm tổn hao công suất phản kháng 13 CHƯƠNG 2: TỔNG QUÁT VỀ THIẾT BỊ FACTS 2.1 Khái niệm thiết bị FACTS 15 2.1.1 Đặt vấn đề 15 2.1.2 Khái niệm 16 2.2 Các điều khiển thiết bị FACTS 17 2.2.1 Ưu điểm thiết bị FACTS 17 2.2.2 Phân loại thiết bị FACTS 17 2.2.2.1 Bộ bù ngang tĩnh SVC ( Static Var Compensator) 18 2.2.2.2 Bộ bù ngang tĩnh STATCOM (Static Synchronous Compensator) 22 2.2.2.3 Bộ bù dọc tĩnh TCSC( Thyristor Controlled Series Capacitor) 23 2.2.2.4 Bộ tích hợp điều khiển luồng công suất UPFC(Unified Power Flow Controlled) 27 CHƯƠNG 3: MÔ HÌNH CÁC BỘ ĐIỀU KHIỂN THIẾT BỊ FACTS 3.1 Giới Thiệu 30 3.2 Các điều khiển dựa Thyristor truyền thống 31 3.2.1Cuộn kháng điều khiển Thyristor (TCR) 31 3.2.2 Bộ bù tĩnh SVC ( Static Var Compensator) 36 function edit_SVC_MIN_Callback(hObject, eventdata, handles) % hObject handle to edit_SVC_MIN (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) % Hints: get(hObject,'String') returns contents of edit_SVC_MIN as text % str2double(get(hObject,'String')) returns contents of edit_SVC_MIN as a double % - Executes during object creation, after setting all properties function edit_SVC_MIN_CreateFcn(hObject, eventdata, handles) % hObject handle to edit_SVC_MIN (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles empty - handles not created until after all CreateFcns called % Hint: edit controls usually have a white background on Windows % See ISPC and COMPUTER if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor')) set(hObject,'BackgroundColor','white'); end function edit_SVC_MAX_Callback(hObject, eventdata, handles) % hObject handle to edit_SVC_MAX (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) % Hints: get(hObject,'String') returns contents of edit_SVC_MAX as text % str2double(get(hObject,'String')) returns contents of edit_SVC_MAX as a double % - Executes during object creation, after setting all properties function edit_SVC_MAX_CreateFcn(hObject, eventdata, handles) % hObject handle to edit_SVC_MAX (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles empty - handles not created until after all CreateFcns called % Hint: edit controls usually have a white background on Windows % See ISPC and COMPUTER if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor')) set(hObject,'BackgroundColor','white'); end function edit_NUMBER_OF_TCSC_Callback(hObject, eventdata, handles) % hObject handle to edit_NUMBER_OF_TCSC (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) % Hints: get(hObject,'String') returns contents of edit_NUMBER_OF_TCSC as text % str2double(get(hObject,'String')) returns contents of edit_NUMBER_OF_TCSC as a double % - Executes during object creation, after setting all properties function edit_NUMBER_OF_TCSC_CreateFcn(hObject, eventdata, handles) % hObject handle to edit_NUMBER_OF_TCSC (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles empty - handles not created until after all CreateFcns called % Hint: edit controls usually have a white background on Windows % See ISPC and COMPUTER if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor')) set(hObject,'BackgroundColor','white'); end function edit_TCSC_MIN_Callback(hObject, eventdata, handles) % hObject handle to edit_TCSC_MIN (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) % Hints: get(hObject,'String') returns contents of edit_TCSC_MIN as text % str2double(get(hObject,'String')) returns contents of edit_TCSC_MIN as a double % - Executes during object creation, after setting all properties function edit_TCSC_MIN_CreateFcn(hObject, eventdata, handles) % hObject handle to edit_TCSC_MIN (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles empty - handles not created until after all CreateFcns called % Hint: edit controls usually have a white background on Windows % See ISPC and COMPUTER if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor')) set(hObject,'BackgroundColor','white'); end function edit_TCSC_MAX_Callback(hObject, eventdata, handles) % hObject handle to edit_TCSC_MAX (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) % Hints: get(hObject,'String') returns contents of edit_TCSC_MAX as text % str2double(get(hObject,'String')) returns contents of edit_TCSC_MAX as a double % - Executes during object creation, after setting all properties function edit_TCSC_MAX_CreateFcn(hObject, eventdata, handles) % hObject handle to edit_TCSC_MAX (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles empty - handles not created until after all CreateFcns called % Hint: edit controls usually have a white background on Windows % See ISPC and COMPUTER if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor')) set(hObject,'BackgroundColor','white'); end function edit_ANPHA_1_Callback(hObject, eventdata, handles) % hObject handle to edit_ANPHA_1 (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) % Hints: get(hObject,'String') returns contents of edit_ANPHA_1 as text % str2double(get(hObject,'String')) returns contents of edit_ANPHA_1 as a double % - Executes during object creation, after setting all properties function edit_ANPHA_1_CreateFcn(hObject, eventdata, handles) % hObject handle to edit_ANPHA_1 (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles empty - handles not created until after all CreateFcns called % Hint: edit controls usually have a white background on Windows % See ISPC and COMPUTER if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor')) set(hObject,'BackgroundColor','white'); end function edit_ANPHA_2_Callback(hObject, eventdata, handles) % hObject handle to edit_ANPHA_2 (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) % Hints: get(hObject,'String') returns contents of edit_ANPHA_2 as text % str2double(get(hObject,'String')) returns contents of edit_ANPHA_2 as a double % - Executes during object creation, after setting all properties function edit_ANPHA_2_CreateFcn(hObject, eventdata, handles) % hObject handle to edit_ANPHA_2 (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles empty - handles not created until after all CreateFcns called % Hint: edit controls usually have a white background on Windows % See ISPC and COMPUTER if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor')) set(hObject,'BackgroundColor','white'); end function edit_LOCATION_OF_SVC_Callback(hObject, eventdata, handles) % hObject handle to edit_LOCATION_OF_SVC (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) % Hints: get(hObject,'String') returns contents of edit_LOCATION_OF_SVC as text % str2double(get(hObject,'String')) returns contents of edit_LOCATION_OF_SVC as a double % - Executes during object creation, after setting all properties function edit_LOCATION_OF_SVC_CreateFcn(hObject, eventdata, handles) % hObject handle to edit_LOCATION_OF_SVC (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles empty - handles not created until after all CreateFcns called % Hint: edit controls usually have a white background on Windows % See ISPC and COMPUTER if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor')) set(hObject,'BackgroundColor','white'); end function edit_RATE_OF_SVC_Callback(hObject, eventdata, handles) % hObject handle to edit_RATE_OF_SVC (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) % Hints: get(hObject,'String') returns contents of edit_RATE_OF_SVC as text % str2double(get(hObject,'String')) returns contents of edit_RATE_OF_SVC as a double % - Executes during object creation, after setting all properties function edit_RATE_OF_SVC_CreateFcn(hObject, eventdata, handles) % hObject handle to edit_RATE_OF_SVC (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles empty - handles not created until after all CreateFcns called % Hint: edit controls usually have a white background on Windows % See ISPC and COMPUTER if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor')) set(hObject,'BackgroundColor','white'); end function edit_LOCATION_OF_TCSC_Callback(hObject, eventdata, handles) % hObject handle to edit_LOCATION_OF_TCSC (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) % Hints: get(hObject,'String') returns contents of edit_LOCATION_OF_TCSC as text % str2double(get(hObject,'String')) returns contents of edit_LOCATION_OF_TCSC as a double % - Executes during object creation, after setting all properties function edit_LOCATION_OF_TCSC_CreateFcn(hObject, eventdata, handles) % hObject handle to edit_LOCATION_OF_TCSC (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles empty - handles not created until after all CreateFcns called % Hint: edit controls usually have a white background on Windows % See ISPC and COMPUTER if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor')) set(hObject,'BackgroundColor','white'); end function edit_RATE_OF_TCSC_Callback(hObject, eventdata, handles) % hObject handle to edit_RATE_OF_TCSC (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) % Hints: get(hObject,'String') returns contents of edit_RATE_OF_TCSC as text % str2double(get(hObject,'String')) returns contents of edit_RATE_OF_TCSC as a double % - Executes during object creation, after setting all properties function edit_RATE_OF_TCSC_CreateFcn(hObject, eventdata, handles) % hObject handle to edit_RATE_OF_TCSC (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles empty - handles not created until after all CreateFcns called % Hint: edit controls usually have a white background on Windows % See ISPC and COMPUTER if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor')) set(hObject,'BackgroundColor','white'); end function edit_X_OF_TCSC_Callback(hObject, eventdata, handles) % hObject handle to edit_X_OF_TCSC (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) % Hints: get(hObject,'String') returns contents of edit_X_OF_TCSC as text % str2double(get(hObject,'String')) returns contents of edit_X_OF_TCSC as a double % - Executes during object creation, after setting all properties function edit_X_OF_TCSC_CreateFcn(hObject, eventdata, handles) % hObject handle to edit_X_OF_TCSC (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles empty - handles not created until after all CreateFcns called % Hint: edit controls usually have a white background on Windows % See ISPC and COMPUTER if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor')) set(hObject,'BackgroundColor','white'); end function edit_QLOSS_Callback(hObject, eventdata, handles) % hObject handle to edit_QLOSS (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) % Hints: get(hObject,'String') returns contents of edit_QLOSS as text % str2double(get(hObject,'String')) returns contents of edit_QLOSS as a double % - Executes during object creation, after setting all properties function edit_QLOSS_CreateFcn(hObject, eventdata, handles) % hObject handle to edit_QLOSS (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles empty - handles not created until after all CreateFcns called % Hint: edit controls usually have a white background on Windows % See ISPC and COMPUTER if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor')) set(hObject,'BackgroundColor','white'); end function edit_COST_Callback(hObject, eventdata, handles) % hObject handle to edit_COST (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) % Hints: get(hObject,'String') returns contents of edit_COST as text % str2double(get(hObject,'String')) returns contents of edit_COST as a double % - Executes during object creation, after setting all properties function edit_COST_CreateFcn(hObject, eventdata, handles) % hObject handle to edit_COST (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles empty - handles not created until after all CreateFcns called % Hint: edit controls usually have a white background on Windows % See ISPC and COMPUTER if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor')) set(hObject,'BackgroundColor','white'); end function edit_VMIN_Callback(hObject, eventdata, handles) % hObject handle to edit_VMIN (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) % Hints: get(hObject,'String') returns contents of edit_VMIN as text % str2double(get(hObject,'String')) returns contents of edit_VMIN as a double % - Executes during object creation, after setting all properties function edit_VMIN_CreateFcn(hObject, eventdata, handles) % hObject handle to edit_VMIN (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles empty - handles not created until after all CreateFcns called % Hint: edit controls usually have a white background on Windows % See ISPC and COMPUTER if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor')) set(hObject,'BackgroundColor','white'); end function edit_QLOSS_INITIAL_Callback(hObject, eventdata, handles) % hObject handle to edit_QLOSS_INITIAL (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) % Hints: get(hObject,'String') returns contents of edit_QLOSS_INITIAL as text % str2double(get(hObject,'String')) returns contents of edit_QLOSS_INITIAL as a double % - Executes during object creation, after setting all properties function edit_QLOSS_INITIAL_CreateFcn(hObject, eventdata, handles) % hObject handle to edit_QLOSS_INITIAL (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles empty - handles not created until after all CreateFcns called % Hint: edit controls usually have a white background on Windows % See ISPC and COMPUTER if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor')) set(hObject,'BackgroundColor','white'); end function edit_VMIN_INITIAL_Callback(hObject, eventdata, handles) % hObject handle to edit_VMIN_INITIAL (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) % Hints: get(hObject,'String') returns contents of edit_VMIN_INITIAL as text % str2double(get(hObject,'String')) returns contents of edit_VMIN_INITIAL as a double % - Executes during object creation, after setting all properties function edit_VMIN_INITIAL_CreateFcn(hObject, eventdata, handles) % hObject handle to edit_VMIN_INITIAL (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles empty - handles not created until after all CreateFcns called % Hint: edit controls usually have a white background on Windows % See ISPC and COMPUTER if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor')) set(hObject,'BackgroundColor','white'); end % - Executes on button press in pushbutton_START function pushbutton_START_Callback(hObject, eventdata, handles) % hObject handle to pushbutton_START (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) tcscsvc(handles); % - Executes on button press in pushbutton_RESET function pushbutton_RESET_Callback(hObject, eventdata, handles) % hObject handle to pushbutton_RESET (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) set(handles.edit_QLOSS,'string','MVAr'); set(handles.edit_COST,'string','$'); set(handles.edit_RATE_OF_TCSC,'string','MVAr'); set(handles.edit_LOCATION_OF_TCSC,'string','Branch'); set(handles.edit_X_OF_TCSC,'string','K'); set(handles.edit_LOCATION_OF_SVC,'string','Bus'); set(handles.edit_RATE_OF_SVC,'string','MVar'); set(handles.edit_QLOSS_INITIAL,'string','Initial Qloss'); set(handles.edit_VMIN_INITIAL,'string','Initial Vmin'); set(handles.edit_VMIN,'string','New Vmin'); set(handles.edit_OPTIMAL_INDEX,'string','Index'); function edit_OPTIMAL_INDEX_Callback(hObject, eventdata, handles) % hObject handle to edit_OPTIMAL_INDEX (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) % Hints: get(hObject,'String') returns contents of edit_OPTIMAL_INDEX as text % str2double(get(hObject,'String')) returns contents of edit_OPTIMAL_INDEX as a double % - Executes during object creation, after setting all properties function edit_OPTIMAL_INDEX_CreateFcn(hObject, eventdata, handles) % hObject handle to edit_OPTIMAL_INDEX (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles empty - handles not created until after all CreateFcns called % Hint: edit controls usually have a white background on Windows % See ISPC and COMPUTER if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor')) set(hObject,'BackgroundColor','white'); end function [Qloss,rateTCSC]=main_tcscsvc(XY,busdata2,linedata2,old) busdata=busdata2; linedata=linedata2; lfybus lfnewton lineflow Ploss=deltaP;Qloss=deltaQ; new=thien; nhanh=[linedata2(branch_k,1),linedata2(branch_k,2)]; [m,n]=size(new); for i=1:m tachnhanh=[new(i,1),new(i,2)]; if nhanh==tachnhanh rateTCSC=abs(old(i,5))- abs(new(i,5)); break; end end function tcscsvc(handles) clear all; clc %doc du lieu ban dau busdata0; linedata0; %goi chuong trinh tinh toan popsize=get(handles.edit_POPSIZE,'string'); popsize=str2num(popsize); maxit=get(handles.edit_ITERATION,'string'); maxit=str2num(maxit); mutrate=get(handles.edit_MUTATION_RATE,'string'); mutrate=str2num(mutrate); selection=get(handles.edit_SELECTION,'string'); selection=str2num(selection); anpha1=get(handles.edit_ANPHA_1,'string'); anpha1=str2num(anpha1); anpha2=get(handles.edit_ANPHA_2,'string'); anpha2=str2num(anpha2); npar=4; varlo=[1 -0.8 -100]; %[min position TCSC, rate TCSC, position_SVC, varhi=[41 0.2 30 100]; busdata=busdata0; linedata=linedata0; lfybus lfnewton lineflow busout_plotV; vmin_Initial = vmin; old_Voltage = thuan(:,2); Old_Qloss=deltaQ; old=thien; [bestcost,bestsolution,globalbestcost,rate_TCSC,new_Qloss,cost] = GA(npar,varhi, varlo,popsize,maxit,busdata0,linedata0,old,mutrate,selection,anpha1, anpha2); figure(123); plot(globalbestcost,'b'); grid on title('\fontsize{12}\bf Optimal location of FACTS devices'); xlabel('\fontsize{12}\bf Iteration');ylabel('\fontsize{12}\bf Optimal Index'); legend('\fontsize{10}\bf Genetic Algorithim',1); set(handles.edit_OPTIMAL_INDEX,'string',num2str(bestcost)); set(handles.edit_QLOSS,'string',num2str(new_Qloss)); set(handles.edit_COST,'string',num2str(cost)); set(handles.edit_RATE_OF_TCSC,'string',num2str(-rate_TCSC)); set(handles.edit_LOCATION_OF_TCSC,'string',num2str(bestsolution(1,1))); set(handles.edit_X_OF_TCSC,'string',num2str(bestsolution(1,2))); set(handles.edit_LOCATION_OF_SVC,'string',num2str(bestsolution(1,3))); set(handles.edit_RATE_OF_SVC,'string',num2str(bestsolution(1,4))); set(handles.edit_QLOSS_INITIAL,'string',num2str(Old_Qloss)); set(handles.edit_VMIN_INITIAL,'string',num2str(vmin_Initial)); [V_plot,Vmin]=take_new_vmin(bestsolution,busdata0,linedata0);%chuong trinh lay vmin va ve V sau uu set(handles.edit_VMIN,'string',num2str(Vmin)); figure(456); plot(V_plot,'-r'); title('\fontsize{12}\bf Compare voltages after optimization'); xlabel('\fontsize{12}\bf Bus');ylabel('\fontsize{12}\bf P.U.'); legend('\fontsize{10}\bf Optimal',1); hold on; plot(old_Voltage,'-b'); grid on; function [bestcost,bestsolution,globalbestcost,rate_TCSC,new_Qloss,money ]=GA(npar,varhi, varlo, popsize,maxit) % _ ff='testfunction'; % objective function Nt=npar; % keep=floor(selection*popsize); % # (lam tron so nguyen nho nhat) nmut = ceil((popsize-1)*Nt*mutrate); % (lam tron ve so nguyen lon hon) M=ceil((popsize-keep)/2); % so luong sinh san % _ iga=0; % Nvar = size(FieldD,2); aux = rand(popsize,npar); m=[-1 1]*FieldD; ublb=ones(popsize,1)*m; lb=ones(popsize,1)*FieldD(1,:); par=ublb.*aux+lb; [cost,TCSC]=feval(ff,par,busdata1,linedata1,old,anpha1, anpha2); % [cost,ind]=sort(cost); % cost in element sap xep theo thu tu tang dan, % _ while iga