Bối cảnh và vấn đề nghiên cứu

Việt Nam nằm trong khu vực tâm sét châu Á, đặc trưng bởi tần suất và mật độ hoạt động dông sét trung bình hàng năm ở mức cao trên thế giới. Sự phóng điện trong khí quyển không chỉ đe dọa trực tiếp đến tính mạng con người và kết cấu công trình xây dựng mà còn tạo ra các xung quá điện áp lan truyền phá hủy các trang thiết bị kỹ thuật. Trong tiến trình công nghiệp hóa và hiện đại hóa, mức độ thiệt hại do sét gây ra có xu hướng gia tăng rõ rệt khi mật độ sử dụng các thiết bị điện tử, viễn thông và hệ thống điều khiển tự động nhạy cảm ngày càng dày đặc. Việc nghiên cứu, đánh giá và xây dựng giải pháp bảo vệ chống sét đồng bộ cho các công trình công cộng, trường học là một yêu cầu kỹ thuật cấp thiết.

Đề tài "Nghiên cứu giải pháp chống sét hợp lý cho công trình xây dựng khu vực Cơ sở II – ĐHGTVT" tập trung giải quyết 03 nhiệm vụ và mục tiêu nghiên cứu cụ thể:

  1. Hệ thống hóa cơ sở lý thuyết về hiện tượng dông sét, cơ chế phóng điện và các nguyên lý chống sét trực tiếp, chống sét lan truyền cho công trình xây dựng.
  2. Tiến hành khảo sát, đo đạc thông số và đánh giá thực trạng kỹ thuật của hệ thống chống sét hiện hữu tại các khu nhà thuộc Trường Đại học Giao thông Vận tải Cơ sở II.
  3. Ứng dụng công cụ tính toán mô phỏng và đề xuất phương án thiết kế hệ thống chống sét kỹ thuật hợp lý, an toàn cho khu giảng đường của nhà trường.

Đối tượng nghiên cứu của đề tài là các quy luật vật lý của dòng điện sét, phương pháp tính toán bán kính bảo vệ, điện trở tiếp đất và các giải pháp công nghệ chống sét công trình. Phạm vi nghiên cứu về mặt không gian tập trung vào các khối công trình thuộc khuôn viên Trường Đại học Giao thông Vận tải Cơ sở II (bao gồm khu nhà C2, Thư viện, E7, E10, D1 và khu giảng đường).

Cơ sở lý thuyết và phương pháp

Khung lý thuyết và tiêu chuẩn áp dụng

Đề tài xây dựng các mô hình tính toán dựa trên hệ thống quy chuẩn, tiêu chuẩn kỹ thuật chuyên ngành trong nước và quốc tế:

  • TCXD 46:1984: Chống sét cho công trình xây dựng – Tiêu chuẩn thiết kế thi công.
  • TCXD 46:2007: Chống sét cho công trình xây dựng – Hướng dẫn thiết kế, kiểm tra và bảo trì hệ thống.
  • NFC 17-102 (Pháp): Tiêu chuẩn thiết kế, đánh giá cấp bảo vệ và kiểm định hệ thống chống sét phát tia tiên đạo sớm (ESE - Early Streamer Emission).
  • TCN 68-141:1995: Tiêu chuẩn ngành về tiếp đất cho các công trình viễn thông, quy định yêu cầu kỹ thuật và phương pháp tính toán điện trở nối đất.
  • NFPA 780: Tiêu chuẩn lắp đặt hệ thống chống sét của Hoa Kỳ (tham chiếu phương pháp chống sét cổ điển Franklin).

Phương pháp nghiên cứu và thu thập dữ liệu

  • Phương pháp phân tích - tổng hợp lý thuyết: Khái quát hóa các mô hình thu sét cổ điển (Franklin), lồng Faraday, công nghệ phân tán điện tích, công nghệ phát tia tiên đạo sớm và hệ thống phòng chống tích cực qua thiết bị cảnh báo sớm (EFM, ESID).
  • Phương pháp khảo sát và đo đạc thực nghiệm: Khảo sát hiện trạng kết cấu kim thu sét, tuyến dây dẫn sét và thực hiện đo đạc giá trị điện trở tiếp đất $R_đ$ bằng thiết bị chuyên dụng tại các điểm đo của từng khối nhà.
  • Phương pháp đo điện trở suất của đất: Ứng dụng phương pháp 4 cực đối xứng Wenner ($A, B, M, N$) với công thức tính điện trở suất $\rho = 2\pi a R$, hỗ trợ khảo sát địa chất công trình ở độ sâu từ 15 m đến 20 m phục vụ thiết kế bãi tiếp địa.
  • Phương pháp tính toán và lập trình tự động hóa: Sử dụng phần mềm MATLAB xây dựng giao diện tính toán tham số bảo vệ, bán kính quét và mô hình hóa hệ thống nối đất có hoặc không có hóa chất giảm điện trở đất GEM (Grounding Enhancement Materials).

