Tổng quan nghiên cứu

Việt Nam nằm trọn trong tâm dông châu Á, là một trong những khu vực có hoạt động dông sét mạnh mẽ và phức tạp nhất thế giới. Theo thống kê khí tượng toàn cầu, mỗi năm Trái Đất ghi nhận khoảng 3 tỷ cú sét đánh, cướp đi sinh mạng của khoảng 5000 người và gây thiệt hại hàng nghìn tỷ đồng cho các ngành điện lực, hàng không, viễn thông và dầu khí. Tại Việt Nam, nhiều khu vực như Hải Dương, Đông Anh (Hà Nội) hay khu vực Đồng bằng sông Cửu Long thường xuyên chịu tổn thất nghiêm trọng về người và thiết bị viễn thông nhạy cảm do sét đánh trực tiếp và cảm ứng điện từ.

Vấn đề nghiên cứu cốt lõi đặt ra là sự phân bố và biến thiên của cấu trúc điện tích trong mây dông diễn biến rất nhanh theo không gian và thời gian, gây khó khăn lớn cho công tác dự báo. Mục tiêu cụ thể của luận văn là xây dựng mô hình tính toán, mô phỏng cấu trúc điện trường tĩnh trong mây dông dựa trên dữ liệu quan trắc bề mặt, từ đó ứng dụng thử nghiệm đánh giá định lượng cho các cơn dông tại khu vực Hà Nội.

Phạm vi nghiên cứu tập trung vào không gian tầng đối lưu có độ cao từ 0 đến 10 km, bán kính quét 40 km quanh Hà Nội, với ca nghiên cứu điển hình là cơn dông ngày 22 tháng 3 năm 2022 kéo dài 86 phút (từ 18 giờ 35 phút đến 20 giờ 00 phút). Luận văn mang ý nghĩa khoa học và thực tiễn sâu sắc khi cung cấp thuật toán số hóa giúp xác định vị trí đám mây tích điện với chỉ số trùng lặp không gian đạt tới 80,88% so với ảnh radar thời tiết, mở ra khả năng cảnh báo sớm dông sét trước từ 30 đến 40 phút.

Cơ sở lý thuyết và phương pháp nghiên cứu

Khung lý thuyết áp dụng

Nghiên cứu ứng dụng nền tảng của Điện động lực học cổ điển kết hợp Định luật Cu-lông (Coulomb) và phương pháp điện tích ảnh để giải bài toán phân bố điện trường tĩnh do các khối điện tích trong mây dông gây ra trên bề mặt đất phẳng dẫn điện. Theo mô hình mạch điện toàn cầu, khí quyển hoạt động như một tụ điện khổng lồ với lớp điện tích mặt đất mang giá trị khoảng -5x10^5 C và điện trường thời tiết tốt tại mặt đất duy trì mức ổn định khoảng -100 V/m đến -150 V/m.

Cấu trúc điện tích mây dông được luận văn mô hình hóa qua hai cấu hình chính: mô hình lưỡng cực thẳng đứng và mô hình tam cực điện tích (được phát triển từ các nghiên cứu kinh điển của Simpson, Scrase và Krehbiel). Trong cấu trúc tam cực, không gian mây dông gồm ba khối điện tích:

  1. Khối điện tích dương trên cao đặt tại độ cao 6 km với độ lớn +40 C, bán kính phân bố 5 km (nơi nhiệt độ dao động từ -10°C đến -25°C).
  2. Khối điện tích âm chính đặt tại độ cao 4 km với độ lớn -40 C, bán kính 5 km.
  3. Khối điện tích dương nhỏ vùng chân mây đặt tại độ cao 300 m với độ lớn +2 C, bán kính 1 km, hình thành do các dòng phóng điện điểm corona từ các vật thể nhọn trên mặt đất.

Các khái niệm then chốt được phân tích bao gồm: điện trường phương thẳng đứng (Ez), mật độ điện tích không gian, dòng điện dịch chuyển Maxwell và dòng phóng điện điểm corona với mật độ dòng trung bình từ 50 đến 100 nA/m².

Phương pháp nghiên cứu

Nguồn dữ liệu thực nghiệm được thu thập đồng bộ từ mạng lưới quan trắc gồm:

  • Ba trạm đo điện trường tĩnh khí quyển sử dụng thiết bị cảm biến xoay EFM-100 đặt tại Phú Thụy (Gia Lâm), Phú Xuyên và Chương Mỹ thuộc Hà Nội.
  • Hệ thống radar thời tiết Doppler băng X (bước sóng từ 2,5 cm đến 4 cm, tần số hoạt động 8 đến 12 GHz, bán kính quan sát hiệu dụng 60 km) đặt tại Học viện Nông nghiệp Việt Nam.

