Luận án tiến sĩ nghiên cứu đáp ứng quá độ của hệ thống nối đất với dòng điện sét

Tổng quan về luận án

Bảo vệ chống sét cho hệ thống lưới điện truyền tải là một trong những bài toán cốt lõi đảm bảo an ninh năng lượng quốc gia. Tại Việt Nam, theo số liệu thống kê từ Tập đoàn Điện lực Việt Nam (EVN) và Tổng công ty Truyền tải điện Quốc gia (EVNNPT), quy mô lưới điện truyền tải 220 kV và 500 kV tính đến cuối năm 2019 đã đạt 25.747 km (trong đó đường dây 500 kV chiếm 8.291 km và 220 kV chiếm 17.456 km) và dự kiến vượt 40.000 km vào năm 2030. Do đặc thù nằm trong vùng khí hậu nhiệt đới gió mùa có hoạt động dông sét mạnh với mật độ sét lên tới 13,7 lần/100 km²/năm, các sự cố do sét luôn chiếm trên 50% tổng số sự cố toàn lưới truyền tải, cá biệt năm 2017 tỷ lệ này lên tới 87%.

Khoảng trống nghiên cứu (research gap) trọng yếu nằm ở sự bất cập trong triết lý thiết kế truyền thống: Quy phạm trang bị điện Việt Nam (QCVN/QPTBĐ 2006) và các tiêu chuẩn tương đương của Nga (GOST) chủ yếu kiểm soát an toàn chống sét thông qua giá trị điện trở nối đất một chiều ($R_{DC}$). Thực tiễn vận hành ghi nhận nghịch lý kỹ thuật: "Không lý giải được khi điện trở một chiều nhỏ, nhỏ hơn giá trị quy định nhiều lần nhưng vẫn xảy ra phóng điện do sét". Nguyên nhân cốt lõi là dòng điện sét mang bản chất xung dốc cao (thời gian đầu sóng $0,2 - 1,2,\mu\text{s}$, dải tần tương đương lên đến $1,\text{MHz}$), khiến phản ứng của hệ thống nối đất bị chi phối hoàn toàn bởi đặc tính quá độ (transient impulse impedance), hiện tượng cảm kháng ($\omega L$), dung kháng ($\omega C$), quá trình ion hóa đất phi tuyến và hiện tượng truyền sóng trên thân cột điện.

Luận án tiến sĩ kỹ thuật điện của NCS. Nguyễn Xuân Phúc (người hướng dẫn: TS. Phạm Hồng Thịnh, PGS. Trần Văn Tớp; Trường Đại học Bách Khoa Hà Nội, 2021, Mã số: 9520201) với đề tài "Nghiên cứu đáp ứng quá độ của hệ thống nối đất đối với dòng điện sét trên đường dây truyền tải Việt Nam" đã giải quyết triệt để khoảng trống này. Nghiên cứu đặt ra 4 câu hỏi nghiên cứu (RQs) và 4 giả thuyết khoa học (Hs) chính:

• RQ1 & H1: Đáp ứng tần số cao làm biến đổi tổng trở nối đất như thế nào so với $R_{DC}$? Giả thuyết: Ở dải tần sét ($>100,\text{kHz}$), thành phần điện cảm làm tổng trở xung vượt xa $R_{DC}$, đặc biệt đối với các cấu hình có tính khép kín.
• RQ2 & H2: Kích thước hình học và tương tác giữa các điện cực có tồn tại ngưỡng bão hòa tản sét? Giả thuyết: Tồn tại chiều dài hiệu dụng ($l_{eff}$) và hiệu ứng che chắn hỗ cảm giữa các cọc/tia làm suy giảm hiệu quả tản sét của phần tử ở trung tâm.
• RQ3 & H3: Hiện tượng ion hóa đất dưới xung dòng biên độ lớn ($30,\text{kA}$) bù trừ thành phần điện cảm ra sao? Giả thuyết: Điện trường vượt ngưỡng tới hạn $E_c$ tạo ra kênh dẫn phóng điện xung quanh điện cực, làm giảm mạnh tổng trở quá độ cục bộ.
• RQ4 & H4: Sự phối hợp sóng giữa thân cột điện và hệ thống nối đất ảnh hưởng thế nào đến quá điện áp trên chuỗi cách điện? Giả thuyết: Mô hình điện trở tập trung bỏ qua tổng trở sóng của cột ($Z_T$) gây đánh giá thấp điện áp đỉnh cột lên tới 15%.

Phạm vi nghiên cứu bao quát hơn 12.000 vị trí cột thực tế thuộc Công ty Truyền tải điện 1 (PTC1) và tuyến đường dây 220 kV Hòa Bình – Tây Hà Nội, cung cấp luận cứ khoa học định lượng để chuyển dịch từ thiết kế $R_{DC}$ sang thiết kế đáp ứng tổng trở xung quá độ toàn diện.

Literature Review và Positioning

Lịch sử phát triển lý thuyết nối đất khởi đầu từ các giải pháp giải tích tĩnh điện và dừng của Dwight (1936), Sunde (1949), Schwarz (1954) và Tagg (1964), vốn là nền tảng xây dựng nên các công thức tính toán điện trở nối đất một chiều đang áp dụng trong IEEE Std 80-2013 và IEEE Std 142-2007. Khi tiếp cận dòng xung sét tần số cao, trường phái nghiên cứu quốc tế phân hóa thành 3 luồng phương pháp luận chính:

                  ┌─────────────────────────────────────────────────────────┐
                  │ CÁC PHƯƠNG PHÁP TÍNH TOÁN NỐI ĐẤT QUÁ ĐỘ DÒNG SÉT       │
                  └────────────────────────────┬────────────────────────────┘
                                               │
         ┌─────────────────────────────────────┼─────────────────────────────────────┐
         │                                     │                                     │
         ▼                                     ▼                                     ▼
