Page 128 - Tạp chí Giao thông Vận Tải - Số Tết Dương Lịch

Số 12/2024 (748) KHOA HỌC CÔNG NGHỆ
Tập
4
6
Tập 64
TT Nội dung/Yếu tố Tham số
I.2 Bước thời gian 30 giây
Trong vòng 15
ngày (01/2020,
I.3 Thời gian bắt đầu
01/2022,
4/2022, 7/2022,
10/2022,
02/2023,
I.4 Thời gian kết thúc
6/2023 và Hình 3.4: Trường tổng hợp quỹ đạo di chuyển của các thành phần
12/2023)
II Các mô hình được lựa chọn mô phỏng hạt dầu tràn giai đoạn từ 26/02/2023 đến 11/3/2023
Vào tháng 02/2023 xuất hiện vệt dầu không rõ nguồn
II.1 Mô hình thủy động lực gốc tại khu vực bến cảng Lạch Huyện. Sau đó, dưới tác
động của triều xuống, vệt dầu có xu hướng chảy ra phía
II.3 Mô hình truyền tải dầu tràn
ngoài khu vực cầu số 3, 4. Sang các ngày tiếp theo, vệt dầu
Bảng 3.2. Thông tin lựa chọn của mô hình thủy động lực có xu hướng chảy theo hướng Bắc Tây Bắc tiếp cận khu vực
cảng nhà máy Vinfast. Sau đó, vệt dầu chảy theo hướng
TT Nội dung Yếu tố Tham số chính là Đông Nam và Tây Bắc.
- Theo thời gian là tính toán bậc cao; 3.3. Kết quả mô phỏng tràn dầu bằng phần mềm
Phương - Theo không gian là tính toán bậc cao; tự xây dựng
trình nước - Bước thời gian nhỏ nhất là 0,01s;
Giải pháp nông - Bước thời gian lớn nhất là 30,0s; 3.3.1. Chuyển đổi miền tính và lưới tính
1 công - Điều kiện hội tụ (CFL) là 0,85. Sử dụng công cụ của phần mềm MIKE để phục vụ
nghệ công việc nghiên cứu thuận lợi. Công cụ chuyển đổi này
Phương - Bước thời gian nhỏ nhất là 0,01s;
trình truyền sẽ thực hiện định dạng tệp có đuôi “.dfsu” sang định dạng
tải khuếch - Bước thời gian lớn nhất là 30,0s; tệp có đuôi “.dfs2”.
- Điều kiện hội tụ (CFL) là 0,85.
tán
- Độ sâu khô là 0,005 m;
Điều kiện tiêu chuẩn khô
2 - Độ sâu ngập là 0,010 m;
và ướt
- Độ sâu ướt là 0,025 m.
Nhớt rối Theo công thức Smagorisky với giá trị
ngang hàng số là 0,28.
3 Nhớt rối 2
Tham số rối - Tham số nhỏ nhất là 1,0e-3 m /s;
2
- Tham số lớn nhất là 1,0e+10 m /s.
Theo số Manning với giá trị hằng số là
4 Ma sát đáy
31,5 m /s.
1/3
5 Lực Coriolis - Biến đổi theo không gian của miền tính.
- Biến đổi theo thời gian và đồng nhất Hình 3.5: Miền tính, lưới tính và trường địa hình khu vực nghiên cứu
theo không gian của miền tính; 3.3.2. Chuyển đổi dữ liệu thủy động lực
- Hệ số ma sát nhỏ nhất 0,001 với gió Tương tự như chuyển đổi tệp địa hình, tệp chứa dữ
6 Lực gió bề mặt
nhỏ hơn 2,5 m/s; liệu về kết quả tính toán mô phỏng thủy động lực cũng
- Hệ số ma sát lớn nhất nhất 0,0032 với
gió bằng hoặc lớn hơn 25 m/s. được tiến hành. Các kết quả được thể hiện trình bày như
hình dưới đây.
- Mực nước là 0,0 m;
7 Điều kiện ban đầu
- Dòng chảy (u, v) là 0,0 m/s.
- Biên mực
nước ở các
biên hở; - Với biên sông là Trạm Thủy văn Cửa
- Biến đổi Cấm và Trạm Thủy văn Do Nghi;
Điều kiện
8 ban đầu là - Với biên biển là dự báo mực nước
biên
3.600s; theo hằng số điều hòa của phần mềm
- Có sự tham MIKE sẵn có.
gia của lực
Coriolis; gió.
Hình 3.6: Trường độ sâu khu vực nghiên cứu
- Theo số liệu trường từng giờ là mực 3.3.3. Đặc trưng trường hợp tràn dầu thực tế tháng
nước, dòng chảy;
10 Đầu ra 02/2023
- Theo số liệu điểm trạm cụ thể từng
giờ là mực nước, dòng chảy.
3.2.2. Kết quả tính toán mô hình tràn dầu MIKE 21/3 OS
Từ kết quả tính toán mô phỏng thủy động lực đã được
trình bày, kết quả tính toán mô phỏng tràn dầu với kịch
bản tràn dầu cụ thể xảy ra vào tháng 02/2023 ở khu vực
nghiên cứu. Hình 3.4 thể hiện tổng hợp quỹ đạo di chuyển
của các thành phần hạt dầu từ 26/02/2023 đến 11/3/2023.
Dưới tác động của yếu tố gió, dòng chảy và sóng, vệt dầu
chủ yếu tập trung tại khu vực luồng Nam Triệu, luồng Lạch
Huyện và khu vực vịnh Hải Phòng. Hình 3.7: Vệt dầu tràn từ mô hình MIKE và mô hình xây dựng ngày
26/02/2023 (trái) và ngày 04/3/2023 (phải)
127