Cỡ mẫu nghiên cứu bao gồm chuỗi số liệu đo liên tục từng giây trong suốt 86 phút của cơn dông ngày 22 tháng 3 năm 2022, kết hợp tập dữ liệu đối soát 521 ngày quan trắc khí quyển tại trạm Phú Thụy. Phương pháp chọn mẫu tập trung vào các đợt dông đối lưu nhiệt điển hình có sự nhiễu loạn mạnh của điện trường mặt đất vượt ngưỡng 1000 V/m.

Lý do lựa chọn phương pháp phân tích số trên phần mềm MATLAB kết hợp thuật toán tối ưu hóa bình phương tối thiểu sai số là vì phương pháp này cho phép xử lý hiệu quả bài toán ngược phi tuyến: từ các giá trị điện trường Ez đo được tại 3 tọa độ mặt đất, thuật toán quét trên không gian lưới chia nhỏ 1x1 km với độ phân giải ô khối 100x100x100 m để tìm ra tọa độ không gian (x, y, z) và độ lớn điện tích (từ 0 C đến 80 C) có độ sai lệch nhỏ nhất giữa tính toán lý thuyết và thực đo.

Kết quả nghiên cứu và thảo luận

Những phát hiện chính

Quá trình mô phỏng và phân tích dữ liệu thực nghiệm đã mang lại 4 phát hiện quan trọng:

Thứ nhất, vị trí đám mây điện tích tính toán bằng mô hình quang học - tĩnh điện trùng khớp cao với ảnh phản hồi vô tuyến của radar thời tiết băng X. Chỉ số trùng lặp không gian giữa tọa độ đám mây tích điện tính toán và vùng mây gây mưa có cường độ giáng thủy trên 2 mm/h đạt mức 80,88%, và đạt tới 89,88% đối với vùng mây gây mưa trên 1 mm/h.

Thứ hai, biên độ biến thiên điện trường phụ thuộc chặt chẽ vào khoảng cách từ tâm dông đến trạm quan trắc. Tại trạm Phú Thụy (nằm ngay dưới vùng mây dông hoạt động), cường độ điện trường Ez dao động dữ dội từ -8 kV/m đến +6 kV/m. Trong khi đó, tại trạm Phú Xuyên (cách rìa cơn dông khoảng 25 km) và Chương Mỹ (cách rìa cơn dông khoảng 20 km), giá trị Ez ghi nhận được lần lượt từ -2 kV/m đến +2 kV/m và từ -4 kV/m đến +2 kV/m, cao hơn từ 15 đến 40 lần so với điện trường thời tiết đẹp.

Thứ ba, sự tiến triển của tổng lượng điện tích mây dông chia thành hai pha rõ rệt trong suốt chu kỳ 86 phút:

  • Pha 1 (35 phút đầu): Dòng thăng phát triển với vận tốc lên tới 30 m/s đưa các hạt nước chậm đông lên độ cao 10 km, làm tăng nhanh ma sát và tích lũy điện tích, điện trường biến thiên theo mô hình lưỡng cực dương.
  • Pha 2 (sau 35 phút đến khi kết thúc): Các hạt giáng thủy rơi xuống kéo theo dòng khí giáng, phóng điện corona mặt đất hình thành lớp ion dương dày đặc dưới chân mây, chuyển dịch điện trường sang cấu trúc tam cực cho đến khi điện tích tiêu tán hoàn toàn.

Thứ tư, khối tâm điện tích có xu hướng dịch chuyển rõ ràng từ phía Đông Hà Nội ở khoảng cách hoành độ x từ 30 đến 40 km lùi dần về khu vực hoành độ x từ 15 đến 25 km trước khi tan rã.

Thảo luận kết quả

Nguyên nhân của sự chuyển pha cấu trúc điện trường từ lưỡng cực sang tam cực bắt nguồn từ sự hình thành lớp điện tích đảo cực dưới chân mây do hiệu ứng corona từ thảm thực vật và các công trình cao tầng khi điện trường mây dông vượt ngưỡng đánh thủng cục bộ.

Khi so sánh với các nghiên cứu của Zeng (năm 2013) tại khu vực Đông Á hay các công trình quan trắc của Nguyễn Xuân Anh và cộng sự, kết quả của luận văn hoàn toàn tương thích về mặt vật lý. Thời gian điện trường bề mặt biến đổi đột biến trên 1000 V/m xuất hiện trước cú sét đánh xuống đất đầu tiên từ 30 đến 40 phút, tương đồng với tỷ lệ cảnh báo đúng đạt 86,99% của phương pháp cảnh báo hai vùng (vùng trung tâm 10 km và vùng mở rộng 30 km).