┌──────────────────┐                 ┌──────────────────┐                 ┌──────────────────┐
│  Mạch tập trung  │                 │  Đường dây dài   │                 │   Trường điện từ │
│ (Circuit Models) │                 │(Distributed TLM) │                 │(Full-wave Field) │
├──────────────────┤                 ├──────────────────┤                 ├──────────────────┤
│• Meliopoulos1987 │                 │• Sunde (1940)    │                 │• Grcev (1992)    │
│• Geri (1999)     │                 │• Bewley (1934)   │                 │• Rachidi (2011)  │
│• Đơn giản, nhanh │                 │• Lorentzou (2003)│                 │• Phương pháp MoM │
│• Bỏ qua sóng &   │                 │• Thottappillil   │                 │  FDTD, và FEM    │
│  hỗ cảm không gian│                │• Giới hạn cấu    │                 │• Giải Maxwell đầy│
│                  │                 │  hình phức tạp   │                 │  đủ, độ chính xác│
│                  │                 │                  │                 │  cao nhất        │
└──────────────────┘                 └──────────────────┘                 └──────────────────┘

1. Phương pháp mô hình mạch tập trung và tham số rải (Circuit & Network Approach): Điển hình bởi Meliopoulos (1987), Geri (1999), Ramamoorty (1989). Phương pháp này rời rạc hóa hệ thống nối đất thành mạng lưới $R-L-C$, thuận tiện tích hợp vào các phần mềm mô phỏng EMTP/ATP nhưng gặp sai số lớn khi mô tả tương tác điện từ không gian 3 chiều và bức xạ sóng tần số cao.
2. Phương pháp đường dây truyền tải (Transmission Line Modeling - TLM): Khởi xướng từ Sunde (1940), Bewley (1934) và phát triển bởi Thottappillil (2001, 2005, 2007), Lorentzou et al. (2003). TLM phân đoạn điện cực thành các đoạn vi phân có tính đến độ trễ truyền sóng. Tuy nhiên, TLM chỉ áp dụng hiệu quả cho các cấu hình dây thẳng đơn giản (tia counterpoise), không xử lý trọn vẹn được góc uốn cong, mối nối phức tạp và sự không đồng nhất đa chiều của địa tầng.
3. Phương pháp số phân tích trường điện từ sóng đầy đủ (Full-wave Electromagnetic Field Approach): Sử dụng các phương pháp số giải hệ phương trình vi tích phân Maxwell như Phương pháp Mô-men (Method of Moments - MoM) bởi Grcev (1992, 1996, 2008), Rachidi (2007, 2011), Jinliang He (2005, 2019); Phương pháp Sai phân hữu hạn miền thời gian (FDTD) bởi Ametani (2006), Tatematsu (2015); và Phương pháp Phần tử hữu hạn (Finite Element Method - FEM) bởi Meunier (1996), Hatziargyriou (1999), Sadiku (2001), CIGRE WG C4.401 (Brochure 83, 2013).

Tranh luận học thuật cốt lõi (Scientific Debates):

• Quan điểm 1 (Trường phái quy phạm truyền thống - IEEE Std 81, GOST): Cho rằng duy trì $R_{DC} \le 10-30,\Omega$ là đủ để đảm bảo an toàn tản sét cho cột đường dây, việc tăng chiều dài tia và bổ sung cọc luôn tỷ lệ thuận với khả năng hạ thấp trở kháng.
• Quan điểm 2 (Trường phái sóng điện từ hiện đại - Grcev, Visacro 2006, Darveniza 2005): Chứng minh rằng ở tần số cao, cảm kháng tự cảm và hỗ cảm chặn dòng sét thâm nhập sâu vào các phân đoạn xa của điện cực, tạo ra hiện tượng "chiều dài hiệu dụng", khiến việc kéo dài tia vượt ngưỡng trở nên vô nghĩa về mặt kinh tế - kỹ thuật.

Định vị học thuật của luận án: Luận án tiến sĩ của Nguyễn Xuân Phúc định vị ở đỉnh cao của trường phái phân tích trường điện từ bằng cách tiên phong ứng dụng 3D-FEM kết hợp thuật toán lặp phi tuyến động mô phỏng trọn vẹn hệ tương tác đa vật lý: Xung sét $\rightarrow$ Thân cột thép 3D $\rightarrow$ Hệ thống nối đất phức hợp $\rightarrow$ Môi trường đất phân tầng và phi tuyến ion hóa. Công trình thiết lập bước tiến vượt bậc khi đối sánh trực tiếp với 2 nghiên cứu thực nghiệm kinh điển quốc tế:

• Kiểm chứng mô hình điện cực cọc trong đất hai tầng của Ahmeda et al. [88]: Trị số FEM tính toán đạt $75,9,\Omega$, hoàn toàn khớp với thực nghiệm ($76,\Omega$) và lời giải giải tích Tagg ($75,8,\Omega$).
• Kiểm chứng mô hình đáp ứng quá độ tia nối đất của Harid et al. [89] (xung $4,9/13,\mu\text{s}$, biên độ $2,4,\text{A}$) và phản ứng quá độ cột truyền tải của Yamada et al. [90] (xung dốc $1/70,\mu\text{s}$, cột cao $120,\text{m}$), khẳng định tính chuẩn xác tuyệt đối của mô hình số.

Đóng góp lý thuyết và khung phân tích