Thiết kế và triển khai công cụ tính toán mô phỏng

Đề tài đã thiết kế công cụ phần mềm tính toán chống sét trên nền tảng MATLAB 6 thông qua môi trường lập trình giao diện người dùng GUIDE (Graphical User Interface Development Environment).

Kiến trúc và chức năng của phần mềm

Phần mềm bao gồm các mô-đun chức năng chính:

  • Mô-đun xác định cấp bảo vệ: Tiếp nhận các thông số hình học công trình (chiều dài $L$, chiều rộng $W$, chiều cao $H$), mật độ sét phóng điện $N_{gmax}$, và các hệ số vị trí, môi trường ($C_1, C_2, C_3, C_4, C_5$). Chương trình tự động tính diện tích thu sét tương đương $A_c$, tần số sét đánh trực tiếp $N_d$, tần số chấp nhận $N_c$ và hiệu quả bảo vệ $E$ để đưa ra cấp bảo vệ (Cấp I, II, III).

  • Mô-đun tính toán kim thu sét phát tia tiên đạo: Xác định bán kính bảo vệ $R_p$ theo độ cao kim $h$, độ lợi thời gian $\Delta T$ và khoảng cách tiên đạo $\Delta L$ dựa trên công thức tiêu chuẩn NFC 17-102:

    $$R_p = \sqrt{h(2D - h) + \Delta L(2D + \Delta L)} \quad (với\ h \ge 5m)$$

  • Mô-đun tính toán điện trở bãi tiếp đất: Lập trình thuật toán xác định điện trở của cọc tiếp đất đơn, thanh dẫn nằm ngang, vòng xuyến, tấm tròn, hệ thống nhiều tia hình sao và hệ nhiều cọc nối song song có xét đến hệ số sử dụng $\eta_1, \eta_2$.

  • Mô-đun tính toán lượng hóa chất GEM: Tính toán số bao hóa chất GEM cần thiết ($n = V / 0,014$) dựa trên thể tích rãnh/hố tiếp địa nhằm hạ điện trở suất của các vùng đất có tính dẫn điện kém (cát, sỏi, đá).

Giao diện và cơ chế xử lý dữ liệu

Giao diện được xây dựng từ các thành phần điều khiển giao tiếp: Push Button, Edit Text, Static Text, và Axes. Cơ chế kích hoạt tính toán được lập trình thông qua hàm gọi lại Callback (như buttonCalculate), kết hợp các lệnh truy xuất dữ liệu get, set và hàm chuyển đổi kiểu dữ liệu str2num, num2str.

Nội dung chính theo từng chương

Chương 1: Lý thuyết chống sét cho công trình xây dựng