Về mặt trình bày dữ liệu, các kết quả tính toán có thể được trực quan hóa hiệu quả thông qua:

  • Biểu đồ đường thời gian biểu diễn biến thiên điện trường Ez thực đo và mô phỏng tại 3 trạm trắc nghiệm.
  • Bản đồ mặt bằng không gian 2D và 3D thể hiện các điểm điện tích Cu-lông phân tầng theo thang màu, lồng ghép trực tiếp trên bản đồ phản hồi vô tuyến radar phản ánh cường độ mưa (từ 0 mm/h đến trên 2 mm/h).
  • Bảng thống kê đối chiếu hệ số tương quan sai số và chỉ số trùng lặp không gian qua từng giai đoạn phát triển của cơn dông.

Đề xuất và khuyến nghị

Dựa trên các kết luận khoa học đạt được, luận văn đưa ra 4 khuyến nghị và giải pháp hành động cụ thể:

  1. Mở rộng mạng lưới trạm đo điện trường mặt đất EFM-100: Bổ sung từ 5 đến 8 trạm quan trắc điện trường tĩnh phân bố đều trong bán kính 50 km quanh các trung tâm đô thị lớn và khu công nghiệp trọng điểm. Mục tiêu là nâng cao độ chính xác định vị tâm điện tích thêm 15% đến 20%, triển khai trong giai đoạn 2024-2025 do Viện Vật lý Địa cầu chủ trì.
  2. Tích hợp thuật toán nghịch đảo điện tích vào hệ thống radar và vệ tinh viễn thám: Xây dựng module phần mềm tự động tính toán vị trí khối điện tích theo thời gian thực kết hợp ảnh mây vệ tinh Himawari và radar băng X. Mục tiêu rút ngắn độ trễ tính toán xuống dưới 15 giây/chu kỳ quét, nâng thời gian cảnh báo sét sớm lên 35 đến 45 phút với tỷ lệ cảnh báo khống dưới 18%, thực hiện trong lộ trình 2024-2026 bởi Tổng cục Khí tượng Thủy văn.
  3. Chuẩn hóa quy trình vận hành an toàn sét hai vùng cho hạ tầng trọng yếu: Ban hành quy chuẩn vận hành ngắt điện tự động và bảo vệ quá áp cho các trạm biến áp 220 kV - 500 kV, sân bay quốc tế và tháp truyền thông dựa trên ngưỡng cảnh báo điện trường Ez vượt mức 1,5 kV/m ở vùng cảnh báo 10 km. Mục tiêu giảm thiểu 90% nguy cơ cháy nổ thiết bị viễn thông, thực hiện trong vòng 12 tháng tới do Tập đoàn Điện lực Việt Nam và Cục Hàng không phối hợp ban hành.
  4. Ứng dụng thuật toán học máy để lọc nhiễu phóng điện cục bộ: Huấn luyện các mô hình trí tuệ nhân tạo trên tập dữ liệu hơn 500 ngày quan trắc để tự động hiệu chỉnh độ méo điện trường do địa hình và thảm thực vật gây ra, nâng độ tin cậy của thuật toán lên trên 92%, giao cho các nhóm nghiên cứu viện - trường thực hiện từ năm 2025.

Đối tượng nên tham khảo luận văn