┌──────────────────────────────────────────────────────────────────────────┐
│              KHUNG PHÂN TÍCH QUÁ ĐỘ ĐIỆN TỪ ĐA TẦNG (FEM-HYBRID)        │
└────────────────────────────────────┬─────────────────────────────────────┘
                                     │
         ┌───────────────────────────┼───────────────────────────┐
         │                           │                           │
         ▼                           ▼                           ▼
┌──────────────────┐       ┌───────────────────┐       ┌───────────────────┐
│ MIỀN TẦN SỐ      │       │ MIỀN THỜI GIAN    │       │ TƯƠNG TÁC CỘT     │
│ (1 Hz - 1 MHz)   │       │ (Xung 1,2/50 μs)  │       │ & NỐI ĐẤT TOÀN DIỆN│
├──────────────────┤       ├───────────────────┤       ├───────────────────┤
│• 3 dải tần số:   │       │• Chiều dài hiệu   │       │• Tổng trở sóng cột│
│  Thuần trở,      │       │  dụng l_eff       │       │  Z_T (500kV vs    │
│  Cảm kháng,      │       │• Thuật toán lặp   │       │  220kV)           │
│  Dung kháng      │       │  ion hóa đất:     │       │• Sóng phản xạ chân│
│• Hiệu ứng vòng   │       │  E >= E_c         │       │  cột              │
│  kín (Quấn vòng  │       │• Bán kính ion hóa │       │• Điện áp chuỗi    │
│  tăng 50% Z)     │       │  r_i(t)           │       │  cách điện        │
└──────────────────┘       └───────────────────┘       └───────────────────┘

Đóng góp cho lý thuyết

Luận án đã mở rộng và tái cấu trúc hệ thống lý thuyết tiếp địa chống sét qua các đóng góp nền tảng:

1. Bác bỏ giả định an toàn của mô hình tĩnh điện: Chứng minh khái niệm "cảm giác an toàn giả tạo" thông qua bằng chứng định lượng: Một hệ thống có $R_{DC}$ đạt chuẩn ($10,\Omega$) vẫn có thể tạo ra tổng trở xung $Z_p > 35,\Omega$ dưới tác động của xung sét đầu sóng dốc, dẫn đến quá điện áp vượt ngưỡng cách điện $U_{50%}$ gây sự cố phóng điện ngược (back-flashover).
2. Mô hình hóa lý thuyết đáp ứng 3 vùng tần số của điện cực tiếp địa:
   • Vùng tần số thấp ($f < 10,\text{kHz}$): Trở kháng không đổi, bằng $R_{DC}$, quyết định bởi độ dẫn điện của đất $\sigma$.
   • Vùng tần số trung bình ($10,\text{kHz} \le f \le 100,\text{kHz}$): Xuất hiện tính chất cảm kháng rõ rệt ($\omega L$), góc pha tổng trở chuyển dịch dương ($+10^\circ \div +45^\circ$), tổng trở tăng dần theo tần số.
   • Vùng tần số cao ($f > 100,\text{kHz} - 1,\text{MHz}$): Hiện tượng dung kháng ($\omega \varepsilon$) bắt đầu phát huy tác dụng làm suy giảm tổng trở ở đất có điện trở suất cao, tuy nhiên đối với cấu hình tia dài hoặc quấn vòng, cảm kháng vẫn áp đảo hoàn toàn.
3. Mô hình toán học quá trình phóng điện ion hóa phi tuyến trong đất: Thiết lập mối quan hệ phụ thuộc không gian - thời gian của mật độ dòng rò $J(t)$, điện trường cục bộ $E(r,t)$ và bán kính ion hóa $r_{ion}(t)$ dựa trên tiêu chuẩn điện trường đánh thủng tới hạn $E_c = 241\cdot \rho^{0,1},(\text{kV/m})$ (theo mô hình Oettle, Darveniza và CIGRE).

Khung phân tích độc đáo

Khung phân tích của luận án tích hợp đồng thời 3 hệ lý thuyết: Lý thuyết điện từ trường vi phân Maxwell, Lý thuyết đường dây truyền sóng phi đồng nhất (Non-uniform Transmission Lines) và Lý thuyết nhiệt - điện ly trong môi trường bán dẫn xốp (đất).

Mô hình thiết lập điều kiện biên không giới hạn (Unbounded domain) thông qua lớp phần tử vô hạn (Infinite Element Domain) dày $2,\text{m}$ bao bọc biên không gian vật lý, kết hợp cổng kích thích tập trung (Lumped Port) bơm trực tiếp dòng xung $i(t)$ dạng Heidler hoặc xung chuẩn $1,2/50,\mu\text{s}$: $$i(t) = \frac{I_p}{\eta} \cdot \frac{(t/\tau_1)^n}{1 + (t/\tau_1)^n} \cdot \exp(-t/\tau_2)$$

Điều kiện giới hạn biên (Boundary conditions) được xác định rõ: Môi trường đất đồng nhất hoặc phân 2 lớp địa điện ($\rho_1, h_1; \rho_2$), độ từ thẩm tương đối $\mu_r = 1$, hằng số điện môi tương đối $\varepsilon_r \in [5, 80]$, điện trở suất $\rho \in [100, 5000],\Omega\cdot\text{m}$.

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

Thiết kế nghiên cứu

Nghiên cứu tuân thủ chặt chẽ thế giới quan thực chứng (positivism) kết hợp mô hình hóa toán - lý chính xác cao (computational electromagnetics). Thiết kế nghiên cứu đa tầng kết nối từ cấu trúc vi mô phân đoạn điện cực đến vĩ mô cột tháp truyền tải 500 kV/220 kV:

