ĐẠI HỌC QUỐC GIA HÀ NỘI
TRƯỜNG ĐẠI HỌC KHOA HỌC TỰ NHIÊN
---------------------
Đào Thị Hồng Vân
NGHIÊN CỨU KHẢ NĂNG ỨNG DỤNG MÔ HÌNH WRF-CHEM
VÀO KHU VỰC VIỆT NAM
Chuyên ngành: Khí tượng và khí hậu học
Mã số: 62.44.87
LUẬN VĂN THẠC SĨ KHOA HỌC
NGƯỜI HƯỚNG DẪN KHOA HỌC:
GS. TS. Phan Văn Tân
Hà Nội - 2013
2
LỜI CẢM ƠN
Trong thời gian qua, việc nghiên cứu các nội dung luận văn luôn được thuận
lợi và đúng hướng ngoài sự nỗ lực của cá nhân, học viên còn nhận được sự chỉ bảo,
giúp đỡ nhiệt tình và sâu sắc của thầy cô, bạn bè và đồng nghiệp.
Học viên xin gửi lời cảm ơn đến thầy giáo hướng dẫn, GS.TS. Phan Văn Tân
đã những chỉ bảo, định hướng quan trọng trong toàn bộ quá trình thực hiện. Học
viên xin bày tỏ lòng biết ơn chân thành tới thầy, người đã hết lòng quan tâm cũng
như kiên trì giúp đỡ từng bước nghiên cứu của học viên.
Học viên cũng xin gửi lời cảm ơn đến các Thầy, Cô trong Bộ môn khí tượng
nói riêng Khoa Khí tượng Thủy Văn Hải dương học i chung đã truyền đạt
các kiến thức bổ ích cho học viên trong toàn bộ khóa học.
Ngoài ra, học viên cũng xin chân thành cảm ơn tập thể Trung tâm Nghiên
cứu Quan trắc hình hóa Môi trường đã luôn giúp đỡ, động viên tạo điều
kiện để học viên hoàn thành luận văn, hơn thế nữa còn mang lại một môi trường
làm việc hiếm có cho học viên.
Học viên
Đào Thị Hồng Vân
3
MỤC LỤC
DANH MỤC HÌNH .................................................................................................... 4
Chương 1 TNG QUAN .......................................................................................... 10
1.1 Các nghiên cứu trên thế giới ............................................................................... 10
1.2 Các nghiên cứu trong nước ................................................................................. 21
Chương 2 PHƯƠNG PHÁP NGHIÊN CỨU ........................................................... 24
2.1 Sơ lược về mô hình WRF/Chem ......................................................................... 24
2.1.1 Mô hình WRF .................................................................................................. 24
2.1.2 Mô đun CHEM ................................................................................................ 27
2.2 Thiết kế thí nghiệm ............................................................................................. 29
2.2.1 Miền tính và thời gian thí nghiệm .................................................................... 29
2.2.2 Các thí nghiệm ................................................................................................. 30
2.2.3 Nguồn số liệu ................................................................................................... 33
2.4 Tạo bộ số liệu phát thải cho WRF/Chem ............................................................ 34
Chương 3 KT QU VÀ NHẬN XÉT .................................................................... 36
3.1 Đánh giá kết quả của thí nghiệm I ...................................................................... 36
3.1.1 Hoàn lưu, nhiệt độ và lượng mưa từ đầu ra của WRF/Chem .......................... 36
3.1.2 Mô phỏng nồng độ bụi từ WRF/Chem ............................................................ 43
3.2 Đánh giá kết quả của thí nghiệm II ..................................................................... 49
3.2.1 Trường nhiệt độ và lượng mưa với các tùy chọn của WRF/Chem .................. 49
3.2.2 Mô phỏng các chất phát thải từ WRF/Chem ................................................... 58
KT LUN ............................................................................................................... 65
TÀI LIỆU THAM KHO ........................................................................................ 66
PH LC .................................................................................................................. 69
4
DANH MỤC HÌNH
Hình 1.1 So sánh sự thay đổi lượng mưa của các tháng 6,7,8 trong thí nghim A
(hình trên) khi tính đến ảnh hưởng của BC thí nghiệm B (hình dưới) khi
không tính đến ảnh hưởng ca BC (Menon ccs., 2002) ........................................... 12
Hình 1.2 Sự thay đổi lượng phát thải NOx Bc Mỹ, Châu Âu Châu Á từ năm
1970 đến 2000 (Akimoto, 2003) ............................................................................... 13
Hình 2.1. Cấu trúc tng quan của mô hình WRF ..................................................... 25
Hình 2.2. Sơ đồ h thống mô hình WRF/Chem phiên bản 3.4. ................................ 28
Hình 2.3 Miền tính của WRF trong các thí nghiệm. Độ phân gii ngang 30 km ..... 29
Hình 3.1 Trường lượng mưa trung bình mực t đầu ra ca WRF_DUST (trái)
WRF_NOCHEM (phi) của các ngày 02 04/01/2006 (t trên xuống dưới) ...... 37
Hình 3.2 Trường nhiệt độ không khí trung bình mực 2m t đầu ra ca WRF_DUST
(trái) và WRF_NOCHEM (phi) của ngày 02 và 04/01/2006 (từ trên xuống dưới) 38
Hình 3.3 Trường nhiệt độ không khí trung bình mực 2m t đầu ra ca WRF_DUST
(trái), số liu APHRODITE (giữa) và hiệu giữa chúng (phải) của các ngày t 01 đến
04/01/2006 (t trên xuống dưới) ............................................................................... 39
Hình 3.4 Trường lượng mưa trung bình mực t đầu ra ca WRF_DUST (trái), số
liu APHRODITE (giữa) và hiệu gia chúng (phải) của các ngày t 01 đến
04/01/2006 (t trên xuống dưới) ............................................................................... 41
Hình 3.5 Trường độ cao địa thế v trường gió tại mc 850 mb t đầu ra ca
WRF/Chem (bên trái) so sánh với NNRP (bên phải) t 01 đến 04/01/2006 ............ 43
Hình 3.6 Phân bố bi loi 1 (DUST_01) trường gió tại mực 850 mb lúc 00, 06,
12, 18h t 02/01/2006 đến 04/01/2006 mô phỏng bi WRF/Chem ......................... 46
Hình 3.7 Mt cắt kinh hướng phân bố bi loi 1 (DUST_01) nhân với gió kinh
ớng (trung bình từ 102 đến 110E) lúc 00, 12h từ 01/01/2006 đến 04/01/2006 .... 48
Hình 3.8 Trường nhiệt độ mc 2m t đầu ra ca WRF với các tùy chọn hóa học
300, 301, 11 (t trái qua phải) tr đi WRF_NOCHEM của các ngày t 01 đến
04/01/2006 (t trên xuống dưới) ............................................................................... 50
5
Hình 3.9 Trường lượng mưa từ đầu ra ca WRF với các tùy chọn hóa học 300, 301,
11 (t trái qua phải) tr đi WRF_NOCHEM của các ngày từ 01 đến 04/01/2006 (t
trên xuống dưới) ........................................................................................................ 52
Hình 3.10 Trường lượng mưa trung bình mực t đầu ra ca WRF_C300,
WRF_C301, WRF_C011 (t trái qua phải) tr đi số liu APHRODITE ca các
ngày từ 01 đến 04/01/2006 (t trên xuống dưới) ...................................................... 54
Hình 3.11 Trường nhiệt độ mc 2m t đầu ra ca WRF_C300, WRF_C301,
WRF_C011 (t trái qua phải) tr đi số liu APHRODITE của các ngày từ 01 đến
04/01/2006 (t trên xuống dưới) ............................................................................... 56
Hình 3.12 Profile nhiệt độ t đầu ra ca WRF_C300, WRF_C301, WRF_C011 so
sánh với WRF_NOCHEM của các ngày từ 01 đến 04/01/2006 ............................... 57
Hình 3.13 Profile của t s xáo trộn hơi nước t đầu ra ca WRF_C300,
WRF_C301, WRF_C011 so sánh với WRF_NOCHEM của các ngày từ 01 đến
04/01/2006 ................................................................................................................ 58
Hình 3.14 Profile ca bi PM2.5 t đầu ra ca WRF_C300, WRF_C301,
WRF_C011 của các ngày từ 01 đến 04/01/2006 ...................................................... 59
Hình 3.15 Profile ca bi PM10 t đầu ra ca WRF_C300, WRF_C301,
WRF_C011 của các ngày từ 01 đến 04/01/2006 ...................................................... 60
Hình 3.16 Profile ca nồng độ SO2 t đầu ra ca WRF_C300, WRF_C301,
WRF_C011 của các ngày từ 01 đến 04/01/2006 ...................................................... 61
nh 3.17 Phân bố ca nồng độ PM2.5 t đầu ra ca WRF_C300, WRF_C301,
WRF_C011 (trái qua phải) mc 1000 (trên) và 850 mb (dưới) ngày 04/01/2006 ... 62
Hình 3.18 Phân bố ca nồng độ PM10 t đầu ra ca WRF_C300, WRF_C301,
WRF_C011 (trái qua phải) mc 1000 (trên) và 850 mb (dưới) ngày 04/01/2006 ... 63
Hình 3.19 Phân bố ca nồng độ SO2 t đầu ra ca WRF_C300, WRF_C301,
WRF_C011 (trái qua phải) mc 850 (trên) và 500 mb (dưới) ngày 04/01/2006 ..... 64
Hình P.1 Trường độ cao địa thế v và trường gió tại mc 500 mb t đầu ra ca WRF
(bên trái) so sánh với s liệu NNRP (bên phải) t 01 đến 04/01/2006 ..................... 70