Nghiên cứu phương pháp dự báo các đại lượng dịch chuyển đất đá và biến dạng bề mặt phù hợp với điều kiện khai thác hầm lò ở Việt Nam

MỞ ĐẦU 1. Tính cấp thiết của đề tài Kết quả của công tác khai thác khoáng sản bằng phương pháp hầm lò là lấy đi một khối lượng lớn đất đá và khoáng sản ra khỏi lòng đất, tạo ra các khoảng trống lớn làm mất trạng thái cân bằng ứng suất nguyên thủy trong lòng đất. Dưới tác động của trọng lực, đất đá phía trên có xu hướng chuyển dịch lấp đầy khoảng trống để tạo lập trạng thái cân bằng mới. Quá trình dịch chuyển lan truyền từ vách qua các lớp đất đá lên đến bề mặt đất gây ra sự biến dạng khối đá và bề mặt địa hình [14, 59, 73] là mối đe dọa nguy hiểm đối với các công trình nằm trong địa tầng đất đá bao quanh khoảng trống khai thác và các công trình tự nhiên, nhân tạo trên bề mặt đất. Đặc biệt, khi khai thác tiến hành dưới các đối tượng chứa nước tự nhiên và nhân tạo, có thể gây ra sự cố bục nước nghiêm trọng làm đình trệ sản xuất, gây thiệt hại về tài sản và nguy hiểm đến tính mạng con người. Trên thế giới, ngay từ đầu thế kỷ 18, tại một số nước Châu Âu như Đức, Bỉ, Pháp, ... việc khai thác than hầm lò đã gây ra những tác động tiêu cực tới môi trường tự nhiên và các công trình trong lòng đất cũng như trên bề mặt, gây nhiều tổn thất nặng nề về người và của cải [74]. Tại Việt Nam, những thiệt hại do dịch chuyển biến dạng bề mặt mỏ đã xảy ra khá phổ biến trong quá khứ cũng như hiện tại. Dưới đấy là một số ví dụ điển hình: - Ở mỏ Mạo Khê trạm quạt ở mức 142 m bị hỏng, phải sửa chữa nhiều lần rất tốn kém, bị bục nước vào lò vỉa 9 qua giếng khoan VIa, lò xuyên vỉa 56 bị biến dạng, bục cát và nước vào lò qua phay A (năm 1995) gây tai nạn chết người; - Mỏ Mông Dương có khu trung tâm vỉa G9 bị bục nước từ lò khai thác cũ với lưu lượng có lúc 200m3/giờ, dịch chuyển đất đá mỏ làm nghiêng cột điện cao thế 110 Kv, gây biến dạng giếng đứng, nứt tường nhà dân [64]; - Hiện tượng sụt lún gây xuống cấp nghiêm trọng đường ô tô lên mỏ Đèo Nai [21]; - Dịch chuyển bề mặt làm hư hại hệ thống đường ray ở khu vực Công ty Than Dương Huy; Mỏ Tây Khe Sim bị bục nước vào tháng 8 năm 2008, ... Cùng với việc gia tăng quy mô và sản lượng khai thác, hậu quả của quá trình dịch chuyển biến dạng mỏ tác động ngày càng lớn với phạm vi rộng hơn, tần suất thường xuyên hơn và theo đó mức độ thiệt hại về người và của cũng lớn hơn. Theo Quyết định 403/QĐ-TTg về việc phê duyệt điều chỉnh Quy hoạch phát triển ngành than Việt Nam đến năm 2020 xét đến năm 2030 [28], các phương pháp khai thác hầm lò ngày càng chiếm ưu thế cả về số lượng mỏ và sản lượng. Bể than Đồng bằng Sông Hồng đang bước vào giai đoạn nghiên cứu khai thác thử nghiệm [15], theo dự báo, sẽ được đưa vào kế hoạch khai thác trong tương lai gần [4]. Dù công nghệ khai thác bể than Đồng Bằng Sông Hồng bằng các phương pháp nào thì ảnh hưởng của nó đối với quá trình sụt lún bề mặt sẽ xảy ra và kèm theo đó là những tác động tiêu cực tới hiện trạng bề mặt. Do vậy, nhu cầu dự báo quy luật sụt lún bề mặt nhằm ngăn ngừa và giảm thiểu các tác động tiêu cực chắc chắn đã, đang và sẽ là vấn đề rất cấp thiết. Quy luật, đặc điểm và cường độ dịch chuyển, biến dạng bề mặt được đặc trưng bởi các các đại lượng, bao gồm: dịch chuyển đứng (lún), dịch chuyển ngang, biến dạng ngang, độ nghiêng, độ cong địa hình. Các đại lượng này là các chỉ số thể hiện mức độ dịch chuyển đất đá và biến dạng bề mặt trong các điều kiện kỹ thuật và công nghệ khai thác khác nhau. Các đại lượng dịch chuyển và biến dạng trên có thể được xác định bằng nhiều phương pháp khác nhau như quan trắc thực địa, nghiên cứu trong phòng thí nghiệm hoặc trên các mô hình toán học, vật lý và địa cơ học. Quy luật, cường độ và quá trình dịch chuyển đất đá và bề mặt chịu sự chi phối của đặc điểm địa hình, địa chất, thế nằm và tính chất cơ lý của các lớp địa tầng và công nghệ khai thác khoáng sản. Do đó, đối với từng vùng mỏ khác nhau, hoặc từng khu vực của một mỏ, cần xây dựng các phương pháp dự báo thích hợp, qua đó cho phép xác định sớm các quy luật, đặc điểm định tính và định lượng của các đại lượng biến dạng làm cơ sở cho việc lựa chọn các phương pháp công nghệ và kỹ thuật khai thác phù hợp, các phương pháp xây dựng, gia cố công trình bề mặt thích ứng, … nhằm mục đích ngăn ngừa và giảm thiểu ảnh hưởng tiêu cực của quá trình khai thác mỏ đối với các công trình. Bên cạnh đó, các kết quả dự báo sớm trị số các đại lượng dịch chuyển, biến dạng theo thời gian sẽ góp phần nâng cao độ tin cậy các kết quả đánh giá tác động

BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC MỎ - ĐỊA CHẤT NGUYỄN QUỐC LONG NGHIÊN CỨU PHƯƠNG PHÁP DỰ BÁO CÁC ĐẠI LƯỢNG DỊCH CHUYỂN ĐẤT ĐÁ VÀ BIẾN DẠNG BỀ MẶT PHÙ HỢP VỚI ĐIỀU KIỆN KHAI THÁC HẦM LÒ Ở VIỆT NAM Ngành: Kỹ thuật Trắc địa - Bản đồ Mã số: 9520503 LUẬN ÁN TIẾN SĨ KỸ THUẬT NGƯỜI HƯỚNG DẪN KHOA HỌC: 1.GS.TS VÕ CHÍ MỸ GVC.TS VƯƠNG TRỌNG KHA Hà Nội - 2019 i MỤC LỤC LỜI CAM ĐOAN LỜI CẢM ƠN MỤC LỤC I DANH MỤC CÁC CHỮ VIẾT TẮT IV DANH MỤC BẢNG BIỂU V DANH MỤC CÁC HÌNH VẼ VII MỞ ĐẦU Chương Tổng quan nghiên cứu dự báo dịch chuyển biến dạng ảnh hưởng khai thác mỏ hầm lò .7 1.1 Dịch chuyển biến dạng bề mặt mỏ khai thác hầm lò 1.1.1 Bản chất toán dịch chuyển đất đá bề mặt khai thác hầm lò 1.1.2 Bồn dịch chuyển 1.1.3 Các thơng số góc trình dịch chuyển 10 1.1.4 Các đại lượng dịch chuyển biến dạng 11 1.1.5 Các yếu tố ảnh hưởng tới dịch chuyển biến dạng bề mặt 11 1.2 Tổng quan cơng trình nghiên cứu dịch chuyển, biến dạng bề mặt mỏ 16 1.2.1 Phân nhóm phương pháp dự báo dịch chuyển, biến dạng mỏ 16 1.2.2 Các kết nghiên cứu tiêu biểu giới 20 1.2.3 Tình hình nghiên cứu nước 26 1.3 Nhận xét chương 30 Chương Cơ sở khoa học lựa chọn phương pháp dự báo dịch chuyển biến dạng bề mặt mỏ hầm lò vùng Quảng Ninh 32 2.1 Đặc điểm địa chất-mỏ bể than Quảng Ninh 32 2.1.1 Đặc điểm tính chất lý đất đá 32 2.1.2 Đặc điểm cấu tạo vỉa than 34 2.1.3 Đặc điểm cơng nghệ khai thác hầm lò vùng Quảng Ninh 38 2.1.4 Đặc điểm lò chợ 40 2.2 Đặc điểm liệu quan trắc dịch chuyển biến dạng 40 ii 2.3 Lựa chọn phương pháp dự báo phù hợp với điều kiện khai thác mỏ hầm lò Việt Nam 42 2.3.1 Tiêu chuẩn lựa chọn phương pháp dự báo 42 2.3.2 Phân tích lựa chọn phương pháp dự báo phù hợp 42 2.4 Nhận xét chương 45 Chương Nghiên cứu mơ hình Asadi dự báo dịch chuyển đứng khai thác hầm lò vỉa dốc Việt Nam 46 3.1 Bồn dịch chuyển trường hợp khai thác vỉa dốc 46 3.2 Mơ hình dự báo dịch chuyển đứng cho vỉa nghiêng 47 3.3 Triển khai mơ hình dự báo Asadi vào thực tiễn Việt nam 49 3.3.1 Xác định dịch chuyển đứng cực đại 49 3.3.2 Xác định góc biên 50 3.3.3 Xác định góc lún cực đại 52 3.3.4 Bán bồn dịch chuyển 53 3.3.5 Xác định hệ số f, g, p, q từ số liệu quan trắc 53 3.3.6 Đánh giá độ xác kết dự báo 56 3.4 Xây dựng chương trình xác định hệ số mơ hình 57 3.5 Xây dựng mơ hình dự báo dịch chuyển đứng khai thác than hầm lò Quảng Ninh, Việt Nam 58 3.5.1 Khu vực nghiên cứu (mỏ Thống Nhất) 58 3.5.2 Khu vực nghiên cứu (mỏ Mông Dương) 62 3.5.3 Xây dựng mơ hình dự báo dịch chuyển đứng sở liệu quan trắc điều kiện địa chất - khai thác mỏ Thống Nhất 66 3.5.4 Áp dụng mơ hình dự báo điều kiện mỏ than Mông Dương 75 3.6 Nhận xét chương 88 Chương Nghiên cứu ứng dụng mạng nơ-ron nhân tạo dự báo dịch chuyển biến dạng bề mặt mỏ theo thời gian 89 4.1 Lý sử dụng mạng nơ-ron nhân tạo 89 4.2 Thông tin mạng nơ-ron nhân tạo 89 4.2.1 Nơ-ron sinh học 89 4.2.2 Nơ-ron nhân tạo 90 iii 4.2.3 Hàm xử lý 92 4.2.4 Mơ hình mạng nơ-ron 93 4.2.5 Huấn luyện mạng 94 4.3 Mạng nơ-ron truyền thẳng nhiều lớp 95 4.4 Thuật toán lan truyền ngược (Back-Propagation) 97 4.5 Lựa chọn yếu tố đầu vào cho mạng nơ-ron dự báo dịch chuyển đứng bề mặt theo thời gian 98 4.6 Dự báo dịch chuyển đứng theo thời gian mạng nơ-ron nhân tạo 100 4.6.1 Lựa chọn cấu hình mạng 100 4.6.2 Huấn luyện mạng 102 4.6.3 Kết dự báo đánh giá độ xác 108 4.7 Nhận xét chương 137 KẾT LUẬN VÀ KIẾN NGHỊ .138 DANH MỤC CÁC CƠNG TRÌNH ĐÃ CƠNG BỐ LIÊN QUAN ĐẾN LUẬN ÁN CỦA NGHIÊN CỨU SINH 140 TÀI LIỆU THAM KHẢO .143 PHỤ LỤC .152 iv DANH MỤC CÁC CHỮ VIẾT TẮT ANN Artificial Neutral Network - Mạng nơ-ron nhân tạo BP Back-Propagation - Thuật toán lan truyền ngược DCBD Dịch chuyển biến dạng DT Dynamic Trough - Bồn động DCĐ Dịch chuyển đứng (độ lún) ĐTM Đánh giá tác động môi trường FA Feedback architecture - Cấu trúc phản hồi, mạng phản hồi FC Feedback connections - Kết nối ngược FW Feed-forward - Truyền thẳng, mạng truyền thẳng KHCNM Viện Khoa học Công nghệ Mỏ (Vinacomin) LMS Least mean square - Nguyên lý bình phương tối thiểu MAE Mean Absolute Error - Sai số tuyệt đối trung bình MSE Mean Square Error - Sai số trung bình MFN Multilayer Feedforward Networks - Mạng truyền thẳng đa lớp RMS Root Mean Square Error - Sai số trung phương ST Static Trough - Bồn tĩnh SW Synaptic weight - Trọng số liên kết TNMT Tài nguyên Môi trường TDKT Thăm dò khai thác VINACOMIN Tập đồn cơng nghiệp Than - Khoáng sản Việt Nam VNIMI Viện nghiên cứu Trắc địa mỏ Địa mỏ Liên bang (Nga) v DANH MỤC BẢNG BIỂU Bảng 1.1 Hàm độ lún mẫu chuẩn 29 Bảng 2.1 Tổng hợp tính chất lý đá số mỏ hầm lò 33 Bảng 2.2 Bảng tổng hợp đặc điểm vỉa than mỏ Nam Mẫu 35 Bảng 2.3 Bảng tổng hợp đặc điểm vỉa than mỏ Hà Lầm 36 Bảng 2.4 Bảng tổng hợp đặc điểm vỉa than mỏ Mông Dương 36 Bảng 2.5 Bảng tổng hợp đặc điểm vỉa than mỏ Thống Nhất 36 Bảng 2.6 Đặc điểm lò chợ khai thác 40 Bảng 2.7: Đặc điểm tuyến quan trắc mỏ 41 Bảng 3.1 Phân loại nhóm mỏ theo độ cứng đất đá 51 Bảng 3.2 Góc biên theo hướng đường phương hướng ngược dốc 52 Bảng 3.3 Góc biên theo hướng xuôi dốc 52 Bảng 3.4 Xác định hệ số K1 52 Bảng 3.5 Số liệu quan trắc tuyến D mỏ Thống Nhất 59 Bảng 3.6 Số liệu quan trắc tuyến P mỏ Thống Nhất 60 Bảng 3.7 Số liệu quan trắc tuyến D mỏ Mông Dương 63 Bảng 3.8 Số liệu quan trắc tuyến P mỏ Mông Dương 64 Bảng 3.9 Các tham số dịch động - địa chất - khai thác 66 Bảng 3.10 Các hệ số mơ hình dự báo lún theo hướng dốc vỉa 66 Bảng 3.11 Đánh giá độ xác kết quan trắc dự báo 67 Bảng 3.12 Độ xác dự báo 68 Bảng 3.13 Các tham số dịch động - địa chất - khai thác 71 Bảng 3.14 Các hệ số mơ hình dự báo 71 Bảng 3.15 Độ xác mơ hình dự báo 72 Bảng 3.16 Kết dự báo độ xác dự báo tuyến P mỏ Thống Nhất 73 vi Bảng 3.17 Các tham số dịch động - địa chất - khai thác tuyến D mỏ Mông Dương 75 Bảng 3.18 Kết dự báo độ xác dự báo tuyến D mỏ Mông Dương 76 Bảng 3.19 Kết dự báo độ xác dự báo tuyến D mỏ Mông Dương (phương pháp VNIMI) 78 Bảng 3.20 Kết dự báo độ xác dự báo tuyến D mỏ Mông Dương (phương pháp KHCNM) 79 Bảng 3.21 Các tham số dịch động - địa chất - khai thác tuyến P mỏ Mông Dương 82 Bảng 3.22 Kết dự báo độ xác dự báo tuyến P mỏ Mông Dương 82 Bảng 3.23 Kết dự báo độ xác dự báo tuyến P mỏ Mông Dương (phương pháp VNIMI) 84 Bảng 3.24 Kết dự báo độ xác dự báo tuyến P mỏ Mông Dương (phương pháp Viện KHCNM) 86 Bảng 4.1 Khả biểu diễn mạng có số lớp ẩn khác 100 Bảng 4.2 Số liệu tập huấn mạng kiểm tra kết dự báo 103 Bảng 4.3 So sánh kết dự báo mạng nơ-ron thực tế trường hợp 108 Bảng 4.4 So sánh kết dự báo mạng nơ-ron thực tế trường hợp 111 Bảng 4.5 So sánh kết dự báo mạng nơ-ron thực tế trường hợp 114 Bảng 4.6 So sánh kết dự báo mạng nơ-ron thực tế trường hợp 120 Bảng 4.7 So sánh kết dự báo mạng nơ-ron thực tế trường hợp 128 vii DANH MỤC CÁC HÌNH VẼ Hình 1.1 Bồn dịch chuyển khai thác hầm lò Hình 1.2 Các vùng ảnh hường bồn dịch chuyển Hình 1.3 Các thơng số bồn dịch chuyển 10 Hình 1.4 Tính chất phân bố dịch chuyển biến dạng mặt đất 11 Hình 1.5 Sự phân bố biến dạng khai thác vỉa có độ dốc khác 12 Hình 1.6 Mối quan hệ độ sâu khai thác dịch chuyển đứng 13 Hình 1.7 Ảnh hưởng chiều rộng lò chợ đại lượng dịch chuyển đứng 14 Hình 1.8 Quan hệ dịch chuyển đứng bề mặt tiến độ khai thác 15 Hình 1.9 Các giai đoạn trình dịch chuyển đứng theo thời gian 16 Hình 1.10 Hàm mặt cắt dịch chuyển đứng 17 Hình 1.11 Mơ hình vật lý làm từ cát khơ 18 Hình 2.1 Đặc điểm độ dốc số vỉa than mỏ Mơng Dương 35 Hình 2.2 Biểu đồ phân bố trữ lượng dự án mỏ hầm lò 38 Hình 2.3 Hệ thống khai thác cột dài theo phương 39 Hình 2.4 Công nghệ chống giữ cột thủy lực đơn 39 Hình 2.5 Tuyến quan trắc V2 chu kỳ mỏ than Mạo Khê 42 Hình 3.1 Bồn dịch chuyển trường hợp khai thác vỉa 46 Hình 3.2 Độ lệch tâm bồn dịch chuyển khai thác vỉa nghiêng 47 Hình 3.3 Các tham số mơ hình dự báo dịch chuyển 48 Hình 3.4 Sơ đồ thuật giải tính hệ số mơ hình dự báo 57 Hình 3.5 Chương trình tính hệ số mơ hình dự báo 58 Hình 3.6 Mặt cắt tuyến quan trắc D mỏ Thống Nhất 60 Hình 3.7 Mặt cắt tuyến quan trắc P mỏ Thống Nhất 61 Hình 3.8 Mặt cắt tuyến quan trắc D mỏ Mông Dương 64 viii Hình 3.9 Mặt cắt tuyến quan trắc P mỏ Mông Dương 65 Hình 3.10 Biểu đồ so sánh đánh giá độ xác hệ số mơ hình dự báo tuyến D 68 Hình 3.11 So sánh đường cong lún thực tế dự báo tuyến D mỏ Thống Nhất 69 Hình 3.12 Tương quan kết quan trắc dự báo tuyến D mỏ Thống Nhất 70 Hình 3.13 Biểu đồ so sánh đường cong lún đánh giá độ xác hệ số mơ hình dự báo tuyến P mỏ Thống Nhất 72 Hình 3.14 So sánh đường cong lún thực tế dự báo tuyến P mỏ Thống Nhất 74 Hình 3.15 Tương quan độ lún quan trắc dự báo tuyến P mỏ Thống Nhất 74 Hình 3.16 So sánh đường cong lún thực tế dự báo tuyến D mỏ Mơng Dương 77 Hình 3.17 So sánh đường cong lún thực tế dự báo tuyến D mỏ Mông Dương (phương pháp VNIMI) 79 Hình 3.18 So sánh đường cong lún thực tế dự báo tuyến D mỏ Mông Dương (phương pháp KHCNM) 81 Hình 3.19 So sánh đường cong lún thực tế dự báo tuyến P mỏ Mơng Dương 83 Hình 3.20 So sánh đường cong lún thực tế dự báo tuyến P mỏ Mông Dương (phương pháp VNIMI) 85 Hình 3.21 So sánh đường cong dịch chuyển đứng thực tế kết dự báo tính từ phương pháp KHCNM 87 Hình 4.1 Mơ hình Nơ-ron sinh học 90 Hình 4.2 Nơ-ron nhân tạo 91 Hình 4.3 Đồ thị hàm Sigmoid 93 Hình 4.4 Kiến trúc mạng nơ-ron truyền thẳng 94 Hình 4.5 Kiến trúc mạng nơ-ron phản hồi 94 Hình 4.6 Quy tình huấn luyện mạng có giám sát 95 Hình 4.7: Mạng nơ-ron truyền thẳng đa lớp tổng quát 96 Hình 4.8 Các yếu số đầu vào mạng nơ-ron 99 Hình 4.9 Mơ hình mạng nơ-ron dự báo dịch chuyển đứng theo thời gian 102 ix Hình 4.10: Định dạng tệp huấn luyện mạng 106 Hình 4.11 Đường cong lỗi huấn luyện mạng 106 Hình 4.12 Biểu đồ sai số huấn luyện mạng 107 Hình 4.13 Biểu đồ tương quan huấn luyện mạng 107 Hình 4.14 So sánh đường cong dự báo thực tế chu kỳ 12 (trường hợp 1) 110 Hình 4.15 Tương quan giá trị dự báo thực tế chu kỳ 12 (trường hợp 1) 110 Hình 4.16 So sánh đường cong lún dự báo thực tế chu kỳ 11 (trường hợp 2) 112 Hình 4.17 Tương quan giá trị dự báo thực tế chu kỳ 11 (trường hợp 2) 113 Hình 4.18 So sánh đường cong lún dự báo thực tế chu kỳ 12 (trường hợp 2) 113 Hình 4.19 Tương quan giá trị dự báo thực tế chu kỳ 12 (trường hợp 2) 114 Hình 4.20 So sánh đường cong dự báo thực tế chu kỳ 10 (trường hợp 3) 117 Hình 4.21 Tương quan giá trị dự báo thực tế chu kỳ 10 (trường hợp 3) 118 Hình 4.22 So sánh đường cong dự báo thực tế chu kỳ 11 (trường hợp 3) 118 Hình 4.23 Tương quan giá trị dự báo thực tế chu kỳ 11 (trường hợp 3) 119 Hình 4.24 So sánh đường cong dự báo thực tế chu kỳ 12 (trường hợp 3) 119 Hình 4.25 Tương quan giá trị dự báo thực tế chu kỳ 12 (trường hợp 3) 120 Hình 4.26 So sánh đường cong dự báo thực tế chu kỳ (trường hợp 4) 124 Hình 4.27 Tương quan giá trị dự báo thực tế chu kỳ (trường hợp 4) 124 Hình 4.28 So sánh đường cong dự báo thực tế chu kỳ 10 (trường hợp 4) 125 Hình 4.29 Tương quan giá trị dự báo thực tế chu kỳ 10 (trường hợp 4) 125 Hình 4.30 So sánh đường cong dự báo thực tế chu kỳ 11 (trường hợp 4) 126 Hình 4.31 Tương quan giá trị dự báo thực tế chu kỳ 11 (trường hợp 4) 126 Hình 4.32 So sánh đường cong dự báo thực tế chu kỳ 12 (trường hợp 4) 127 Hình 4.33 Tương quan giá trị dự báo thực tế chu kỳ 12 (trường hợp 4) 127 Hình 4.34 So sánh đường cong dự báo thực tế chu kỳ (trường hợp 5) 132 158   -40.000 -1.300 -14.000 -1.377 0.000 -1.386 32.000 -1.076 49.000 -0.802 60.000 -0.630 93.000 -0.243 108.000 -0.142 123.000 -0.078 156.000 -0.017 -Vong lap: Ma tran A 0.0259 0.0996 0.1165 0.1281 0.1088 0.0894 0.0592 0.0210 0.0023 0 0 0 0 -0.1363 -0.2083 -0.3274 -0.4017 -0.3885 -0.3196 -0.1587 -0.0267 0 0 0 0 0 0 0 0 0.0616 0.0991 0.1122 0.0956 0.0730 0.0505 0.0165 0 0 0 0 -0.4320 -0.5081 -0.4747 -0.2189 -0.1189 -0.0531 Ma tran N 0.0643 -0.1805 0 -0.1805 0.6096 0 0 0.0435 -0.1625 0 -0.1625 0.7350 159   Ma tran N-1 91.9887 27.2452 0 27.2452 9.7100 0 0 132.2891 29.2414 So hieu chinh df: -1.4847 dg: -0.4385 dp: -0.0599 dq: -0.0101 He so cua mo hinh f: 5.4651 g: 2.6082 p: 4.4887 q: 1.8158 Danh gia chinh xac RMS: 0.080 MAE: 0.060 r: 0.987 Do lun tinh tu -268.000 -190.000 -171.000 -141.000 -106.000 -90.000 -69.000 -40.000 -14.000 0.000 32.000 49.000 60.000 93.000 mo hinh -0.006 -0.149 -0.255 -0.498 -0.852 -1.009 -1.183 -1.334 -1.383 -1.386 -1.076 -0.801 -0.628 -0.240 0 29.2414 7.8241 160   108.000 -0.139 123.000 -0.075 156.000 -0.016 -Vong lap: Ma tran A 0.0059 0.0609 0.0790 0.0933 0.0758 0.0586 0.0343 0.0093 0.0006 0 0 0 0 -0.1144 -0.1940 -0.3274 -0.3844 -0.3495 -0.2547 -0.0971 -0.0101 0 0 0 0 0 0 0 0 0.0606 0.0978 0.1108 0.0937 0.0711 0.0488 0.0156 0 0 0 0 -0.4312 -0.5086 -0.4751 -0.2176 -0.1174 -0.0520 Ma tran N 0.0291 -0.1121 0 -0.1121 0.5022 0 0 0.0420 -0.1598 0 -0.1598 0.7340 Ma tran N-1 243.7985 54.4274 0 54.4274 14.1420 0 0 139.0448 30.2685 0 30.2685 7.9515 161   So hieu chinh df: -0.7774 dg: -0.1085 dp: -0.0078 dq: -0.0013 He so cua mo hinh f: 6.2425 g: 2.7167 p: 4.4965 q: 1.8171 Danh gia chinh xac RMS: 0.078 MAE: 0.060 r: 0.988 Do lun tinh tu mo hinh -268.000 -0.003 -190.000 -0.119 -171.000 -0.220 -141.000 -0.466 -106.000 -0.839 -90.000 -1.004 -69.000 -1.185 -40.000 -1.338 -14.000 -1.383 0.000 -1.386 32.000 -1.076 49.000 -0.801 60.000 -0.628 93.000 -0.239 108.000 -0.138 123.000 -0.075 156.000 -0.015 162   Vong lap: Ma tran A 0.0027 0.0469 0.0648 0.0814 0.0675 0.0518 0.0297 0.0076 0.0005 0 0 0 0 -0.1007 -0.1818 -0.3262 -0.3908 -0.3530 -0.2516 -0.0905 -0.0084 0 0 0 0 0 0 0 0 0.0605 0.0977 0.1106 0.0935 0.0709 0.0486 0.0155 0 0 0 0 -0.4310 -0.5086 -0.4751 -0.2174 -0.1172 -0.0519 Ma tran N 0.0212 -0.0959 0 -0.0959 0.4985 0 0 0.0418 -0.1594 0 -0.1594 0.7339 Ma tran N-1 361.4707 69.5278 0 69.5278 15.3795 0 0 139.9360 30.4032 0 30.4032 7.9682 So hieu chinh df: -0.2216 dg: -0.0312 dp: -0.0009 dq: -0.0002 163   He so cua mo hinh f: 6.4641 g: 2.7479 p: 4.4974 q: 1.8173 Danh gia chinh xac RMS: 0.078 MAE: 0.061 r: 0.988 Do lun tinh tu -268.000 -190.000 -171.000 -141.000 -106.000 -90.000 -69.000 -40.000 -14.000 0.000 32.000 49.000 60.000 93.000 108.000 123.000 156.000 mo hinh -0.002 -0.112 -0.211 -0.458 -0.836 -1.004 -1.187 -1.339 -1.383 -1.386 -1.076 -0.801 -0.628 -0.239 -0.138 -0.075 -0.015 164   XAC DINH HE SO MO HINH TUYEN P MO THONG NHAT So lieu quan trac -186.000 0.000 -149.000 -0.123 -102.000 -0.241 -76.000 -0.513 -61.000 -0.810 0.000 -1.418 32.000 -1.325 76.000 -1.112 105.000 -0.783 130.000 -0.547 146.000 -0.345 198.000 -0.014 207.000 0.000 He so khoi tinh f0: 11.5271 g0: 2.7129 p0: 1.0662 q0: 1.4755 -Vong lap: Ma tran A 0.0000 0.0014 0.0290 0.0452 0.0393 0 0 0 -0.0036 -0.2010 -0.4668 -0.5057 0 0 0 0 0 0.0843 0.2535 0.3521 0.4173 0.4480 0.4892 0 0 -0.1678 -0.2708 -0.2548 -0.2070 -0.1668 -0.0232 165   Ma tran N 0.0044 -0.0469 0 0.4882 -0.0469 0.5141 0 0 1.0479 -0.3450 0 -0.3450 0.2376 Ma tran N-1 5.9338 0.5409 0 0.5409 0.0513 0 0 0.0018 0.0027 0 0.0027 0.0081 So hieu chinh df: 6.6675 dg: 0.5707 dp: -1.0658 dq: -0.7796 He so cua mo hinh f: 4.8595 g: 2.1423 p: 2.1321 q: 2.2551 Danh gia chinh xac RMS: 0.247 MAE: 0.159 r: 0.898 Do lun tinh tu -186.000 -149.000 -102.000 -76.000 -61.000 0.000 mo hinh -0.000 -0.003 -0.148 -0.513 -0.810 -1.418 166   32.000 -1.325 76.000 -1.112 105.000 -0.958 130.000 -0.829 146.000 -0.750 198.000 -0.522 207.000 -0.488 -Vong lap: Ma tran A 0.0110 0.0430 0.1023 0.1020 0.0833 0 0 0 0 -0.0463 -0.2988 -0.4438 -0.4514 0 0 0 0 0 0 0.0204 0.1185 0.1934 0.2354 0.2446 0.1864 0.1682 0 0 -0.0812 -0.2531 -0.2799 -0.2334 -0.1820 -0.0177 Ma tran N 0.0298 -0.1155 0 -0.1155 0.4921 0 0 0.2301 -0.1886 0 -0.1886 0.2370 Ma tran N-1 369.6823 86.7328 0 86.7328 22.3809 0 0 12.4873 9.9360 0 9.9360 12.1259 167   So hieu chinh df: 0.0052 dg: 0.3279 dp: -0.9921 dq: -0.3318 He so cua mo hinh f: 4.8543 g: 1.8143 p: 3.1242 q: 2.5869 Danh gia chinh xac RMS: 0.047 MAE: 0.040 r: 0.996 Do lun tinh tu mo hinh -186.000 -0.011 -149.000 -0.055 -102.000 -0.277 -76.000 -0.545 -61.000 -0.746 0.000 -1.418 32.000 -1.383 76.000 -1.122 105.000 -0.827 130.000 -0.555 146.000 -0.400 198.000 -0.088 207.000 -0.062 -Vong lap: Ma tran A 0.0111 0.0369 -0.0397 0 0 168   0.0932 0.1074 0.0987 0 0 0 -0.2718 -0.4664 -0.5341 0 0 0 0 0 0.0110 0.0840 0.1428 0.1666 0.1620 0.0780 0.0624 0 -0.0644 -0.2630 -0.3028 -0.2422 -0.1767 -0.0108 Ma tran N 0.0314 -0.1296 0 -0.1296 0.5784 0 0 0.0916 -0.1359 0 -0.1359 0.2550 Ma tran N-1 416.9181 93.4259 0 93.4259 22.6646 0 0 52.1834 27.8045 0 27.8045 18.7363 So hieu chinh df: -0.3444 dg: -0.0983 dp: -0.5471 dq: -0.1657 169   He so cua mo hinh f: 5.1987 g: 1.9127 p: 3.6713 q: 2.7527 Danh gia chinh xac RMS: 0.039 MAE: 0.032 r: 0.997 Do lun tinh tu mo hinh -186.000 -0.008 -149.000 -0.047 -102.000 -0.273 -76.000 -0.555 -61.000 -0.766 0.000 -1.418 32.000 -1.388 76.000 -1.124 105.000 -0.805 130.000 -0.511 146.000 -0.348 198.000 -0.055 207.000 -0.036 -Vong lap: Ma tran A 0.0078 0.0309 0.0865 0.1001 0.0908 0 0 -0.0357 -0.2702 -0.4660 -0.5261 0 0 0 0 0.0081 0.0712 0.1242 0 0 -0.0558 -0.2621 -0.3095 170   0 0 0 0 0.1421 0.1332 0.0487 0.0361 -0.2426 -0.1707 -0.0080 Ma tran N 0.0268 -0.1189 0 -0.1189 0.5682 0 0 0.0622 -0.1151 0 -0.1151 0.2556 Ma tran N-1 529.0434 110.7086 110.7086 24.9272 0 0 0 97.1567 43.7629 0 43.7629 23.6242 So hieu chinh df: -0.1114 dg: -0.0197 dp: -0.0221 dq: 0.0110 He so cua mo hinh f: 5.3101 g: 1.9324 p: 3.6934 q: 2.7417 Danh gia chinh xac RMS: 0.038 MAE: 0.031 r: 0.997 Do lun tinh tu mo hinh -186.000 -0.007 -149.000 -0.045 171   -102.000 -0.269 -76.000 -0.553 -61.000 -0.766 0.000 -1.418 32.000 -1.387 76.000 -1.119 105.000 -0.798 130.000 -0.505 146.000 -0.343 198.000 -0.054 207.000 -0.035 -Vong lap: Ma tran A 0.0070 0.0291 0.0842 0.0981 0.0888 0 0 0 0 -0.0342 -0.2686 -0.4661 -0.5259 0 0 0 0 0 0 0.0083 0.0717 0.1242 0.1412 0.1318 0.0477 0.0353 0 0 -0.0572 -0.2655 -0.3113 -0.2425 -0.1700 -0.0078 Ma tran N 0.0255 -0.1160 0 -0.1160 0.5671 0 0 0.0615 -0.1152 0 -0.1152 0.2585 172   Ma tran N-1 571.4828 116.9248 116.9248 25.6861 0 0 So hieu chinh df: -0.0380 dg: -0.0074 dp: 0.0074 dq: 0.0038 He so cua mo hinh f: 5.3481 g: 1.9398 p: 3.6860 q: 2.7378 Danh gia chinh xac RMS: 0.038 MAE: 0.031 r: 0.997 Do lun tinh tu -186.000 -149.000 -102.000 -76.000 -61.000 0.000 32.000 76.000 105.000 130.000 146.000 198.000 207.000 mo hinh -0.007 -0.044 -0.268 -0.553 -0.767 -1.418 -1.387 -1.119 -0.798 -0.506 -0.344 -0.054 -0.036 0 98.8231 44.0488 0 44.0488 23.5031 ... CHUYỂN BIẾN DẠNG DO ẢNH HƯỞNG CỦA KHAI THÁC MỎ HẦM LÒ 1.1 Dịch chuyển biến dạng bề mặt mỏ khai thác hầm lò 1.1.1 Bản chất toán dịch chuyển đất đá bề mặt khai thác hầm lò Việc đào lò chuẩn bị khai thác. .. án tiến sĩ Nghiên cứu phương pháp dự báo đại lượng dịch chuyển đất đá biến dạng bề mặt phù hợp với điều kiện khai thác hầm lò Việt Nam lựa chọn xuất phát từ nhu cầu thực tế, đáp ứng yêu cầu thực... chọn phương pháp dự báo dịch chuyển biến dạng bề mặt mỏ hầm lò bể than Quảng Ninh Chương Nghiên cứu mơ hình Asadi dự báo dịch chuyển đứng khai thác hầm lò vỉa dốc Việt Nam Chương Nghiên cứu ứng