- Triển vọng loại A là các đới chưa được tiến hành kiểm tra mặt đất nhưng có kết quả tốt trong các phương pháp phân tích nhận dạng.
- Triển vọng loại B là các đới chưa được tiến hành kiểm tra mặt đất và các kết quả phân tích nhận dạng theo các phương pháp đã tiến hành khơng hồn tồn trùng nhau.
Kết quả là đã khoanh định được 9 đới triển vọng khống sản. Trong đó có 3 đới triển vọng loại A1 là những đới triển vọng đã được kiểm tra, 4 đới triển vọng loại A là những đới chưa được kiểm tra, 2 đới triển vọng loại B.
Khoáng sản trong vùng biểu hiện khá phong phú, đa dạng, đáng chú ý là vàng, thiếc, volfram, molipden, đất hiếm, trong đó nổi bật nhất là vàng. Vàng được phát hiện nhiều nơi trên nhiều đối tượng địa chất với nhiều thành hệ khác nhau, được biểu hiện chủ yếu bằng các đới dị thường mang bản chất Kali. Khoáng sản đáng chú ý thứ 2 là volfram trong các đới biến đổi như greizen hóa, các mạch pecmatit được biểu hiện bằng các đới dị thường bản chất thori và hỗn hợp.
KẾT LUẬN
Trên cơ sở thu thập tài liệu, tìm hiểu các phương pháp phân tích tổ hợp tài liệu địa vật lý đặc biệt là tài liệu từ - phổ gamma hàng không trong nước và trên thế giới, học viên đã tiến hành nghiên cứu ứng dụng phương pháp phân tích hệ số tương quan trong xử lý phân tích tài liệu phổ gamma hàng khơng.
Với việc hồn thành các mục tiêu, nhiệm vụ đặt ra luận văn đã đạt được một số kết quả chính như sau:
- Nghiên cứu phương pháp hệ số tương quan từ đó đưa ra khả năng ứng dụng của phương pháp trong xử lý phân tích số liệu phổ gamma hàng khơng nhằm góp phần nâng cao hiệu quả trong xử lý tài liệu.
- Thành lập một chương trình tự động xử lý số liệu phổ gamma hàng không theo phương pháp hệ số tương quan.
- Thực hiện phân tích thử nghiệm chương trình với số liệu giả định và số liệu thực tế vùng Đông tỉnh Đak Lak. Từ các kết quả đạt được khẳng định khả năng ứng dụng của phương pháp hệ số tương quan góp phần nâng cao hiệu quả xử lý tài liệu phổ gamma hàng không
Từ các kết quả đạt được học viên đã xử dụng phương pháp, kết hợp với các tài liệu thu thập được từ đó thành lập một sơ đồ dự báo triển vọng khoáng sản của khu vực rộng khoảng 87km2 thuộc phần Đông tỉnh Đak Lak.
Các kết quả nghiên cứu đạt được đã góp phần nói lên tính đúng đắn và độ tin cậy của phương pháp cũng như ý nghĩa thực tiễn và phạm vi ứng dụng của các phương pháp được lựa chọn. Sơ đồ được thành lập là những kết quả mới, khách quan góp phần làm sáng tỏ thêm về đặc điểm khoáng sản khu vực nghiên cứu, đồng thời mở ra hướng nghiên cứu hoàn thiện tiếp theo.
TÀI LIỆU THAM KHẢO Tiếng Việt
1. Trần Trọng Hòa (2005), Nghiên cứu điều kiện thành tạo và quy luật phân bố khoáng sản quý hiếm liên quan đến hoạt động magma khu vực miền Trung và Tây Nguyên, Đề tài độc lập cấp Nhà nước.
2. Đặng Mai (2004), Toán ứng dụng trong địa chất, Nxb Đại học Quốc gia Hà Nội,
Hà Nội.
3. Tăng Mười, Võ Thanh Quỳnh (1998), Ứng dụng phương pháp phổ gamma hàng
khơng trong tìm kiếm Uran và các khống sản có ích khác liên quan với phóng xạ, Báo cáo tại Hội nghị khoa học vật lý hạt nhân ứng dụng, Viện Khoa học Việt
Nam, Hà Nội
4. Lê Khánh Phồn (2004), Thăm dị phóng xạ, NXB giao thông vận tải.
5. Lê Khánh Phồn (10/1996), Phương pháp xử lý luận giải tài liệu phổ gamma mặt
đất mơ phỏng theo địa hố, Báo cáo khoa học tại Hội nghị khoa học lần thứ 12,
Trường Đại học Mỏ - Địa chất.
6. Võ Thanh Quỳnh (chủ biên), Nguyễn Xuân Sơn, Nguyễn Văn Phùng (1995),
Báo cáo kết quả bay đo từ - phổ gamma tỉ lệ 1:25.000 vùng Tuy Hòa, Lưu trữ
Cục Địa Chất và Khoáng Sản Việt Nam.
7. Võ Thanh Quỳnh, (2007), “Một cách tiếp cận mới giải quyết bài toán nhận dạng trong xử lý, phân tích tài liệu địa vật lý”, TC Địa chất, A/302:76-80, Hà Nội.
8. Võ Thanh Quỳnh (02/2008), “phương pháp đánh giá phân loại cụm dị thường trong xử lý phân - tích tài liệu phổ gamma hàng không”, TC Địa chất, A/304:70-75, Hà Nội.
9. Võ Thanh Quỳnh, Nguyễn Xuân Ngoan, Vũ Tuấn Hùng (2002), Thành lập bộ
chương trình xử lý tài liệu địa vật lý máy bay, Lưu trữ Cục Địa chất và Khoáng
sản Việt Nam, Hà Nội.
10. Tống Duy Thanh, Vũ Khúc (chủ biên) (2005), Các phân vị địa tầng Việt Nam, Nhà xuất bản Đại học Quốc gia Hà Nội.
11. Nguyễn Tài Thinh, Võ Thanh Quỳnh và n.n.k (1997), Áp dụng các phương pháp
mới trong xử lý, phân tích, biểu diễn và lưu trữ các tài liệu địa vật lý, Lưu trữ
Cục Địa chất Việt Nam, Hà Nội.
12. Cao Hào Thi (2008), Giáo trình xác suất thống kê. TP. Hồ Chí Minh
13. Cao Đình Triều (2005), Trường Địa vật lý và cấu trúc thạch quyển lãnh thổ Việt
Nam, Nhà xuất bản Khoa học và kỹ thuật.
14. Phạm Năng Vũ (2002), Bài giảng cơ sở lý thuyết xử lý số liệu địa vật lý, Trường Đại học Mỏ Địa chất.
PHỤ LỤC
- Phụ lục 1: Các bản vẽ, tài liệu thu thập vùng Đông Đak Lak - Phụ lục 2: Mã nguồn chương trình tính hệ số tương quan
Phụ lục 2: Mã nguồn chương trình tính hệ số tương quan #include<stdio.h> #include<conio.h> #include<math.h> #include<stdlib.h> #include <iostream> int main(){ int n,m,i,j,j5,o,i1,i2,i3,i4,i5,k1,k2,k3,k4,k5,k[20] ; double tong2; float tb1, tb2, tb3, p[20000][5], cs[2000][5], tq[20000][5], xiyi1,xiyi2,xiyi3, nxy1, nxy2, nxy3, a, b, c, d, e, f, g, tq1, tq2, tq3, h[20][2000], h1[20][2000], tx[20], tx1[20], minx, maxx, miny, maxy, chuyenx, chuyeny, x, y, y1, y2, y3, y4, y5, y6, y7, y8, y9, xy, xy2, khoang[20][2], max, min , tong, p1[20][2000],;
char filename[20],filename2[20];
FILE *file, *output, *file2;
printf ("\n\n\n\t\t\tCHUONG TRINH PHAN TICH HE SO TUONG
QUAN
\n\t\t************************************************** **\n\n\n\t--- NHAP DU LIEU --- \n\n\tNhap vao file du lieu \n\n\t luu y!!! File du lieu la file tex (***.txt) phai nam trong cung \n\tmot foder voi chuong trinh\n\tCau truc cua file du lieu la:\n\t\t- Dong dau tien la tong so diem du lieu\n\t\t- Moi dong tiep theo la so lieu cua 1 diem bao gom:\n 1.Toa do x\t 2.Toa do y\t3. Ham luong U \t4.Ham luong Th\t5.Ham luong K\n\n\n**Vi
du:\n\t3564\n25.3\t37.4\t3.443\t43.552\t0.233\n....\t... .\t....\t....\t....\n ");
fflush(stdin);
thoi la 20.000 diem du lieu:\n\n\t\tNhap ten file du lieu:\t ");
scanf("%s",filename);
file= fopen(filename,"r"); fscanf(file,"%d",&m);
printf ("\n\n\n\n\n\n\n\n\n\n\n\n\n\n\t\t Tong so diem
du lieu M = %d \n\n\n\t Press any key to continue.....",m); getch(); j=0; for (j=0;j<m ;j++) { i=0; printf ("\n\n"); for (i=0;i<5;i++) { fscanf(file,"%f",&a); p[j][i] = a; printf ("%13.4f ",p[j][i]); } printf ("\n\n"); } fclose(file); printf ("\n\n\n\n\n\n\n\n\n\n\n\t\t\t------- DA NHAP XONG DU LIEU -------\n\n\n\t Press any key to continue.....");
getch();
printf ("\n\n\n\t ------ CHON PHUONG AN XU LY SO LIEU
------ \n\n\n\t De tinh he so tuong quan cua 1 cum di thuong chon 1\n\n\tDe xu ly theo toan dien tich chon 2 ");
scanf("%d",&i);
//Tinh he so tuong quan theo 1 cum di thuong ******
if(i==1)
{
printf ("\n\n\n\n\n\n\n\n\n\n\n\t\t\t------- Chon phuong an 1: Ting tuong quan cua 1 cum di thuong ------- \n\n\n\t Press any key to continue.....");
getch();
//tinh gia tri trung binh cua ham luong U, Th, K tb1=0; tb2=0; tb3=0; for(j= 0;j<m;j++) { tb1= tb1+p[j][2]; tb2= tb2+p[j][3]; tb3= tb3+p[j][4]; } tb1 = tb1/m ; tb2 = tb2/m ; tb3 = tb3/m ;
printf ("\n\n\n\t Gia tri trung binh: \n\n\tU =
%f\tTh = %f\tK = %f \n\t ",tb1,tb2,tb3);
// Tinh he so tuong quan ******
i1=0; xiyi1=0; xiyi2=0; xiyi3=0;tq1=0;tq2=0;tq3=0;a=0;b=0;c=0;
nxy1=tb1*tb2*m; nxy2=tb1*tb3*m; nxy3=tb2*tb3*m;
for (i1=0 ; i1<m ; i1++) {
tq2=tq2+ (p[i1][2]-tb1)*(p[i1][4]-tb3); tq3=tq3+ (p[i1][3]-tb2)*(p[i1][4]-tb3); }
for (i1=0 ; i1<m ; i1++) {
a= a + pow((p[i1][2]-tb1),2); b= b + pow((p[i1][3]-tb2),2); c= c + pow((p[i1][3]-tb2),2); }
printf ("\n\n\n\t tu : \n\n\tU/Th = %f\t U/K =
%f\tTh/K = %f \n\n\t Mau a = %f \t b =%f \t c = %f ",tq1,tq2,tq3,a,b,c); getch(); y1=a*b; y1=sqrt(y1); y2=a*c; y2=sqrt(y2); y3=b*c; y3=sqrt(y3); tq1= tq1/y1; tq2= tq2/y2; tq3= tq3/y3;
printf ("\n\n\n\t He so tuong quan: \n\n\tU/Th = %f\t
U/K = %f\tTh/K = %f \n\t ",tq1,tq2,tq3);
getch();
printf ("\n\n\n\n\n\t\t\t------- DA TINH XONG HE SO
TUONG QUAN-------\n\n\n\t\t Press any key to QUIT.....");
}
// Phan 2: Tinh he so tuong quan tren toan dien tich
if (i==2)
printf ("\n\n\n\n\n\t\t\t------- Chon phuong an 2:
tinh tuong quan tren toan dien tich -------\n\n\n\t Press any key to continue.....");
printf ("\n\n\tSo lieu cua khu vuc nghien cuu da duoc
nhap xong\n\tNhap vao do dai canh cua o vuong cua so quet\n\n\t");
printf ("\n\n\t Cua so quet rong =\t "); scanf ("%f",&f);
printf ("\n\n\tBan muon lam voi O vuong co chieu rong
la: %15.2f\n\n\n\n\t Nhap vao ten file xuat du lieu\n\n\t ",f);
scanf("%s",filename);
output = fopen(filename,"r");
/*tim gia tri lon nhat va nho nhat cua toa do x, y*/
minx = maxx = p[0][0]; miny = maxy = p[0][1]; for (j=1; j<m ;j++) { if (minx>p[j][0]) minx=p[j][0]; if (maxx<p[j][0]) maxx=p[j][0]; if (miny>p[j][0]) miny=p[j][0]; if (maxy<p[j][0]) maxy=p[j][0]; }
printf ("\n\tMax X = %6.6f\n\tMinX= %6.6f\n\tMaxY=
%6.6f\n\tMinY= %6.6f\n\t",maxx,minx,maxy,miny);
// tinh toa do tam cua cua so quet
y8=minx; y9=miny; k1=0;k2=0;
while (y8 < (maxx-f*0.3))
{ k1++;
y8=y8+f/8;
while (y9 < (maxy-f*0.3)) {
k2++;
y9=y9+f/8; tq[k2][0]=y8;
tq[k2][1]=y9;
//vong lap nhat du lieu cho cua so quet va tinh he
so tuong quan
for (i1=0; i1 < m; i1++) {
i2=0;
y6 = p[i1][0]; y7 = p[i1][1];
if(y6 > (y8 - f/2) && y6 < (y8 + f/2)&& y7 > (y9 - f/2) && y7 < (y9 + f/2))
{
printf ("\n\t x = %f \t y= %f\n\t", p[i1][0], p[i1][1]);
for (i3=0; i3 < 5; i3++)
{cs[i2][i3]=p[i1][i3];}
//nhap vao du lieu cho cua so quet
i2++; }
//tinh he so tuong quan cua cua so quet tb1=0; tb2=0; tb3=0;
for(j5= 0;j<i2;j++) {
tb1= tb1+cs[j5][2]; tb2= tb2+cs[j5][3]; tb3= tb3+cs[j5][4]; } tb1 = tb1/i2 ; tb2 = tb2/i2 ; tb3 = tb3/i2 ;
j5=0; xiyi1=0; xiyi2=0; xiyi3=0; tq1=0; tq2=0; tq3=0; a=0 ; b=0; c=0; nxy1=tb1*tb2*i2; nxy2=tb1*tb3*i2; nxy3=tb2*tb3*i2; for (j5=0 ; j5<i2 ; j5++) { tq1=tq1+ (cs[j5][2]-tb1)*(cs[j5][3]-tb2); tq2=tq2+ (cs[j5][2]-tb1)*(cs[j5][4]-tb3); tq3=tq3+ (cs[j5][3]-tb2)*(cs[j5][4]-tb3); } for (j5=0 ; j5<i2 ; j5++) { a= a + pow((cs[j5][2]-tb1),2); b= b + pow((cs[j5][3]-tb2),2); c= c + pow((cs[j5][3]-tb2),2); } y1=a*b; y1=sqrt(y1); y2=a*c; y2=sqrt(y2); y3=b*c; y3=sqrt(y3); tq1= tq1/y1;
tq3= tq3/y3;
printf ("\n\n\n\t He so tuong quan: \n\n\tU/Th
= %f\t U/K = %f\tTh/K = %f \n\t ",tq1,tq2,tq3); } tq[k2][2]=tq1; tq[k2][3]=tq2; tq[k2][4]=tq3; fprintf (output, "\n\t %f\t%f\t%f\t%f\t%f", &tq[k2][0], &tq[k2][1], &tq[k2][2], &tq[k2][3], &tq[k2][4]);
} }
printf ("\n\t So luong diem tren luoi = %d
\t\n\t",k1*k2);
fclose(file);
printf ("\n\n\n\t\t------- DA TINH XONG -------\n\n\n\t
Press any key to QUIT.....");
getch(); }
if (i!=1& i!=2)
printf ("\n\n\n\n\n\t\t\t------- Chon sai phuong an ---
----\n\n\n\t Press any key to QUIT.....");