Bài viết sử dụng mô hình địa cơ nhằm phân tích, dự báo các yếu tố ảnh hưởng đến dịch chuyển biến dạng đất đá trong mối tương quan với góc dốc vỉa, độ sâu và hệ thống khai thác.
66 Tạp chí Khoa học Kỹ thuật Mỏ - Địa chất Tập 58, Kỳ (2017) 66-72 Nghiên cứu ảnh hưởng góc dốc vỉa độ sâu khai thác đến dịch chuyển biến dạng bề mặt mơ hình địa Phạm Văn Chung 1,*, Nguyễn Quốc Long 1, Vương Trọng Kha 1, Nguyễn Quang Phích Khoa Khoa Trắc địa - Bản đồ Quản lý đất đai, Trường Đại học Mỏ - Địa chất, Việt Nam Xây dựng, Trường Đại học Mỏ - Địa chất, Việt Nam THƠNG TIN BÀI BÁO TĨM TẮT Q trình: Nhận 15/3/2017 Chấp nhận 25/7/2017 Đăng online 31/8/2017 Bài báo sử dụng mơ hình địa nhằm phân tích, dự báo yếu tố ảnh hưởng đến dịch chuyển biến dạng đất đá mối tương quan với góc dốc vỉa, độ sâu hệ thống khai thác Mô hình số địa cho phép nghiên cứu tốn có tính đến đặc điểm phân lớp (khơng đồng nhất, không đẳng hướng khối đá) Các quy luật nhận phân bố ứng suất, biến dạng, phát triển vùng phá hủy, cường độ chúng phụ thuộc vào đặc điểm cấu trúc, quy luật thứ tự xếp lớp, hình dạng tiết diện ngang cơng trình mỏ Từ kết thu nhận phần mềm FLAC 2D cho phép dự báo quy luật dịch chuyển xảy khai thác than hầm lò Việt Nam Từ khóa: Dự báo dịch chuyển biến dạng Khai thác hầm lò Mơ hình địa © 2017 Trường Đại học Mỏ - Địa chất Tất quyền bảo đảm Mở đầu Công tác khai thác mỏ gây tác động lớn đến bề mặt đất lấy phần vật chất từ phần vỏ trái đất, dẫn đến thay đổi trạng thái vật chất trạng thái cấu trúc ban đầu khối đá làm cân ứng suất nén ban đầu địa tầng chứa than Hậu khối đá tự nhiên rơi vào trạng thái “mất cân bằng” ứng suất “tự điều chỉnh” để xác lập trạng thái cân mới, trạng thái cân học có ý nghĩa đặc biệt (Phạm Văn Chung, 2010) Quá trình diễn dạng phân bố lại trạng thái ứng suất kéo theo tượng dịch chuyển biến dạng đất đá có tính chất đa dạng, phức tạp phụ thuộc vào nhiều yếu tố địa chất, khai thác, lý đá… _ *Tác giả liên hệ E-mail: phamvanchung@humg.edu.vn Trong trình biến đổi trên, khối đất đá tiến đến trạng thái ổn định hoàn toàn, trạng thái ổn định mức độ khác Trạng thái ổn định cân ứng suất dịch chuyển không đồng gây biến dạng Khoảng trống đất đá than khấu làm cân lực nên lớp đất đá vách (nóc lò) uốn võng sập đổ Khoảng sập đổ lớn hay nhỏ tùy thuộc vào kích thước khoảng trống khai thác Sự sập đổ lớp đất đá vách tượng dịch chuyển biến dạng đất đá mỏ Hiện tượng giảm dần đất đá lấp đầy khoảng trống khai thác đất đá vùng trụ vỉa ổn định lan truyền lên phía chí lên tới mặt đất đất đá yếu độ sâu khai thác không lớn Hậu lan truyền gây bồn dịch chuyển mặt đất xuất tầng bậc, ảnh hưởng đến cơng trình tùy theo cường độ biến dạng thể qua tượng Phạm Văn Chung nnk./Tạp chí Khoa học Kỹ thuật Mỏ - Địa chất 58 (4), 66-72 dịch chuyển phá hủy (nén vỡ đá, sạt lở, sập lở…) Cho đến nay, lĩnh vực học đá, khai thác mỏ, có nhiều phương pháp ứng dụng để nghiên cứu dịch chuyển biến dạng đất đá có nhiều kết nghiên cứu công bố giới Mơ hình nghiên cứu dự báo chiều cao vùng phá hủy, tách lớp nhà khoa học nghiên cứu, điển tác giả (Kratzsch, 1983; Mindlin, 1939) Tuy nhiên nước ta, lĩnh vực khai thác mỏ vấn đến ý Từ năm 2011 đến có số cơng trình nghiên cứu dịch chuyển biến dạng mơ hình vật liệu tương đương cơng bố nhiều tài liệu khác (Nguyễn Anh Tuấn nnk, 2011; Nguyễn Anh Tuấn nnk, 2012; Nguyễn Văn Sỹ nnk, 2012; Nguyễn Tam Sơn, 2012; Trần Trung Anh, 2012) Tuy nhiên, qua tài liệu cho thấy việc nghiên cứu nhiều công sức, thực mơ hình nhỏ, nên hệ số tương đương lớn Nghiên cứu mang tính lý thuyết, khơng kể đến cơng trình cơng bố nhóm thực đề tài này, để dự báo lún sụt số chuyên gia thực (Nguyễn Anh Tuấn nnk, 2011) sử dụng chương trình Phase2 phân tích lún sụt trình biến đổi học khai thác hỗn hợp hầm lò, lộ thiên Trong phạm vi viết tác giả phân tích, dự báo kết dịch chuyển, biến dạng mơ hình khai thác lò chợ dài với góc dốc vỉa khác mơ hình khai thác lò chợ dài mức độ sâu khác thông qua phần mềm FLAC 2D Phần mềm FLAC 2D FLAC (Fast Lagrangian Analysis of Continua) phát triển trì ITASCA Các phần mềm FLAC dựa phương pháp sai phân hữu hạn hàm hiện, qua cho phép mơ tính chất vật liệu đất, đá kết cấu xây dựng để giải loạt vấn đề phức tạp địa kỹ thuật, cơng trình dân dụng khai thác mỏ Cơ sở để giải toán tổng quát trình biến đổi học khối đá phương pháp giải tích phương pháp số lý thuyết động lực học Điều kiện cân động học trạng thái đơn trục hệ có bậc tự có dạng: PM+PC+PK=P(t) (1) Trong đó: PM lực quán tính; PC lực cản nhớt; PK lực kháng đàn hồi; P(t)- véc tơ ngoại 67 lực cưỡng có dạng hàm số thời gian Ngoại lực P(t) tác dụng lên vật thể bao gồm loại lực khác gây vật thể nội lực dạng: lực kháng đàn hồi, có hướng tác dụng ngược với hướng chuyển vị, lực cản nhớt (tắt dần), ngược với tốc độ chuyển vị, lực quán tính, hay lực cản gia tốc Như phương trình (1) phương trình chuyển động biểu diễn cân tất lực Theo nguyên lý Dalambe (khối lượng m gây lực quán tính, tỷ lệ với gia tốc hướng ngược với gia tốc), biểu diễn ngoại lực dạng: (2) 𝑃(𝑡) = −𝑚𝑢𝑔̈ (𝑡) Trong (2): 𝑢𝑔̈ (𝑡)là gia tốc dịch chuyển khối đá Với 𝑃𝑀 = 𝑚 𝑢̈ , 𝑃𝑀 = 𝑚 𝑢̇ , 𝑃𝐾 = 𝑘 𝑢, phương trình (2) có dạng: (3) 𝑚 𝑢̈ + 𝑚 𝑢̇ + 𝑘 𝑢 = −𝑚𝑢𝑔̈ (𝑡) Cơ sở lý thuyết phương trình tổng quát để giải tốn nghiên cứu q trình biến đổi học khối đá phương pháp giải tích tham khảo (Mindlin, 1939), để giải phương pháp số tham khảo (Mindlin, 1948) Xuất phát từ mục tiêu tốn xây dựng mơ hình địa dự báo tai biến địa chất - kỹ thuật chủ yếu phương pháp số sử dụng phần mềm FLAC2D, nên trình bày khái quát sở lý thuyết phương pháp Trong trường hợp tổng quát để giải toán phương pháp số, phương trình động lực học hệ học viết dạng ma trận là: 𝜕𝑢2 𝜕𝑢 [𝑀] { } + [𝐶] { } + [𝐾]{𝑢} = {𝑃(𝑡)} (4) 𝜕𝑡 𝜕𝑡 Trong đó: [M], [C], [K] ma trận khối lượng, ma trận cản ma trận độ cứng; {u}, {∂u / ∂t}, {∂u2 / ∂t2} véc tơ dịch chuyển (chuyển vị), véc tơ tốc độ dịch chuyển véc tơ gia tốc dịch chuyển {P(t)} véc tơ lực (tải trọng) t thời gian Trường hợp tốn tĩnh học, {∂u/ ∂t} = nên phương trình (4) có dạng đơn giản: [𝐾]{𝑢} = {𝑃} (5) Khi gặp toán động lực học, thiết phải giải phương trình vi phân (5) xét đến tất thành phần lực cản (cản nhớt) lực khối lượng (quán tính) 68 Phạm Văn Chung nnk./Tạp chí Khoa học Kỹ thuật Mỏ - Địa chất 58 (4), 66-72 Với tải trọng động đất gây ra, viết (5) dạng : 𝜕𝑢2 }+ 𝜕𝑡 2 𝜕𝑢 𝑔𝑥 [𝑀] { 𝜕𝑢 [𝐶] { } + [𝐾]{𝑢} = 𝜕𝑡 𝜕𝑢2 𝑔𝑦 𝜕𝑢2 𝑔𝑧 (6) −[𝑀𝑥 ] { 𝜕𝑡 } − [𝑀𝑦 ] { 𝜕𝑡 } − { 𝜕𝑡 } Trong đó: [Mx] - Là ma trận cột khối lượng theo phương x [My]- Ma trận cột khối lượng theo phương y [Mz]- Ma trận cột khối lượng theo phương z 𝜕𝑢2𝑔𝑥 𝜕𝑢2𝑔𝑦 𝜕𝑢2𝑔𝑧 (𝑡), (𝑡), = 𝑢 ̈ = 𝑢 ̈ { } { } { } 𝑔𝑥 𝑔𝑦 𝜕𝑡 𝜕𝑡 𝜕𝑡 = 𝑢𝑔𝑧 ̈ (𝑡) Tương ứng với thành phần nằm ngang (x), thẳng đứng (y) nằm ngang (z) gia tốc có tác động động đất Phân tích dịch chuyển biến dạng với mơ hình khối đá đồng Những kết xây dựng mô hình số địa nghiên cứu dịch chuyển biến dạng đất đá kết hợp với việc ứng dụng phần mềm FLAC 2D đưa đến khả áp dụng phát triển công tác dự báo dịch chuyển biến dạng khai thác hầm lò Mỏ than Mạo Khê phía Nam bể than Quảng Ninh trục đường quốc lộ 18 từ Hà Nội Hòn Gai cách thị trấn Đông Triều 10 km phía Đơng Vị trí mỏ than Mạo Khê thuận tiện cho giao thông đường thuỷ, đường sắt đường bộ, cầu nối Hà Nội thành phố Hạ Long Mỏ than Mạo Khê có trữ lượng than lớn, điều kiện kiến tạo vỉa không phức tạp, khu vực tập trung đơng dân cư có nhiều tiềm phát triển các mặt kinh tế, văn hóa, xã hội Hiện mỏ Mạo Khê tập trung khai thác than vỉa: vỉa 3, 4, 5, 6, 7, 8, 9, vỉa 9b Mở vỉa vào ruộng than lò dọc vỉa đá vào cánh vỉa Phân chia ruộng than thành khu vực thuận tiện cho việc vận tải than điều khiển áp lực mỏ Sử dụng công nghệ khai thác cột dài theo phương, chống lò gỗ giá thuỷ lực Điều khiển đá vách phá hỏa toàn phần Các điều kiện địa chất nằm vỉa 9b thể Bảng Thực tiễn cho thấy, có nhiều yếu tố ảnh hưởng đến dịch chuyển, biến dạng đất đá, báo tính hai trường hợp áp dụng vỉa 9b mỏ than Mạo Khê: + Khai thác lò chợ dài với góc dốc vỉa thay đổi; + Khai thác lò chợ dài mức độ sâu khác Bảng Điều kiện địa chất vỉa 9b STT Vỉa than Các thông số lò Đơn vị chợ 9b Mức khai thác m -80 -25 Chiều dày vỉa m 2.5 Góc dốc vỉa độ 25-27 Chiều dày đất phủ m 10 Chiều dài lò chợ theo m 110 hướng dốc Chiều dài lò chợ theo m 600 phương Chiều sâu trung bình m 380 - 400 lò chợ 3.1 Mơ hình khai thác lò chợ dài với góc dốc vỉa thay đổi Mặt cắt địa chất tuyến IV mỏ than Mạo Khê, địa tầng đất đá bao gồm lớp đá cát kết, bột kết, sét kết than nằm nghiêng với tham số học thể Bảng Miền khảo sát có kích thước 200m x 100m (rộng x cao) Vì nhằm rút nhận định mang tính chất định tính, nên để đơn giản chúng tơi chọn mơ hình cho khối đá, với giả thiết coi môi trường đất đá mỏ đồng đẳng hướng Bảng Các tham số học sử dụng cho mơ hình Mật độ Mơ đun Mơ đun Góc ma Lực dính ρ nén thể trượt G sát kết C (k/m3) tích K (Pa) (Pa) φ (Pa) 9 2200 0,788.10 0,484.10 30 2,51.106 Lò khai thác khảo sát với thông số: chiều dài 80 m, chiều cao m Hệ thống khai thác lò chợ dài, phá hỏa tồn phần Các lò chợ khảo sát có góc dốc thay đổi Do tiến hành cố định điểm lò chợ, nên thay đổi góc dốc vỉa, điểm ranh giới lò chợ gần mặt đất Trên hình thể tập hợp kết phân tích độ lún mặt đất, biểu thị qua biến hình lưới sai phân ảnh hưởng khai thác vỉa than quặng với góc dốc vỉa khác Phạm Văn Chung nnk./Tạp chí Khoa học Kỹ thuật Mỏ - Địa chất 58 (4), 66-72 Hình Sơ đồ tốn khai thác lò chợ dài 69 Từ kết thực nghiệm cho phép rút kết luận: - Khi độ sâu khai thác gần mặt đất, mức độ lún sụt tăng; - Khi góc dốc vỉa tăng hình dạng bề mặt bồn dịch chuyển dần tính đối xứng; điểm lún cực đại chuyển dịch phía ranh giới lò chợ đầu thân quặng - Vùng phá hủy đất đá phát triển chủ yếu phía lò chợ, có chiều cao hình dạng phụ thuộc vào nằm vỉa khoảng cách tương đối khu vực khoảng trống khai thác đến mặt đất; Trong trường hợp này, bề mặt đất xuất vùng phá hủy lân cận lộ vỉa trình lún; biến dạng tổng thể gây phá hủy kéo, xuất vùng biên miền nghiên cứu Vị trí quy mô hai vùng thay đổi tùy thuộc độ dốc vỉa (thân quặng) khoảng cách vùng khai thác đến mặt đất 3.2 Mơ hình khai thác lò chợ dài với độ sâu khác Cũng với liệu vào trường hợp 1, trường hợp này, xây dựng mơ hình khai thác lò chợ dài, nằm ngang có chiều dài theo dốc 80m, độ sâu khác từ 45m đến 300m, tính từ mặt đất đến hay vách lò khai thác Các kết phân tích cho phép nhận đường độ lún, tương ứng với độ sâu Hình Bảng Số liệu độ dịch chuyển độ lún sụt lớn Hình a) Sự phụ thuộc biên độ lún hình dạng bồn dịch chuyển; b) Sự phát triển vùng phá hủy (màu tím) vào góc dốc lò chợ Độ sâu khai Dịch chuyển Độ lún lớn thác (m) lớn (m) (m) 45 1.786 0.390 75 1.815 0.270 90 1.824 0.245 105 1.844 0.230 120 1.865 0.220 135 1.883 0.210 150 1.907 0.204 165 1.922 0.203 180 1.940 0.212 210 1.957 0.216 300 2.017 0.219 Tính chất học prop dens 2200 bulk 0.788e9 shear 0.484e9 fric 30 coh 2.51e6 70 Phạm Văn Chung nnk./Tạp chí Khoa học Kỹ thuật Mỏ - Địa chất 58 (4), 66-72 Hình Đường cong lún chịu ảnh hưởng độ sâu khai thác, với mơ hình khối đá đồng Hình Sự phát triển vùng phá hủy xung quanh lò khai thác, phụ thuộc độ sâu khai thác, với mơ hình khối đá đồng Phạm Văn Chung nnk./Tạp chí Khoa học Kỹ thuật Mỏ - Địa chất 58 (4), 66-72 Kết phân tích cho thấy: Lún sụt mặt đất ban đầu giảm dần theo độ sâu kể từ độ sâu định lại tăng theo độ sâu; Khi khai thác gần mặt đất đường cong lún có giá trị biến thiên lớn, nên thường có dạng nhọn Khai thác xuống sâu, vùng lún sụt mặt đất phát triển rộng thêm nên độ cong giảm dần, khả gây nghiêng, trượt, sụp đổ cơng trình kiên trúc giảm Trên Bảng số liệu độ dịch chuyển lớn điểm lò khai thác độ lún lớn mặt đất Từ Bảng thấy rằng, khai thác xuống sâu, độ dịch chuyển lớn tăng Với liệu đầu vào mơ hình này, độ lún giảm dần chiều sâu khai thác lớn đạt cực tiểu độ sâu khoảng 165 m Khai thác xuống sâu độ lún có gia tăng khơng đáng kể Trên Hình hình ảnh xuất vùng phá hủy khối đá xung quanh khu vực lò khai thác, mức khai thác sâu khác Từ kết mô thu nhận thấy quy luật lún sụt phát triển vùng phá hủy, khai thác độ sâu khác nhau: Ở độ sâu 45 m, 75 m: vùng phá hủy kéo tập trung lò khai thác lan đến mặt đất; phía lò xuất vùng phá hủy kéo, phủ kín phần đá trụ vỉa; sát đỉnh hai đầu lò chợ xuất vùng phá hủy cắt Trong trường hợp khai thác khơng gây tai biến lún, mà gây tai biến trụt lở Ở độ sâu 105 m chiều cao vùng phá hủy kéo cắt lò tăng lên, phá hủy kéo lò giảm dần, thay vào xuất vị trí khối đá bị phá hủy bề mặt với xu lò khai thác sâu vùng phá hủy bề mặt có khoảng cách xa Ở độ sâu 120 m 135 m, vùng phá hủy kéo phía lò đặc biệt phá hủy cắt tăng, nhiên không xuất phá hủy mặt đất Kết luận Từ phân tích kết thực nghiệm nêu cho thấy sử dụng mơ hình địa để dự báo định tính quy luật dịch chuyển, biến dạng đất đá mơ hình tốn học sơ đồ đường cong lún (Kratzsch, 1983; Whittaker Reddish, 1989) 71 Mơ hình địa cho phép giải tốn có tính đến đặc điểm biến động điều kiện địa chất - khai thác mỏ góc dốc vỉa, độ sâu khai thác Các quy luật nhận phân bố ứng suất, biến dạng, phát triển vùng phá hủy, cường độ chúng phụ thuộc rõ vào đặc điểm Nói chung, biến động điều kiện địa chất tính chất lý đá nên khơng thể có nhận định tổng qt quy luật biến đổi học Để nhận hình ảnh xác quy luật biến đổi học khối đá với cấu trúc địa chất phức tạp, thiết phải phân tích cụ thể cho trường hợp, mơ hình xác lập cho khối đá tương ứng với cấu trúc phân lớp, tham số học lớp Với phương thức phân tích khả phát triển “điều kiện địa chất” cụ thể dẫn đến “sự cố, tai nạn”, thơng qua xác định loại “tai biến địa chất” ngẫu nhiên xuất yếu tố nhân sinh Tài liệu tham khảo Phạm Văn Chung, 2010 Nghiên cứu xác định thông số dịch chuyển biến dạng bề mặt đất điều kiện địa chất đặc biệt khai thác hầm lò bể than Quảng Ninh, Báo cáo đề tài Bộ Công Thương Kratzsch, H., 1983 Mining Subsidence Engineering, (Springer-Verlag, Berlin Heidelberg New York) Mindlin, R D., 1939 Stress distribution around a tunnel Transactions American Society of Civil Engineers, Proceedings of the American Society of Civil Engineers, 2082 Mindlin, R D., 1948 Stress distribution around a hole near the edge of a plate under tension Proceedings of the Society for Experimental Stress Analysis 5(2), 56-68 Singh, M M., 1985 “Review of Coal Mine Subsidence Control Measures,” Transctions SME-AIME 278, Littleton, CO, 1988-1992 Singh M.M Chapter 10.6 Mine Subsidence Nguyễn Anh Tuấn, Lê Đức Nguyên, Nguyễn Văn Sỹ, 2011 Đánh giá kết nghiên cứu ảnh hưởng dịch động sơ đồ khai thác hỗn hợp hầm lò-lộ thiên mơ hình vật liệu 72 Phạm Văn Chung nnk./Tạp chí Khoa học Kỹ thuật Mỏ - Địa chất 58 (4), 66-72 tương đương Thông tin Khoa học công nghệ Mỏ 12 Nguyễn Anh Tuấn, Lê Đức Nguyên, Lê Quang Phục, 2011 Nghiên cứu ảnh hưởng dịch động sơ đồ khai thác hỗn hợp hầm lò-lộ thiên mơ hình số Thông tin Khoa học công nghệ Mỏ 11 Nguyễn Anh Tuấn, Đào Hồng Quảng, Lê Đức Nguyên, 2012 Nghiên cứu dịch chuyển biến dạng đất đá mỏ khai thác than hầm lò vùng Quảng Ninh mơ hình vật liệu tương đương Thông tin khoa học công nghệ mỏ 10 18-23 Nguyễn Văn Sỹ, Trần Văn Thanh, Lê Tiến Dũng, 2012 Xác định thông số dịch chuyển đất đá ảnh hưởng khai thác hầm lò đến bề mặt khống sản than Mơng Dương-Khe Chàm mơ hình tương đương Tuyển tập báo cáo Hội nghị Khoa học lần thứ 20, Đại học Mỏ - Địa chất, Hà Nội, 15/11/2012 Nguyễn Tam Sơn, 2012 Báo cáo tổng kết đề tài “Nghiên cứu xác định thông số dịch chuyển, biến dạng đất đá khai thác vỉa dày phương pháp hầm lò mơ hình vật liệu tương đương Viện KHCN mỏ 2012 Trần Trung Anh, 2012 Quan trắc dịch chuyển đất đá biến dạng bề mặt mơ hình vật liệu tương đương cơng nghệ đo ảnh Tạp chí KHKT Mỏ- Địa chất 39, 7/2012, (Chuyên đề Ttrắc địa mỏ), 5-10 ABSTRACT Research on the influence of the dip and the depth of excavation seam on surface movement and deformation using geomacherical digital models Chung Van Pham 1,*, Long Quoc Nguyen 1, Kha Trong Vuong 1, Phich Quang Nguyen Faculty of Geomatics and Land Administration, Hanoi University of Mining and Geology, Vietnam 2Faculty of Civil Engineering, Hanoi University of Mining and Geology, Vietnam Geomachenical models are able to solve the problems which take into asccount layers’ features such as heterogeneity and anisotropy of rock The obtained characteristics like stress distribution, deformation, the progress of damaged regions, and their intensity depend on structure, arrangement of layers, and shape of objects This paper presents the application of geomachenical models on analyzing and predicting movement and deformation of rock considering the relationship with slope angle and depth of seams, as well as excavation methods From the result achieved from the FLAC 2D software, it is confirmed that geomachenical models can be used to predict the law of movement and deformation in underground mining in Vietnam Keywords: Subsidence, Displacement, Deformation, Underground Mining, Geomachenical models ... với độ sâu Hình Bảng Số liệu độ dịch chuyển độ lún sụt lớn Hình a) Sự phụ thuộc biên độ lún hình dạng bồn dịch chuyển; b) Sự phát triển vùng phá hủy (màu tím) vào góc dốc lò chợ Độ sâu khai Dịch. .. đổi góc dốc vỉa, điểm ranh giới lò chợ gần mặt đất Trên hình thể tập hợp kết phân tích độ lún mặt đất, biểu thị qua biến hình lưới sai phân ảnh hưởng khai thác vỉa than quặng với góc dốc vỉa. .. trình biến đổi học khai thác hỗn hợp hầm lò, lộ thiên Trong phạm vi viết tác giả phân tích, dự báo kết dịch chuyển, biến dạng mô hình khai thác lò chợ dài với góc dốc vỉa khác mơ hình khai thác