Khảo sát và thiết kế hệ truyền động cho lò quay

MỤC LỤC DANH MỤC HÌNH VẼ Hình 1.1 Sơ đồ dây chuyền công nghệ sản xuất xi măng Hình 2.2 Xu hướng kinh tế hệ truyền động 14 Hình 2.2 Sơ đồ thay động không đồng .16 Hình 2.3 Đặc tính động không đồng 18 Hình 2.4 Chuyển đổi tương đương cuộn dây động không đồng 21 Hình 2.5 Mặt cát ngang động KĐB với hệ tọa độ gắn Stator Rotor 22 Hình 2.6 Dòng điện stator biểu diễn hai trục tọa độ gắn với stator rotor 28 Hình 2.7 Biên độ vector không gian góc quy trê hệ tọa độ khác 29 Hình 2.8 Sơ đồ nguyên lý đặc tính thay đổi điện trở phụ mạch rotor .35 Hình 2.9 Họ đặc tính điều chỉnh điện áp 38 Hình 2.10 Sơ đồ cấu trúc hệ điều khiển vectơ 42 Hình 2.11 Quỹ đạo từ thông stator tác dụng vector điện áp đóng ngắt 44 Hình 2.12 Sơ đồ khối DTC .46 Hình 3.1 Cấu trúc biến tần .50 Hình 3.2 Mạch điện thay pha động không đồng 53 Hình 3.3 Sơ đồ nghịch lưu diode kẹp điện áp mức .54 Hình 3.4 Sơ đồ mạch lực biến tần 63 Hình 3.5 Các đâu vào biến tân .64 Hình 3.7 Panel vận hành biến tần .65 Hình 3.8 Chế độ điều khiển chỗ/từ xa 66 Hình 4.1 Vòng điều chỉnh tốc độ 66 Hình 4.1 Mô hình mô ước lượng từ thông stator mô men điện từ 68 Hinh 4.3 Mô hình khối điều khiển DTC 70 Hình 4.4 Sơ đồ matlab khối bảng chọn vector 71 Hình 4.5 Sơ đồ matlab khối khối tính mô men 72 Hình 4.6 Sơ đồ matlab khối tính thành phần từ thông stator hệ tọa độ dq .72 Hình 4.7 Sơ đồ matlab khối tính điện áp stator hệ tọa độ cố định 73 Hình 4.8 Sơ đồ matlab khối tính dòng điện stator hệ tọa độ cố định .73 Hình 4.9 Quỹ đạo từ thông tải Mc=Mdm=8108 Nm 73 Hình 4.10 Đồ thị mô men tải Mc= Mdm =8108 Nm 74 Hình 4.11 Quỹ đạo từ thông có tải Mc= 1000 Nm 74 Hình 4.12 Đồ thị mô men tải Mc=1000 N 75 Lời nói đầu LỜI NÓI ĐẦU Hiện động xoay chiều ngày ứng dụng phổ biến, đặc biệt động xoay chiều rotor lồng sóc Do phương pháp điều khiển động xoay chiều cải thiện Phương phương điều khiển trực tiếp mô men (DTC) phương pháp điều khiển vị trí vector từ thông để điều khiển mô men động Sự kết hợp dùng biến tần điều khiển động rotor lồng sóc mang lại hiệu cao mặt kinh tế kỹ thuật Chính việc thiết kế biến tần ứng dụng phương pháp điều khiển động cho động xoay chiều có ý nghĩa thiết thực cao Trong trình thực tập nhà máy xi măng Vicem Hoàng Thạch em tìm hiểu nghiên cứu công nghệ sản xuất cách vận hành nhà máy Dây chuyền dây chuyền sử dụng động xoay chiều không đồng rotor lồng sóc đem lại hiệu suất cao trình sản xuất nhà máy Vì trước trường em làm đồ án tốt đề tài: “Khảo sát thiết kế hệ truyền động cho lò quay” Đồ án em gồm chương : - Chương I Tổng quan công nghệ xi măng sản xuất xi măng lò quay Chương II Phân tích chọn lựa phương án truyền động cho lò quay Chương III Giới thiệu biến tần tính chọn van mạch lực biến tần Chương IV Tổng hợp mô hệ thống Em xin chân thành cảm ơn thầy Võ Minh Chính thầy môn Tự động hóa xí nghiệp công nghiệp – Viện Điện – Đại học Bách Khoa Hà Nội kỹ sư nhà máy xi măng Hoàng Thạch ân cần bảo hướng dẫn em hoàn thành đồ án tốt nghiệp Hà Nội, ngày tháng 06 năm 2014 Sinh viên thực Nguyễn Xuân Hải Chương Tổng quan công nghệ sản xuất xi măng lò quay Chương TỔNG QUAN CÔNG NGHỆ SẢN XUẤT XI MĂNG VÀ LÒ QUAY 1.1 Giới thiệu chung xi măng Xi măng loại vật liệu xây dựng, chất kết dính xây dựng mà nhà khoa học tìm vào cuối kỉ thứ 19 sản xuất trước tiên số nước tư Đan Mạch, Anh, Pháp, Mỹ… Đầu kỷ 20, xi măng thực trở thành nhu cầu thiếu xây dựng phát triển kinh tế Xi măng xuất hầu hết khắp thị trường giới Xi măng dùng chủ yếu ngày xi măng Poocland Vậy xi măng Poocland gì? Xi măng Poocland chất kết dính chịu nước, sản xuất từ nguyên liệu đá vôi đá sét nung đến nhiệt độ 14600C sau làm nguội nghiền mịn Ngoài nguyên liệu đá vôi đá sét dùng thêm phụ gia hoạt tính xỉ pyrit, cát thạch anh pha trộn theo tỷ lệ định Thành phần hoá học có oxít quan trọng : S - SiO2=21- 25 % Al- Al2O3 =4-7 % F - Fe2O3 = 2,5- % C- CaO = 63- 67 % 1.2 Sơ đồ công nghệ Chương Tổng quan công nghệ sản xuất xi măng lò quay Hình 1.1 Sơ đồ công nghệ sản xuất xi măng Chương Tổng quan công nghệ sản xuất xi măng lò quay 1.3 Các công đoạn sản xuất xi măng 1.3.1 Chuẩn bị nguyên liệu Đá vôi: Đá vôi khai thác phương pháp khoan nổ, cắt tầng theo quy trình quy hoạch khai thác, sau đá vôi xúc vận chuyển tới máy đập búa thiết bị vận chuyển có trọng tải lớn, đá vôi đập nhỏ thành đá dăm cỡ 25 x 25 vận chuyển băng tải kho đồng sơ rải thành đống riêng biệt, đống khoảng 15.000 Đá sét: Đá sét khai thác phương pháp cày ủi khoan nổ mìn bốc xúc vận chuyển thiết bị vận tải có trọng tải lớn máy đập búa Đá sét đập máy đập búa xuống kích thước 75 mm (đập lần 1) đập máy cán trục xuống kích thước 25 mm (đập lần 2) Sau đập đá sét vận chuyển rải thành đống riêng biệt kho đồng sơ bộ, đống khoảng 6.600 Phụ gia điều chỉnh:Để đảm bảo chất lượng Clinker, Công ty kiểm soát trình gia công chế biến hỗn hợp phối liệu theo Modun, hệ số xác định Do đá vôi đá sét có nguyên liệu điều chỉnh quặng sắt (giàu hàm lượng ô xít Fe2O3), quặng bôxit (giàu hàm lượng ô xít Al2O3) đá Silíc ( giàu hàm lượng SiO2) 1.3.2 Nghiền nguyên liệu Đá vôi, đá sét phụ gia điều chỉnh cấp vào máy nghiền qua hệ thống cân DOSIMAT cân băng điện tử Máy nghiền nguyên liệu sử dụng hệ thống nghiền bi sấy nghiền liên hợp có phân ly trung gian Các điều khiển tự động khống chế tỷ lệ % đá vôi, đá sét, bô xít quặng sắt cấp vào nghiền điều khiển máy tính điện tử thông qua số liệu phân tích hệ thống QCX, đảm bảo khống chế hệ số chế tạo theo yêu cầu Bột liệu sau máy nghiền vận chuyển đến xilô đồng nhất, hệ thống gầu nâng, máng khí động 1.3.3 Lò nung Dây chuyền I xi măng Hoàng Thạch dây chuyền sản xuất xi măng lò quay, phương pháp khô, chu trình kín, có hệ thống trao đổi nhiệt tầng(Cyclon) hệ thống làm nguội kiểu hành tinh gồm 10 lò Nhiên liệu hỗn hợp gồm 85% than cám 15% dầu MFO, công ty cải tạo lại vòi phun đốt 100% than cám 3, dầu nặng MFO Chương Tổng quan công nghệ sản xuất xi măng lò quay dùng cho sấy lò sử dụng nghiền than gặp cố thiếu than mịn Dây chuyền I xi măng Hoàng Thạch từ khâu nguyên liệu đến nghiền, đóng bao xuất xi măng tự động hoàn toàn Dây chuyền II III xi măng Hoàng Thạch dây chuyền sản xuất xi măng lò quay, phương pháp khô, chu trình kín, có hệ thống tiền nung(Canciner) tiêu hao nhiệt lượng thấ p 715 kcal/kg clinker, làm nguội kiểu ghi, tăng hiệu làm mát, chất lượng sản phẩm tốt, dễ nghiền Hệ thống điều khiển tự động hoàn toàn, đại công nghệ PJC Master Piece ABB 1.3.4 Nghiền xi măng Clinker từ xilô, thạch cao phụ gia từ kho chứa tổng hợp vận chuyển lên két máy nghiền hệ thống băng tải gầu nâng, từ két máy nghiền Clinker, thạch cao, phụ gia cấp vào máy nghiền định lượng hệ thống cân DOSIMAS Máy nghiền xi măng làm việc theo chu trình kín (có phân ly trung gian) Xi măng khỏi máy nghiền độ mịn đạt 3.200 cm2/g, vận chuyển tới xilô chứa xi măng bột hệ thống băng tải, máng khí động 1.3.5 Đóng bao xi măng Từ đáy silo chứa, qua hệ thống cửa tháo liệu xi măng vận chuyển tới két chứa máy đóng bao, phận xuất xi măng rời đường bộ, bao xi măng sau đóng xong qua hệ thống băng tải vận chuyển đến máng xuất đường bộ, đường sắt đường thuỷ 1.4 Tổng quan công nghệ lò quay Công nghệ sản xuất ximăng đòi hỏi hệ thống lò làm việc liên tục, qua lò bán thành phẩm clinker lò đóng vai trò quan trọng dây chuyền sản xuất xi măng Hiện thực tế có nhiều loại lò sử dụng công nghiệp: Lò đứng, lò phòng, lò quay, lò bể, lò hầm, lò tầng sôi, lò điện, lò đúc Mỗi loại lò có ưu điểm hạn chế khác nhau, tùy theo lĩnh vực yêu cầu xản suất mà người ta sử dụng loại lò định Trong công nghiệp sản xuất xi măng thường sử dụng loại lò đứng lò quay Chương Tổng quan công nghệ sản xuất xi măng lò quay Lò đứng có ưu điểm cấu tạo đơn giản, khí lò thải có nhiệt độ thấp tỉ lệ tiêu thụ nhiên liệu thấp, hệ số sử dụng nhiệt cao, làm việc liên tục Trong lò đứng, vật liệu chuyển động từ xuống sản phẩm cháy từ lên ngược chiều Trước lò đứng làm việc gián đoạn, nghĩa vật liệu đổ đầy lò đứng yên Ngọn lửa sản phẩm chảy từ lên đạt yêu cầu tắt lửa tháo sản phẩm Như lò làm việc có suất hiệu suất thấp Hiện đại đa số lò đứng cố gắng tiến lên liên tục có trang bị thiết bị tháo nạp sản phẩm nhiên liệu khí hóa tự động hóa, đồng thời việc cung cấp không khí để đốt nhiên liệu dùng quạt hút thổi có áp suất cao Do công suất lò đứng đưa lên mức cao Hiện lò đứng thường dùng để nung vôi, nung ximăng, nung sản phẩm bán sản phẩm ngành vật liệu chịu lửa, nung gạch, sấy đất sét, cát,… 1.4.1 Giới thiệu lò quay Lò quay ống hình trụ bên vật liệu chịu lửa, vỏ thép Lò quay tròn xung quanh đường tim Vật liệu sản phẩm cháy chuyển động ngược chiều Do lò quay nên vật liệu sáo trộn tiếp xúc với sản phẩm cháy tạo điều kiện trao đổi nhiệt tốt, nâng cao hiệu làm việc Lò quay loại lò làm việc liên tục, thường sử dụng nhiều để nung ximăng, nung vôi, nung vật liệu chịu lửa, sấy đất sét, cát,… ngành khác sử dụng nhiều Đặc điểm lò quay lò đặt nghiêng bệ đỡ với lăn đỡ (số bệ đỡ phụ thuộc kích thước kết cấu lò) dẫn động hệ thống vành lò lắp cố định vỏ lò động lò Các vành băng đa lắp lỏng vỏ lò nhờ guốc lò, có nhiệm vụ đỡ thân lò giúp trình quay lò dễ dàng Để lò vị trí làm việc ổn định, hệ thống lăn đẩy thủy lực với bơm dầu thủy lực lắp đặt bệ đỡ gần đầu vào lò Quá trình cung cấp nhiên liệu điều chỉnh trình đốt lò quay thực nhờ vòi đốt (vòi đốt đa kênh) 1.4.2 Cấu tạo hoạt động lò quay a) Cấu tạo lò quay Lò quay hay ống lò bốn thành phần hệ thống lò bao gồm: cyclon tháp trao đổi nhiệt, buồng phân hủy, ống lò bột làm nguội Lò quay ống trụ dài Lò Chương Tổng quan công nghệ sản xuất xi măng lò quay có vỏ làm thép, có lớp gạch chịu lửa bên lớp lót Đi kèm với lò phận: • Chi đốt, đèn đốt, miệng lửa, buồng đốt để đốt cháy nhiên liệu • Buồng làm việc, phận chủ yếu để chứa vật liệu cần gia công nhiệt, tất yêu cầu thực • Các thiết bị nạp vật liệu tháo sản phẩm • Các thiết bị hút thổi khí • Các kênh, ống dẫn khí, ống khói • Các thiết bị lợi dụng nhiệt khí thải • Các thiết bị khử bụi • Các thiết bị động lực làm chuyển động lò phận lò • Các thiết bị giúp cho việc vận hành kiểm tra lò Lò thiết kế cho phép trình nung nóng nguyên liệu xảy nguyên liệu chuyển xuống đầu ra, đồng thời khí cháy kiểm soát chuyển lên lò vào ống dẫn đứng Để tránh cho trình nung khỏi xâm nhập khí giả làm tăng nhiệt lượng tiêu thụ lò bố trí lắp đặt thêm thiết bị kín đầu vào đầu lò b) Hoạt động lò quay Tất chất cháy đưa vào lò đốt nóng vòi đốt Quạt gió (quạt ID) đặt sau tháp trao đổi nhiệt tạo ra, điều chỉnh điều tiết lượng khí cần thiế t lò qua cyclon tháp trao đổi nhiệt Lưu lượng khí tháp trao đổi nhiệt tổng của: chất khí cháy nhiên liệu, khí thừa lò để đảm bảo nhiên liệu cháy hết, khí thoát từ nguyên liệu trình nung nóng canxi hoá, khí giả thâm nhập vào hệ thống qua hệ thống làm kín đầu vào lò rò rỉ tháp trao đổi nhiệt • Khí: Khí để đốt cháy lò bao gồm gió đưa vào qua vòi đốt, gió lò khí giả vào lò sức hút quạt ID tạo Trước vào lò nung gió trao đổi nhiệt với clinker nóng rời khỏi lò qua máy làm nguội kiểu hành tinh kiểu ghi Trong suốt trình nhiệt độ khí tăng lên tới khoảng 800-1000oC 10 Chương Giới thiệu biến tần tính chọn van mạch lực 3.5.2 Đầu vào biến tần Biến tần ACS 800 1060_7SR có + đầu vào tương tự : AI1 (0-10V) AI2 (0-20 mA AI3(0-20 mA) + đầu tương tự : AO1 (0-20 mA) AO2 (0-20 mA) + đầu vào số : DI1-DI6 (24VDC) Hình 3.5 Các đâu vào biến tân68 72 Chương Giới thiệu biến tần tính chọn van mạch lực 3.5.3 Panel vận hành + Màn hình hiển thị LCD dòng với 20 kí tự + Ngôn ngữ lựa chọn cài đặt ta lựa chon tham số 99.01 Tại người dung lựa chọn ngôn ngữ thích hợp để cài đặt,hiện tai biến tần ACS 800 1060_7SR có 14 ngôn ngữ tiếng Anh,Anh Mĩ,Đức, Italia, Tây Ban Nha,Bồ Đào Nha,Hà Lan, Pháp, Đan Mạch, Phần Lan, Thụy Điển, Séc, Nga, Polish + Có mode vận hành : Hiển thị tín hiệu thực (phím ACT) Tham số (phím PAR) Chức (Phím FUNC) Lựa chọn hệ TĐ (phím DRIVE) +Chức phím khác 1.Khởi động 2.Dừng 3.Đặt lượng chủ đạo 4.Chiều quay thuận 5.Chiều quay ngược 6.Xóa lỗi 7.ĐK chỗ/từ xa Hình 3.7 Panel vận hành biến tần 73 Chương Giới thiệu biến tần tính chọn van mạch lực 3.5.4 Chế độ điều khiển chỗ/từ xa Hình 3.8 Chế độ điều khiển chỗ/từ xa Biến tần ACS 800 1060_7SR điều khiển panen hiển thị chỗ, điều khiển thông qua panen ngoài, thông qua đầu cực profiledbus.Khi lựa chọn phương thức điều khiển ta cần sét bit trạng thái đánh dấu ô tròn Nếu để trống ta lựa chọn phương thức điều khiển thông qua đầu cực profiledbus, kí hiệu L điều khiển panen hiển thị chỗ kí tự R tức điều khiển thông qua đầu cực profiledbus 3.5.5 Lựa chon phương pháp điều khiển động Ta lựa chọ tham số 99.04 MOTOR CTRL MODE chọn phương pháp điều khiển 99.04 MOTOR CTRL MODE = 0: DTC Điều khiển trực tiếp mô men = 1: SCALAR Điều khiển vô hướng Phương pháp điều khiển Vô hướng thích hợp với hệ thống : 74 Chương Giới thiệu biến tần tính chọn van mạch lực • • • • • Truyền động nhiều động với loại động khác Khi dòng định mức động nhỏ 1/16 lần dòng định mức BBT Dùng để test tải động Phương pháp điều khiển trực tiếp mô men thích hợp với hệ thống : Khi làm việc cần mô men lớn 75 Chương Tổng hợp mô hệ thống Chương TỔNG HỢP VÀ MÔ PHỎNG HỆ THỐNG 4.1 Tổng hợp vòng điều chỉnh tốc độ Hình 4.1 Vòng điều chỉnh tốc độ Ta coi DTC khâu Deadbeat tác động nhanh, tổng hợp điều chỉnh tốc độ, ta coi khâu DTC khâu quán tính bậc với số thời gian Tm Mô DTC với tín hiệu đặt dạng bậc thang chế độ không tải, ta thấy mô men đáp ứng sau khoảng 1-2 ms tùy thuộc tốc độ đặt Mặc dù tổng hợp điều chỉnh ta phải tìm cách bù ảnh hưởng T m, nhiên Tm nhỏ gây điều chỉnh lớn với mô men, cho Tm =2ms để có đáp ứng vừa phải Do đó, ta có hàm truyền xấp xỉ khâu DTC: GDTC = 1 = + sT ω + 0.002 s 76 Chương Tổng hợp mô hệ thống Phản hồi tốc độ dùng máy phát tốc có hàm truyền coi khâu quán tính bậc có T số thời gian nhỏ ω = 0.5ms , hệ số khếch đại kω = Khâu phát tốc có hàm truyền: Gw = 1 = + sTω + 0.0005s Đối tượng vòng vòng điều chỉnh tốc độ không tải G0 = ⇒ G0 = 1/ j s ( + sTm ) ( + sTω ) 1/ j s 1 + ( Tm + Tω ) s + Tm * Tω * s    Tm * Tω = 0.002*0.0005 = 10 − ≈ Ta có , bỏ qua thành phần Tm * Tω *s2 hàm truyền đối tượng T = Tm + T ω =0.0025 ta đối tượng vòng điều chỉnh tốc độ thành khâu Đặt tích phân quán tính có hàm truyền xấp xỉ: G0 = Hàm truyền đạt kín có dạng Fk = k s (1 + sT1 ) ω + as = ω * + as + bs + cs3 Chỉnh định tham số điều khiển theo tiêu chuẩn tối ưu đối xứng ta có a = 4T , b = 8T , c = 8T 77 Chương Tổng hợp mô hệ thống   Rω = k p  + ÷  TI s  Ta có điều khiển J kp = 2T , TI = 4T Bộ lọc hàm truyền Floc = 1 + 4sT 4.2 Mô hệ thống Hệ thống mô với mô hình động không đồng matlab, với thông số tính mục 3.2.1 Hình 4.2 Mô hình mô ước lượng từ thông stator mô men điện từ 78 Chương Tổng hợp mô hệ thống 4.2.1 Khối điều khiển Hinh 4.3 Mô hình khối điều khiển DTC 4.2.2 Khối chọn Mục đích lựa chọn vector đóng cắt tối ưu phát xung mở van đóng ngắt mạch nghịch lưu 79 Chương Tổng hợp mô hệ thống Hình 4.4 Sơ đồ matlab khối bảng chọn vector 80 Chương Tổng hợp mô hệ thống 4.2.3 Khối tính mô men te = P(ψ sD isQ −ψ sQ isD ) Theo công thức (2.71) ta có Hình 4.5 Sơ đồ matlab khối khối tính mô men 4.2.4 Khối tính thành phần từ thông stator hệ tọa độ dq Hình 4.6 Sơ đồ matlab khối tính thành phần từ thông stator hệ tọa độ dq 4.2.5 Khối tính điện áp dòng điện stator hệ tọa độ cố định Theo công thức (2.43) (2.15) 81 Chương Tổng hợp mô hệ thống Hình 4.7 Sơ đồ matlab khối tính điện áp stator hệ tọa độ cố định Hình 4.8 matlab dòng hệ định Sơ đồ khối tính điện stator tọa độ cố 4.3 Kết mô Hình 4.9 Quỹ đạo từ thông tải Mc=Mdm=8108 Nm 4.10 Đồ thị mô men tải Mc=Mdm=8108 Nm 82 Chương Tổng hợp mô hệ thống Hình 4.11 Quỹ đạo từ thông có tải Mc= 1000 Nm Hình 4.12 Đồ thị mô men tải Mc=1000 Nm 83 Chương Tổng hợp mô hệ thống Nhận xét : Đặc tính từ thông ổn định giá trị 1.6 Wb suốt trình hoạt động dù có thay đổi tải Mc= 1000 Nm , Mc= Mdm =8108 Nm Đường dặc tính mô men tác động nhanh, sau thời gian 1,1s tốc độ ổn định mômen tải trở gẫn giá trị đặt, có dao động nhẹ, bám theo giá tri đặt Đường đặc tính mô men chưa thật ổn định mà dao động nhẹ quanh giá trị đặt phép ước lượng độ phân giải chưa cao 84 Chương Tổng hợp mô hệ thống KẾT LUẬN Đồ án thực đầy đủ nội dung thầy hướng dẫn giao Phần tính toán mô đáp ưng yêu cầu hệ thống Tuy nhiên độ xác chưa tuyệt đối độ xác phụ thuộc vào độ xác phép ước lượng tham số mô men, từ thông, vị trí Việc nâng cao độ phân giải khả tính toán vi xử lý làm tăng độ xác phép ước lượng, qua cải thiện tính xác phương pháp Do giới hạn thời gian nên em chưa thực phương pháp cải tiển để nâng cao chất lượng điều khiển 85 Chương Tổng hợp mô hệ thống TÀI LIỆU THAM KHẢO [1] Bùi Quốc Khánh , Trần Trọng Minh: Đảm bảo chất lượng dòng đầu vào Matric Converter điều kiện điện áp lưới cân bằng; Tạp chí khoa học công nghệ trường ĐHKT, số 57, (2006) [2] Nguyễn Phùng Quang , Andreas, Truyền động điện thông minh, NXB Khoa học Kỹ Thuật, 2004 [3] Nguyễn Phùng Quang, Điều khiển tự động truyền động điện xoay chiều ba pha , NXB Giáo Dục 1998 [4] Võ Minh Chính, Phạm Quốc Hải, Trần Trọng Minh,Điện tử công suất , NXB Khoa Học Kỹ Thuật [5] Bùi Quốc Khánh ,Nguyễn Văn Liễn , Cơ sở truyền động điện, NXB Khoa Học Kỹ Thuật 2007 [6] http://www.ximanghoangthach.com/home.aspx, truy cập cuối ngày 20/5/2014 86