Thực trạng ô nhiễm asen trong nước giếng khoan dùng cho ăn uống và nguy cơ sức khỏe của người dân xã chuyên ngoại, duy tiên hà nam

BỘ GIÁO DỤC VÀ ĐÀO TẠO BỘ Y TẾ TRƢỜNG ĐẠI HỌC Y TẾ CÔNG CỘNG BÙI HUY TÙNG H P THỰC TRẠNG Ô NHIỄM ASEN TRONG NƢỚC GIẾNG KHOAN DÙNG CHO ĂN UỐNG VÀ NGUY CƠ SỨC KHỎE CỦA NGƢỜI DÂN XÃ CHUYÊN NGOẠI, DUY TIÊN, HÀ NAM NĂM 2011 H U LUẬN VĂN THẠC SĨ Y TẾ CÔNG CỘNG Mã số: 60.72.76 HÀ NỘI - 2011 BỘ GIÁO DỤC VÀ ĐÀO TẠO BỘ Y TẾ TRƢỜNG ĐẠI HỌC Y TẾ CÔNG CỘNG BÙI HUY TÙNG H P THỰC TRẠNG Ô NHIỄM ASEN TRONG NƢỚC GIẾNG KHOAN DÙNG CHO ĂN UỐNG VÀ NGUY CƠ SỨC KHỎE CỦA NGƢỜI DÂN XÃ CHUYÊN NGOẠI, DUY TIÊN, HÀ NAM NĂM 2011 U H LUẬN VĂN THẠC SĨ Y TẾ CÔNG CỘNG Mã số: 60.72.76 Hƣớng dẫn khoa học: TS Nguyễn Việt Hùng HÀ NỘI - 2011 LỜI CẢM ƠN Tôi xin trân trọng bày tỏ lòng biết ơn tới: Người hướng dẫn khoa học TS Nguyễn Việt Hùng tận tình dạy dỗ tơi q trình học tập nghiên cứu, người thầy trực tiếp hướng dẫn thực đề tài ThS Trần Thị Tuyết Hạnh, giảng viên Bộ môn Sức khỏe môi trường, Đại học Y tế công cộng hướng dẫn nhiều chuyên mơn suốt q trình làm luận văn H P Viện Y học lao động Vệ sinh môi trường - Bộ Y tế, Trung tâm Nghiên cứu Công nghệ Môi trường Phát triển bền vững (CETASD), Trường Đại học Khoa học tự nhiên - Đại học Quốc gia Hà Nội giúp đỡ tơi q trình xét nghiệm đảm bảo chất lượng xét nghiệm mẫu sinh học Cục Quản lý Môi trường y tế - Bộ Y tế, Trung tâm Y tế dự phòng tỉnh Hà U Nam, Trung tâm Y tế dự phòng huyện Duy Tiên, Ủy ban nhân dân xã Chuyên Ngoại, Trạm Y tế xã Chuyên Ngoại, hộ gia đình xã Chuyên Ngoại giúp H đỡ tạo điều kiện cho tơi q trình xác định vấn đề nghiên cứu triển khai nghiên cứu địa bàn xã Chuyên Ngoại, huyện Duy Tiên, tỉnh Hà Nam Trung tâm quốc gia lực nghiên cứu Bắc - Nam Thụy Sĩ (NCCR N-S) hỗ trợ kinh phí tạo điều kiện cho thực đề tài nghiên cứu Tôi xin chân thành cảm ơn Trường Cao đẳng Y tế Hà Nội tạo điều kiện cho tơi học, cảm ơn gia đình, bạn bè, đồng nghiệp giúp đỡ, động viên tinh thần để hoàn thành luận văn Hà Nội, ngày 29 tháng năm 2011 Bùi Huy Tùng i DANH MỤC CÁC CHỮ VIẾT TẮT As : Asen AAS : Atomic Absorption Spectrophotometry (Quang phổ hấp thụ nguyên tử) BYT : Bộ Y tế CETASD : Center for Environmental Technology and Sustainable Development (Trung tâm Nghiên cứu Công nghệ Môi trường Phát triển bền vững) H P CFS : Cancer Slope Factor (Hệ số ước tính nguy ung thư) EHRA : Environmental Health Risk Assessment (Đánh giá nguy sức khỏe mơi trường) HGĐ : Hộ gia đình LD50 : Lethal dose 50% (Liều gây chết nửa) LOAEL : Lowest Observed Adverse Effect Level U (Liều thấp gây đáp ứng có hại quan sát được) MRL : Minimum Risk Level (Mức nguy tối thiểu) H NCCR N-S : National Centre of Competence Research North - South (Trung tâm quốc gia lực nghiên cứu Bắc - Nam Thụy Sĩ) NOAEL : No Observed Adverse Effect Level (Liều khơng gây đáp ứng có hại quan sát được) ppb : Part per billion (Phần tỉ, ppb = 10-3 mg/l = µg/l) QCVN : Quy chuẩn Việt Nam SKMT : Sức khỏe môi trường SMR : Standardized Mortality Ratios (Tỷ lệ tử vọng chuẩn hóa) TCCP : Tiêu chuẩn cho phép TCVN : Tiêu chuẩn Việt Nam TCVS : Tiêu chuẩn vệ sinh TDI : Tolerable Daily Intake (Mức tiêu thụ hàng ngày chịu đựng được) ii TTYTDP : Trung tâm Y tế dự phòng UNICEF : United Nations International Children’s Fund (Quỹ Nhi đồng Liên hợp quốc) USEPA : United States Environmental Protection Agency (Cục Bảo vệ Môi trường Mỹ) WHO : World Health Organization (Tổ chức Y tế Thế giới) YHLĐ : Y học lao động H P H U iii DANH MỤC CÁC BẢNG Bảng 1.1 Cơ sở giá trị hƣớng dẫn cho asen 24 Bảng 1.2 Quy ƣớc yếu tố không chắn ƣớc tính quần thể ngƣời25 Bảng 1.3 MRL asen [30] 26 Bảng 3.1 Số thành viên HGĐ số thành viên thƣờng xuyên ăn uống nhà 45 Bảng 3.2 Giới tính ngƣời lấy tóc 45 Bảng 3.3 Điều kiện kinh tế hộ gia đình 46 Bảng 3.4 Nguồn nƣớc sử dụng cho mục đích ăn uống 47 H P Bảng 3.5 Tuổi trung bình giếng khoan bể nƣớc mƣa 47 Bảng 3.6 Thể tích bể chứa nƣớc mƣa nƣớc giếng khoan 47 Bảng 3.7 Kết xét nghiệm định lƣợng asen mẫu nƣớc 49 Bảng 3.8 Chênh lệch hàm lƣợng asen trƣớc sau lọc, hiệu lọc khử asen bể lọc HGĐ 49 U Bảng 3.9 Hàm lƣợng asen trung bình nƣớc giếng khoan trƣớc lọc thôn 50 H Bảng 3.10 Hàm lƣợng asen trung bình nƣớc giếng khoan trƣớc lọc vị trí ngồi đê sông Hồng 50 Bảng 3.11 Tỷ lệ mẫu nƣớc giếng khoan đạt TCVS nƣớc sinh hoạt 51 Bảng 3.12 Tổng thời gian sử dụng nƣớc giếng khoan cho ăn uống 54 Bảng 3.13 Liều ƣớc lƣợng asen đƣa vào thể ngày qua đƣờng ăn uống với ngƣời trƣởng thành 55 Bảng 3.14 Một số đặc điểm bể lọc nƣớc giếng khoan HGĐ 57 Bảng 3.15 Các vật liệu đƣợc HGĐ sử dụng để lọc nƣớc giếng khoan 58 Bảng 3.16 Bề dày lớp lọc thời gian thay/ rửa vật liệu lọc 59 Bảng 3.17 Mối liên quan hàm lƣợng asen nƣớc giếng khoan trƣớc lọc với độ sâu giếng khoan 59 Bảng 3.18 Mối liên quan hiệu lọc với việc làm giàn mƣa 59 iv Bảng 3.19 Phân tích hồi quy tuyến tính hàm lƣợng asen nƣớc giếng khoan trƣớc lọc với tuổi giếng khoan 60 Bảng 3.20 Phân tích hồi quy tuyến tính biến: Hiệu lọc, bề dày lớp lọc, thời gian thay/ rửa vật liệu lọc 60 Bảng 3.21 Hàm lƣợng asen tóc ngƣời dân sử dụng nƣớc giếng khoan bị ô nhiễm asen cho ăn uống 61 Bảng 3.22 Tỷ lệ mẫu tóc đạt tiêu chuẩn cho phép tiêu chuẩn chẩn đoán nhiễm độc 61 Bảng 3.23 Mối liên quan hàm lƣợng asen tóc với việc sử dụng nƣớc H P giếng khoan bị ô nhiễm asen 62 Bảng 3.24 Phân tích hồi quy tuyến tính biến: Hàm lƣợng asen tóc, hàm lƣợng asen nƣớc giếng khoan sau lọc, tổng thời gian sử dụng nƣớc giếng khoan cho ăn uống tuổi ngƣời lấy tóc 62 Bảng 3.25 Mô tả hàm số với thông tin đánh giá phơi nhiễm 65 Bảng 3.26 Nguy ung thƣ ƣớc tính 66 H U v DANH MỤC CÁC BIỂU ĐỒ Biểu đồ 1.1 Tỷ lệ giếng khoan bị ô nhiễm asen số tỉnh Việt Nam 12 Biểu đồ 3.1 Số lƣợng thành viên hộ gia đình 45 Biểu đồ 3.2 Tuổi ngƣời lấy tóc 45 Biểu đồ 3.3 Trình độ học vấn ngƣời trả lời vấn 46 Biểu đồ 3.4 Mục đích sử dụng nguồn nƣớc mƣa nƣớc giếng khoan 47 Biểu đồ 3.5 Năm xây dựng giếng khoan bể nƣớc mƣa 47 Biểu đồ 3.6 Tỷ lệ mẫu nƣớc giếng khoan bị ô nhiễm asen 50 Biểu đồ 3.7 Tỷ lệ HGĐ sử dụng nƣớc giếng khoan cho ăn uống theo H P tháng năm 54 Biểu đồ 3.8 Số tháng sử dụng nƣớc giếng khoan cho ăn uống năm 54 Biểu đồ 3.9 So sánh liều ƣớc lƣợng asen đƣa vào thể ngày với mức TDI WHO 55 Biểu đồ 3.10 Tƣơng quan hàm lƣợng asen tóc hàm lƣợng asen U nƣớc sau lọc 63 Biểu đồ 3.11 Nguy ung thƣ xác suất xảy 68 H vi DANH MỤC CÁC HÌNH Hình 1: Khung đánh giá nguy sức khỏe môi trƣờng Hình Các đƣờng thâm nhập asen vào thể ngƣời [11] 17 Hình Vị trí địa lý xã Chuyên Ngoại, Duy Tiên, Hà Nam 29 H P H U vii MỤC LỤC DANH MỤC CÁC CHỮ VIẾT TẮT i DANH MỤC CÁC BẢNG iii DANH MỤC CÁC BIỂU ĐỒ v DANH MỤC CÁC HÌNH vi MỤC LỤC vii TÓM TẮT NGHIÊN CỨU x ĐẶT VẤN ĐỀ MỤC TIÊU NGHIÊN CỨU H P Chƣơng TỔNG QUAN TÀI LIỆU Phƣơng pháp đánh giá nguy sức khỏe môi trƣờng 1.1 Khung lý thuyết đánh giá nguy sức khỏe môi trường 1.2 Áp dụng phương pháp đánh giá nguy sức khỏe môi trường độc chất asen .7 U Vấn đề sử dụng nguồn nƣớc ngầm bị ô nhiễm asen 2.1 Thực trạng ô nhiễm asen nước ngầm giới H 2.2 Thực trạng ô nhiễm asen nước ngầm Việt Nam .11 2.3 Người dân Hà Nam vấn đề sử dụng nước giếng khoan bị ô nhiễm asen 14 Xác định yếu tố nguy 15 3.1 Đặc tính lý hóa tồn asen môi trường .15 3.2 Đường phơi nhiễm, tích lũy đào thải asen thể .17 3.3 Ảnh hưởng sức khỏe cấp tính asen .18 3.4 Ảnh hưởng sức khỏe mạn tính asen 19 Đánh giá mối quan hệ Liều - Đáp ứng 23 Chƣơng PHƢƠNG PHÁP NGHIÊN CỨU 28 Đối tƣợng nghiên cứu 28  Không dùng nước mưa Q13 Ông/bà dùng nước mưa vào mục đích gì? (nhiều lựa chọn)  Ăn uống  Giặt giũ  Tắm gội   Tưới Rửa rau thực phẩm khác  Mục đích khác: Q14 Ơng/bà có bể lọc nước giếng khoan khơng?   Có Khơng (=> CHUYỂN CÂU) Q15 Bể lọc xây vào năm nào? _ H P Q16 Thể tích bể chứa nước giếng khoan sau lọc khối? m³ Q17 Ơng/bà có làm giàn mưa khơng?   Có Khơng Q18 Ơng/bà có để lắng nước giếng khoan khơng?  U  Có Khơng Q19 Bể lọc gia đình gồm thành phần nào? Độ dày thành phần này? H  Cát (nói chung)  Cát vàng  Cát đen  Sỏi  Than dày _ cm  Đá dày _ cm  Khác: dày _ cm dày _ cm dày _ cm dày _ cm dày _ cm Q20 Ông/bà xây bể lọc theo hướng dẫn từ đâu?  Tự thiết kế  Theo truyền thông, hướng dẫn y tế   Làm theo hộ khác Khác: _ Q21 Từ xây bể lọc, ông/bà thay lớp lọc chưa?   Có Khơng (=> CHUYỂN CÂU) Q22 Lần thay lớp lọc gần cách tháng? _ Q23 Ơng/bà làm với lớp cát cũ thay lớp lọc?   Đãi rửa cát dùng lại Bỏ bỏ đâu (ghi rõ): _ H P Q24 Ơng/bà có bể chứa nước mưa khơng?   Có Khơng (=> CHUYỂN CÂU) Q25 Bể chứa nước mưa xây vào năm nào? _ Q26 Thể tích bể chứa nước mưa gia đình khối? _ m³ H U PHỤ LỤC XÁC ĐỊNH ASEN TRONG NƢỚC BẰNG KIT THỬ NHANH CỦA HÃNG MERCK Hố chất dụng cụ Bộ kít thử gồm có: - 01 lọ dung dịch hóa chất số (asen-1) - 01 lọ hóa chất số (asen-2) - 01 lọ hóa chất số (asen-3) H P - 01 bình phản ứng có nắp xốy - 01 thìa đong hóa chất màu đỏ - 01 thìa đong hóa chất màu xanh - 01 bảng màu chuẩn để xác định hàm lượng với thang (mg/l): (0 - 0,005 - 0,01 0,025 - 0,05 - 0,1 - 0,25 - 0,5) Các hố chất khác khơng bao gồm kít thử: dung dịch asenic chuẩn U CertiPUR@, 1000 mg/l asen, mã đặt hàng 1.19773 Các bƣớc tiến hành phân tích H Mẫu nước lấy vào chai nhựa sạch, ghi nhãn theo quy định phải phân tích vịng sau lấy mẫu Bước 1: Lấy vào bình phản ứng 60ml mẫu cần đo (lấy mẫu đến vạch chai) Bước 2: Cho thêm vào bình phản ứng 02 giọt tác nhân (asen-1) lắc Bước 3: Cho thêm vào thìa đỏ (bằng miệng thìa) tác nhân (asen-2) lắc cho tan hoàn toàn Bước 4: Cho thêm vào thìa xanh (bằng miệng thìa) tác nhân (asen-3) đóng nắp bình phản ứng lại Bước 5: Lật gá đỡ màu đen nắp bình phản ứng cho chấm trắng quay phía người ngồi thao tác Lấy que thử hộp đưa vào bình phản ứng qua khe hẹp nắp bình (phần phản ứng đưa vào trước) Đưa vào đến vạch mầu xanh que thử ngang bề mặt nắp bình gập gá đỡ màu đen lại Que thử gá đỡ giữ chặt nguyên vị trí Thao tác từ mở gá đỡ đến cho que thử vào bình phản ứng phải thật nhanh để tránh sai lệch kết ảnh hưởng đến sức khoẻ người phân tích Bước 6: Để bình phản ứng vòng 20 phút lắc lần vịng 20 phút Chú ý: tránh khơng để que thử tiếp xúc với dung dịch nước bên Bước 7: Sau 20 phút lấy que thử ra, so sánh với bảng màu kèm theo để xác định hàm lượng mg/l asen3+/5+ tương ứng bảng màu Nếu màu không tương đương hồn tồn với vạch bảng ta đọc ước lượng giá trị trung bình hai giá trị hai vạch màu màu que thử H P H U PHỤ LỤC PHƢƠNG PHÁP PHÂN TÍCH ASEN TRONG DUNG DỊCH BẰNG QUANG PHỔ HẤP THỤ NGUYÊN TỬ Trang thiết bị dụng cụ cần thiết - Máy quang phổ hấp thụ nguyên tử (ASS) Model AA-6800 ghép nối thiết bị tạo khí hydrua HVG-1, đèn catơt rỗng asen cuvet thạch anh, hãng sản Shimadzu - Lò vi sóng hãng National Model NN-C988W (năng lượng tối đa 900W) H P kèm theo ống Teflon MPR 600/10M (10 chiếc) - Tủ sấy, tủ hút - Bếp điện điều nhiệt, nhiệt kế - Các bình Kendan - Pipet thủy tinh 5, 10mL; micropipet 100, 200, 1000, 5000 µL - Các bình định mức nhựa P.P 10, 25, 100 mL U - Các ống nghiệm nhựa 15, 50 mL (Greiner), giá để ống nghiệm - Các lọ nhựa PP đựng dung dịch chuẩn Hóa chất H - Dung dịch chuẩn gốc As (III) g/l, Fluka - Muối As(V) Na2HAsO4.7H2O, pA, Fluka - Các axit HCl 36%, HNO3 65%, H2SO4 98% dung dịch H2O2 30%, pA, Merck - Các hóa chất NaI, NaOH, NaBH4, axit amido sulfonic HSO3NH2, axit L (+) ascorbic, pA, Merck - Các dung dịch chuẩn kiểm chứng ICP multi-element (As 102 mg/L); As(V) 1g/l pha từ muối Na2HAsO4, Merck - Nước cất hai lần - Mẫu tóc kiểm chứng NCS ZC 81002 (GBW 09101), Trung tâm Phân tích Quốc gia Trung Quốc - Khí acetilen 99%, argon 99,9% - Dầu gội trung tính (có thể dùng cồn/ete) Đƣờng chuẩn Dựng đường chuẩn As với nồng độ 0; 1; 2,5; 5; 7,5; 10 µg/L (khoảng tuyến tính asen xác định thiết bị HVG-AAS 6800 - 10 µg/L) Đường chuẩn đường bậc với hệ số tương quan R ≥ 0,995 * Chuẩn bị dung dich chuẩn As (III) mg/L Pha loãng 0,1 ml As (III) chuẩn (1000 mg/l) thành 100 ml HCl 2%  dung dịch As chuẩn mg/l (để tuần) H P * Chuẩn bị dung dịch As (III) 10, 25, 50, 75, 100 µg/L Lấy thể tích 0,1; 0,25; 0,5; 0,75; ml dung dịch As chuẩn 1mg/l vào bình định mức 10ml, định mức tới vạch nước cất * Đường chuẩn: - Lấy 5ml dung dịch As chuẩn (0, 10, 25, 50, 75, 100 µg/L) vào ống U nghiệm chia vạch 50ml, thêm 9ml HCl 8M - Định mức đến 50ml axit HCl 2% (0,55M)  Nồng độ HCl cuối là: 1,8 M H Chuẩn bị dụng cụ kiểm chứng * Dung dịch chuẩn ICP: dung dịch chuẩn gồm 32 nguyên tố khác biết trước nồng độ Chúng sử dụng dung dịch để kiểm tra dộ xác việc chuẩn bị dung dịch chuẩn, đường chuẩn sử dụng để định lượng kết phân tích As dung dịch ICP ngoại suy từ đường chuẩn nằm khoảng 80 - 120% * Chuẩn bị dung dịch chuẩn ICP 2,5 µg As/L - Pha loãng 0,25 ml dung dịch ICP As 102 mg/l thành 25 ml nước cất Dung dịch chuẩn 1,02 mg/L - Sau lấy 0,25 ml dung dịch ICP 1,02 mgAs/L cho vào bình 10 ml, định mức H2O  Dung dịch chuẩn ICP 25,5 µg/L - Lấy ml dung dịch ICP 25,5 µg As/L, thêm 9ml HCl 8M, pha loãng thành 50ml HCl 2%  dung dịch chuẩn ICP 2,55 µg As/L Các thông số đo As thiết bị HVG-AAS - Vạch đo: 193,7 nm - Khe đo: 0,5nm - Cường độ dòng catot rỗng: 10 mA (Imax= 40 mA) - Chiều cao Burner head: 16mm - Dạng lửa: Không khí - C2H2 - Tốc độ khí cháy C2H2: 2,0 L/phút H P - Áp suất đầu vào thiết bị AAS C2H2: - Áp suất đầu vào thiết bị AAS khơng khí: - Tốc độ khí agon: 70 mL/phút H U 0,045 Mpa (0,45 atm) 0,25 Mpa (2,5 atm) PHỤ LỤC KẾT QUẢ XÉT NGHIỆM HÀM LƢỢNG ASEN TRONG NƢỚC GIẾNG KHOAN VÀ TRONG TÓC XÃ CHUYÊN NGOẠI, HUYỆN DUY TIÊN, TỈNH HÀ NAM Mã phiếu 001 002 003 004 005 006 007 008 009 010 011 012 013 014 015 016 017 018 019 020 021 022 023 024 025 026 027 028 029 030 031 032 033 Mẫu nƣớc trƣớc lọc (ppb) 10 135 323 328 417 144 262 213 335 296 422 197 86 53 111 323 240 429 295 319 338 271 321 146 362 340 216 289 332 88 114 289 Mẫu tóc (mg/kg) H P U H Mẫu nƣớc sau lọc (ppb) 0 11 14 12 42 15 12 28 20 11 27 20 10 0 12 21 18 17 38 10 37 28 16 0,04 0,04 0,63 0,02 0,13 0,17 0,08 0,50 0,03 0,63 0,04 0,04 0,17 0,25 0,58 0,58 0,38 0,29 0,17 0,17 0,58 0,21 0,54 0,17 0,33 0,21 0,04 0,04 0,21 1,00 0,63 0,29 0,13 034 035 036 037 038 039 040 041 042 043 044 045 046 047 048 049 050 051 052 053 054 055 056 057 058 059 060 061 062 063 064 065 066 067 068 069 070 071 072 073 074 075 407 398 266 283 349 446 82 291 306 36 383 213 336 384 471 336 259 333 447 474 484 325 406 287 276 329 273 181 453 389 353 261 254 197 579 315 213 476 504 138 204 341 H P U H 32 16 11 38 20 40 25 12 13 33 42 32 19 99 107 21 32 14 21 60 82 36 14 19 163 33 68 16 23 40 32 20 39 33 23 0,08 0,04 0,04 0,50 0,29 0,25 0,17 0,21 0,33 0,02 0,25 0,63 0,63 0,13 0,33 1,46 0,21 0,21 0,21 0,38 0,17 0,17 0,17 1,58 0,50 0,58 0,08 0,13 0,25 0,54 0,13 0,63 0,29 0,08 0,42 0,21 0,25 0,02 0,08 0,25 0,13 0,33 076 077 078 079 080 081 082 083 084 085 086 087 088 089 090 091 092 093 094 095 096 097 098 099 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 438 426 325 316 318 355 475 207 215 186 214 72 213 242 206 195 280 231 280 296 453 412 460 324 502 494 319 410 310 83 291 156 349 241 398 276 403 248 414 310 354 238 H P U H 21 13 30 23 18 12 33 14 38 58 12 12 44 54 33 15 19 64 42 13 23 15 20 29 66 18 25 34 146 15 15 28 28 13 47 11 26 0,13 0,42 0,29 0,08 0,25 0,04 0,63 0,83 0,46 0,42 0,79 1,79 0,29 0,08 0,93 0,38 0,38 0,21 0,21 0,01 0,75 0,02 0,25 0,46 0,02 0,17 0,17 0,29 1,50 0,63 0,04 0,01 0,25 0,04 0,42 0,21 0,42 0,63 0,25 0,17 0,13 0,21 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 214 187 183 199 188 171 177 354 217 323 304 333 417 374 400 332 298 571 259 255 89 164 267 520 365 71 391 520 417 430 561 392 507 H P U H 46 50 44 54 57 39 26 20 16 33 15 25 41 31 59 38 53 23 17 33 13 17 25 52 48 27 24 15 0,67 0,08 0,08 0,54 0,29 0,42 0,13 0,38 0,02 0,21 0,13 0,63 0,17 0,54 0,63 0,46 0,21 0,29 0,25 0,42 0,08 0,21 0,13 0,21 0,08 0,54 0,02 0,04 0,42 0,50 0,25 0,21 0,02 PHỤ LỤC QUY TRÌNH TÍNH TỐN NGUY CƠ THEO PHƢƠNG PHÁP XÁC SUẤT Bƣớc Mô tả phân bố hàm lƣợng asen mẫu phân tích phơi nhiễm Để thực phương pháp xác suất, thông tin phơi nhiễm hàm lượng asen mẫu phân tích sử dụng dạng hàm phân bố thay giá trị trung bình Sau xác định hàm lượng asen mẫu phân tích, thời gian H P phơi nhiễm (số ngày phơi nhiễm năm, số năm phơi nhiễm) tương ứng với mẫu môi trường, tập hợp dãy số liệu chuyển thành cột phần mềm Microsoft Excel 2007 Chạy phần mềm @Risk 5.5 với lệnh @Risk/Model/Fit distribution to data cho cột số liệu Phần mềm @risk mô tả cho biết hàm phân bố tối ưu U cho dãy số liệu đưa vào (được mơ hình bên dưới) H Hình Hàm phân bố đƣợc hiển thị @Risk 5.5 cho dãy số liệu Hàm lượng asen sau lọc (đơn vị: mg/l) H P Hình Hàm phân bố đƣợc hiển thị @Risk 5.5 cho dãy số liệu Số ngày phơi nhiễm năm (đơn vị: ngày) U H Hình Hàm phân bố đƣợc hiển thị @Risk 5.5 cho dãy số liệu Số năm phơi nhiễm (đơn vị: năm) Bƣớc Tính tốn nguy Các kết bước copy vào Microsofr Excel 2007 thao tác Write to cell Các kết hiển thị ô nhỏ Excel dạng trị số trung bình Sau lập cơng thức tính nguy (đã trình bày Chương Phương pháp nghiên cứu; Phần Phương pháp xử lý phân tích số liệu) ta tính nguy tương ứng từ số liệu Nguy dạng hàm số, hiển thị Excel giá trị trung bình H P U Hình Thiết lập cơng thức tính nguy Microsoft Excel 2007 tích hợp @Risk 5.5 H Bƣớc Chạy phân bố nguy theo phƣơng pháp xác suất phần mềm @Risk 5.5 Tại ô kết vừa thu Bước 2, chạy lệnh @Risk/Simulation/Start (Mô theo Monte Carlo với 10.000 lần hoán vị) ta phân bố xác suất nguy mơ tả theo hình bên dưới: H P Hình Phân bố xác suất xảy nguy theo @Risk 5.5 Hình cho ta biết giá trị là: nguy xác suất xảy Với giá trị U nguy cụ thể trục hoành ta có xác suất xảy tương ứng trục tung Đồ thị cho biết tính biến đổi nguy giới hạn đường cong nguy Xác suất xảy trường hợp dao động từ đến Nguy H phần diện tích phía đường cong có tổng diện tích