Chương này tổng kết các nguyên lý khoa học và thực tiễn về hiện tượng dông sét:

  • Bản chất vật lý: Sét là sự phóng điện trong khí quyển giữa các đám mây dông tích điện hoặc giữa mây và đất. Vận tốc dòng ion đạt khoảng 36.000 dặm/giờ, nhiệt độ dòng plasma lên tới 30.000 °C (gấp nhiều lần nhiệt độ làm nóng chảy cát silica để tạo thành đá sét fulgurite). Hàng năm có khoảng 16 triệu cơn dông trên toàn cầu, trung bình khoảng 100 cú sét đánh xuống đất mỗi giây.
  • Phân loại tác động: Sét đánh trực tiếp (nguy cơ phá hủy cơ học và hỏa hoạn nghiêm trọng nhất); sét đánh gián tiếp (xung điện truyền qua mạng lưới dây dẫn điện, dây thông tin); sét cảm ứng (bao gồm cảm ứng tĩnh điện và cảm ứng điện từ gây quá điện áp xung phá hỏng linh kiện điện tử).
  • Lợi ích tự nhiên: Cung cấp nguồn nhiệt ban đầu trong lịch sử văn minh; tạo khí quyển ozon ($2O_2 \rightarrow O_3 + [O]$); cố định nitơ thành hợp chất $NO$ giúp cải tạo độ phì của đất; hỗ trợ nghiên cứu địa chất, khoáng sản và mạch nước ngầm.
  • Hệ thống hóa các phương pháp bảo vệ:
    1. Phương pháp cổ điển (Kim Franklin): Tạo điểm đón sét trên cao và dẫn xuống đất, góc bảo vệ quy chuẩn từ 45° đến 60° theo TCXD 46:2007.
    2. Phương pháp kim phát tia tiên đạo sớm (ESE): Tích lũy năng lượng điện trường để phát tia tiên đạo đi lên sớm hơn, mở rộng vùng thu sét (Collection Volume).
    3. Phương pháp phân tán điện tích: Sử dụng các điện cực đa điểm nhọn để ion hóa không khí, giảm mật độ tích tụ điện trường trên mái công trình.
    4. Phương pháp lồng Faraday: Bao bọc công trình bằng lưới kim loại liên tục nhằm triệt tiêu điện trường bên trong (thường dùng cho kho vũ khí, công trình hạt nhân, trung tâm dữ liệu).
    5. Phương pháp hút sét bằng laser: Định hướng dòng phóng điện bằng chùm tia laser ion hóa (nghiên cứu thực nghiệm tại Nhật Bản năm 1997 thu được 2 lần sét đánh).
    6. Phương pháp phòng ngừa tích cực: Cảnh báo sớm dông sét qua cảm biến điện trường và định vị phóng điện (như hệ thống EFM, ESID tại trạm Phú Thụy, Gia Lâm).

Chương 2: Mô phỏng và tính toán chống sét bằng phần mềm MATLAB

Chương 2 trình bày toàn bộ cơ sở toán học, thuật toán và quy trình xây dựng công cụ tính toán tự động:

  • Cơ sở xác định cấp bảo vệ: Triển khai thuật toán theo NFC 17-102. Xác định diện tích thu sét $A_c = LW + 6H(L+W) + 9\pi H^2$, tần số sét đánh trực tiếp $N_d = N_{gmax} \cdot A_c \cdot C_1 \cdot 10^{-6}$, và tần số chấp nhận $N_c = 5,5 \cdot 10^{-3} / (C_2 C_3 C_4 C_5)$.

  • Thuật toán tính toán tiếp địa: Chuyển hóa các công thức giải tích của TCN 68-141:1995 và TCXD 46:1984 thành mã lệnh xử lý. Tính toán điện trở cọc đứng ($l \le 3$ m, $d \le 6$ cm), thanh nằm ngang ($l \le 12$ m, $h \ge 0,6$ m), và hệ cọc nối song song:

    $$R_{đ} = \frac{R_0}{n \cdot \eta_1}$$

  • Quy trình vận hành phần mềm: Hướng dẫn thao tác nhập số liệu hình học, chọn hệ số công trình, tra cứu mật độ sét theo địa bàn và xuất kết quả tính toán cấp bảo vệ, bán kính kim thu và cấu hình bãi cọc tiếp địa.

Chương 3: Hiện trạng hệ thống chống sét của các công trình ở ĐH GTVT Cơ sở II

Chương này công bố số liệu khảo sát thực tế và đánh giá an toàn hệ thống chống sét tại các khối nhà thuộc khuôn viên trường:

