Đề bài đồ án : Các tham số của mô hình bài toán thuận như sau: Giả sử ta khảo sát trên một mạng lưới có diện tích 310x310 m2, các nút lưới theo phương trục x (Bắc), và phương trục y (Đông) cách đều nhau một khoảng là 10m. Trong diện tích khảo sát có hai đối tượng dạng cầu bị từ hóa hoàn toàn do cảm ứng dưới tác dụng của trường từ T0 = 44000 g. Quả cầu thứ nhất nằm ở vị trí x1 = 150 m, y1 = 150 m, độ sâu z1 = 200 m, có bán kính R1 = 80 m và độ từ cảm c1 = 90000x106 theo đơn vị CGS. Quả cầu thứ hai nằm ở vị trí x2 = 250 m, y2 = 250 m, độ sâu z2 = 100 m, có bán kính R2 = 40 m và độ từ cảm c2 = 45000x106 theo đơn vị CGS. Vùng nghiên cứu được giả thiết có góc nghiêng I = 150 và góc lệch D = 00. Hình vẽ mô tả đề bài:
Phần I: Lời Mở Đầu Giới thiệu phương pháp thăm dò từ Thăm dò từ phương pháp địa vật lý áp dụng để giải nhiệm vụ địa chất sở nghiên cứu đặc điểm trường từ đất mức độ nhiễm từ khác loại đất đá quặng tạo 2.Tính cấp thiết toán thăm dò từ Trong số phương pháp Địa vật lý phương pháp từ có nhiều nét giống với phương pháp trọng lực Những nét riêng phương pháp từ đặc điểm trường từ trường hai cực, có phương chiều thay đổi Cũng phương pháp trường khác, hạn chế đáng kể phương pháp từ tính đa trị giải thích tài liệu Trong thăm dò từ, để có sở bố trí tuyến điểm đo, hoạch định kỹ thuật đo ghi, chỉnh lý số liệu lựa chọn phép xử lý giải thích phù hợp trường dị thường liên quan đến đối tượng địa chất mà ta muốn phát nghiên cứu, phải dự đoán hay hình dung nét phân bố trường gây gán cho chúng đặc điểm hình học vật lý theo thông tin ban đầu mà ta có Bài toán suy trường biết nguồn gọi toán thuận Do thực tế đối tượng địa chất đa dạng phức tạp nên việc xác định trường cho đối tượng cụ thể khó khăn nhiều trường hợp làm Giải pháp thường sử dụng xấp xỉ trường đối tượng địa chất cụ thể hay nhiều nguồn đơn giản mà biểu thức tính trường biết Một nguồn đơn giản điển hình cầu Trong đề tài trình bày phần trường dị thường từ lý thuyết vật thể dạng cầu gây nên hình dung trực quan phân bố trường từ với mặt quan sát mặt đất Z = Từ số liệu phân bố trường Z = 0, để làm mờ hiệu ứng nguồn nằm nông từ khuếch đại hiệu ứng nguồn nằm sâu , tiến hành nâng trường thực phép xử lý nhằm xác định phân bố trường mặt phẳng song song với mặt quan sát cao Z = h (xa nguồn hơn) Trong phần thứ hai đề tài, thuật toán tính phân phố trường nâng trường xây dựng đồ đẳng trị dị thường từ mức Đề đồ án : Các tham số mô hình toán thuận sau: Giả sử ta khảo sát mạng lưới có diện tích 310x310 m 2, nút lưới theo phương trục x (Bắc), phương trục y (Đông) cách khoảng 10m Trong diện tích khảo sát có hai đối tượng dạng cầu bị từ hóa hoàn toàn cảm ứng tác dụng trường từ T0 = 44000 γ Quả cầu thứ nằm vị trí x = 150 m, y1 = 150 m, độ sâu z1 = 200 m, có bán kính R = 80 m độ từ cảm χ1 = 90000x10-6 theo đơn vị CGS Quả cầu thứ hai nằm vị trí x = 250 m, y2 = 250 m, độ sâu z2 = 100 m, có bán kính R = 40 m độ từ cảm χ2 = 45000x10-6 theo đơn vị CGS Vùng nghiên cứu giả thiết có góc nghiêng I = 150 góc lệch D = 00 Hình vẽ mô tả đề bài: Phần II: Giải đề tài A Nhiệm vụ Tính trường dị thường từ hai cầu gây mặt đất: ∆T(0) Tính trường dị thường từ hai cầu gây mặt z = -0.65m (mặt cao mặt đất 0.65m): ∆T(0.65) Tính đạo hàm thẳng đứng bậc trường dị thường từ hai cầu gây mặt đất So sánh giá trị gradient thẳng đứng trường dị thường từ hai cầu gây cách tính trực tiếp từ số liệu thu mục 1,2 theo công thức cho với đạo hàm bậc mục Gradient = [ ∆T(0) - ∆T(0.65) ]/0.65 Trong trường hợp máy đo trực tiếp gradient thẳng đứng có cần thiết phải tính đạo hàm thẳng đứng bậc không? B Cách tính Cách tính phép biến đổi trường a Công thức tính trường nguồn dạng cầu - Nguồn cho thể khoáng hóa dạng cầu, bị từ hóa hoàn toàn cảm ứng tác dụng trường từ T0, coi vật có từ hóa - Vì trường đạo hàm trường thế, mà có tính chất cộng nên suy trường phải có tính chất cộng Tính: Trường cầu 1: ΔT1 Trường cầu 2: ΔT2 Sau tính tổng trường: ΔT = ΔT1+ ΔT2 Trong đó: - Giả sử có vật từ hóa dạng cầu, bán kính R, cường độ từ hóa J không đổi; góc nghiêng J so với mặt phẳng nằm ngang ghi I, góc lệch J lên mặt phẳng ngang so với trục Ox D - Khi giá trị tuyệt đối ΔH nhỏ nhiều lần so với giá trị tuyệt đối trường địa từ bình thường T0 vùng nghiên cứu, nghĩa : |ΔH | f(x,y): hàm ifft2, kết hợp hàm real Cách tính nâng trường - Nâng trường phép xử lý nhằm xác định phân bố trường mặt phẳng song song với mặt phẳng quan sát cao (Z= 0,65 ( m)) Ở toán này, hiệu ứng nguồn nằm nông bị mờ, khuếch đại hiệu ứng nguồn nằm sâu, ta thấy rõ hiệu ứng nguồn nằm sâu - Có phương pháp nâng trường như: dựa vào khai triển thành chuỗi Poa-xông, phương pháp dựa vào biến đổi Furie - Trong phương pháp biến đổi Furie phương pháp hay áp dụng vì: − Có công thức tính dạng phổ đơn giản thực dể dàng − Dưới dạng phổ phép lọc trường có tính chọn lọc cao − Thông tin không bị trình biến đổi Công thức tính trường theo biến đổi Furie sau: C Kết Bài Để tính dị thường hai cầu gây mặt đất Ta viết code : m=0:10:310; %gia tri truc x n=0:10:310; %gia tri truc y x1=150-m; %toa x qua y1=150-n; %toa y qua z1=200; %gia tri z1 x2=250-m; % toa x qua y2=250-n ; %toa y qua z2=100; % gia tri z2 I=(15*pi)/180; % goc nghieng D=0*pi; % goc lech R1=80; % ban kinh qua cau R2=40; % ban kinh qua cau J1=90000*44000*10^(-6); %cuong tu hoa qua cau J2=44000*45000*10^(-6); %cuong tu hoa qua cau kq=[]; T=zeros(length(m),length(n)); for k=1:length(m) for l=1:length(n) G1=(4/3)*pi*((R1/z1)^3)*J1*((1+((x1(k)/z1)^2)+((y1(l)/z1)^2))^(-5/2)); H1=((cos(I))^2)*((cos(D))^2)*((2*(x1(k)/z1)^2) -((y1(l)/z1)^2) -1); M1=(cos(I))^2*((sin(D))^2)*(2*((y1(l)/z1)^2) -((x1(k)/z1)^2) -1); N1=((sin(I))^2)*(2-((x1(k)/z1)^2) -((y1(l)/z1)^2)); L1=2*((cos(I))^2)*(cos(D))*(sin(D))*((3*x1(k)*y1(l))/(z1^2)); V1=(2*(cos(I))*(sin(I))*(cos(D))*(3*x1(k)/z1)); P1=(2*(cos(I))*(sin(I))*(sin(D))*(3*y1(l)/z1)); T1(k,l)=G1*(H1+M1+N1+L1+V1+P1); % cong thuc tinh truong tu qua G2=(4/3)*pi*((R2/z2)^3)*J2*((1+((x2(k)/z2)^2)+((y2(l)/z2)^2))^(-5/2)); H2=((cos(I))^2)*((cos(D))^2)*((2*(x2(k)/z2)^2) -((y2(l)/z2)^2) -1); M2=(cos(I))^2*((sin(D))^2)*(2*((y2(l)/z2)^2) -((x2(k)/z2)^2) -1); N2=((sin(I))^2)*(2-((x2(k)/z2)^2) -((y2(l)/z2)^2)); L2=2*((cos(I))^2)*(cos(D))*(sin(D))*((3*x2(k)*y2(l))/(z2^2)); V2=(2*(cos(I))*(sin(I))*(cos(D))*(3*x2(k)/z2)); P2=(2*(cos(I))*(sin(I))*(sin(D))*(3*y2(l)/z2)); T2(k,l)=G2*(H2+M2+N2+L2+V2+P2); % cong thuc tinh truong tu qua T(k,l)=T1(k,l)+T2(k,l); % tong truong kq=[kq;n(l) m(k) T(k,l)] end end save 'truongtinhtai0.dat' kq -ascii; Sau cho chạy matlab Rồi từ ta lấy số liệu đưa vào phần mềm surfer để vẽ đồ đẳng trị Kết thu sau: Hình 1: Bản đồ đẳng trị dị thường từ hai khối cầu gây mặt đất Mô tả đồ: Nhìn vào đồ đẳng trị ta thấy : - Giá trị nằm khoảng: (- 1100 : 600) nT - Gồm phần dị thường âm biến đổi từ:-1100 đến 0(nT),nằm hầu hết khu vực phần phía đồ đẳng trị khu vực, chạy dọc theo hướng Đông –Tây; - Và phần có dị thường dương dao động từ: đến 250(nT), nằm phía đồ đẳng trị khu vực, tức chủ yếu phân bố phía Nam đồ; - Ảnh hưởng vị trí cầu đến dị thường: dị thường cầu gây bị dịch chuyển thay đổi vị trí Kích thước cầu lớn dị thường thu rõ hơn, ngược lại Hai cầu càn nằm nông dị thường bật, ngược lại sâu Bài 2: Tính trường dị thường từ hai cầu gây mặt z = -0.65m (mặt cao mặt đất 0.65m): ∆T(0.65) ? Ta dùng code sau : m=0:10:310; %gia tri truc x n=0:10:310; %gia tri truc y x1=150-m; %toa x qua y1=150-n ; %toa y qua z1=200.65; %gia tri z1 x2=250-m; %toa x qua y2=250-n; %toa y qua z2=100.65 ; % gia tri z2 I=(15*pi)/180 % goc nghieng D=0 % goc lech R1=80 % ban kinh qua cau R2=40 % ban kinh qua cau J1=44000*90000*10^(-6) %cuong tu hoa qua cau J2=44000*45000*10^(-6) %cuong tu hoa qua cau kq=[] T=zeros(length(m),length(n)); for k=1:length(m) for l=1:length(n) G1=(4/3)*pi*((R1/z1)^3)*J1*((1+((x1(k)/z1)^2)+((y1(l)/z1)^2))^(-5/2)); H1=((cos(I))^2)*((cos(D))^2)*((2*(x1(k)/z1)^2)-((y1(l)/z1)^2)-1); M1=(cos(I))^2*((sin(D))^2)*(2*((y1(l)/z1)^2)-((x1(k)/z1)^2)-1); N1=((sin(I))^2)*(2-((x1(k)/z1)^2)-((y1(l)/z1)^2)); L1=2*((cos(I))^2)*(cos(D))*(sin(D))*((3*x1(k)*y1(l))/(z1^2)); V1=(2*(cos(I))*(sin(I))*(cos(D))*(3*x1(k)/z1)); P1=(2*(cos(I))*(sin(I))*(sin(D))*(3*y1(l)/z1)); T1(k,l)=G1*(H1+M1+N1+L1+V1+P1); % cong thuc tinh truong tu qua G2=(4/3)*pi*((R2/z2)^3)*J2*((1+((x2(k)/z2)^2)+((y2(l)/z2)^2))^(-5/2)); H2=((cos(I))^2)*((cos(D))^2)*((2*(x2(k)/z2)^2)-((y2(l)/z2)^2)-1); M2=(cos(I))^2*((sin(D))^2)*(2*((y2(l)/z2)^2)-((x2(k)/z2)^2)-1); N2=((sin(I))^2)*(2-((x2(k)/z2)^2)-((y2(l)/z2)^2)); L2=2*((cos(I))^2)*(cos(D))*(sin(D))*((3*x2(k)*y2(l))/(z2^2)); V2=(2*(cos(I))*(sin(I))*(cos(D))*(3*x2(k)/z2)); P2=(2*(cos(I))*(sin(I))*(sin(D))*(3*y2(l)/z2)); T2(k,l)=G2*(H2+M2+N2+L2+V2+P2); % cong thuc tinh truong tu qua T(k,l)=T1(k,l)+T2(k,l); % tong truong kq=[kq;n(l) m(k) T(k,l)] end end save 'nangtruongtai065.dat' kq -ascii; Hình 2: Bản đồ đẳng trị dị thường từ khối cầu Z=0.65 (m) Từ hình ta thấy: Phần dị thường từ thay đổi từ -1050 đến 450(nT): - Gồm phần dị thường âm biến đổi từ:-1050 đến 0(nT),nằm hầu hết khu vực phần phía đồ đẳng trị khu vực, chạy dọc theo hướng Đông Tây - Và phần có dị thường dương dao động từ: đến 450(nT), nằm phía đồ đẳng trị khu vực, tức chủ yếu phân bố phía Nam đồ - Ảnh hưởng vị trí cầu đến dị thường: dị thường cầu gây bị dịch chuyển thay đổi vị trí Kích thước cầu lớn dị thường thu rõ hơn, ngược lại Hai cầu càn nằm nông dị thường bật ngược lại sâu → So sánh giá trị dị thường từ cầu gây mặt đất nâng trường lên độ cao 0.65m: nâng trường lên độ cao 0.65m giá trị dị thường thu thấp mặt đất xa nguồn nâng lên độ cao không lớn nên rõ đồ Khi nâng trường lên ta có thê thấy cầu lớn nông nên có hình dạng bật lên so với dị thường đo mặt đất cầu nhỏ nằm sâu nên bị mờ 10 Bài 3: Tính đạo hàm thẳng đứng bậc trường dị thường từ hai cầu gây mặt đất ? 3.1 Cơ sở phương pháp a Nội dung toán Giả sử trường dị thường quan sát ∆T(x, y, 0) Đạo hàm theo phương thẳng đứng phép xử lý nhằm thu gradien thẳng đứng hay đạo hàm bậc n theo phương thẳng đứng từ phân bố trường dị thường quan sát b Tầm quan trọng Tài liệu đạo hàm theo phương thẳng đứng thể rõ ràng hiệu ứng nguồn nông (càng gần mặt phẳng quan sát trường) bậc đạo hàm cao, sử dụng để tách hiệu ứng gộp, khoanh vị trí nguồn dị thường nằm nông, xóa nhòa dị tường ứng với nguồn nằm sâu Ngoài ra, sử dụng để xác định tham số nguồn số phương pháp giải thích định lượng c Phương pháp tính đạo hàm theo phương thẳng đứng dựa biến đổi Fourier Phương pháp quy trình tính đạo hàm thẳng đứng bậc n dựa biến đổi Fourier tương tự cách tính nâng trường hạ trường, cần thay đặc trưng tần số H(ω) phép biến đổi thành: H(ω) = tính đạo hàm thẳng đứng bậc n mặt quan sát H(ω) = tính đạo hàm thẳng đứng bậc n độ cao h 3.2 Cách tính kết Nhập code vào phần mềm Matlab chạy số liệu đồ đẳng trị gradient thẳng đứng khối cầu vị trí Z=0 m=0:10:310; %gia tri truc x n=0:10:310; %gia tri truc y x1=150-m y1=150-n z1=200 %gia tri z1 x2=250-m y2=250-n z2=100 % gia tri z2 I=(15*pi)/180 % goc nghieng D=0*pi % goc lech R1=80 % ban kinh qua cau R2=40 % ban kinh qua cau J1=44000*90000*10^(-6) %cuong tu hoa qua cau 11 J2=44000*45000*10^(-6) %cuong tu hoa qua cau T=zeros(length(m),length(n)); for k=1:length(m) for l=1:length(n) G1=(4/3)*pi*((R1/z1)^3)*J1*((1+((x1(k)/z1)^2)+((y1(l)/z1)^2))^(-5/2)); H1=((cos(I))^2)*((cos(D))^2)*((2*(x1(k)/z1)^2)-((y1(l)/z1)^2)-1); M1=(cos(I))^2*((sin(D))^2)*(2*((y1(l)/z1)^2)-((x1(k)/z1)^2)-1); N1=((sin(I))^2)*(2-((x1(k)/z1)^2)-((y1(l)/z1)^2)); L1=2*((cos(I))^2)*(cos(D))*(sin(D))*((3*x1(k)*y1(l))/(z1^2)); V1=(2*(cos(I))*(sin(I))*(cos(D))*(3*x1(k)/z1)); P1=(2*(cos(I))*(sin(I))*(sin(D))*(3*y1(l)/z1)); T1(k,l)=G1*(H1+M1+N1+L1+V1+P1); G2=(4/3)*pi*((R2/z2)^3)*J2*((1+((x2(k)/z2)^2)+((y2(l)/z2)^2))^(-5/2)); H2=((cos(I))^2)*((cos(D))^2)*((2*(x2(k)/z2)^2)-((y2(l)/z2)^2)-1); M2=(cos(I))^2*((sin(D))^2)*(2*((y2(l)/z2)^2)-((x2(k)/z2)^2)-1); N2=((sin(I))^2)*(2-((x2(k)/z2)^2)-((y2(l)/z2)^2)); L2=2*((cos(I))^2)*(cos(D))*(sin(D))*((3*x2(k)*y2(l))/(z2^2)); V2=(2*(cos(I))*(sin(I))*(cos(D))*(3*x2(k)/z2)); P2=(2*(cos(I))*(sin(I))*(sin(D))*(3*y2(l)/z2)); T2(k,l)=G2*(H2+M2+N2+L2+V2+P2); T(k,l)=T1(k,l)+T2(k,l); % tong truong c?a qua W1(k,l)=(2*pi*k)/(31*10); W2(k,l)=(2*pi*l)/(31*10); W3(k,l)=(((W1(k,l))^2)+((W2(k,l))^2))^(1/2); % ham loc cua end end T_fou=fft2(T); %chuyen sang tan so, T_dh=(T_fou).*(W3(k,l)); T_h= ifft2(T_dh); % furie nguoc T_dhtd=real(T_h); % chon phan thuc kq=[] for k=1:length(m) for l=1:length(n) kq=[kq;n(l) m(k) T_dhtd(k,l) ] end 12 end save 'gradientinhtai0.dat' kq -ascii; Hình 3: Bản đồ đạo hàm bậc mặt đất Mô tả đồ: Nhìn vào đồ đẳng trị ta nhận thấy: - Giá trị ∆� (0) nằm khoảng : (- 950 - 450) nT - Phần giá trị dị thường âm biến đổi từ: -950 đến (nT), nằm hầu hết khu vực phần phía đồ, vùng trung tâm kéo dài có giá trị dị thường nhỏ lớn dần xung quanh Vùng kéo dài gần chạy dọc theo hướng ĐB-TN Phần có giá trị dị thường dương dao động từ: đến 450 (nT), nằm phía đồ chạy khu vực phía Nam có giá trị dị thường cao giảm dần xung quanh - Bài 4: So sánh giá trị gradient thẳng đứng trường dị thường từ hai cầu gây cách tính trực tiếp từ số liệu thu mục 1,2 theo công thức cho với đạo hàm bậc mục Gradient = [ ∆T(0) - ∆T(0.65) ]/0.65 Trong trường hợp máy đo trực tiếp gradient thẳng đứng có cần thiết phải tính đạo hàm thẳng đứng bậc không? 13 4.1 Tính gradient m=0:10:310; %gia tri truc x n=0:10:310; %gia tri truc y x1=160-m y1=80-n z11=110 %gia tri z1 z12=110.65 x2=180m y2=220-n z21=80 z22=80.65 % gia tri z2 I=(15*pi)/180 % goc nghieng D=0*pi % goc lech R1=80 % ban kinh qua cau R2=40 % ban kinh qua cau J1=44000*70000*10^(-6) %cuong tu hoa qua cau J2=44000*90000*10^(-6) %cuong tu hoa qua cau kq=[] T=zeros(length(m),leng th(n)); for k=1:length(m) for l=1:length(n) G1=(4/3)*pi*((R1/z11)^3)*J1*((1+((x1(k)/z11)^2)+((y1(l)/z11)^2))^(5/2)); H1=((cos(I))^2)*((cos(D))^2)*((2*(x1(k)/z11)^2)-((y1(l)/z11)^2)-1); M1=(cos(I))^2*((sin(D))^2)*(2*((y1(l)/z11)^2)-((x1(k)/z11)^2)-1); N1=((sin(I))^2)*(2-((x1(k)/z11)^2)-((y1(l)/z11)^2)); L1=2*((cos(I))^2)*(cos(D))*(sin(D))*((3*x1(k)*y1(l))/(z11^2)); V1=(2*(cos(I))*(sin(I))*(cos(D))*(3*x1(k)/z11)); P1=(2*(cos(I))*(sin(I))*(sin(D))*(3*y1(l)/z11)); T1(k,l)=G1*(H1+M1+N1+L1+V1+P1); G2=(4/3)*pi*((R2/z21)^3)*J2*((1+((x2(k)/z21)^2)+((y2(l)/z21)^2))^(5/2)); H2=((cos(I))^2)*((cos(D))^2)*((2*(x2(k)/z21)^2)-((y2(l)/z21)^2)-1); M2=(cos(I))^2*((sin(D))^2)*(2*((y2(l)/z21)^2)-((x2(k)/z21)^2)-1); N2=((sin(I))^2)*(2-((x2(k)/z21)^2)-((y2(l)/z21)^2)); L2=2*((cos(I))^2)*(cos(D))*(sin(D))*((3*x2(k)*y2(l))/(z21^2)); 14 V2=(2*(cos(I))*(sin(I))*(cos(D))*(3*x2(k)/z21)); P2=(2*(cos(I))*(sin(I))*(sin(D))*(3*y2(l)/z21)); T2(k,l)=G2*(H2+M2+N2+L2+V2+P2); T(k,l)=T1(k,l)+T2(k,l); % truong tinh G12=(4/3)*pi*((R1/z12)^3)*J1*((1+((x1(k)/z12)^2)+((y1(l)/z12)^2))^(5/2)); H12=((cos(I))^2)*((cos(D))^2)*((2*(x1(k)/z12)^2)-((y1(l)/z12)^2)-1); M12=(cos(I))^2*((sin(D))^2)*(2*((y1(l)/z12)^2)-((x1(k)/z12)^2)-1); N12=((sin(I))^2)*(2-((x1(k)/z12)^2)-((y1(l)/z12)^2)); L12=2*((cos(I))^2)*(cos(D))*(sin(D))*((3*x1(k)*y1(l))/(z12^2)); V12=(2*(cos(I))*(sin(I))*(cos(D))*(3*x1(k)/z12)); P12=(2*(cos(I))*(sin(I))*(sin(D))*(3*y1(l)/z12)); T12(k,l)=G12*(H12+M12+N12+L12+V12+P12); G22=(4/3)*pi*((R2/z22)^3)*J2*((1+((x2(k)/z22)^2)+((y2(l)/z22)^2))^(5/2)); H22=((cos(I))^2)*((cos(D))^2)*((2*(x2(k)/z22)^2)-((y2(l)/z22)^2)-1); M22=(cos(I))^2*((sin(D))^2)*(2*((y2(l)/z22)^2)-((x2(k)/z22)^2)-1); N22=((sin(I))^2)*(2-((x2(k)/z22)^2)-((y2(l)/z22)^2)); L22=2*((cos(I))^2)*(cos(D))*(sin(D))*((3*x2(k)*y2(l))/(z22^2)); V22=(2*(cos(I))*(sin(I))*(cos(D))*(3*x2(k)/z22)); P22=(2*(cos(I))*(sin(I))*(sin(D))*(3*y2(l)/z22)); T22(k,l)=G22*(H22+M22+N22+L22+V22+P22); V(k,l)=T12(k,l)+T22(k,l); %truong P(k,l)= (T(k,l)-V(k,l))./0.65; %gradien kq=[kq;n(l) m(k) P(k,l)] end end save 'gradiendo.dat' kq -ascii; 15 Hình Bản đồ đẳng trị gradient 4.2 So sánh giá trị gradient thẳng đứng trường dị thường từ hai cầu gây (hình 4) với đạo hàm bậc (hình 3) : - Ta thấy đồ giống hình dáng lại khác nhiều giá trị : giá trị đạo hàm bậc mặt đất gấp khoảng 50 lần giá trị gradient thẳng đứng khối cầu gây - Mỗi đồ có vùng dị thường vùng dị thường âm vùng dị thường dương, dị thường âm phân bố phía dị thường dương phân bố phía : + Bản đồ hình phần dị thường âm có biên độ dao động từ -950 – (nT) hình dao động khoảng từ -18 đến (nT) + Dị thường dương hình dao động khoảng từ đến 450 (nT) hình lại dao động khoảng từ đến (nT) - Trong trường hợp máy đo trực tiếp gradient thẳng đứng cần phải tính đạo hàm thẳng đứng bậc 1: Tài liệu đạo hàm theo phương thẳng đứng thể rõ ràng hiệu ứng nguồn nông (càng gần mặt phẳng quan sát trường) bậc đạo hàm cao, sử dụng để tách hiệu ứng gộp, khoanh vị trí nguồn dị thường nằm nông, xóa nhòa dị tường ứng với nguồn nằm sâu Ngoài ra, sử dụng để xác định tham số nguồn số phương pháp giải thích định lượng 4.3 Kết luận 16 Qua kết thu nhận thấy hai đồ tính trường dị thường từ đồ gradient tính toán có hình dạng giống là: dị thường âm nằm phía dị thường dương nằm phía Khi nâng trường từ độ cao h=0 lên độ cao h=0.65(m) hiệu ứng nguồn gây dị thường từ bị mờ phản ánh đối tượng nông có diện tích hẹp Qua đồ gradient tính toán phản ánh rõ đối tượng nằm bên mặt đất dị thường từ mà đối tượng gây PHẦN III: LỜI CẢM ƠN Trong trình nghiên cứu chúng em áp dụng kiến thức Thăm dò từ sở giáo trình “Phương Pháp Từ” KS Bùi Thế Bình với số tài liệu tham khảo khác Nhờ có đề tài giúp chúng em củng cố kiến thức thiếu sót học cách sử dung ngôn ngữ lập trình ứng dụng toán Thăm dò từ Với cố gắng chúng em chuyển đến bạn đọc nhìn chi tiết đề tài Tuy nhiên, trình độ chuyên môn hạn chế trình nghiên cứu tránh khỏi thiếu sót Chúng em mong nhận nhận xét, đóng góp ý kiến từ phía thầy cô toàn thể bạn để đề tài hoàn thiện Đồ án có giúp đỡ nhiệt tình thầy cô giáo môn địa vật lý đặc biệt giúp đỡ thầy giáo KS Phạm Ngọc Kiên chúng em cám ơn giúp đỡ tận tình thầy cô Chúng em xin chân thành cảm ơn! Hà Nội, ngày 26-05-2016 17 [...]... thấy hai bản đồ tính trường dị thường từ và bản đồ gradient tính toán có hình dạng giống nhau là: các dị thường âm nằm ở phía trên và dị thường dương nằm ở phía dưới Khi nâng trường từ độ cao h=0 lên độ cao h=0.65(m) thì hiệu ứng nguồn gây dị thường từ bị mờ đi phản ánh các đối tượng nông và có diện tích hẹp Qua bản đồ gradient tính toán thì phản ánh rõ đối tượng nằm bên dưới mặt đất và dị thường từ... 15 Hình 4 Bản đồ đẳng trị gradient 4.2 So sánh giá trị gradient thẳng đứng của trường dị thường từ do hai quả cầu gây ra (hình 4) với đạo hàm bậc nhất (hình 3) : - Ta thấy 2 bản đồ khá giống nhau về hình dáng nhưng lại khác nhau rất nhiều về giá trị : giá trị của đạo hàm bậc nhất tại mặt đất gấp khoảng 50 lần giá trị gradient thẳng đứng của 2 khối cầu gây ra - Mỗi bản đồ đều có 2 vùng dị thường chính... đó là vùng dị thường âm và vùng dị thường dương, dị thường âm phân bố phía trên và dị thường dương phân bố phía dưới : + Bản đồ ở hình 3 phần dị thường âm có biên độ dao động từ -950 – 0 (nT) còn ở hình 4 thì nó chỉ dao động trong khoảng từ -18 đến 0 (nT) + Dị thường dương ở hình 3 dao động trong khoảng từ 0 đến 450 (nT) còn hình 4 thì lại chỉ dao động trong khoảng từ 0 đến 9 (nT) - Trong trường hợp... save 'gradientinhtai0.dat' kq -ascii; Hình 3: Bản đồ đạo hàm bậc nhất tại mặt đất Mô tả bản đồ: Nhìn vào bản đồ đẳng trị ta nhận thấy: - Giá trị ∆� (0) nằm trong khoảng : (- 950 - 450) nT - Phần trên là giá trị dị thường âm biến đổi từ: -950 đến 0 (nT), nằm ở hầu hết khu vực ở phần phía trên của bản đồ, vùng trung tâm kéo dài sẽ có giá trị dị thường nhỏ nhất và lớn dần ra xung quanh Vùng kéo dài gần... n dựa trên biến đổi Fourier là tương tự như cách tính nâng trường và hạ trường, chỉ cần thay đặc trưng tần số H(ω) của phép biến đổi thành: H(ω) = khi tính đạo hàm thẳng đứng bậc n tại mặt quan sát H(ω) = khi tính đạo hàm thẳng đứng bậc n tại độ cao h 3.2 Cách tính và kết quả Nhập code dưới đây vào phần mềm Matlab sẽ chạy ra được số liệu bản đồ đẳng trị gradient thẳng đứng của 2 khối cầu ở vị trí Z=0... quanh Vùng kéo dài gần như chạy dọc theo hướng ĐB-TN Phần dưới có giá trị dị thường dương dao động từ: 0 đến 450 (nT), nằm ở phía dưới của bản đồ và chạy khu vực phía Nam sẽ có giá trị dị thường cao nhất rồi giảm dần ra xung quanh - Bài 4: So sánh giá trị gradient thẳng đứng của trường dị thường từ do hai quả cầu gây ra bằng cách tính trực tiếp từ các số liệu thu được ở mục 1,2 theo công thức cho dưới...Bài 3: Tính đạo hàm thẳng đứng bậc nhất của trường dị thường từ do hai quả cầu gây ra tại mặt đất ? 3.1 Cơ sở phương pháp a Nội dung bài toán Giả sử trường dị thường quan sát được là ∆T(x, y, 0) Đạo hàm theo phương thẳng đứng là phép xử lý nhằm thu được gradien thẳng đứng hay đạo hàm bậc n theo phương thẳng đứng từ sự phân bố trường dị thường quan sát được b Tầm quan trọng... hạn chế trong quá trình nghiên cứu không thể tránh khỏi những thiếu sót Chúng em rất mong nhận được sự nhận xét, đóng góp ý kiến từ phía các thầy cô và toàn thể các bạn để đề tài trên được hoàn thiện hơn Đồ án có sự giúp đỡ nhiệt tình của các thầy cô giáo trong bộ môn địa vật lý và đặc biệt là sự giúp đỡ của thầy giáo KS Phạm Ngọc Kiên chúng em rất cám ơn sự giúp đỡ tận tình của thầy cô Chúng em xin... mặt phẳng quan sát trường) khi bậc đạo hàm càng cao, do vậy nó được sử dụng để tách hiệu ứng gộp, khoanh vị trí nguồn dị thường nằm nông, xóa nhòa những dị tường ứng với nguồn nằm sâu Ngoài ra, nó còn được sử dụng để xác định tham số nguồn trong một số phương pháp giải thích định lượng c Phương pháp tính đạo hàm theo phương thẳng đứng dựa trên biến đổi Fourier Phương pháp và quy trình tính đạo hàm thẳng... máy có thể đo trực tiếp gradient thẳng đứng thì vẫn cần phải tính đạo hàm thẳng đứng bậc 1: Tài liệu đạo hàm theo phương thẳng đứng thể hiện rõ ràng hơn hiệu ứng của nguồn càng nông (càng gần mặt phẳng quan sát trường) khi bậc đạo hàm càng cao, do vậy nó được sử dụng để tách hiệu ứng gộp, khoanh vị trí nguồn dị thường nằm nông, xóa nhòa những dị tường ứng với nguồn nằm sâu Ngoài ra, nó còn được sử dụng