Luận văn là tài liệu tham khảo chuyên sâu và hữu ích cho 4 nhóm đối tượng trọng tâm:

  • Nhóm nghiên cứu và học viên cao học chuyên ngành Vật lý lý thuyết, Vật lý Khí quyển và Địa vật lý: Nắm bắt phương pháp giải bài toán ngược tĩnh điện, kỹ thuật xử lý dữ liệu cảm biến xoay EFM-100 và mô hình hóa cấu trúc tam cực điện tích trong khí quyển nhiệt đới.
  • Kỹ sư thiết kế và vận hành hệ thống lưới điện truyền tải: Sử dụng các số liệu thực nghiệm về biên độ biến thiên điện trường (từ -8 kV/m đến +6 kV/m) để tính toán dung lượng chống sét van, thiết kế hệ thống tiếp địa và cài đặt rơ-le bảo vệ cho các trạm biến áp 110 kV đến 500 kV.
  • Chuyên gia điều hành bay và quản lý an toàn công trình hàng không, bưu chính viễn thông: Áp dụng thuật toán cảnh báo sớm hai vùng để điều phối lịch cất hạ cánh trong điều kiện thời tiết nguy hiểm và bảo vệ an toàn cho hệ thống radar hàng không, trung tâm dữ liệu.
  • Cán bộ quản lý rủi ro thiên tai tại các địa phương: Xây dựng bản đồ phân vùng nguy cơ dông sét cục bộ, ban hành khuyến cáo an toàn lao động ngoài trời cho người dân tại các điểm nóng dông sét như Hải Dương hay khu vực Đồng bằng sông Cửu Long.

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

  1. Mô hình cấu trúc tam cực giải thích hiện tượng điện trường mây dông ưu việt hơn mô hình lưỡng cực như thế nào? Mô hình tam cực bổ sung khối điện tích dương nhỏ (khoảng +2 C) ở độ cao 300 m gần mặt đất bên cạnh hai khối tích điện chính (+40 C ở 6 km và -40 C ở 4 km). Cấu trúc này giải thích chính xác hiện tượng điện trường mặt đất đảo dấu từ dương sang âm tại trạm Phú Thụy sau 35 phút hoạt động của cơn dông.

  2. Cảm biến điện trường tĩnh EFM-100 hoạt động theo nguyên lý đo lường nào? Thiết bị sử dụng cánh quay rotor kim loại quay liên tục để che chắn và phơi bày luân phiên các bản cực stator tĩnh trước điện trường khí quyển. Dòng điện cảm ứng xoay chiều sinh ra trên stator tỷ lệ thuận với cường độ điện trường ngoài, được bộ khuếch đại thuật toán xử lý để đo chính xác điện trường từ -8 kV/m đến +6 kV/m.

  3. Dữ liệu radar thời tiết băng X đóng vai trò gì trong việc kiểm chứng mô hình tính toán? Radar băng X sử dụng bước sóng ngắn từ 2,5 đến 4 cm, tần số 8-12 GHz, cực kỳ nhạy trong việc phát hiện các hạt nước nhỏ và mưa phùn trong bán kính 60 km. Vị trí đám mây điện tích tính từ mô hình toán cho thấy sự trùng khớp không gian tới 80,88% với các tâm phản hồi mưa trên 2 mm/h của radar.

  4. Sự biến thiên của điện trường mặt đất giúp cảnh báo sớm sét trước bao nhiêu phút? Trong điều kiện thực tế, khi mây dông đối lưu hình thành, cường độ điện trường bề mặt bắt đầu biến đổi mạnh và vượt ngưỡng 1000 V/m đến 1500 V/m trước khi xuất hiện cú phóng điện xuống đất đầu tiên từ 30 đến 40 phút, tạo khoảng thời gian vàng để kích hoạt các biện pháp phòng ngừa.

  5. Tại sao các trạm đo ở xa cơn dông như Phú Xuyên hay Chương Mỹ vẫn ghi nhận điện trường bất thường? Dù nằm cách rìa tâm dông từ 20 đến 25 km, các trạm này vẫn ghi nhận điện trường dao động từ -4 kV/m đến +2 kV/m (vượt xa mức thời tiết đẹp -100 V/m) do đường sức điện trường từ các khối tích điện quy mô lớn trên độ cao 6 km lan tỏa rộng trong không gian và tương tác với các đám mây đối lưu lân cận.

Kết luận

  • Xây dựng thành công thuật toán số hóa mô phỏng phân bố không gian điện tích 3D trong mây dông dựa trên nền tảng Định luật Cu-lông và phương pháp điện tích ảnh.
  • Xác thực độ chính xác của mô hình qua cơn dông ngày 22 tháng 3 năm 2022 tại Hà Nội, đạt tỷ lệ trùng khớp không gian 80,88% so với dữ liệu radar thời tiết băng X.
  • Định lượng chi tiết sự phân hóa điện trường cực đại tại trạm tâm dông Phú Thụy (-8 kV/m đến +6 kV/m) và các trạm ngoại vi (-4 kV/m đến +2 kV/m), làm rõ cơ chế chuyển pha lưỡng cực - tam cực.
  • Thiết lập cơ sở khoa học vững chắc phục vụ bài toán cảnh báo sớm sét đánh trước từ 30 đến 40 phút với lộ trình tích hợp hệ thống số liệu đa nguồn giai đoạn 2024-2026.
  • Đề nghị các cơ quan khí tượng, ngành điện lực và các đơn vị quản lý công trình khẩn trương phối hợp mở rộng mạng lưới trạm đo EFM-100 để hiện đại hóa công tác phòng chống thiên tai sét tại Việt Nam.