Tên công trình Kích thước ($L \times W \times H$) Thiết bị thu sét hiện trạng Dây dẫn sét Điện trở tiếp đất đo đạc ($R_đ$) Đánh giá kỹ thuật và đề xuất
Khu nhà C2 50m $\times$ 10m $\times$ 20m Kim chủ động PREVECTRON S6 (S2.60) Cáp đồng trần trong ống nhựa, đi thẳng, liên tục $R_đ = 1,81\ \Omega$ (Đo 3 lần: 1,80; 1,82; 1,81) Hệ thống đạt chuẩn an toàn, kim thu đủ bán kính bao phủ toàn bộ công trình và khu lân cận.
Khu nhà Thư viện 25m $\times$ 25m $\times$ 13,5m Kim thu bị động (cổ điển) Thép xây dựng Ø8, hàn sơ sài, rỉ sét, uốn cong nhiều Cọc 1: $1,60\ \Omega$
Cọc 2: $2,46\ \Omega$
Cọc 3: $6,22\ \Omega$
Cọc 4: $3,73\ \Omega$		 Không an toàn; công trình nằm trọn trong vùng bảo vệ của kim C2. Kiến nghị tháo dỡ để tránh hiện tượng "gọi sét" hoặc phải thiết kế lại.
Khu nhà E7 55m $\times$ 13m $\times$ 10m Cột thu sét thụ động cao 3m trên mái Thép Ø8; khoảng cách 2 dây là 55m (>20m quy định); dây bị đứt Cọc 1: $10,8\ \Omega$ (không đạt)
Cọc 2: $8,80\ \Omega$		 Nguy hiểm do dây dẫn đứt làm kim mất nối đất (thành bẫy sét); cọc tiếp đất không đạt yêu cầu ($>10\ \Omega$). Kiến nghị gỡ bỏ và thiết kế lại.
Khu nhà E10 56m $\times$ 11m $\times$ 10m Sử dụng nhiều kim thu bị động Thép Ø8; khoảng cách 2 dây là 28m (>20m quy định); thiếu số lượng dây Cọc 1: $3,06\ \Omega$
Cọc 2: $5,24\ \Omega$
Cọc 3: $8,80\ \Omega$		 Bố trí quá nhiều kim gây lãng phí; số lượng dây dẫn xuống không đủ mật độ quy chuẩn. Cần thiết kế lại.
Khu nhà D1 22m $\times$ 11m $\times$ 8m Kim chủ động PDC-E15 Cáp đồng trần trong ống nhựa, đạt chuẩn $R_đ = 9,23\ \Omega$ (trung bình dưới $10\ \Omega$) Hệ thống hoạt động an toàn, đảm bảo yêu cầu kỹ thuật.
Khu nhà E4, E5, E6 Khối giảng đường Sử dụng 01 cột chống sét thụ động Thép dẫn sét thông thường Nối đất theo kết cấu hiện hữu Cần tính toán nâng cấp đồng bộ theo tiêu chuẩn hiện hành.

Chương 4: Tính toán chống sét cho khu giảng đường trường Đại học GTVT Cơ sở II

Chương cuối cùng trình bày phương án kỹ thuật thiết kế hoàn chỉnh cho khu giảng đường:

  • Lựa chọn phương án công nghệ: Đề xuất thay thế hệ thống kim cổ điển phân tán bằng giải pháp kim thu sét phát tia tiên đạo sớm công nghệ chủ động PREVECTRON kết hợp thiết bị đếm số lần sét đánh.
  • Tính toán vùng phủ bảo vệ: Xác định cao độ lắp đặt kim thu trên đỉnh mái công trình để bán kính bảo vệ $R_p$ bao trùm toàn bộ các điểm nhọn, góc cạnh của công trình và các khối nhà lân cận.
  • Thiết kế hệ thống tiếp địa: Sử dụng cọc thép bọc đồng kết hợp thanh đồng liên kết dẹt chôn sâu dưới đất $\ge 0,6$ m; bổ sung giải pháp sử dụng hợp chất giảm điện trở đất GEM tại các vị trí nền đất có điện trở suất cao nhằm duy trì trị số tiếp đất tổng hợp luôn đạt $R_đ < 10\ \Omega$.
  • Quy trình thi công và kiểm tra: Hướng dẫn chi tiết kỹ thuật liên kết đẳng thế, lắp đặt hộp kiểm tra điện trở cách ly phục vụ công tác đo kiểm định kỳ hàng năm.

Kết quả và đóng góp

