Trong bài báo này, trình bày các kết quả nghiên cứu xác định một số tham số hoạt động dông sét từ chuỗi số liệu vị trí phóng điện mô phỏng dựa theo phương pháp của Finke cho ba loại dông khác nhau (tương tự như cách phân chia của Peckham).
33(2)[CĐ], 134-141 Tạp chí CÁC KHOA HỌC VỀ TRÁI ĐẤT 6-2011 XÁC ĐỊNH MỘT SỐ THAM SỐ HOẠT ĐỘNG DÔNG SÉT TỪ CHUỖI SỐ LIỆU MƠ PHỎNG HỒNG HẢI SƠN, NGUYỄN XUÂN ANH, LÊ VIỆT HUY, PHẠM XUÂN THÀNH E-mail: hhson@igp-vast.vn Viện Vật lý Địa cầu - Viện Khoa học Công nghệ Việt Nam Ngày nhận bài: 10-4-2011 Mở đầu Việt Nam nằm trong ba tâm dông sét hoạt động mạnh giới, [3, 8] Thiệt hại dông sét gây ngày tăng với q trình cơng nghiệp hố, đại hố đất nước [1, 2] Để giảm thiểu thiệt hại dơng sét gây ra, địi hỏi phải có kiến thức, hiểu biết hoạt động dông sét khu vực, nắm bắt thông số hoạt động dông sét Để nghiên cứu hoạt động dông sét khu vực, số tác giả thường dựa số liệu định vị sét độ phân giải cao để phân loại mây dông nghiên cứu đặc trưng mây dông [4-6] Theo Finke (1999), xác định tham số hoạt động dông sét như: tốc độ lan truyền, thời gian tồn kích thước đám mây dơng từ số liệu định vị phóng điện thông qua hàm tự tương quan chúng Tuy nhiên, để xác định tham số hoạt động dơng sét từ số liệu định vị phóng điện địi hỏi chuỗi số liệu phải có độ phân giải cao Trong tương lai, hệ thống định vị phóng điện có độ phân giải cao thiết lập Việt Nam Hiện tại, để phục vụ mục đích nghiên cứu ứng dụng chuỗi số liệu cho cảnh báo hoạt động dông sét, tạo chuỗi số liệu mô với đặc trưng tương tự chuỗi số liệu thực nghiên cứu chuỗi số liệu Trong báo này, chúng tơi trình bày kết nghiên cứu xác định số tham số hoạt động dông sét từ chuỗi số liệu vị trí phóng điện mơ dựa theo phương pháp Finke cho ba loại dông khác (tương tự cách phân chia Peckham) Tiến trình thực nghiên cứu chia thành bước: (1) Phân loại dông thành 134 dạng khác dựa tham số chúng (tham số giả định); (2) Mơ vị trí tia sét từ tham số giả định dạng mây dơng; (3) Tính hệ số tương quan không gian, thời gian từ chuỗi số liệu mơ từ tính số tham số hoạt động dông sét; (4) Mô tả kết tính tốn so sánh giá trị tham số hoạt động dơng sét tính với tham số giả định ban đầu để đánh giá hiệu phương pháp tính tốn Các bước từ (1) đến (3) trình bày phần phương pháp Bước (4) trình bày phần kết phân tích Phương pháp 2.1 Phân loại mây dông Để nghiên cứu hoạt động dông sét khu vực, số tác giả phân loại mây dông thành dạng khác Theo Peckham [7], mây dơng chia thành ba loại (I, II, IIII) dựa tham số quan sát thực tế trình bày bảng 1; Mark (2005) [6], từ chuỗi số liệu định vị sét khảo sát vùng trung tâm phía đơng nam nước Mỹ với độ xác thiết bị đo đến 0,5km, phân chia mây dông thành dạng: dông địa phương, dông đa địa phương, dơng theo vùng: dạng elíp, dạng tuyến tính, dạng không xác định dông xẩy diện rộng Hagen (1999) [5], dựa chuỗi số liệu định vị sét khảo sát vùng phía nam nước Đức với độ xác thiết bị đo đến 1,0km để phân chia mây dông thành ba loại khác nhau: loại I có thời gian kéo dài nhỏ giờ, loại II có chiều dài cực tiểu 50km chiều rộng cực tiểu 3km; loại III có chiều rộng lớn 100km Bảng Phân loại số dạng mây dông theo tham số, [7] (Par: tham số đo đạc, Std: độ lệch chuẩn) Mây dông loại I Các tham số Thời gian dơng (phút) Diện tích (km ) Số lần phóng điện xuống đất -1 -1 Mật độ phóng điện xuống đất (km phút ) -1 Tốc độ phóng điện trung bình xuống đất (phút ) -1 Tốc độ phóng điện cực đại xuống đất (phút ) Mây dông loại II Mây dông loại III Par Std Par Std Par Std 41 16 77 26 130 51 103 63 256 154 900 841 73 69 270 261 887 720 0,018 0,011 0,015 0,008 0,010 0,006 1,7 1,2 3,4 2,3 6,8 4,7 3,7 2,6 7,3 4,4 14 Dựa theo cách phân chia Peckham, chia mây dông thành ba loại với số tham số tương tự Peckham bổ sung thêm số tham số khác Đây tham số hoạt động dông sét giả định 2.2 Tạo chuỗi số liệu mô hoạt động dông sét Ứng với dạng mây bảng 2, tạo chuỗi số liệu (chuỗi số liệu mô phỏng) dựa đặc trưng giả định dạng mây Chuỗi số liệu mô bao gồm đặc trưng sau: thời gian xẩy sét đánh, kinh độ vỹ độ vị trí sét đánh Đối với dạng mây có ổ dơng (mây dơng loại I) q trình mơ tiến hành sau: (1) Xác định vị trí ban đầu thời gian tia sét làm tâm đám mây phóng điện; (2) Căn vào diện tích bao phủ độ lệch chuẩn vị trí phóng điện giả định, sử dụng hàm phân bố chuẩn hoá (phân bố Gaussian) để tạo tập hợp ngẫu nhiên cú phóng điện phân bố quanh vị trí tâm đám mây phóng điện; (3) Căn vào tốc độ di chuyển bước thời gian giả định, di chuyển tọa độ tâm đám mây phóng điện thực lại bước (1) (2) Quá trình từ (1) đến (3) lặp lại hết thời gian kéo dài dông Kết thu chuỗi số liệu mơ vị trí thời gian phóng điện cho mây dơng loại I Đối với mây dơng loại II, q trình tạo chuỗi số liệu mơ tương tự mây dông loại I, nhiên mây dơng loại II có hai ổ dông di chuyển đồng thời khoảng thời gian số tham số giới hạn khác (bảng 2) Đối với mây dơng loại III, q trình tạo chuỗi số liệu mô tương tự mây dông loại I, nhiên tham số độ lệch chuẩn vị trí phóng điện theo kinh độ vỹ độ khơng nên ổ mây dơng có dạng hình elíp, kéo dài theo vỹ hướng Bảng Các tham số đầu vào lựa chọn để tạo chuỗi số liệu mô dạng mây dông I, II III Các tham số Mây dông loại I Mây dông loại II Mây dông loại III Số ổ dông Khoảng cách ổ dông (km) 40 Thời gian dông (phút) 45 75 120 Bước thời gian chia khoảng (phút) 15 15 20 Tổng số phóng điện 75 270 900 Diện tích bao phủ (km ) 110 250 1200 Độ lệch chuẩn vị trí phóng điện theo kinh độ (km) 1,5 2,5 5,0 Độ lệch chuẩn vị trí phóng điện theo vỹ độ (km) 1,5 2,5 10,0 Tốc độ di chuyển (km/h) 15,0 20,0 30,0 2.3 Phương pháp xác định số tham số hoạt động dông sét từ chuỗi số liệu mô Để xác định số tham số hoạt động dông sét theo số dạng mây dông từ chuỗi số liệu mô phỏng, áp dụng phương pháp phân tích tương quan khơng gian thời gian Finke [4] Phân bố sét mô tả hàm mật độ thực nghiệm phụ thuộc vào thời gian không gian: N f ( t, r ) = ∑ δ t- t i δ r - ri i ( ) ( ) (1) 135 Trong δ: hàm Dirac, N: tổng số sét đánh, ri : khoảng cách từ vị trí sét đánh thứ i đến gốc tọa độ ti : thời gian xẩy cú sét thứ i 12 ⎛ ( ρ - cτ )2 ⎞ ⎛ ρ ⎞ ⎜⎟ Brt ( τ,ρ ) ≈ N T ⎜ exp ⎟ ⎜ 4σ ⎟ ⎝ 4πσ cτ ⎠ ⎝ ⎠ (6) Hàm tự tương quan thực nghiệm B(τ, ρ), khoảng không gian thời gian, xác định sau: - Bán kính mây dơng (d) tính qua độ lệch chuẩn (d = 4*σ), với giá trị τ = 0, σ xác định qua biểu thức sau: B (τ , ρ ) = ∫ΔT dt ∫F d r f ( t + τ, r + ρ ) f ( t, r ) Brt ( τ,ρ ) ≈ N 2T 2σ (2) Tích phân tính khoảng khơng gian F khoảng thời gian ΔT = (T1, T2) τ: khoảng thời gian, ρ: véc tơ khoảng cách Thay hàm phân bố sét (1) vào phương trình (2) nhận được: ( ) ( ) N B ( τ, ρ ) = ∑ δ ⎡ τ − t i - t j ⎤ δ ⎡ρ - ri - rj ⎤ ⎦⎥ ⎢⎣ ⎦⎥ i, j ⎣⎢ (3) Hàm B(τ, ρ) biểu diễn phân bố mật độ khoảng không gian (ri – rj) khoảng thời gian (ti – tj) Để đơn giản hố, chúng tơi đưa vào hai hàm tích phân riêng hàm tương quan B(τ, ρ): - Tương quan bán kính-thời gian, tích phân B tất hướng (Φ): B rt ( τ, ρ ) = ρ ∫02π B ⎣⎡τ, ρ = (ρ, Φ )⎦⎤ dΦ (4) Trong đó, (ρ, Φ) thành phần cực véc tơ chuyển vị ρ - Tương quan không gian khoảng thời gian (τ1, τ2): τ B xy ( ξ, η ) = ρ ∫τ B ⎡⎣ τ, ρ = ( ξ, η )⎤⎦ dτ (5) Trong đó, (ξ, η) tọa độ đề ρ Phân tích hai hàm (4), (5) xác định số đặc tính phân bố sét di chuyển mây dông Qua số bước biến đổi, xác định công thức biểu diễn quan hệ tham số hoạt động dơng sét với hàm tương quan bán kính thời gian, hàm tương quan không gian: - Tốc độ di chuyển (c) mây dơng xác định qua biểu thức sau: 136 ( ) -1 ⎛ ρ + c2 τ ⎞ ρexp ⎜ ⎟ 4σ ⎠ ⎝ (7) - Tốc độ di chuyển (c) hướng di chuyển trung bình mây dơng xác định cx cy qua biểu thức sau: 2 -1 Bxy ( ξ, η ) = N Tδτ 4πσ × ( ) 2⎞ ⎛ ξ - c x τ0 ) + η - c y τ0 ⎟ ( ⎜ ×exp ⎜ ⎟ 4σ ⎜ ⎟ ⎝ ⎠ ( ) (8) Kết tính tốn phân tích 3.1 Kết tính tốn phân tích cho mây dông loại I Chuỗi số liệu mô mây dơng loại I, biểu diễn hình 1, giao diện phần mềm sở hệ thống định vị sét Việt Nam (phiên 1.2) Sự thay đổi khơng gian vị trí phóng điện tương ứng với thay đổi thời gian, độ tập trung vị trí phóng điện có dạng gần trịn (quy luật Gaussian) Sử dụng chuỗi số liệu mô mây dơng loại I này, chúng tơi tính tốn hàm tương quan (trình bày mục 2.3): tương quan bán kính-thời gian, tương quan khơng gian Kết tính tốn trình bày hình 2, Căn theo hai hàm tương quan xác định số tham số hoạt động dông sét mây dông loại I, bao gồm: tốc độ di chuyển trung tâm mây dơng: 15,566km/h (tính theo hàm tương quan bán kính - thời gian, lệch 0,566km/h so với tốc độ giả định ban đầu: 15,0km/h) 14,142km/h (tính theo hàm tương quan không gian, lệch 0,858km/h so với tốc độ giả định ban đầu: 15,0km/h), đường kính mây dơng: 6,197km (lệch 0,197km so với đường kính mây dơng mô phỏng: 6km), thời gian dông: 40,42 phút (lệch 4,58 phút so với thời gian giả định: 45 phút) Hình Biến đổi theo không gian thời gian vị trí sét đánh mây dơng loại I, theo chuỗi số liệu mơ Hình Tương quan bán kính-thời gian mây dông loại I, theo chuỗi số liệu mô Để thấy rõ kết tính tốn so sánh trên, xem xét hình hình Trong hình 2, ta thấy góc nghiêng vùng hàm tương quan bán kính - thời gian, bước thời gian trễ lớn (tương ứng với tốc độ di chuyển trung tâm đám mây) có giá trị xấp xỉ xung quanh giá trị 15,0km/h Sử dụng giá trị hàm tương quan bán kính - thời gian thời điểm thời gian trễ 0, xác định độ lệch chuẩn (1,549), qua xác định đường kính mây Hình Tương quan không gian mây dông loại I, theo chuỗi số liệu mô dông qua ước lượng (d = 4*độ lệch chuẩn, [4]), đường kính mây dơng xác định gần với giá trị mô Đỉnh cực đại hàm tương quan khơng gian (hình 3), vị trí (6, 6), có khoảng cách so với gốc (0, 0) 8,485km, khoảng cách tương ứng với thời gian di chuyển mây dông khoảng thời gian 0,6 giờ; khoảng giờ, mây dông di chuyển 14,142km, gần với giá trị tốc độ giả định 15,0km/h Xu dịch chuyển trung tâm hàm tương quan khơng gian 137 có xu từ trái qua phải, qua ta thấy hướng di chuyển mây dơng có xu từ tây qua đơng 3.2 Kết tính tốn phân tích cho mây dông loại II Chuỗi số liệu mô mây dông loại II, biểu diễn hình 4, cho thấy thay đổi khơng gian vị trí phóng điện tương ứng với thay đổi thời gian hai ổ dơng, độ tập trung vị trí phóng điện có dạng gần trịn hai ổ dông (quy luật Gaussian) Chuỗi số liệu mô mây dơng loại II chúng tơi sử dụng để tính hàm tương quan; kết tính tốn trình bày hình hình Hình Biến đổi theo khơng gian thời gian vị trí sét đánh mây dông loại II, theo chuỗi số liệu mơ Hình Tương quan bán kính-thời gian mây dông loại II, theo chuỗi số liệu mô 138 Hình Tương quan khơng gian mây dơng loại II, theo chuỗi số liệu mô Căn theo hai hàm tương quan xác định số tham số hoạt động dông sét mây dông loại II, bao gồm: tốc độ di chuyển trung tâm mây dơng: 22,724km/h (tính theo hàm tương quan bán kính - thời gian, lệch 2,724km/h so với tốc độ giả định: 20,0km/h) 19,799km/h (tính theo hàm tương quan không gian, lệch 0,201km/h so với tốc độ giả định: 20,0km/h), đường kính mây dơng: 10,238km (lệch 1,762km so với đường kính mây dơng mơ phỏng: 12,0km), khoảng cách hai ổ dông: 40km, thời gian dông 73,64 phút (lệch 1,36 phút so với thời gian giả định: 75 phút) 40km, giá trị khoảng cách hai ổ dông Đỉnh cực đại (B11+B22), hàm tương quan không gian (hình 6), vị trí (14, 14), có khoảng cách so với gốc (0, 0) 19,799km, khoảng giờ, mây dông di chuyển 19,799km, gần với giá trị tốc độ giả định 20km/h; khoảng cách đỉnh cực đại (B11+B22) B21 (B11+B22) B12 khoảng 40km, qua ta xác định khoảng cách hai ổ dông Xu dịch chuyển trung tâm hàm tương quan không gian (B11+B22) có xu từ trái qua phải, qua ta thấy hướng di chuyển mây dơng có xu từ tây qua đơng Xem hình 5, ta thấy góc nghiêng vùng hàm tương quan bán kính-thời gian (ở vùng B11+B22), bước thời gian trễ lớn (tương ứng tốc độ di chuyển trung tâm đám mây) có giá trị xấp xỉ xung quanh giá trị 20,0km/h Sử dụng giá trị hàm tương quan bán kính-thời gian thời điểm thời gian trễ 0, xác định độ lệch chuẩn (2,560), qua xác định đường kính mây dông gần với giá trị mô Khoảng cách hai đỉnh cực đại hàm tương quan bán kính-thời gian có khoảng cách 3.3 Kết tính tốn phân tích cho mây dơng loại III Chuỗi số liệu mô mây dông loại III, biểu diễn hình 7, cho thấy thay đổi khơng gian vị trí phóng điện tương ứng với thay đổi thời gian đám mây dông ổ dơng, độ tập trung vị trí phóng điện có dạng hình elíp (quy luật Gaussian) Chuỗi số liệu mô mây dông loại III sử dụng để tính hàm tương quan; kết tính tốn trình bày hình hình Hình Biến đổi theo khơng gian thời gian vị trí sét đánh mây dơng loại III, theo chuỗi số liệu mô Căn theo hai hàm tương quan xác định số tham số hoạt động dông sét mây dông loại III, bao gồm: tốc độ di chuyển trung tâm mây dơng: 30,328km/h (tính theo 139 hàm tương quan bán kính - thời gian, lệch 0,328km/h so với tốc độ giả định: 30km/h) 28,284km/h (tính theo hàm tương quan không gian, lệch 1,716km/h so với tốc độ giả định: 30km/h), đường kính mây dơng 40,251km (lệch 0,251km so với đường kính mây dơng mơ phỏng: 40km), thời gian dông 117,52 phút (lệch 2,48 phút so với thời gian giả định: 120 phút) Hình Tương quan bán kính-thời gian mây dông loại III, theo chuỗi số liệu mô Hình Tương quan khơng gian mây dơng loại III, theo chuỗi số liệu mơ Xem hình 8, ta thấy góc nghiêng vùng hàm tương quan bán kính-thời gian, bước thời gian trễ lớn (tương ứng tốc độ di chuyển trung tâm đám mây) có giá trị xấp xỉ xung quanh giá trị 30km/h Sử dụng giá trị hàm tương quan bán kính-thời gian thời điểm thời gian trễ 0, xác định độ lệch chuẩn (10,063), qua xác định đường kính mây dơng gần với giá trị mô số tham số hoạt động dông sét từ chuỗi số liệu mô theo phương pháp Finke Từ kết trình bày trên, rút số kết luận sau: Hàm tương quan khơng gian (hình 9), vị trí (20, 20), có khoảng cách so với gốc (0, 0) 28,284km, khoảng giờ, mây dông di chuyển 28,284km, gần với giá trị tốc độ giả định 30km/h Xu dịch chuyển trung tâm hàm tương quan khơng gian có xu từ trái qua phải, qua ta thấy hướng di chuyển mây dơng có xu từ tây qua đông Kết luận Mục tiêu báo xác định số tham số hoạt động dông sét từ chuỗi số liệu định vị phóng điện mơ Để thực mục tiêu trên, tiến hành xây dựng chuỗi số liệu định vị phóng điện mô từ thông số hoạt động dông sét giả định cho ba dạng mây dông, tương tự cách phân chia Peckham Sau đó, xác định 140 - Các tham số hoạt động dông sét xác định từ chuỗi số liệu mô cho dạng mây dông khác (I, II, III), thông qua phương pháp Finke, chênh lệch không nhiều so với tham số giả định ban đầu sử dụng cho mô Chênh lệch tốc độ dạng mây I: dVI = 0,566 km/h, dạng mây II: dVII = 2,724km/h, dạng mây III: dVIII = 0,328km/h Chênh lệch bán kính mây dơng: dRI = 0,197km, dRII = 1,762km, dRIII = 0,251km) Chênh lệch thời gian dông: dTI = 4,58 phút, dTII = 1,36 phút, dTIII = 2,48 phút) Từ hàm tương quan không gian hàm tương quan bán kính thời gian xác định số ổ dông, khoảng cách ổ dông, giá trị giống với giá trị giả định ban đầu dùng cho mô - Kết báo sở để áp dụng phương pháp Finke, cho việc tính tốn, xác định số tham số hoạt động dông sét như: tốc độ di chuyển mây dơng, bán kính mây dơng, thời gian dông, số ổ dông, khoảng cách ổ dơng, từ chuỗi số liệu định vị phóng điện quan sát thực tế, độ xác số liệu bảo đảm sử dụng tham số cho nghiên cứu cảnh báo hoạt động dơng sét Lời cảm ơn: Cơng trình nhận hỗ trợ kinh phí Nhiệm vụ Khoa học Công nghệ cấp thiết thực địa phương năm 2010: “Nghiên cứu đề xuất giải pháp phòng chống sét địa bàn tỉnh Quảng Nam”, đề tài nghiên cứu bản: “Nghiên cứu khả dự báo ngày bắt đầu gió mùa mùa hè - mùa mưa khu vực Nam Bộ”, thuộc Quỹ phát triển Khoa học Công nghệ Quốc gia (NAFOSTED) TÀI LIỆU DẪN [1] Nguyễn Xuân Anh, Lê Việt Huy, Hoàng Hải Sơn, 2010: Về số kết nghiên cứu dông sét phòng chống sét gần Viện Vật lý Địa cầu Tuyển tập báo cáo Hội nghị Khoa học kỷ niệm 35 năm Viện Khoa học Công nghệ Việt Nam 1975-2010, Tiểu ban Khoa học Trái Đất, 295-299 [2] Nguyễn Xuân Anh, Lê Việt Huy, Hoàng Hải Sơn, 2008: Nghiên cứu hoạt động dông sét giải pháp phịng chống sét Việt Nam Tuyển tập Các cơng trình nghiên cứu Vật lý Địa cầu, 289-306 [3] Nguyễn Xuân Anh, Lê Việt Huy, Hoàng Hải Sơn, 2007: Một số kết nghiên cứu hoạt động dông sét qua mạng trạm định vị sét Việt Nam Tuyển tập báo cáo Hội nghị Khoa học kỹ thuật Địa Vật lý lần thứ 5, 11-20 [4] Finke U., 1999: Space-Time Correlations of Lightning Distributions Monthly Weather Review, Volume 127, pp 1850-1861 [5] Hagen M., et al 1999: Motion charateristcs of thunderstorms in southern Germany Meteorol Appl 6, 227-239 [6] Mark S M and Charles E K., 2005: Spatial and Temporal Patterns of Thunderstorm Events that Produce Cloud-to-Ground Lightning in the Interior Southeastern United States Monthly Weather Review, Volume 127, pp 1417-1433 [7] Peckham D W., et al, 1984: Lightning Phenomenology in the Tampa Bay Area J Geophys Res., 89:11,789-11, 805 [8] Hồng Hải Sơn, 2006: Nghiên cứu hoạt động dơng sét thiết bị định vị phóng điện Luận văn Thạc sỹ Khoa học, Trường Đại học Khoa học Tự nhiên, ĐHQG Hà Nội, 99tr SUMMARY Determination of lightning activity parameters from lightning simulation data In this paper, lightning activity parameters could be listed, including: moving speed of thundercloud, radius of thundercloud, storm duration, number of cell thundercloud, and distance between the cell of thundercloud for three difference types of thundercloud is determined using lightning simulation data The first of all, thundercloud is divided into three different forms (similarity the way division of Peckham) which is based on their presumed parameters Next, determination of lightning simulation data based on the parameters assumed above Finally, using the method of correlation analysis space and time of Finke and lightning simulation data of three types of thundercloud above to determine the radial-temporal correlation and space correlation function Besides, two correlation functions was used to identify the lightning activity parameters As a result, values of lightning activity parameters and assumed initial value is not much difference The difference in speed of cloud form I (dVI = 0.566 km/h), cloud form II (dVII = 2.724 km/h), and cloud form III (dVIII = 0.328 km/h) The difference in radius of cloud (dRI = 0.197 km, dRII = 1.762 km, dRIII = 0.251 km) The difference in storm duration (dTI = 4.42 minutes, dTII = 1.36 minutes, dTIII = 2.48 minutes) The number of cell thundercloud and the distance between the cell of thundercloud identified with the same initial value assumption 141 ... 20,0 30,0 2.3 Phương pháp xác định số tham số hoạt động dông sét từ chuỗi số liệu mô Để xác định số tham số hoạt động dông sét theo số dạng mây dông từ chuỗi số liệu mô phỏng, chúng tơi áp dụng... báo xác định số tham số hoạt động dông sét từ chuỗi số liệu định vị phóng điện mơ Để thực mục tiêu trên, tiến hành xây dựng chuỗi số liệu định vị phóng điện mô từ thông số hoạt động dông sét. .. Peckham, chia mây dông thành ba loại với số tham số tương tự Peckham bổ sung thêm số tham số khác Đây tham số hoạt động dông sét giả định 2.2 Tạo chuỗi số liệu mô hoạt động dông sét Ứng với dạng