Tiểu luận môn phương pháp nghiên cứu khoa học mô hình toán học cải tiến đánh giá lợi ích hợp tác quốc tế trong việc giảm

Trường Đai Học Bách Khoa Tp Hồ Chí Minh Khoa Khoa Học & Kỹ Thuật Máy Tính MƠN PHƯƠNG PHÁP NGHIÊN CỨU KHOA HỌC Báo cáo tập Mơ Hình Tốn Học Cải Tiến Đánh Giá Lợi Ích Hợp Tác Quốc Tế Trong Việc Giảm Thiểu Khí C02 CBGD: TS Nguyễn Văn Minh Mẫn Nhóm 5: Phạm Phi Long - 7140245 Lê Trương Trọng Duy - 7140227 Trịnh Thị Hồng Nhung - 7140672 Trường Đại Học Bách Khoa Tp.Hồ Chí Minh Khoa Khoa Học Kỹ Thuật Máy Tính Mục lục Giới thiệu Biến đổi khí hậu 2.1 Khái niệm 2.2 Nguyên nhân 2.3 Hậu biến đổi khí hậu 3 Mô 3.1 3.2 3.3 3.4 4 7 8 Mơ hình tốn học cải tiến 4.1 Tiếp cận trực tiếp 4.2 Tiếp cận gián tiếp 9 Mơ 5.1 5.2 5.3 hình hóa tốn: Mơ hình EFOM MARKAL quốc gia Mơ hình EFOM-MARKAL đa quốc gia Thuế suất thống thị trường cho giấy phép phát thải Thanh toán chuyển khoản 3.4.1 Phương pháp thị trường cho phép phát thải 3.4.2 Phương pháp dựa giá trị Shapley nước 10 Mơ hình hợp tác quốc tế 10 Đánh giá hợp tác 11 Chia sẻ cổ tức hợp tác 12 Kết luận PPNCKH 13 Trang 1/14 Mơ Hình Tốn Học Cải Tiến Đánh Giá Lợi Ích Hợp Tác Quốc Tế Trong Việc Giảm Thiểu Khí C02 Ngày tháng năm 2015 Tóm tắt nội dung Hiện tình hình biến đổi khí hậu diễn biến vơ phức tạp mà ngun nhân gây khí nhà kính đặc biệt khí CO2 thải mơi trường với mức độ ngày nghiêm trọng Nhằm khuyến khích quốc gia tích cực tham gia giảm phát xạ khí CO2 Bài báo đề xuất phương pháp ghép mơ hình lượng quốc gia nhằm đánh giá lợi ích có quốc gia tham gia xử lý khí thải CO2 Mơ hình đánh giá dựa tham gia nước Châu Âu: Hà Lan, Bỉ, Thụy sĩ, Đức Keywords: Kiểm sốt lượng khí thải CO2 , mơ hình lượng , kỹ thuật giải toán học Trường Đại Học Bách Khoa Tp.Hồ Chí Minh Khoa Khoa Học Kỹ Thuật Máy Tính Giới thiệu Mục đích báo gồm phần: • Đề xuất phương pháp để ghép mơ hình lượng để đánh giá lợi ích hợp tác việc giảm khí thải gây nhiễm mơi trường • Mơ hình tốn học kỹ thuật phân tích giải mơ hình lượng đa quốc gia • Mơ hợp tác bốn nước châu Âu Bỉ, Đức, Hà Lan Thụy Sĩ việc kiềm chế lượng khí thải CO2 họ Biến đổi khí hậu Trước vào nội dung báo Ta tìm hiểu số khái niệm ảnh hưởng biến đổi khí hậu 2.1 Khái niệm Biến đổi khí hậu biến đổi khác thường điều kiện khí hậu bầu khí Trái Đất kéo theo tác động tiêu cực lên nhiều phần Trái Đất 2.2 Nguyên nhân Một nguyên nhân dẫn đến tượng biến đổi khí hậu gia tăng phát xạ khí nhà kính (Hình1):CO2 , CH4 , N2 O, HF Cs, P F Cs SF6 , đặc biệt khí C02 từ người PPNCKH Trang 3/14 Trường Đại Học Bách Khoa Tp.Hồ Chí Minh Khoa Khoa Học Kỹ Thuật Máy Tính Hình 1: Khí nhà kính 2.3 Hậu biến đổi khí hậu • Mưa acid, thủng tần Ozon, cháy rừng, lũ lụt, hạn hán, sa mạc hóa, băng tan • Hệ sinh thái bị phá hủy • Mất đa dạng sinh học • Thiệt hại kinh tế • Dịch bệnh Mơ hình hóa tốn: Trong phần này,chúng ta xem xét mơ hình EFOM MARKAL quốc gia Sau đó, mơ hình kết hợp EFOM MARKAL đa quốc gia, thấy mơ hình có liên quan bối cảnh thuế suất thống chung quốc gia thị trường giấy phép phát thải cân Cuối cùng, đề xuất hai phương pháp để tính tốn chuyển khoản khoản chi phí cần thiết (trong phương pháp tính thuế) để tạo cơng lợi ích nước hợp tác với 3.1 Mơ hình EFOM MARKAL quốc gia EFOM MARKAL mơ hình bottom-up hay cịn gọi tổng hợp quy nạp Mơ hình mô thị trường cạnh tranh hãng lượng công nghệ lượng PPNCKH Trang 4/14 Trường Đại Học Bách Khoa Tp.Hồ Chí Minh Khoa Khoa Học Kỹ Thuật Máy Tính quốc gia khu vực định Họ dùng số lượng lớn thiết bị hỗ trợ người sử dụng, để đáp ứng dịch vụ lượng cạnh tranh sản xuất hãng tiếp tế lượng Trong MARKAL có hoạt động, có liên quan với cơng nghệ giai đoạn: đầu tư, lực sản xuất Những hoạt động ảnh hưởng từ đầu tư, bảo trì hoạt động chi phí tương ứng Trong EFOM, hệ thống lượng mô tả mạng lưới theo định hướng, liên kết đại diện cho dịng lượng hình thành, nút đại diện cho công nghệ lượng Các hoạt động sau gia kết dòng lượng liên kết gia tăng lực cần thiết để thực dòng chảy Việc phát thải CO2 tỷ lệ thuận với việc tiêu thụ nhiên liệu hóa thạch mức độ hoạt động công nghệ mà nhiên liệu sử dụng EFOM MARKAL mô lâu dài tiến hóa hệ thống lượng để đáp ứng cách ngoại sinh mức quy định dịch vụ lượng Những yêu cầu đưa điều khoản nhu cầu kinh tế-xã hội khác Giả định nhu cầu kinh tế-xã hội tỷ lệ nghịch với đàn hồi việc cân giá lượng Từ quan điểm này, EFOM MARKAL phần đề xuất lượng cân thị trường,thu thông qua việc tối ưu hóa sử dụng lượng Để tinh sử dụng mơ MARKAL-MACRO, nơi mà nhu cầu dịch vụ lượng biến nội sinh xác định kinh tế vĩ mơ EFOM MARKAL mơ hình quy hoạch tuyến tính đa chiều (LP) có cơng thức viết sau: min{cT x : Ex ≤ e, P x ≤ p}, x Trong đó: • x vector hoạt động cho hệ thống lượng quốc gia • c vector chi phí cho hoạt động • E hệ số phát thải CO2 ma trận • e vector mục tiêu phát thải CO2 • P mơ tả ràng buộc ma trận hệ thống lượng • p Right-hand-side vector.[3] Trong EFOM, biến số dòng chảy lượng di chuyển qua hệ thống lực cần thiết để vận chuyển dòng chảy Trong MARKAL,các biến mô tả chủ yếu công nghệ.Hơn nội dung chi tiết ràng buộc P x ≤ p khác hai mơ hình, đặc biệt việc xây dựng chế cân PPNCKH Trang 5/14 Trường Đại Học Bách Khoa Tp.Hồ Chí Minh Khoa Khoa Học Kỹ Thuật Máy Tính 3.2 Mơ hình EFOM-MARKAL đa quốc gia Bắt đầu từ khác mơ hình EFOM mơ hình MARKAL quốc gia, chúng tơi thiết kế mơ hình EFOM - MARKAL đa quốc gia lớn Nó liên kết mơ hình quốc gia thơng qua tập ràng buộc tồn cầu khí thải CO2, cho phép nghiên cứu hợp tác quốc tế để hài hòa nỗ lực giảm khí thải.Dạng tuyến tính mơ hình EFOM MARKAL cho N quốc gia biểu đồ hóa sau : N cTr xr (1a) cTr xr (1b) Pr xr ≤ pr = N (1c) r=1 N r=1 Trong • r=1 N số quốc gia tham gia mơ hình • cr vector chi phí hoạt động hệ thống lượng • xr vector cấp độ • Er hệ số phát thải CO2 ma trận • Pr ràng buộc ma trận mô tả hệ thống lượng • pr liên quan đến vector bên phải • e vector mục tiêu tổng thể cho lượng CO2 thải N nước với Ràng buộc (1c) ràng buộc cục bộ, chúng tương ứng với mơ hình EFOM MARKAL cho N quốc gia khơng hạn chế khí thải có khối cấu trúc chéo Ràng buộc (1b) gọi ràng buộc tồn cục hạn chế, chúng liên kết nước N Họ đảm bảo tổng lượng khí thải CO2 nước N giảm xuống mức e Chú ý hạn chế toàn cục mở rộng để đối phó với việc nhiễm khơng khí xun biên giới, trao đổi lượng N quốc gia [1] Hàm mục tiêu vấn đề đưa (1a) Nó tương ứng để giảm thiểu tổng chi phí hệ thống cho N nước Các chi phí chuyển đổi trước với loại tiền tệ tăng thêm lên Bài báo giả định kinh tế N nước tương đối có trạng thái phát triển Để đưa vào tài khoản kinh tế quốc gia phát triển khác nhau, chi phí quốc gia có trọng lượng khác nhau, sử dụng tỷ suất điện mua khách Cuối cùng, lưu ý tổng giá trị hàm mục tiêu vơ nghĩa, mơ hình PPNCKH Trang 6/14 Trường Đại Học Bách Khoa Tp.Hồ Chí Minh Khoa Khoa Học Kỹ Thuật Máy Tính quốc gia khơng thiết chiếm chi phí cố định Tuy nhiên, biến thiên hàm mục tiêu có giá trị hai tình khác (ví dụ, có khơng có phối hợp sách) có ý nghĩa đại diện cho tiết kiệm thêm chi phí phát sinh 3.3 Thuế suất thống thị trường cho giấy phép phát thải Các mô hình đa quốc gia EFOM-MARKAL (1) cho phép xem xét nhóm quốc gia, mà họ cam kết giảm lượng khí thải CO2 họ theo tỷ lệ định Để thực cam kết này,mỗi nước r = 1, ,N đạt mục tiêu CO2 tương ứng er Có nhiều cách hiệu cho nước tồn đạt hàm mục tiêu tổng thể N e= er cách chia sẻ giảm lượng khí thải để chi phí xử lý cận biên nước r=1 Điều thu cách áp đặt khoản thuế CO2 thống nước Chúng tơi tính thuế cách giải (1) Thuế giá trị đưa biến tối ưu kép kết hợp với ràng buộc (1b) Thật vậy, đặt [1] vấn đề (1) có: N (cr − ErT v ∗ )T xr :Pr xr ≤ pz (xr )N r=1 Trong trường hợp không suy biến, giống nguyên thủy giải pháp tối ưu kép Nhận thấy vấn đề sau tương ứng với tình hình thuế thống v> áp đặt r nước khí thải Er xr Thơng qua thuế này, nước r có biên độ tương đối cao chi phí xử lý - kiềm chế phát thải cá nhân họ Phát nhiều so với mục tiêu ban đầu họ er Sau quốc gia mua (tương ứng, bán) - chi phí xử lý cận biên (tương ứng, thấp) cao - giấy phép phát thải thị trường hình thành nước N 3.4 Thanh toán chuyển khoản Giả sử nhóm quốc gia N cam kết kiềm chế lượng khí thải CO2 họ lượng khí định (giảm đến mức e = N r=1 er )bẳng cách sử dụng loại thuế thống Cách hiệu nước tự quy định mức giảm, mang lại lợi ích tài gọi "cổ tức hợp tác" Tuy nhiên, theo thuế suất thống vậy, số nước r số N thải nhiều mục tiêu ban đầu er nước lại phải bồi thường cách giảm nhiều khí thải lại Để đảm bảo lợi ích cho nước tham gia giảm lượng khí thải phương pháp địi hỏi phân chia cơng khoản tốn chuyển giao nước Việc tính tốn thể qua công thức sau: PPNCKH Trang 7/14 Trường Đại Học Bách Khoa Tp.Hồ Chí Minh Khoa Khoa Học Kỹ Thuật Máy Tính 3.4.1 Phương pháp thị trường cho phép phát thải Có cơng thức là: t∗r = v ∗T (er − Er x∗r ) (2) • t∗r (t∗r ≥) khoản nhận • t∗r (t∗r ≤) khoản trả • v ∗ x∗r vector tối ưu 3.4.2 Phương pháp dựa giá trị Shapley Đây khái niệm lý thuyết phân bổ lợi ích hợp tác đáp ứng cơng cho tất nước tham gia Công thức: t∗r = v(r) − cTr (x∗r − xG r ) (3) • t∗r (t∗r ≥) khoản nhận • t∗r (t∗r ≤) khoản trả • v(r) hàm tính tốn khả giảm lợi nhuận tham gia hợp tác • x∗r vector tối ưu • G thể số liên minh nước Lưu ý sử dụng công thức thời gian tính tốn tăng lên nhanh chóng số lượng nước liên minh tăng Hai phương pháp cách tính khác nhau, tùy vào trường hợp mà ta sử dụng Đối với phương pháp thị trường cho phép tất nước tham gia khơng có việc khơng tham gia mà có lợi ích, cịn phương pháp tính tốn theo giá trị Shapley có quyền khỏi liên minh sau tính tốn khoản lợi nhuận chi trả cần thiết Mơ hình tốn học cải tiến Mơ hình EFOM-MARKAL đưa vấn đề tốn qui hoạch tuyến tính lớn Cụ thể với mơ hình tính tốn cho nước nêu trên, ma trận tạo đạt đến 51,900 hàng 54,800 cột Vì độ phức tạp tăng lên nhanh chóng Để giải tốn, tác giả phân tích theo hướng tiếp cận PPNCKH Trang 8/14 Trường Đại Học Bách Khoa Tp.Hồ Chí Minh Khoa Khoa Học Kỹ Thuật Máy Tính 4.1 Tiếp cận trực tiếp Trong hướng tiếp cận này, phương pháp chủ yếu sử dụng "interior point" "simplex" Theo đánh giá, phương pháp interior point dù có ưu điểm so với simplex, có nhu cầu cần nhớ vượt khả cung cấp máy tính đơn lẻ Để vượt qua trở ngại này, số tác giả khác sử dụng biện pháp song song hóa Tuy nhiên, vào thời điểm năm 1998, phần cứng phần mềm hỗ trợ cho giải pháp song song hóa chưa phổ biến 4.2 Tiếp cận gián tiếp Phương pháp chọn để sử dụng "phân rã" Ý tưởng chủ đạo mơ hình hóa tốn ban đầu vốn chứa nhiều biến thành tốn có nghiệm Bài tốn đảm bảo tính "lồi" khơng khả phân Những kĩ thuật đặc biệt dùng để giải khó khăn vấn đề tốn học Ngun tắc "phân rã" giải thích với ý nghĩa kinh tế sau: vấn đề có người điều phối quốc tế N đại diện quốc gia Người điều phối có nhiệm vụ quản lý tài nguyên chung (ví dụ lượng khí thải tối đa phép phát xạ khí quyển) cố gắng để thỏa mãn ràng buộc Mỗi quốc gia xem xét mơ hình EFOM MARKAL tìm cách tối thiểu hóa chi phí hệ thống Để đạc mức tối ưu toàn cục, cần phải thiết lập giao tiếp mơ hình điều chỉnh đơn vị điều phối quốc gia Điều thực cách thiết lập tập mức giá đơn vị điều phối quy định ứng với quyền sử dụng tài ngun chung mơ hình phân rã theo hướng chi phí, hay mức cung cấp tài nguyên quốc gia, mơ hình phân rã theo tài nguyên Tại bước, quan điều phối tạo tập giá quyền phối phối tài ngun chung tùy theo mơ hình phân rã chọn gửi thông tin đến quốc gia Sau đó, quốc gia sử dụng thơng tin để xác định hành động tối ưu cho quốc gia họ, gửi lại cho quan điều phối Sự trao đổi thông tin tiếp diễn đạt cân bằng, nghĩa hành động tối ưu quốc gia không cịn thiện cho tồn cục Trong mơ hình phân rã giá, việc sử dụng tập giá v ≥ thiết lập quan điều phối ràng buộc Lagrang, vấn đế mơ hình lại sau Trong (dr )N = -(Cr )N Hàm mục tiêu diễn giải tối đa hóa tổng lợi nhuận dTr xr , trừ chi phí v để sử dụng Er xr tài nguyên chung v xem thuế thống quan điều phối cho phát thải khí CO2 nhằm đảm bảo N quốc gia phát thải nhiều e Do đó, cách phân rã tương ứng với cách tiếp cận thuế thống N min(eT v + v≥0 PPNCKH max{(dr − ErT v)T xr : Pr xr ≤ pr }) r=1 r=1 (4) Trang 9/14 Trường Đại Học Bách Khoa Tp.Hồ Chí Minh Khoa Khoa Học Kỹ Thuật Máy Tính Trong phương pháp phân rã theo tài ngun, quan điều phối phải tính tốn phân bổ tối ưu (e∗r )N = (Er x∗r )N tài nguyên chung N quốc gia, theo tối ưu đạt khii quốc gia tối ưu hàm sau với er = e∗r : gr (er ) = min{cTr xr : Er xr ≤ er , Pr xr ≤ pr } xr (5) Vấn đề tổ chức điều phối tương ứng với công thức sau: N { (er )N N er ≤ e} gr (er ) : r=1 (6) r=1 Sự chia sẻ tài nguyên chung quan điều phối đến quốc gia, lương phát thải khí CO2 phép Phương pháp phân rã theo tài nguyên tương ứng với hướng tiếp cận theo thị trường hạn ngạch khí thải Chức thị trường mô trao đổi thông tin quan điều phối quốc gia Sự tối ưu phân bổ khí thải quốc gia tương ứng chia hạn ngạch điểm cân thị trường Một lợi ích chung mơ hình phân rã giảm thiểu nhu cầu nhớ cách chia toán ban đầu toán nhỏ Bên cạnh đó, phương pháp phân rã cho phép dễ dàng thực song song hóa Mơ nước Để minh họa cho lợi ích thu từ việc hợp tác xử lý khí CO2 , báo mô dựa hợp tác bốn nước châu Âu: Bỉ, Đức, Hà Lan Thụy Sĩ Trong phần mô tả ngắn gọn mô hình quốc gia sử dụng để hạn chế phát xạ khí CO2 , đánh giá lợi ích việc hợp tác, chi trả chuyển giao khí CO2 để đạt công cho nước tham gia 5.1 Mơ hình hợp tác quốc tế Dưới biểu đồ cho thấy lượng khí thải CO2 nước Châu Âu chưa có kết hợp mơ hình lượng Ba nước Bỉ, Hà Lan, Thụy sĩ sử dụng mơ hình MARKAL cịn Đức mơ hình EFOM Quan sát biểu đồ ta thấy đường gạch đứt thể mong đợi tổng lượg khí thải nước giảm 30% năm 2020 Nhưng thực tế nước sử dụng mơ hình riêng tổng lượng khí thải CO2 giảm Đồng thời nhận thấy lượng khí thải CO2 nước Đức cịn cao nhiều so với ba nước lại Để đạt mục tiêu giảm 30% lượng khí thải CO2 tạo cơng lợi ích kinh tế nước mơ hình cải tiến EFOM-MARKAL đời PPNCKH Trang 10/14 Trường Đại Học Bách Khoa Tp.Hồ Chí Minh Khoa Khoa Học Kỹ Thuật Máy Tính Hình 2: Biểu đồ phát xạ khí CO2 5.2 Đánh giá hợp tác Dưới bảng thể lợi nhuận nước có Nhận thấy chi phí chi trả cho việc Hình 3: Bảng giảm thiểu chi phí lợi ích hợp tác thải khí CO2 nước Bỉ, Hà Lan Thụy sĩ chưa hợp tác nhiều sau hợp tác, cịn Đức ngược lại Nhắm khuyến khích Đức tiếp tục tham gia mơ hình ba nước cịn lại đồng ý chia lợi nhuận với Đức, với mục đích bên có lợi cách cho phép mua bán khí thải CO2 Nghĩa theo quy định nước phải tuân theo hạn mức cho phép thải khí CO2 , nước vượt q phải đóng thuế Khi tham gia mơ hình cải tiến Bỉ, Hà Lan Thụy sĩ đạt mức hạn mức cho phép nên tiết kiệm khoản chi phí, nước dùng khoản tiền mua lượng khí thải mà nước Đức thải ra, nhằm chia sẻ lợi ích cho nước Đức tham gia mơ hình PPNCKH Trang 11/14 Trường Đại Học Bách Khoa Tp.Hồ Chí Minh Khoa Khoa Học Kỹ Thuật Máy Tính 5.3 Chia sẻ cổ tức hợp tác Ở phần báo đưa bảng cho thấy lợi ích nước tham gia mơ hình cải tiến EFOM-MARKAL • Đầu tiên quan sát bảng lợi ích tính tốn thị trường cho phép ta thấy Đức nhận từ Bỉ $1,433, từ Hà Lan $4,874 từ Thụy sĩ $2,822 Lúc lợi nhuận Đức thu $2,193 khơng cịn -$6936 Trong cách tính đảm bảo cơng cho Đức Bỉ khơng tính khác Hình 4: Bảng lợi ích hợp tác thị trường cho phép Do cách cách tính khác đưa • Tính tốn lợi ích theo giá trị Shapley Đối với cách tính đảm bảo cơng cho nước tham gia mơ hình liên kết, số lượng nước tham gia tăng lên việc tính tốn khó khăn Hình 5: Bảng lợi ích hợp tác với giá trị Shapley PPNCKH Trang 12/14 Trường Đại Học Bách Khoa Tp.Hồ Chí Minh Khoa Khoa Học Kỹ Thuật Máy Tính Kết luận Hợp tác quốc tế ủng hộ hiệp định UNFCCC lựa chọn hiệu để giảm lượng khí thải CO2 Các mơ hình lượng EFOM MARKAL cơng cụ tốt để đánh giá chi phí xử lý CO2CO2 Và kết hợp mơ hình EFOM-MARKAL (mơ hình xem xét thuế suất thống quốc tế thị trường giấy phép phát thải quốc tế) Đã đánh giá lợi ích hợp tác xử lý chất thải CO2 nước khác cung cấp nhìn sâu sắc tác động hợp tác sách kế hoạch lượng quốc gia Mơ hình EFOM-MARKAL vấn để cần giải lập trình tốn học Để làm vậy, tác giả sử dụng kỹ thuật phân rã để tách mơ hình hợp tác đa quốc gia Sau tiếp tục giải ACCPM, thuật tốn tiên tiến để tối ưu hóa Các ưu điểm phương pháp phân rã giảm nhớ khả tính tốn song song Kết trình bày minh họa cho lợi ích hợp tác quốc tế cắt giảm khí thải CO2 Sự hợp tác mang lại lợi ích cho quốc gia nhờ vào khoản toán bồi thường PPNCKH Trang 13/14 Trường Đại Học Bách Khoa Tp.Hồ Chí Minh Khoa Khoa Học Kỹ Thuật Máy Tính Tài liệu [1] O Bahn, Techniques de decomposition dans l’analyse des choix én- ergétiques et environnementaux, Ph.D thesis, University of Geneva,Geneva (1994) [2] G.B Dantzig, Linear Programming and Extensions(Princeton University Press, Princeton, 1963) [3] Right hand vector http : //www.maths.unsw.edu.au/sites/def ault/f iles/M atlabSelf P aced/lesson6/M atlabLesson6L ine PPNCKH Trang 14/14 ... 12 Kết luận PPNCKH 13 Trang 1/14 Mô Hình Tốn Học Cải Tiến Đánh Giá Lợi Ích Hợp Tác Quốc Tế Trong Việc Giảm Thiểu Khí C02 Ngày tháng năm 2015 Tóm tắt nội dung Hiện tình hình biến đổi... cho lợi ích hợp tác quốc tế cắt giảm khí thải CO2 Sự hợp tác mang lại lợi ích cho quốc gia nhờ vào khoản toán bồi thường PPNCKH Trang 13/14 Trường Đại Học Bách Khoa Tp.Hồ Chí Minh Khoa Khoa Học. .. lên việc tính tốn khó khăn Hình 5: Bảng lợi ích hợp tác với giá trị Shapley PPNCKH Trang 12/14 Trường Đại Học Bách Khoa Tp.Hồ Chí Minh Khoa Khoa Học Kỹ Thuật Máy Tính Kết luận Hợp tác quốc tế