Kết quả nghiên cứu chính

  • Đánh giá toàn diện thực trạng tiếp địa và thu lôi: Số liệu đo đạc thực tế chỉ ra sự chênh lệch lớn về chất lượng giữa các khối nhà. Trong khi khu C2 ($R_đ = 1,81\ \Omega$) và D1 ($R_đ = 9,23\ \Omega$) đạt độ an toàn cao nhờ trang bị kim ESE và cáp đồng trần, thì khu Thư viện, E7 và E10 tồn tại nhiều vi phạm kỹ thuật nghiêm trọng (dây dẫn sét thép Ø8 bị rỉ sét, đứt liên kết tiếp đất, cọc số 1 khu E7 có $R_đ = 10,8\ \Omega$ vượt ngưỡng cho phép).
  • Phát hiện hiện tượng trùng lặp và xung đột vùng bảo vệ: Phân tích hình học cho thấy khu nhà Thư viện nằm hoàn toàn trong bán kính bảo vệ của đầu thu PREVECTRON tại khối nhà C2 lân cận. Việc duy trì hệ thống thu lôi thụ động bị hỏng tiếp địa tại Thư viện tạo ra mối nguy cơ "gọi sét" đánh thẳng vào công trình thay vì phân tán an toàn.
  • Tự động hóa tính toán kỹ thuật: Xây dựng thành công bộ công cụ mô phỏng trên MATLAB 6 cho phép tính nhanh cấp bảo vệ theo NFC 17-102 và mô phỏng cấu hình bãi cọc tiếp địa kèm lượng hóa chất GEM.

Giải pháp và kiến nghị của tác giả

  • Khẩn trương tháo gỡ hoặc thiết kế lại toàn bộ hệ thống kim thu lôi bị đứt dây dẫn tại khu nhà E7 và hệ thống thu sét hư hỏng tại khu Thư viện.
  • Bổ sung số lượng dây dẫn sét cho khu E7 và E10 nhằm đảm bảo khoảng cách giữa các dây xuống không vượt quá 20 m theo tiêu chuẩn TCXD 46:2007.
  • Quy chuẩn hóa việc sử dụng cáp đồng trần đi trong ống bảo vệ phi kim loại thay thế cho việc sử dụng thép xây dựng Ø8 dễ bị oxy hóa ăn mòn.
  • Thiết lập định kỳ kiểm tra, đo đạc điện trở tiếp địa không quá 12 tháng/lần, đặc biệt trước mùa mưa bão.

Hạn chế và hướng nghiên cứu tiếp

  • Hạn chế về dữ liệu khảo sát ngầm: Do các cọc tiếp đất tại một số khối công trình (C2, Thư viện, E7, E10, D1) đã bị phủ kín dưới lớp nền bê tông sân bãi sau quá trình cải tạo, nhóm tác giả không thể đào khảo sát trực tiếp tình trạng ăn mòn kim loại của cọc và mối hàn dưới lòng đất mà chỉ có thể thu thập dữ liệu gián tiếp qua hộp kiểm tra.
  • Hạn chế trong mô phỏng bãi cọc: Thuật toán mô phỏng tiếp địa giả định môi trường đất đồng nhất và bỏ qua điện trở tiếp xúc giữa vật dẫn với hóa chất GEM bao quanh.
  • Hướng nghiên cứu tiếp: Tích hợp thêm các mô hình tính toán chi tiết cho hệ thống chống sét lan truyền trên đường nguồn hạ thế và mạng truyền dẫn số liệu; mở rộng khảo sát điện trở suất đất theo mặt cắt địa tầng chuyên sâu tại nhiều vị trí khác nhau trong khuôn viên trường.

Giá trị tham khảo

Đồ án là tài liệu tham khảo hữu ích cho sinh viên, kỹ sư các chuyên ngành Kỹ thuật Xây dựng, Kỹ thuật Công trình, Điện công nghiệp và Hệ thống kỹ thuật công trình (MEP):

  • Cung cấp tài liệu tổng hợp rõ ràng về quy trình tính toán cấp bảo vệ theo tiêu chuẩn Pháp NFC 17-102 và các tiêu chuẩn Việt Nam liên quan (TCXD 46:2007, TCN 68-141:1995).
  • Cung cấp mô hình giao diện MATLAB mẫu giúp sinh viên ngành kỹ thuật không chuyên về chống sét có thể tự xây dựng công cụ hỗ trợ thiết kế bãi cọc tiếp địa và vùng phủ kim thu lôi.
  • Là căn cứ kỹ thuật phục vụ trực tiếp cho Ban Quản trị thiết bị và Cơ sở vật chất của Trường Đại học Giao thông Vận tải Cơ sở II trong việc cải tạo, nâng cấp mạng lưới chống sét.

