Nâng cao độ chính xác ổn định cho thiết bị quan trắc và cảnh báo phóng xạ môi trường

Nâng cao độ chính xác ổn định cho thiết bị quan trắc và cảnh báo phóng xạ môi trường Nâng cao độ chính xác ổn định cho thiết bị quan trắc và cảnh báo phóng xạ môi trường Nâng cao độ chính xác ổn định cho thiết bị quan trắc và cảnh báo phóng xạ môi trường luận văn tốt nghiệp,luận văn thạc sĩ, luận văn cao học, luận văn đại học, luận án tiến sĩ, đồ án tốt nghiệp luận văn tốt nghiệp,luận văn thạc sĩ, luận văn cao học, luận văn đại học, luận án tiến sĩ, đồ án tốt nghiệp

ĐẠI HỌC QUỐC GIA HÀ NỘI TRƯỜNG ĐẠI HỌC KHOA HỌC TỰ NHIÊN - Nguyễn Thị Minh NÂNG CAO ĐỘ CHÍNH XÁC, ỔN ĐỊNH CHO THIẾT BỊ QUAN TRẮC VÀ CẢNH BÁO PHĨNG XẠ MƠI TRƯỜNG LUẬN VĂN THẠC SĨ KHOA HỌC Hà Nội - 2014 ĐẠI HỌC QUỐC GIA HÀ NỘI TRƯỜNG ĐẠI HỌC KHOA HỌC TỰ NHIÊN - Nguyễn Thị Minh NÂNG CAO ĐỘ CHÍNH XÁC, ỔN ĐỊNH CHO THIẾT BỊ QUAN TRẮC VÀ CẢNH BÁO PHÓNG XẠ MÔI TRƯỜNG Chuyên ngành: Vật lý nguyên tử Mã số: 60440106 LUẬN VĂN THẠC SĨ KHOA HỌC NGƯỜI HƯỚNG DẪN KHOA HỌC: TS Đặng Quang Thiệu Hà Nội - 2014 Luận văn tốt nghiệp LỜI CẢM ƠN Luận văn kết trình học tập suốt hai năm trường Đại học Khoa học tự nhiên – Đại học Quốc gia Hà Nội trình làm luận văn thân Viện Khoa học Kỹ thuật hạt nhân – Viện Năng lượng nguyên tử Việt Nam Lời đầu tiên, em xin gửi lời cảm ơn chân thành đến TS Đặng Quang Thiệu - Trung tâm Chiếu xạ Hà Nội, tận tình giảng dạy, truyền đạt kiến thức chuyên ngành hướng dẫn em hoàn thành luận văn Đồng thời, em xin gửi lời cảm ơn đến Th.S Nguyễn Thị Bảo Mỹ - Viện Khoa học Kỹ thuật hạt nhân, nhiệt tình giúp đỡ bảo thêm kiến thức chuyên ngành suốt trình nghiên cứu đề tài Em xin cảm ơn thầy giảng dạy lớp cao học khóa 2011 - 2013, đặc biệt thầy cô môn Vật lý hạt nhân - Trường Đại học Khoa học tự nhiên tận tình dạy dỗ, tạo điều kiện thuận lợi cho em suốt thời gian học tập trường Em gửi lời cảm ơn đến gia đình, bạn bè bên cạnh em, động viên, giúp em vượt qua khó khăn để hồn thành đề tài Mặc dù nỗ lực cố gắng, chắn luận văn không tránh khỏi thiếu sót Rất mong nhận ý kiến đóng góp, bổ sung thầy bạn bè Hà Nội, tháng 03 năm 2014 Học viên Nguyễn Thị Minh i Nguyễn Thị Minh Luận văn tốt nghiệp MỤC LỤC MỤC LỤC DANH MỤC CÁC BẢNG BIỂU .3 DANH MỤC CÁC HÌNH VẼ MỞ ĐẦU CHƯƠNG 1: TỔNG QUAN VỀ THIẾT BỊ QUAN TRẮC VÀ CẢNH BÁO PHÓNG XẠ MÔI TRƯỜNG 1.1 Nhu cầu thực tiễn 1.2 Cơ sở lý thuyết 10 1.2.1 Hiệu ứng quang điện .10 1.2.2 Hiệu ứng Compton 12 1.2.3 Hiệu ứng tạo cặp electron - positron .13 1.2.4 Tổng hợp hiệu ứng gamma tương tác với vật chất .14 1.2.5 Cấu trúc phổ gamma .15 1.3 Thiết bị quan trắc cảnh báo phóng xạ mơi trường 17 1.3.1 Đầu dị chứa khí .19 a) Buồng ion hóa 22 b) Ống đếm tỉ lệ .22 c) Ống đếm Geiger-Muller (G-M) 22 1.3.2 Đầu dò bán dẫn .23 1.3.3 Đầu dò nhấp nháy 25 CHƯƠNG 2: ĐẦU DÒ NHẤP NHÁY NaI(Tl) 26 2.1 Đầu dò nhấp nháy .26 2.1.1 Chất nhấp nháy vô 26 2.1.2 Chất nhấp nháy hữu 27 2.2 Đầu dò nhấp nháy NaI(Tl) 29 2.2.1 Đặc điểm cấu tạo .29 2.2.2 Sự hình thành xung lối 33 2.3 Sử dụng phương pháp chuyển phổ thành liều với đầu dò NaI(Tl) 36 CHƯƠNG 3: NÂNG CAO ĐỘ CHÍNH XÁC VÀ ỔN ĐỊNH CHO THIẾT BỊ40 Nguyễn Thị Minh Luận văn tốt nghiệp 3.1 Các yếu tố ảnh hưởng tới đầu dò NaI(Tl) 40 3.1.1 Ảnh hưởng việc xác định đỉnh 40 3.1.2 Ảnh hưởng nhiệt độ 40 3.2 Nâng cao độ xác 41 3.2.1 Làm trơn phổ 42 3.2.2 Xác định đỉnh 44 3.2.3 Chuyển phổ thành liều .49 3.3 Nâng cao độ ổn định .51 3.3.1 Ảnh hưởng nhiệt độ đỉnh phổ K-40 51 3.3.2 Bù nhiệt độ cho thiết bị 53 KẾT LUẬN 55 TÀI LIỆU THAM KHẢO 56 PHỤ LỤC 57 Nguyễn Thị Minh Luận văn tốt nghiệp DANH MỤC CÁC BẢNG BIỂU Bảng 2.1: Các hệ số A(K) cho tinh thể NaI(Tl) hình trụ 2.5 cm x 2.5 cm 39 Bảng 3.1: Kết đo suất liều thiết bị vị trí có suất liều chuẩn khác 50 Bảng 3.2: Kết đo đỉnh K-40 nhiệt độ khác Nguyễn Thị Minh 52 Luận văn tốt nghiệp DANH MỤC CÁC HÌNH VẼ Hình 1.1: Mơ hình tán xạ Compton 12 Hình 1.2: Mơ hình tạo cặp electron - positron 13 Hình 1.3: Sự phụ thuộc tiết diện vào lượng [1] 14 Hình 1.4: Cấu trúc phổ lý tưởng tia gamma theo hiệu ứng: a) quang điện; b) tạo cặp; c) Compton, d) phổ thực [1] 15 Hình 1.5: Sơ đồ khối thiết bị quan trắc cảnh báo phóng xạ mơi trường 17 Hình 1.6: Một số thiết bị đo xạ cầm tay 18 Hình 1.7: Hình thiết bị đo xạ đặt trạm quan trắc Nhật Bản 19 Hình 1.8: Các tượng xảy với cặp ion dương điện tử điện cực hình trụ với cường độ xạ hiệu điện thay đổi [9] 20 Hình 1.9: Mơ hình buồng ion hóa 22 Hình 1.10: Đầu dị bề mặt 24 Hình 1.11: Đầu dị Si Ge-photon 24 Hình 2.1: Sự phụ thuộc hệ số suy giảm tinh thể NaI(Tl) vào lượng tia gamma [4] 30 Hình 2.2: Lý thuyết thực nghiệm phổ tia gamma gây tương tác Comton quang điện đầu dị NaI(Tl) [5] 31 Hình 2.3: Biên độ xung theo số tích phân [3, 6] 36 Hình 3.1: Sơ đồ q trình nâng cao độ xác cho thiết bị quan trắc cảnh báo phóng xạ mơi trường 41 Hình 3.2: Sơ đồ q trình xác định đỉnh phương pháp fit Gauss 46 Hình 3.3: Số liệu thu 47 Nguyễn Thị Minh Luận văn tốt nghiệp Hình 3.4: Số liệu thu biểu diễn theo thang log 47 Hình 3.5: Phổ đo sau làm trơn 48 Hình 3.6: Phổ đo sau fit Gauss 49 Hình 3.7: Phổ K-40 nhiệt độ T = 20oC T = 26oC 52 Hình 3.8: Sự phụ thuộc đỉnh K-40 vào nhiệt độ 53 Nguyễn Thị Minh Luận văn tốt nghiệp MỞ ĐẦU Ngày nay, lượng nguyên tử dần trở thành nguồn lượng thay cho nguồn lượng hóa thạch dần cạn kiệt Không vậy, việc sử dụng nguồn lượng hạt nhân cịn cách để giảm thiểu lượng khí thải CO2, nguyên nhân gây hiệu ứng nhà kính nóng lên trái đất Khơng vậy, xạ hạt nhân cịn sử dụng rộng dãi nhiều lĩnh vực khác y tế, khoa học, quân sự… đem lại lợi ích to lớn đời sống Tuy nhiên, lượng nguyên tử tiềm ẩn nguy hiểm người Thực tế, chứng kiến ảnh hưởng nghiêm trọng chúng vụ nổ bom nguyên tử Tokyo Hirosima, vụ nổ nhà máy điện hạt nhân Chernobyl gần nhà máy điện hạt nhân Fukushima Chính vậy, yêu cầu đặt nhà môi trường nhà quản lý phải kiểm soát liều lượng phóng xạ mơi trường để có biện pháp ứng phó kịp thời với cố liên quan tới phóng xạ hạt nhân Các thiết bị quan trắc mơi trường cảnh báo phóng xạ cơng cụ để nhà quản lý nhà mơi trường theo dõi kiểm sốt nhiễm phóng xạ Vì hầu quan tâm tới việc xây dựng trạm quan trắc cảnh báo phóng xạ mơi trường Hiện nay, Nhật Bản có 37 trạm quan trắc, Hàn Quốc có 13 trạm, Ấn Độ có 16 trạm,… Tại Việt Nam, có trạm quan trắc phóng xạ môi trường Viện Khoa học kỹ thuật hạt nhân, Viện Nghiên cứu hạt nhân thuộc Viện Năng lượng nguyên tử, Trung tâm Công nghệ xử lý môi trường (Bộ Quốc phòng) Nguyễn Thị Minh Luận văn tốt nghiệp Qua thực tế hoạt động, trạm quan trắc cảnh báo phóng xạ mơi trường nước ta bộc lộ nhiều hạn chế như: chưa đồng bộ, khả thu thập phân tích chưa đáp ứng đầy đủ tiêu tính liên tục, độ nhạy theo tiêu chuẩn quốc tế Ngoài trạm chưa có chức cảnh báo trực tuyến cố rò rỉ để phục vụ cho việc ứng phó với trường hợp khẩn cấp Nắm vai trò quan trọng trạm quan trắc cảnh báo phóng xạ mơi trường, Viện Năng lượng nguyên tử giao cho Viện Khoa học kỹ thuật hạt nhân nghiên cứu xây dựng thiết bị quan trắc cảnh báo phóng xạ mơi trường Qua trình nghiên cứu triển khai đo đạc thực nghiệm, tốn đặt cho nhà nghiên cứu làm để nâng cao độ xác ổn định cho thiết bị đo đạc điều kiện môi trường khác Luận văn thực với mục tiêu nâng cao độ xác, tính ổn định thiết bị quan trắc cảnh báo phóng xạ mơi trường Trong sử dụng phương pháp JAERI, phương pháp chuyển phổ thành liều dùng hàm G(E), ổn định phổ phương pháp bù nhiệt độ ghim đỉnh K-40 Luận văn gồm chương với nội dung sau: Chương trình bày tổng quan vai trò, nguyên lý cấu tạo thiết bị quan trắc cảnh báo phóng xạ mơi trường Chương nghiên cứu đầu dò nhấp nháy phương pháp đo số liệu sử dụng loại đầu dò Chương đề xuất phương pháp nâng cao độ xác tính ổn định cho thiết bị Nguyễn Thị Minh Luận văn tốt nghiệp Từ Hình 3.3 3.4, ta dễ dàng thấy phổ thu có lẫn nhiễu tương đối nhiều Do ta cần làm trơn để giảm bớt ảnh hưởng nhiễu (Hình 3.5) Tuy nhiên ta thấy rằng, phổ có đỉnh khoảng kênh gần 400 nên đây, ta cần quan tâm vùng phổ Hình 3.5: Phổ đo sau làm trơn Bằng cách sử dụng thuật tốn bình phương cực tiểu để xác định phân bố Gauss đỉnh, ta thu hàm Gauss sau: y = 77 * e   x − 407   −     16     + 86 (3.12) Qua đó, ta dễ dàng xác định đỉnh phổ độ rộng đỉnh phổ theo (3.11) xP = 407 FWHM = 16 * * ln ( ) ≈ 37 Ta dễ dàng tính độ lệch bình phương số liệu so với phân bố Gauss: 48 Nguyễn Thị Minh Luận văn tốt nghiệp s2 = ∑ f ( x) − y i i N = 25.75 Hay nói cách khác, độ lệch trung bình: s ≈ 5.07 Đồng thời ta có hình biểu diễn số liệu đường fit Hình 3.6 Hình 3.6: Phổ đo sau fit Gauss Ta thấy rằng, đường phổ có phân bố gần trùng khít với phân bố Gauss Điều chứng tỏ việc xác định đỉnh phổ phương pháp fit Gauss có độ xác tương đối tốt 3.2.3 Chuyển phổ thành liều Như trình bày chương 2, đầu dò NaI(Tl) cho phép sử dụng phương pháp chuyển phổ thành liều với ưu điểm xét đến lượng xạ thu thông qua hàm G(E) để xác định suất liều phóng xạ Do vậy, hiển nhiên phương pháp góp phần cải thiện độ xác thiết bị đo quan trắc cảnh báo phóng xạ mơi trường 49 Nguyễn Thị Minh Luận văn tốt nghiệp Quá trình chuyển phổ thành liều thực cách tự động dựa vào đỉnh tính sử dụng phương pháp Fit Gauss trình bày áp dụng vào cơng thức (2.12) Tóm lại, việc kết hợp xác định đỉnh fit hàm Gauss sử dụng phương pháp chuyển phổ thành liều ta tăng độ xác thiết bị quan trắc cảnh báo phóng xạ môi trường lên tương đối cao, đáp ứng yêu cầu thiết bị quan trắc cảnh báo phóng xạ mơi trường sử dụng ngồi thực địa Điều chứng minh thông qua số liệu thực nghiệm Bảng 3.1 Bảng 3.1: Kết đo suất liều thiết bị vị trí có suất liều chuẩn khác STT Suất liều Suất liều Độ lệch Suất liều thiết bị ghi nhận (uSv/hr) chuẩn trung bình chuẩn (uSv/hr) (uSv/hr) (uSv) 11.19 10.63 10.60 10.66 10.57 10.40 10.57 0.10 Độ lệch chuẩn (%) 0.95 Sai số chuẩn (%) -5.54 7.77 7.38 7.47 7.32 7.37 7.46 7.40 0.06 0.81 -4.76 5.71 5.40 5.42 5.51 5.41 5.49 5.45 0.05 0.92 -4.55 4.37 4.16 4.13 4.16 4.21 4.12 4.16 0.04 0.96 -4.81 3.45 3.32 3.28 3.31 3.32 3.34 3.31 0.02 0.60 -4.06 2.80 2.63 2.67 2.71 2.66 2.73 2.68 0.04 1.49 -4.29 1.94 1.91 1.95 1.94 1.89 1.90 1.92 0.03 1.56 -1.03 1.66 1.67 1.62 1.62 1.60 1.68 1.64 0.03 1.83 -1.20 1.43 1.43 1.45 1.42 1.41 1.41 1.42 0.02 1.41 -0.70 10 1.24 1.25 1.25 1.23 1.23 1.23 1.24 0.01 0.81 0.00 11 1.06 1.09 1.11 1.11 1.08 1.09 1.10 0.01 0.91 3.77 12 0.86 0.88 0.88 0.86 0.88 0.88 0.88 0.01 1.14 2.33 13 0.78 0.80 0.79 0.80 0.81 0.81 0.80 0.01 1.25 2.56 14 0.70 0.72 0.74 0.71 0.74 0.73 0.73 0.01 1.37 4.29 15 0.63 0.67 0.65 0.70 0.67 0.68 0.67 0.02 2.99 6.35 16 0.58 0.59 0.64 0.61 0.62 0.62 0.62 0.02 3.23 6.90 17 0.17 0.14 0.15 0.13 0.14 0.14 0.14 0.01 7.14 -19.22 50 Nguyễn Thị Minh Luận văn tốt nghiệp 3.3 Nâng cao độ ổn định Như biết, biên độ lối đầu dò phụ thuộc vào nhiều thông số khác như: nhiệt độ, thời gian sử dụng, nuôi cho PMT Đối với thiết bị sử dụng phòng thí nghiệm, ta liên tục thực chuẩn thiết bị với nguồn chuẩn Tuy nhiên, thiết bị sử dụng ngồi mơi trường vấn đề tương đối phức tạp Quá trình nghiên cứu thực nghiệm rằng, nhiệt độ yếu tố ảnh hưởng lớn gây trơi kết đo Do đó, luận án nghiên cứu ảnh hưởng nhiệt độ tới phổ thu được, đặc biệt với đỉnh K-40 đỉnh phóng xạ tự nhiên xuất nhiều mơi trường Qua đó, đưa hệ số bù nhiệt độ thích hợp 3.3.1 Ảnh hưởng nhiệt độ đỉnh phổ K-40 Sự ảnh hưởng nhiệt độ phổ đo đỉnh phổ K40 mơ tả Hình 3.7 Ta thấy rằng, nhiệt độ giảm xuống phổ thu gần khơng thay đổi, nhiên có dịch chuyển đỉnh phổ tương ứng với thay đổi nhiệt độ 51 Nguyễn Thị Minh Luận văn tốt nghiệp Hình 3.7: Phổ K-40 nhiệt độ T = 20oC T = 26oC Ngồi ra, chúng tơi khảo sát vị trí đỉnh phổ K-40 nhiệt độ khác Kết thực nghiệm ghi lại Bảng 3.2 hình 3.8 với giá trị đỉnh xác định phương pháp Fit Gauss trình bày phần 3.1 Bảng 3.2: Kết đo đỉnh K-40 nhiệt độ khác T 16 Kênh 315 18 20 22 24 26 28 30 32 34 317 320 322 325 327 329 333 336 338 52 Nguyễn Thị Minh Luận văn tốt nghiệp Sự phụ thuộc đỉnh K40 vào nhiệt độ 338 336 334 332 330 Kªnh 328 326 324 322 Sè liÖu y=1.3x+293.8 (R=0.9978) 320 318 316 314 16 18 20 22 24 26 28 30 32 34 NhiƯt ®é (C) Hình 3.8: Sự phụ thuộc đỉnh K-40 vào nhiệt độ Từ Hình 3.8, ta thấy rằng, vị trí đỉnh K-40 phụ thuộc vào nhiệt độ tương đối tuyến tính với mức độ lệch khoảng 1.3 kênh/oC 3.3.2 Bù nhiệt độ cho thiết bị Ta thấy nhiệt độ ảnh hưởng tương đối tuyến tính lên giá trị đo Do đó, ta hồn tồn thêm giá trị xác định vào kết đo để bù lại thay đổi nhiệt độ Giá trị bù xác định sau: ∆kênh = k * (T - T0) (3.13) Trong đó, ∆kênh giá trị bù tương ứng với nhiệt độ T, T0 nhiệt độ mà ta xác định vị trí đỉnh K-40 chuẩn thiết bị K hệ số bù xác định độ dịch đỉnh K-40 theo nhiệt độ (1.3 kênh/oC) Bằng cách bù nhiệt độ vậy, ta hồn tồn hạn chế ảnh hưởng nhiệt độ tới kết đo Ngồi ra, để góp phần nâng cao độ xác, ta kết hợp kiểm tra đỉnh K-40 (xuất đất đá ngồi tự nhiên) vị trí sau bù 53 Nguyễn Thị Minh Luận văn tốt nghiệp Phương pháp ổn định bù nhiệt độ góp phần vào kết thực nghiệm đo Bảng 3.1 54 Nguyễn Thị Minh Luận văn tốt nghiệp KẾT LUẬN Trong thời gian nghiên cứu, thực đề tài: “nâng cao độ xác, ổn định cho thiết bị quan trắc cảnh báo phóng xạ mơi trường”, tác giả thực nội dung sau: - Nghiên cứu, tìm hiểu nguyên lý, cấu tạo thiết bị quan trắc cảnh báo phóng xạ mơi trường Đồng thời tìm hiểu khối đầu dò sử dụng NaI(Tl) - Nghiên cứu phương pháp xác định đỉnh phương pháp fit Gauss sử dụng thuật toán chuyển phổ thành liều để nâng cao độ xác cho thiết bị quan trắc cảnh báo phóng xạ mơi trường - Nghiên cứu phụ thuộc vị trí đỉnh phổ vào nhiệt độ, đề xuất sử dụng bù nhiệt độ để nâng cao tính ổn định cho thiết bị quan trắc cảnh báo phóng xạ mơi trường Kết vậy, hồn toàn phù hợp với mục tiêu đề tài Tuy nhiên, thời gian thực đề tài có hạn nên luận văn số hạn chế chưa khảo sát phụ thuộc vào nhiệt độ đỉnh nguồn phóng xạ khác ngồi K-40 chưa thực khảo sát đo suất liều ứng với nhiệt độ khác 55 Nguyễn Thị Minh Luận văn tốt nghiệp TÀI LIỆU THAM KHẢO Tiếng Việt: Nguyễn Đức Hòa (2012), Điện tử hạt nhân, NXB Giáo dục Việt Nam Trần Đại Nghiệp (2007), Giáo trình xử lý xạ sở công nghệ xạ, NXB Đại học Quốc gia Hà Nội Tiếng Anh: Baba H., Fukuchi T., Kurokawa M and Shimoura S (2004), Study of Digital Pulse Shape Analysis for NaI(Tl) Scintillator, CNS Anual Report Hamamatsu, Photomultiplier Tubes handbook Heath R L (1997), Sintillation Spectrometry gamma-ray spectrum catalogue, Gamma-ray spectrometry center Helmuth S (2002), Pulse processing and analysis, IEEE NPSS short course on Radiation detection and measurement Ianakiev K D., Alexandrov B S., Littlewood P B., Browne M C (2009), “Temperature behavior of NaI (Tl) scintillation detectors”, Nuclear Instruments and Methods in Physics Research Section A: Accelerators, Spectrometers, Đầu dòs and Associated Equipment, Vol 607(2) Kamiya K, Sasatani M, “Effects of radiation exposure on human body”, Nihon rinsho Japanese journal of clinical medicine, Vol 70(3), pp 367-374 LoveLock J E (1961), “Ionization Methods for Analysis of Gases and Vapors”, Analytical Chemistry, Vol 33(2), pp 162-178 10 Miglierini M (2004), Đầu dò of Radiation, E Wigner Course on reactor physics experiments 11 Nelson G., Reilly D (1991), “Gamma-Ray Interactions with Matter”, Passive Nondestructive Analysis of Nuclear Material, Los Almos National Laboratory, pp 27-42 12 Reeder P L., Stromswold D C., “Performance of Large NaI(Tl) GammaRay Detectors Over Temperature -50ºC to +60ºC”, Pacific Northwest National Laboratory 56 Nguyễn Thị Minh Luận văn tốt nghiệp PHỤ LỤC 1-Code chương trình tìm đỉnh tự động sử dụng phần mềm MATLAB clear all; close all; clc; [filename, pathname] = uigetfile('*.mca', 'Select a data file file'); fid=fopen([pathname filename]); fgets(fid); fgets(fid); fgets(fid); fgets(fid); fgets(fid); fgets(fid); fgets(fid); fgets(fid); fgets(fid); fgets(fid); fgets(fid); fgets(fid); fgets(fid); fgets(fid); fgets(fid); fgets(fid); fgets(fid); fgets(fid); y=fscanf(fid,'%d'); fclose(fid); plot(y); x=200:500; figure; yy=smooth(y,5); plot(yy,'r'); xlim([250 600]); [xmin ymin]=ginput(1); [xmax ymax]=ginput(1); xmin=round(xmin); xmax=round(xmax); figure; ydata=yy(xmin:xmax); xdata=xmin:xmax; 57 Nguyễn Thị Minh Luận văn tốt nghiệp ymax=max(ydata(:)); plot(xdata,ydata); [a, resnorm]=lsqcurvefit(@(a,xdata) mygausian(a,xdata),[ymax (xmax+xmin)/2 50 0],xdata',ydata); y2=mygausian(a,xdata); hold on; plot(xdata,y2,'g'); Peak=round(a(2)) FWHM=round(a(3)*2*sqrt(2*log(2))) s2=resnorm/(xmax-xmin) hold off; function f=mygausian(a,xdata) f=a(1)*exp(-1/2*((xdata-a(2))/a(3)).^2)+a(4); end 2- Hệ số G(E) tính cho đầu dò NaI(Tl) 2.5 cm x 2.5 cm E(MeV) 0.050 0.055 0.060 0.065 0.070 0.075 0.080 0.085 0.090 0.095 0.100 0.105 0.110 0.115 0.120 0.125 0.130 0.135 0.140 0.145 0.150 0.155 G(E) 6.425E-03 4.589E-03 4.085E-03 3.966E-03 3.936E-03 3.926E-03 3.937E-03 3.883E-03 4.072E-03 4.208E-03 4.390E-03 4.611E-03 4.855E-03 5.145E-03 5.445E-03 5.758E-03 6.081E-03 6.412E-03 6.749E-03 7.090E-03 7.437E-03 7.790E-03 E(MeV) 0.890 0.895 0.900 0.905 0.910 0.915 0.920 0.925 0.930 0.935 0.940 0.945 0.950 0.955 0.960 0.965 0.970 0.975 0.980 0.985 0.990 0.995 G(E) 2.391E-01 2.406E-01 2.421E-01 2.435E-01 2.451E-01 2.465E-01 2.480E-01 2.494E-01 2.508E-01 2.522E-01 2.537E-01 2.551E-01 2.565E-01 2.578E-01 2.592E-01 2.606E-01 2.619E-01 2.633E-01 2.646E-01 2.660E-01 2.673E-01 2.686E-01 58 Nguyễn Thị Minh E(MeV) 2.460 2.470 2.480 2.490 2.500 2.510 2.520 2.530 2.540 2.550 2.560 2.570 2.580 2.590 2.600 2.610 2.620 2.630 2.640 2.650 2.660 2.670 G(E) 4.884E-01 4.895E-01 4.907E-01 4.918E-01 4.930E-01 4.941E-01 4.953E-01 4.964E-01 4.975E-01 4.987E-01 4.998E-01 5.010E-01 5.021E-01 5.032E-01 5.044E-01 5.055E-01 5.067E-01 5.078E-01 5.089E-01 5.101E-01 5.112E-01 5.124E-01 Luận văn tốt nghiệp 0.160 0.165 0.170 0.175 0.180 0.185 0.190 0.195 0.200 0.205 0.210 0.215 0.220 0.225 0.230 0.235 0.240 0.245 0.250 0.255 0.260 0.265 0.270 0.275 0.280 0.285 0.290 0.295 0.300 0.305 0.310 0.315 0.320 0.325 0.330 0.335 0.340 0.345 0.350 0.355 8.152E-03 8.523E-03 8.906E-03 9.353E-03 9.717E-03 1.015E-02 1.061E-02 1.108E-02 1.159E-02 1.212E-02 1.268E-02 1.327E-02 1.390E-02 1.456E-02 1.525E-02 1.598E-02 1.675E-02 1.756E-02 1.810E-02 1.928E-02 2.020E-02 2.116E-02 2.215E-02 2.319E-02 2.426E-02 2.537E-02 2.651E-02 2.770E-02 2.892E-02 3.017E-02 3.146E-02 3.278E-02 3.414E-02 3.553E-02 6.950E-03 3.840E-02 3.988E-02 4.139E-02 4.293E-02 4.450E-02 1.000 1.010 1.020 1.030 1.040 1.050 1.060 1.070 1.080 1.090 1.100 1.110 1.120 1.130 1.140 1.150 1.160 1.170 1.180 1.190 1.200 1.210 1.220 1.230 1.240 1.250 1.260 1.270 1.280 1.290 1.300 1.310 1.320 1.330 1.340 1.350 1.360 1.370 1.380 1.390 2.699E-01 2.725E-01 2.751E-01 2.776E-01 2.801E-01 2.826E-01 2.850E-01 2.874E-01 2.898E-01 2.922E-01 2.945E-01 2.968E-01 2.991E-01 3.013E-01 3.030E-01 3.057E-01 3.079E-01 3.108E-01 3.121E-01 3.142E-01 3.162E-01 3.183E-01 3.203E-01 3.223E-01 3.243E-01 3.262E-01 3.281E-01 3.301E-01 3.319E-01 3.338E-01 3.357E-01 3.375E-01 3.393E-01 3.411E-01 3.429E-01 3.446E-01 3.463E-01 3.481E-01 3.498E-01 3.515E-01 59 Nguyễn Thị Minh 2.680 2.690 2.700 2.710 2.720 2.730 2.740 2.750 2.760 2.770 2.780 2.790 2.800 2.810 2.820 2.830 2.840 2.850 2.860 2.870 2.880 2.890 2.900 2.910 2.920 2.930 2.940 2.950 2.960 2.970 2.980 2.990 3.000 3.050 3.100 3.150 3.200 3.250 3.300 3.350 5.135E-01 5.147E-01 5.158E-01 5.169E-01 5.181E-01 5.192E-01 5.204E-01 5.215E-01 5.227E-01 5.238E-01 5.249E-01 5.261E-01 5.272E-01 5.284E-01 5.295E-01 5.307E-01 5.318E-01 5.329E-01 5.341E-01 5.352E-01 5.364E-01 5.375E-01 5.387E-01 5.398E-01 5.410E-01 5.421E-01 5.433E-01 5.444E-01 5.455E-01 5.467E-01 5.478E-01 5.490E-01 5.501E-01 5.559E-01 5.616E-01 5.673E-01 5.731E-01 5.788E-01 5.845E-01 5.902E-01 Luận văn tốt nghiệp 0.360 0.365 0.370 0.375 0.380 0.385 0.390 0.395 0.400 0.405 0.410 0.415 0.420 0.425 0.430 0.435 0.440 0.445 0.450 0.455 0.460 0.465 0.470 0.475 0.480 0.485 0.490 0.495 0.500 0.505 0.510 0.515 0.520 0.525 0.530 0.535 0.540 0.545 0.550 0.555 4.609E-02 4.771E-02 4.935E-02 5.101E-02 5.270E-02 5.441E-02 5.614E-02 5.789E-02 5.965E-02 6.144E-02 6.324E-02 6.506E-02 6.690E-02 6.875E-02 7.061E-02 7.249E-02 7.438E-02 7.628E-02 7.819E-02 8.012E-02 8.205E-02 8.398E-02 8.593E-02 8.789E-02 8.985E-02 9.182E-02 9.379E-02 9.577E-02 9.775E-02 9.974E-02 1.017E-01 1.037E-01 1.057E-01 1.077E-01 1.097E-01 1.117E-01 1.137E-01 1.157E-01 1.177E-01 1.197E-01 1.400 1.410 1.420 1.430 1.440 1.450 1.460 1.470 1.480 1.490 1.500 1.510 1.520 1.530 1.540 1.550 1.560 1.570 1.580 1.590 1.600 1.610 1.620 1.630 1.640 1.650 1.660 1.670 1.680 1.690 1.700 1.710 1.720 1.730 1.740 1.750 1.760 1.770 1.780 1.790 3.531E-01 3.549E-01 3.564E-01 3.581E-01 3.597E-01 3.613E-01 3.628E-01 3.504E-01 3.660E-01 3.875E-01 3.691E-51 3.706E-01 3.721E-01 3.736E-01 3.751E-01 3.765E-01 3.780E-01 3.755E-01 3.809E-01 3.803E-01 3.838E-01 3.852E-01 3.868E-01 3.880E-01 3.894E-01 3.907E-01 3.921E-01 3.935E-01 3.948E-01 3.962E-01 3.975E-01 3.988E-01 4.002E-01 4.015E-01 4.028E-01 4.041E-01 4.054E-01 4.067E-01 4.080E-01 4.083E-01 60 Nguyễn Thị Minh 3.400 3.450 3.500 3.550 3.600 3.650 3.700 3.750 3.800 3.850 3.900 3.950 4.000 4.050 4.100 4.150 4.200 4.250 4.300 4.350 4.400 4.450 4.500 4.550 4.600 4.650 4.700 4.750 4.800 4.850 4.900 4.950 5.000 5.050 5.100 5.150 5.200 5.250 5.300 5.350 5.950E-01 6.017E-01 6.074E-01 6.131E-01 6.187E-01 6.249E-01 6.300E-01 6.357E-01 6.413E-01 6.469E-01 6.525E-01 6.580E-01 6.635E-01 6.690E-01 6.745E-01 6.800E-01 6.854E-01 6.908E-01 6.962E-01 7.016E-01 7.069E-01 7.122E-01 7.174E-01 7.226E-01 7.278E-01 7.330E-01 7.381E-01 7.432E-01 7.483E-01 7.533E-01 7.583E-01 7.633E-01 7.682E-01 7.731E-01 7.780E-01 7.828E-01 7.876E-01 7.924E-01 7.972E-01 8.019E-01 Luận văn tốt nghiệp 0.560 0.565 0.570 0.575 0.580 0.585 0.590 0.595 0.600 0.605 0.610 0.615 0.620 0.625 0.630 0.635 0.640 0.645 0.650 0.655 0.660 0.665 0.670 0.675 0.680 0.685 0.690 0.695 0.700 0.705 0.710 0.715 0.720 0.725 0.730 0.735 0.740 0.745 0.750 0.755 1.217E-01 1.237E-01 1.257E-01 1.277E-01 1.296E-01 1.316E-01 1.336E-01 1.356E-01 1.375E-01 1.395E-01 1.415E-01 1.434E-01 1.454E-01 1.473E-01 1.493E-01 1.512E-01 1.531E-01 1.550E-01 1.570E-01 1.589E-01 1.608E-01 1.627E-01 1.646E-01 1.664E-01 1.683E-01 1.702E-01 1.720E-01 1.739E-01 1.757E-01 1.776E-01 1.794E-01 1.812E-21 1.830E-01 1.848E-01 1.866E-01 1.884E-01 1.902E-01 1.919E-01 1.937E-01 1.954E-01 1.800 1.810 1.820 1.830 1.840 1.850 1.860 1.870 1.880 1.890 1.900 1.910 1.920 1.930 1.940 1.950 1.960 1.970 1.980 1.990 2.000 2.010 2.020 2.030 2.040 2.050 2.060 2.070 2.080 2.090 2.100 2.110 2.120 2.130 2.140 2.150 2.160 2.170 2.180 2.190 4.108E-01 4.118E-01 4.131E-01 4.143E-01 4.156E-01 4.168E-01 4.181E-01 4.193E-01 4.206E-01 4.218E-01 4.230E-01 4.243E-01 4.255E-01 4.267E-01 4.279E-01 4.291E-01 4.303E-01 4.315E-01 4.327E-01 4.339E-01 4.351E-01 4.363E-01 4.375E-01 4.387E-01 4.399E-01 4.411E-01 4.423E-01 4.434E-01 4.446E-01 4.458E-01 4.470E-01 4.481E-01 4.493E-01 4.505E-01 4.516E-01 4.522E-01 4.540E-01 4.551E-01 4.563E-01 4.574E-01 61 Nguyễn Thị Minh 5.400 5.450 5.500 5.550 5.600 5.650 5.700 5.750 5.800 5.850 5.900 5.950 6.000 6.050 6.100 6.150 6.200 6.250 6.300 6.350 6.400 6.450 6.500 6.550 6.600 6.650 6.700 6.750 6.800 6.850 6.900 6.950 7.000 7.100 7.200 7.300 7.400 7.500 7.600 7.700 8.065E-01 8.112E-01 8.158E-01 8.204E-01 8.249E-01 8.294E-01 8.339E-01 8.384E-01 8.428E-01 8.473E-01 8.516E-01 8.560E-01 8.603E-01 8.646E-01 8.689E-01 8.732E-01 8.779E-01 8.816E-01 8.858E-01 8.907E-01 8.942E-01 8.983E-01 9.024E-01 9.065E-01 9.106E-01 9.147E-01 9.187E-01 9.228E-01 9.263E-01 9.308E-01 9.348E-01 9.388E-01 9.428E-01 9.508E-01 9.587E-01 9.666E-01 9.745E-01 9.824E-01 9.903E-01 9.982E-01 Luận văn tốt nghiệp 0.760 0.765 0.770 0.775 0.780 0.785 0.790 0.795 0.800 0.805 0.810 0.815 0.820 0.825 0.830 0.835 0.840 0.845 0.850 0.855 0.860 0.865 0.870 0.875 0.880 0.885 1.972E-01 1.989E-01 2.006E-01 2.023E-01 2.040E-01 2.057E-01 2.074E-01 2.091E-01 2.107E-01 2.124E-01 2.140E-01 2.157E-01 2.173E-01 2.189E-01 2.205E-01 2.221E-01 2.237E-01 2.253E-01 2.269E-01 2.285E-01 2.300E-01 2.316E-01 2.331E-01 2.346E-01 2.361E-01 2.376E-01 2.200 2.210 2.220 2.230 2.240 2.250 2.260 2.270 2.280 2.290 2.300 2.310 2.320 2.330 2.340 2.350 2.360 2.370 2.380 2.390 2.400 2.410 2.420 2.430 2.440 2.450 4.586E-01 4.597E-01 4.609E-01 4.621E-01 4.632E-01 4.644E-01 4.655E-01 4.667E-01 4.678E-01 4.690E-01 4.701E-01 4.713E-01 4.724E-01 4.735E-01 4.747E-01 4.758E-01 4.770E-01 4.781E-01 4.793E-01 4.804E-01 4.816E-01 4.827E-01 4.838E-01 4.850E-01 4.861E-01 4.873E-01 62 Nguyễn Thị Minh 7.800 7.900 8.000 8.100 8.200 8.300 8.400 8.500 8.600 8.700 8.800 8.900 9.000 9.100 9.200 9.300 9.400 9.500 9.600 9.700 9.800 9.900 10.000 10.100 10.200 10.300 1.006E+00 1.014E+00 1.022E+00 1.030E+00 1.038E+00 1.046E+00 1.054E+00 1.083E+00 1.071E+00 1.079E+00 1.088E+00 1.096E+00 1.105E+00 1.114E+00 1.122E+00 1.131E+00 1.140E+00 1.150E+00 1.159E+00 1.168E+00 1.178E+00 1.181E+00 1.197E+00 1.207E+00 1.217E+00 1.227E+00 ... đặt cho nhà nghiên cứu làm để nâng cao độ xác ổn định cho thiết bị đo đạc điều kiện môi trường khác Luận văn thực với mục tiêu nâng cao độ xác, tính ổn định thiết bị quan trắc cảnh báo phóng xạ. .. Độ rộng đỉnh hấp thụ toàn phần ứng với độ rộng nửa chiều cao đỉnh (FWHM) gọi độ phân giải đầu dò 1.3 Thiết bị quan trắc cảnh báo phóng xạ mơi trường Thiết bị quan trắc cảnh báo phóng xạ mơi trường. .. Tuy nhiên, tốn nâng cao độ xác ổn định cho thiết bị vấn đề nhà khoa học giới quan tâm 1.2 Cơ sở lý thuyết Để hiểu rõ nguyên tắc hoạt động thiết bị quan trắc cảnh báo phóng xạ môi trường, ta cần