┌──────────────────────────────────────────────────────────────────────────┐
│                      MA TRẬN DỮ LIỆU & THIẾT KẾ MÔ PHỎNG                 │
├──────────────────┬───────────────────────────────┬───────────────────────┤
│ CẤU HÌNH ĐIỆN CỰC│ KÍCH THƯỚC HÌNH HỌC KHẢO SÁT  │ ĐIỀU KIỆN ĐỊA CHẤT    │
├──────────────────┼───────────────────────────────┼───────────────────────┤
│ 1. Tia đơn       │ L = 10m, 20m, 30m, 40m, 50m   │ ρ = 1000 - 5000 Ω·m   │
│    (Chôn sâu 0,8m│ Đường kính d = 16 - 20 mm     │ ε_r = 5 - 80          │
├──────────────────┼───────────────────────────────┼───────────────────────┤
│ 2. Nhiều tia     │ 2 tia, 4 tia song song        │ Đất đồng nhất &       │
│    song song     │ Khoảng cách d = 2m, 5m, 10m   │ đất phân 2 tầng       │
├──────────────────┼───────────────────────────────┼───────────────────────┤
│ 3. Cọc - Tia     │ Tia 36m + 12 cọc (cách 3m)    │ Cùng chuẩn hóa về     │
│    kết hợp       │ Tia 40m + 8 cọc (cách 5m)     │ R_DC = 33 Ω           │
│                  │ Tia 44m + 3 cọc (cách 20m)    │                       │
├──────────────────┼───────────────────────────────┼───────────────────────┤
│ 4. Tia quấn vòng │ 1 vòng (8x12m) + tia 15m      │ Khảo sát dòng vòng    │
│    quanh móng    │ 4 vòng đồng tâm (cách 2m)     │ khép kín              │
└──────────────────┴───────────────────────────────┴───────────────────────┘

Quy trình nghiên cứu và công cụ phần mềm

Mô hình hóa được thực hiện bằng cách ghép nối 2 nền tảng tính toán mạnh nhất hiện nay: COMSOL Multiphysics (RF Module / Wave Optics) và MATLAB Engine qua giao thức tự động hóa API LiveLink.

┌────────────────────────────────────────────────────────────────────────┐
│      LƯU ĐỒ THUẬT TOÁN LẶP DỰNG HÌNH & TÍNH TOÁN ION HÓA ĐẤT (MATLAB - COMSOL) │
└───────────────────────────────────┬────────────────────────────────────┘
                                    │
                                    ▼
       ┌──────────────────────────────────────────────────────────┐
       │ 1. Khởi tạo mô hình 3D FEM trong COMSOL:                 │
       │    Geometry, Meshing, Soil parameters (ρ_0, ε_r)         │
       └────────────────────────────┬─────────────────────────────┘
                                    │
                                    ▼
       ┌──────────────────────────────────────────────────────────┐
       │ 2. Thiết lập bước thời gian tm = [0 : 0.05 : 5] μs        │
       │    Bơm dòng kích thích sét I(tm) vào Lumped Port         │
       └────────────────────────────┬─────────────────────────────┘
                                    │
                                    ▼
       ┌──────────────────────────────────────────────────────────┐
       │ 3. Giải hệ phương trình Maxwell trích xuất mật độ dòng:  │
       │    I_k = mphint2(model, 'temw.J · n', 'surface', seg_k)  │
       └────────────────────────────┬─────────────────────────────┘
                                    │
                                    ▼
       ┌──────────────────────────────────────────────────────────┐
       │ 4. Tính toán điện trường bề mặt và bán kính ion hóa:     │
       │    r_ion,k = (ρ · I_k) / (2π · l_k · E_c)                │
       └────────────────────────────┬─────────────────────────────┘
                                    │
                                    ▼
                         /──────────────────────\
                        <  r_ion,k > r_conduct?  >
                         \──────────────────────/
                                    │
                     ┌──────────────┴──────────────┐
                  YES│                             │NO
                     ▼                             ▼
       ┌───────────────────────────┐ ┌───────────────────────────┐
       │ Cập nhật bán kính dẫn mới │ │ Giữ nguyên bán kính       │
       │ ρ_k = 0 (vùng dẫn sét mở  │ │ ban đầu và điện trở       │
       │ rộng)                     │ │ suất nền ρ_0              │
       └─────────────┬─────────────┘ └─────────────┬─────────────┘
                     │                             │
                     └──────────────┬──────────────┘
                                    │
                                    ▼
       ┌──────────────────────────────────────────────────────────┐
       │ 5. Tái cấu trúc ma trận tham số (model.run)              │
       │    Tính điện áp U_0(tm), dòng I_0(tm) => Z_transient(tm) │
       └────────────────────────────┬─────────────────────────────┘
                                    │
                                    ▼
       ┌──────────────────────────────────────────────────────────┐
       │ 6. Tăng tm = tm + dt; Lặp lại cho đến hết chu kỳ sóng sét│
       └──────────────────────────────────────────────────────────┘

Mã nguồn MATLAB tự động hóa quá trình tích phân mật độ dòng điện rò trên từng phân đoạn bề mặt:

% Trích đoạn thuật toán lặp tính toán dòng rò và ion hóa từng phân đoạn
I11 = mphint2(model,'(temw.Jx*nx+temw.Jy*ny+temw.Jz*nz)','surface','selection',seg1_tags);
rion1 = (ro * I11) / (2 * 3.14159 * l1 * Ec);
if rion1 > r_wire
    ro1 = 0; % Điện trở suất vùng ion hóa sụp đổ thành vật dẫn
    model.set('ro1', ro1);
else
    model.set('ro1', ro);
end
model.run; % Tái giải ma trận trường sóng

Kỹ thuật chia lưới (Meshing) và kiểm soát sai số

Để xử lý bài toán đa tỷ lệ (từ bán kính dây thép $r = 8,\text{mm}$ đến không gian đất $50,\text{m} \times 50,\text{m} \times 30,\text{m}$ và chiều cao cột điện $120,\text{m}$), kỹ thuật chia lưới thích ứng cục bộ (Adaptive User-controlled Mesh) được áp dụng:

• Lưới cực mịn (Extremely fine) tại bề mặt dây dẫn và vùng lân cận cọc để nắm bắt gradient điện trường $\nabla E$ cực lớn.
• Lưới thô dần theo hàm mũ khi tiến ra vùng biên xa vô cực (Infinite Element Domain), đảm bảo nghiệm hội tụ với sai số tương đối $< 0,5%$, tối ưu hóa dung lượng RAM và thời gian xử lý ma trận phần tử hữu hạn.

Phát hiện đột phá và implications

┌─────────────────────────────────────────────────────────────────────────┐
│                    TỔNG HỢP 5 PHÁT HIỆN ĐỘT PHÁ CỦA LUẬN ÁN             │
├─────────────────────────────────────────────────────────────────────────┤
│ 1. NỐI ĐẤT QUẤN VÒNG LÀM TĂNG TỔNG TRỞ XUNG LÊN ĐẾN 50% SO VỚI DẠNG CỌC  │
│    Do dòng điện chạy vòng khép kín sinh hỗ cảm ngược cực lớn.           │
├─────────────────────────────────────────────────────────────────────────┤
│ 2. HIỆU ỨNG CHIỀU DÀI HIỆU DỤNG: KÉO DÀI TIA > 30-40m KHÔNG CÒN TÁC DỤNG │
│    Cảm kháng đường dây chặn dòng điện sét lan truyền tới phần đuôi tia.  │
├─────────────────────────────────────────────────────────────────────────┤
│ 3. HIỆU ỨNG CHE CHẮN HỖ CẢM TRONG HỆ CỌC - TIA:                         │
│    Các cọc bố trí ở giữa tia bị vô hiệu hóa; chỉ cọc ở 2 đầu phát huy   │
│    tối đa khả năng tản sét.                                             │
├─────────────────────────────────────────────────────────────────────────┤
│ 4. ĐỘ DỐC ĐẦU SÓNG QUYẾT ĐỊNH QUÁ ĐIỆN ÁP ĐỈNH CỘT:                     │
│    Xung dốc 0,2 μs làm điện áp đỉnh cột tăng gấp 2 lần xung 1,2/50 μs    │
│    mặc dù tổng trở sóng cột Z_T không đổi.                              │
├─────────────────────────────────────────────────────────────────────────┤
│ 5. MÔ HÌNH NỐI ĐẤT TẬP TRUNG GÂY SAI SỐ NGUY HIỂM 15%:                  │
│    Bỏ qua truyền sóng 3D trên thân cột làm đánh giá sai lệch khả năng   │
│    chịu sét thực tế của đường dây 500 kV và 220 kV.                    │
└─────────────────────────────────────────────────────────────────────────┘

Những phát hiện then chốt (Key Empirical Findings)

1. Sự kém hiệu quả của hệ thống nối đất dạng quấn vòng móng: Kết quả mô phỏng miền tần số ($1,\text{MHz}$) và miền thời gian khẳng định hệ thống quấn vòng tạo ra các vòng lặp dòng điện khép kín, làm tích tụ từ thông hỗ cảm cực lớn. Tổng trở xung quá độ của cấu hình quấn vòng cao hơn dạng cọc và dạng tia tới 50% mặc dù chúng có cùng giá trị $R_{DC}$. Đây là phát hiện cảnh báo nghiêm trọng đối với các thiết kế trạm và cột đường dây truyền tải gần khu dân cư đang lạm dụng nối đất quấn vòng.
2. Xác lập giới hạn chiều dài hiệu dụng ($l_{eff}$): Đối với xung sét $1,2/50,\mu\text{s}$ trên nền đất $\rho = 1000,\Omega\cdot\text{m}$, khi chiều dài tia vượt quá $30 - 40,\text{m}$, biên độ điện áp dâng đỉnh ($V_p$) và tổng trở xung ($Z_p$) không hề suy giảm thêm. Việc tăng chiều dài tia lên $50 - 60,\text{m}$ chỉ làm giảm $R_{DC}$ danh nghĩa mà hoàn toàn không mang lại giá trị chống sét quá độ.
3. Hiệu ứng triệt tiêu tản dòng ở hệ thống cọc - tia: Khảo sát hệ thống cọc - tia chuẩn hóa ($R_{DC} = 33,\Omega$) với các cấu hình (Tia $36,\text{m} + 12$ cọc; Tia $40,\text{m} + 8$ cọc; Tia $44,\text{m} + 3$ cọc) chỉ ra rằng: Do hiệu ứng ghép hỗ cảm không gian, các cọc cắm ở vị trí giữa tia chịu sự tích tụ điện trường đẩy lùi dòng rò, mật độ dòng tản vào đất tại các cọc giữa chỉ bằng $15 - 20%$ so với cọc ở đầu tia. Cấu hình phân tán ít cọc nhưng khoảng cách xa ($20,\text{m}$) có hiệu quả tản sét vượt trội so với bố trí cọc dày đặc ($3,\text{m}$).
4. Tác động khuếch đại của sóng phản xạ trên thân cột và thời gian đầu sóng: Khi xét hệ thống liên hợp hoàn chỉnh (Cột truyền tải + Nối đất):
   • Tổng trở sóng của cột đỡ 500 kV ($Z_T = 120 - 150,\Omega$) lớn hơn cột 220 kV ($Z_T = 90 - 110,\Omega$).
   • Dưới cùng một xung dòng sét, điện áp tại xà trên của cột 500 kV cao hơn 1,3 lần so với cột 220 kV, và điện áp dâng chân cột cao hơn 1,1 lần.
   • Dạng xung đầu sóng dốc ($0,2,\mu\text{s}$) làm tăng gấp đôi (200%) điện áp đỉnh cột so với xung tiêu chuẩn ($1,2,\mu\text{s}$), đẩy nguy cơ phóng điện ngược chuỗi cách điện lên mức báo động.
5. Sai số mang tính hệ thống của mô hình điện trở tập trung: So sánh trực tiếp giữa mô hình nối đất tập trung ($R_{DC}$) và mô hình đầy đủ (FEM 3D phân bố) cho thấy: Mô hình tập trung đánh giá thấp điện áp đỉnh xà từ 10% đến 15%. Điều này chứng minh các phần mềm tính toán suất cắt truyền thống đang ước lượng thiếu an toàn đối với các cột cao và điện cực dài.

Bảng tổng hợp so sánh các cấu hình tiếp địa (Chuẩn hóa cùng $R_{DC} = 33,\Omega$, $\rho = 1000,\Omega\cdot\text{m}$)

Khuyến nghị kỹ thuật và lộ trình chuyển đổi (Actionable Roadmap)

• Giải pháp 1 (Tái cấu trúc hình học): Tuyệt đối không áp dụng nối đất quấn vòng khép kín cho mục đích chống sét tại các vị trí cột có nguy cơ cao. Chuyển đổi sang cấu hình phân nhánh tia hình sao (4 tia ngắn $15 - 20,\text{m}$) bố trí góc mở rộng $\ge 60^\circ$ để triệt tiêu hỗ cảm.
• Giải pháp 2 (Tối ưu hóa hệ cọc - tia): Khoảng cách lắp đặt giữa các cọc nối đất dọc theo tia phải đảm bảo tối thiểu $d \ge 10 - 20,\text{m}$ (thay vì $3 - 5,\text{m}$ như thói quen thiết kế cũ). Ưu tiên đóng cọc tại điểm xa nhất của tia để đón đầu sóng phản xạ.
• Giải pháp 3 (Cải cách tiêu chuẩn ngành): Đề xuất EVN và Bộ Công Thương cập nhật Quy phạm trang bị điện: Đối với đường dây cấp điện áp $\ge 220,\text{kV}$, bắt buộc áp dụng tiêu chuẩn tổng trở xung quá độ $Z_p$ (tính toán tại tần số $500,\text{kHz} - 1,\text{MHz}$) thay thế cho chỉ tiêu duy nhất $R_{DC}$.

Limitations và Future Research

Nhằm duy trì chuẩn mực học thuật khách quan, luận án và các phân tích chuyên sâu ghi nhận các giới hạn biên sau:

1. Sự phụ thuộc vào tần số của thông số điện môi đất ($\rho(f), \varepsilon_r(f)$): Mô hình nghiên cứu chính coi $\rho$ và $\varepsilon_r$ là hằng số theo dải tần trong mỗi kịch bản tính toán. Thực tế, các nghiên cứu thực nghiệm của Visacro và Alipio (2012) chỉ ra rằng điện trở suất đất có xu hướng giảm và hằng số điện môi giảm mạnh khi tần số tăng lên hàng MHz.
2. Giả thiết đường dẫn dòng sét đơn lẻ qua cột: Mô hình hệ thống đầy đủ tập trung khảo sát quá trình sét đánh trực tiếp vào đỉnh cột và tản xuống đất qua thân cột duy nhất, chưa tích hợp mô hình phân nhánh dòng sét ngược về 2 khoảng cột lân cận qua dây chống sét (dây chống sét OPGW/phân dòng đa cột).
3. Mô hình ion hóa đất dạng phân đoạn đẳng hướng: Mô hình giả định vùng ion hóa phát triển hình trụ tròn đồng tâm dọc theo dây dẫn, bỏ qua tính chất bất đối xứng vi mô của cấu trúc hạt đất và độ ẩm cục bộ.

Hướng nghiên cứu phát triển tiếp theo (Future Research Agenda):

• Xây dựng mạch tương đương đa cổng phân bố trong miền thời gian (Vector Fitting / Black-box Macromodel) từ ma trận tán xạ FEM để nhúng trực tiếp vào các phần mềm mô phỏng hệ thống lưới điện thời gian thực như EMTP-RV, PSCAD/EMTDC.
• Nghiên cứu mô hình trường điện từ kết hợp dây chống sét đa khoảng cột có xét đến hồ quang phóng điện ngược trên chuỗi cách điện (Flashover dynamics).
• Thực nghiệm đo đạc đáp ứng xung thực địa tại các vị trí cột hiểm trở bằng máy phát xung di động công suất lớn để kiểm chứng mô hình ion hóa đa tầng.

Tác động và ảnh hưởng

┌────────────────────────────────────────────────────────────────────────┐
│               MA TRẬN LAN TỎA TÁC ĐỘNG (IMPACT MATRIX)                 │
├───────────────────┬────────────────────────────────────────────────────┤
│ HỌC THUẬT QUỐC TẾ │ • Trích dẫn trong các kỷ yếu IEEE, CIGRE, ISET     │
│                   │ • Cung cấp bộ dữ liệu benchmark 3D-FEM chuẩn xác   │
├───────────────────┼────────────────────────────────────────────────────┤
│ NGÀNH ĐIỆN LỰC    │ • Tiết kiệm 20 - 30% khối lượng vật liệu thép tiếp │
│ (EVN / EVNNPT)    │   địa thông qua cắt giảm chiều dài tia dư thừa     │
│                   │ • Giảm thiểu suất cắt sự cố do sét trên lưới 500kV │
├───────────────────┼────────────────────────────────────────────────────┤
│ QUẢN LÝ NHÀ NƯỚC  │ • Cơ sở khoa học sửa đổi Quy phạm trang bị điện    │
│                   │ • Đảm bảo an ninh năng lượng và truyền tải liên tục│
└───────────────────┴────────────────────────────────────────────────────┘

• Tác động học thuật: Đóng góp hệ phương pháp luận kết hợp FEM-LiveLink tự động hóa đã được công bố tại các hội nghị khoa học quốc tế uy tín (IEEE ISEIM, RCEEE, ISET) và các tạp chí chuyên ngành (Tạp chí KH&CN Năng lượng, ĐH Điện lực; Tạp chí KH&CN ĐH Công nghiệp Hà Nội). Tạo tiền đề cho các nghiên cứu sinh chuyên ngành Cao áp và Hệ thống điện phát triển các mô hình trường điện từ phức hợp.
• Tác động công nghiệp và kinh tế: Giúp các đơn vị tư vấn thiết kế điện lực (PECC1, PECC2, PECC3, PECC4) tối ưu hóa thiết kế hệ thống nối đất, loại bỏ việc kéo dài vô ích hàng chục mét tia thép mạ kẽm tại các vị trí đồi núi đá hiểm trở, mang lại hiệu quả tiết kiệm hàng chục tỷ đồng chi phí xây lắp hàng năm cho EVNNPT, đồng thời ngăn ngừa sự cố rã lưới do phóng điện lan truyền.
• Lợi ích an toàn xã hội: Làm rõ sự phân bố điện thế tiếp xúc ($U_{tx}$) và điện áp bước ($U_b$) xung quanh móng cột khi có dòng sét đánh, giúp khoanh vùng cảnh báo nguy hiểm chính xác, bảo vệ an toàn tính mạng cho người dân và gia súc tại các hành lang tuyến đi qua nương rẫy, khu dân cư miền núi.

Đối tượng hưởng lợi

• Nghiên cứu sinh & Giảng viên ngành Kỹ thuật điện: Tiếp cận một tài liệu mẫu mực về phương pháp số 3D-FEM ứng dụng trong kỹ thuật điện cao áp; hiểu rõ bản chất vật lý của sự tương tác sóng điện từ - môi trường đất thay vì các công thức giải tích giải trình đơn giản.
• Kỹ sư thiết kế lưới điện truyền tải: Nắm vững cẩm nang định hình cấu hình tiếp địa tối ưu cho từng loại địa hình: Biết cách từ chối thiết kế quấn vòng móng khi chống sét, biết cách tối ưu hóa tỷ lệ chiều dài tia/khoảng cách cọc cho các tuyến đường dây siêu cao áp 500 kV Bắc - Nam.
• Cán bộ quản lý vận hành truyền tải (EVNNPT, các PTC): Có cơ sở lý luận chính xác để giải trình các sự cố phóng điện bí ẩn tại các cột điện có $R_{DC}$ rất thấp; từ đó xây dựng phương án cải tạo chống sét trúng đích (lắp đặt chống sét van đường dây kết hợp cải tạo hình học tiếp địa thay vì chỉ khoan thêm giếng tiếp địa tốn kém).
• Cơ quan ban hành chính sách và quy chuẩn (Cục Điều tiết Điện lực, Viện Năng lượng): Sở hữu cứ liệu khoa học thực chứng để sửa đổi, hoàn thiện Quy chuẩn Kỹ thuật Quốc gia về kỹ thuật điện, tiệm cận các tiêu chuẩn tiên tiến của quốc tế như IEEE Std 1863 và CIGRE.

Câu hỏi chuyên sâu

1. Đóng góp lý thuyết nào của luận án mang tính độc đáo nhất và đã mở rộng lý thuyết nào?

Đóng góp độc đáo nhất là việc lượng hóa bản chất cảm kháng áp đảo của hệ thống tiếp địa quấn vòng khép kín và mở rộng lý thuyết trường điện từ sóng đầy đủ Maxwell vào bài toán hệ thống liên hợp Cột - Nối đất thực tế tại Việt Nam. Luận án đã mở rộng lý thuyết phân tích quá độ của Grcev (1996) và lý thuyết sóng phản xạ cột của Jordan (1934), Chisholm (1983), chứng minh rằng điện trở một chiều $R_{DC}$ truyền thống hoàn toàn mất hiệu lực trong việc đại diện cho khả năng tản sét của các cấu hình điện cực phức tạp.

2. Sự đổi mới về mặt phương pháp nghiên cứu so với các công trình đi trước được thể hiện như thế nào?

So với các nghiên cứu trong nước trước đây (như Nguyễn Xuân Cường 2006 chỉ dùng mô hình đường dây dài; Vũ Phan Tú 2010 chỉ mô phỏng lưới nối đất trạm bằng FEM tĩnh; Đỗ Văn Bản 2008 chỉ xét tia dài đơn thuần), luận án đã tạo ra bước đột phá về công cụ khi kết hợp thuật toán lặp động đa vật lý giữa COMSOL Multiphysics và MATLAB. Mô hình này cho phép tự động tái tính toán bán kính ion hóa $r_{ion}(t)$ và cập nhật ma trận độ dẫn điện $\sigma(r,t)$ của môi trường đất theo từng bước thời gian $0,05,\mu\text{s}$, giải quyết trọn vẹn tính phi tuyến mạnh của đất dưới tác động dòng sét $30,\text{kA}$.

      So sánh tương quan giữa các phương pháp mô hình hóa:
      
      Phương pháp            Độ chính xác     Khả năng giải bài toán phức tạp
      ────────────────────────────────────────────────────────────────────────
      Mạch tập trung (R-L-C)    Thấp           Hạn chế (Bỏ qua tương hỗ 3D)
      Đường dây dài (TLM)       Trung bình     Chỉ áp dụng cho tia thẳng đơn
      3D-FEM của Luận án        Rất cao        Đa vật lý, Đa tầng, Phi tuyến ion hóa
      ────────────────────────────────────────────────────────────────────────

3. Phát hiện bất ngờ nhất (counter-intuitive) có dữ liệu thực nghiệm/mô phỏng chứng minh là gì?

Phát hiện bất ngờ nhất là "Nghịch lý hệ thống nối đất quấn vòng quanh móng cột". Trong thực tế thiết kế, dạng quấn vòng thường được ưu tiên tại các vị trí chật hẹp vì kỹ sư cho rằng việc cuộn dây nhiều vòng sẽ làm tăng diện tích tiếp xúc với đất và giảm điện trở một chiều. Tuy nhiên, dữ liệu mô phỏng trường sóng của luận án chỉ ra rằng ở tần số dòng sét ($1,\text{MHz}$), do hiện tượng cảm ứng điện từ vòng kín, tổng trở quá độ của hệ thống quấn vòng lớn hơn hệ thống dạng cọc tới 50%, biến cấu hình này thành cấu hình tản sét kém an toàn nhất trong tất cả các dạng tiếp địa.

4. Luận án có cung cấp quy trình/giao thức mô phỏng lặp lại (Replication Protocol) hoàn chỉnh không?

Hoàn toàn có. Luận án cung cấp chi tiết tại Phụ lục Chương 2, 3, 4 toàn bộ các bước thiết lập mô hình 3D:

• Định nghĩa hình học cọc đồng $\phi 16,\text{mm}$, dài $2,4,\text{m}$, chôn trong đất 2 tầng ($9,\text{m}$ tầng trên $\rho_1 = 200,\Omega\cdot\text{m}$, tầng dưới $\rho_2 = 30,\Omega\cdot\text{m}$).
• Thiết lập biên xa vô hạn Infinite Element Domain $2,\text{m}$.
• Thiết lập cổng kích thích Lumped Port với dòng $i(t)$.
• Công bố toàn văn đoạn mã kịch bản script MATLAB tương tác tự động với API COMSOL (mphload, mphint2, mphglobal), cho phép bất kỳ nhà nghiên cứu nào cũng có thể tái lập và kiểm chứng kết quả độc lập.

5. Lộ trình nghiên cứu 10 năm (10-Year Research Agenda) tiếp theo được định hình như thế nào?

Lộ trình 10 năm định hình 3 giai đoạn nghiên cứu chiến lược:

1. Giai đoạn 1 (1 - 3 năm): Mở rộng mô hình FEM sang dải tần số biến thiên phụ thuộc của đất $\rho(f)$ và $\varepsilon(f)$, xây dựng phần mềm chuyên dụng tính toán tổng trở xung tiếp địa tích hợp cơ sở dữ liệu địa chất sét Việt Nam.
2. Giai đoạn 2 (4 - 6 năm): Tích hợp toàn diện mô hình phóng điện hồ quang chuỗi cách điện và phân dòng lưới dây chống sét toàn tuyến trong không gian 3D, tối ưu hóa lắp đặt chống sét van đường dây (TLA).
3. Giai đoạn 3 (7 - 10 năm): Ứng dụng trí tuệ nhân tạo (AI/Machine Learning) huấn luyện mô hình thế chân (Surrogate Models) từ dữ liệu Big Data FEM để dự báo thời gian thực rủi ro phóng điện sét trên toàn bộ hệ thống lưới điện truyền tải số quốc gia (Smart Grid).

Kết luận

1. Luận án đã hoàn thành xuất sắc việc phân tích toàn diện đáp ứng quá độ của các hệ thống nối đất đường dây truyền tải 220 kV và 500 kV Việt Nam trong cả miền tần số ($1,\text{Hz} - 1,\text{MHz}$) và miền thời gian ($1,2/50,\mu\text{s}$ và xung dốc $0,2,\mu\text{s}$) bằng phương pháp phần tử hữu hạn (FEM) 3D.
2. Xác lập cơ sở khoa học định lượng chứng minh điện trở một chiều $R_{DC}$ không phản ánh đúng hiệu quả tản dòng điện sét, vạch rõ nguy cơ từ "cảm giác an toàn giả tạo" trong các tiêu chuẩn thiết kế hiện hành.
3. Phát hiện và chứng minh cấu hình nối đất quấn vòng móng có tổng trở xung cao hơn dạng cọc tới 50% do tích tụ cảm kháng vòng kín; đồng thời chỉ rõ hiệu ứng che chắn hỗ cảm làm tê liệt khả năng tản sét của các cọc nối đất bố trí ở giữa tia.
4. Xác định chính xác ngưỡng chiều dài hiệu dụng của tia nối đất ($30 - 40,\text{m}$) ứng với từng loại điện trở suất đất, loại bỏ việc lãng phí vật tư tiếp địa trong xây lắp truyền tải.
5. Tiên phong xây dựng mô hình liên hợp Cột truyền tải - Nối đất 3D hoàn chỉnh, chứng minh mô hình điện trở tập trung truyền thống gây sai số đánh giá thấp điện áp đỉnh cột tới 15%.
6. Đề xuất hệ thống giải pháp kỹ thuật có tính ứng dụng thực tiễn cao cho EVNNPT, mở ra hướng nghiên cứu đột phá trong việc xây dựng tiêu chuẩn thiết kế chống sét quá độ hiện đại cho hệ thống điện Việt Nam và khu vực.