Câu hỏi thường gặp

1. Tại sao đồ án kiến nghị có thể tháo gỡ hệ thống kim thu sét tại khu nhà Thư viện?
Dựa trên kết quả tính toán mô phỏng vùng bảo vệ kết hợp, bán kính quét của đầu thu sét chủ động PREVECTRON lắp đặt tại tòa nhà C2 đã bao phủ toàn bộ diện tích không gian của khu nhà Thư viện. Mặt khác, hệ thống thu lôi hiện hữu của Thư viện sử dụng dây thép Ø8 hàn sơ sài, rỉ sét và nhiều điểm gấp khúc. Nếu tiếp địa hoặc dây dẫn suy giảm chất lượng, các kim thu này sẽ trở thành điểm đón sét nhưng không thể tiêu tán dòng sét an toàn, gây nguy hiểm cho công trình.

2. Tiêu chuẩn nào được tác giả sử dụng để xác định cấp bảo vệ và tính toán bán kính kim ESE?
Tác giả áp dụng tiêu chuẩn chống sét NFC 17-102 của Pháp để tính toán tần số sét đánh trực tiếp $N_d$, tần số chấp nhận được $N_c$, hệ số hiệu quả $E = 1 - N_c / N_d$ và xác định bán kính bảo vệ $R_p$ của kim thu phát tia tiên đạo sớm theo các cấp bảo vệ I, II, III.

3. Phương pháp đo điện trở suất của đất trong đề tài được tiến hành như thế nào?
Đề tài sử dụng phương pháp 4 cực đối xứng Wenner. Bốn điện cực $A, B, M, N$ được cắm thẳng hàng trên mặt đất với khoảng cách đều nhau là $a$. Dòng điện $I$ được phát qua hai cực dòng $A, B$ và hiệu điện thế $V$ được ghi nhận giữa hai cực thế $M, N$. Điện trở suất biểu kiến được xác định theo biểu thức $\rho = 2\pi a (V/I) = 2\pi a R$.

4. Khảo sát thực tế tại khu nhà E7 phát hiện những lỗi kỹ thuật nguy hiểm nào?
Khảo sát chỉ ra khu nhà E7 có cột thu sét cao 3 m trên trần nhưng dây dẫn sét đi xuống đã bị đứt, khiến kim mất liên kết với hệ thống đất (tạo nguy cơ thu hút sét đánh vào công trình nhưng không thoát sét được). Ngoài ra, khoảng cách giữa 2 dây dẫn sét là 55 m (vượt xa giới hạn tối đa 20 m theo quy chuẩn), và cọc tiếp địa số 1 có điện trở đo được là $10,8\ \Omega$ (vượt ngưỡng cho phép $10\ \Omega$).

5. Hóa chất GEM được đưa vào chương trình mô phỏng với quy cách tính toán ra sao?
Hóa chất GEM (Grounding Enhancement Materials) được sử dụng để giảm điện trở tiếp địa cho các vùng đất sỏi đá, bùn lầy hoặc bê tông. Trong mô hình tính toán, mỗi bao hóa chất tạo ra thể tích phủ dẫn điện khoảng $0,014\ m^3$. Số lượng bao GEM cần cung cấp cho hệ thống rãnh cọc tiếp địa được xác định theo công thức $n = V / 0,014$, với $V$ là thể tích không gian xử lý tiếp địa.

Kết luận

Đồ án tốt nghiệp đã hệ thống hóa đầy đủ cơ sở lý thuyết chống sét hiện đại và xây dựng thành công phần mềm tính toán, mô phỏng chuyên dụng trên môi trường MATLAB GUIDE. Qua khảo sát thực tế tại Trường Đại học Giao thông Vận tải Cơ sở II, nghiên cứu đã chỉ rõ các khiếm khuyết kỹ thuật nghiêm trọng tại một số khối nhà và đưa ra giải pháp khắc phục cụ thể. Kết quả tính toán thiết kế cho khu giảng đường bằng kim thu tiên đạo sớm PREVECTRON và bãi tiếp địa đạt chuẩn đóng góp giải pháp an toàn, tối ưu chi phí cho công tác phòng chống thiên tai tại công trình.