CHƯƠNG II: MÔ HÌNH HÓA VÀ CÁC KẾT QUẢ TÍNH TOÁN
II.1. Áp dụng tín hiệu giải tích xác định biên của nguồn gây dị thường từ
II.1.1. Trường hợp nguồn là 1 lăng trụ bị từ hóa
a) Các thông số mô hình
Mô hình đầu tiên chỉ bao gồm một lăng trụ được phân bố rất gần mặt quan sát với các thông số như sau:
- Độ sâu so với đỉnh là 3,5 km.
- Kích thước các cạnh lăng trụ là 30x30x2 km3. - Số điểm quan sát trên tuyến: 128 x 128.
- Khoảng cách giữa các điểm quan sát 0.5/63.5 km - Tọa độ tâm trên mặt phẳng Oxy là (31.5, 31.5) km - Độ từ thiên của Trái đất và vecto từ hóa: -1
- Cường độ từ hóa: 5 - Góc quay: 45
Với các thông số này, theo công thức (2.6), chúng tôi xác định được dị thường từ của lăng trụ với các góc nghiêng I = 90 , I = 30 , I = 45 và I = 1
b) Kết quả tính toán
14 Góc nghiêng từ hóa I = 900
Hình 2.1 - Kết quả tính tín hiệu giải tích của mô hình 1 lăng trụ với I = 90 a) Dị thường từ 2D b) Dị thường từ 3D
c) Biên độ tín hiệu giải tích 2D d) Biên độ tín hiệu giải tích 3D
x (km)
y (km)
0 10 20 30 40 50 60
0 10 20 30 40 50 60
nT
-100 0 100 200 300 400
0 10 20 30 40 50 60
0 10 20 30 40 50 60
nT/km
20 40 60 80 100 120 140 160 180 200 220
0
50
100
150
0 50 100 150
0 50 100 150 200 250
0 20 40 60 80 100 120 140
0 50 100 150 -200 -100 0 100 200 300 400 500
a) b)
c) d)
15 Góc nghiêng từ hóa I = 450
Hình 2.2 - Kết quả tính tín hiệu giải tích của mô hình 1 lăng trụ với I = 45 a) Dị thường từ 2D b) Dị thường từ 3D
c) Biên độ tín hiệu giải tích 2D d) Biên độ tín hiệu giải tích 3D
0 20 40 60 80 100 120 140
0 50 100 150 -400 -200 0 200 400 600 800
x (km)
y (km)
0 10 20 30 40 50 60
0 10 20 30 40 50 60
nT
-300 -200 -100 0 100 200 300 400 500
0 10 20 30 40 50 60
0 10 20 30 40 50 60
nT/km
50 100 150 200
0 20 40 60 80 100 120 140
0 50 100 150
0 50 100 150 200 250 300
a) b)
c) d)
16 Góc nghiêng từ hóa I = 300
x (km)
y (km)
0 10 20 30 40 50 60
0 10 20 30 40 50 60
nT
-300 -200 -100 0 100 200 300 400 500
0 20 40 60 80 100 120 140
0 50 100 150 -400 -200 0 200 400 600
0 10 20 30 40 50 60
0 10 20 30 40 50 60
nT/km
20 40 60 80 100 120 140 160 180 200
0 20 40 60 80 100 120 140
0 50 100 150
0 50 100 150 200 250
Hình 2.3 - Kết quả tính tín hiệu giải tích của mô hình 1 lăng trụ với I = 30 a) Dị thường từ 2D b) Dị thường từ 3D
c) Biên độ tín hiệu giải tích 2D d) Biên độ tín hiệu giải tích 3D
a) b)
c) d)
17 Góc nghiêng từ hóa I = 10
Hình 2.4 - Kết quả tính tín hiệu giải tích của mô hình 1 lăng trụ với I = 1 a) Dị thường từ 2D b) Dị thường từ 3D
c) Biên độ tín hiệu giải tích 2D d) Biên độ tín hiệu giải tích 3D
x (km)
y (km)
0 10 20 30 40 50 60
0 10 20 30 40 50 60
nT
-300 -200 -100 0 100 200
0 20 40 60 80 100 120 140
0 50 100 150 -400 -300 -200 -100 0 100 200 300
0 10 20 30 40 50 60
0 10 20 30 40 50 60
nT/km
20 40 60 80 100 120 140
0 20 40 60 80 100 120 140
0 50 100 150
0 20 40 60 80 100 120 140 160
a) b)
c) d)
18 II.1.2. Mô hình 3 vật thể dạng lăng trụ
II.1.2.1. Trường hợp 3 lăng trụ có độ sâu như nhau a) Các thông số mô hình
Mô hình tiếp theo bao gồm ba lăng trụ được phân bố gần mặt quan sát với các thông số như sau:
- Độ sâu tới mặt trên của cả 3 lăng trụ: 0,5 km
- Kích thước theo trục x và trục y của lăng trụ là (15, 15); (10, 40); (15, 15) .
- Số điểm quan sát trên tuyến: 128 x 128.
- Khoảng cách giữa các điểm quan sát 0.5/63.5 km
- Tọa độ tâm trên mặt phẳng Oxy là (31.5, 12); (31.5, 31.5); (31.5, 51.5) km - Độ từ thiên của Trái đất và vecto từ hóa: -1
- Cường độ từ hóa: 5, 5, 5.
- Góc quay: 45 , 0 45 . b) Kết quả tính toán
Với các thông số đã đưa ra, theo công thức (2.6), chúng tôi xác định được dị thường từ của 3 lăng trụ này lần lượt với các góc nghiêng từ hóa : I = 900 (hình 2.5a,b), I = 300 (hình 2.6a,b) , I = 50 (hình 2.7a,b) sau đó thực hiện việc tính tín hiệu giải tích. Kết quả tính toán tương ứng được đưa ra trên các hình vẽ : hình 2.5c,d;
hình 2.6c,d và hình 2.7c,d.
19 Góc nghiêng từ hóa I = 900
x (km)
y (km)
0 10 20 30 40 50 60
0 10 20 30 40 50 60
nT
-200 -100 0 100 200 300 400 500
0 20 40 60 80 100 120 140
0 100 200 -400 -200 0 200 400 600
x (km)
y (km)
0 10 20 30 40 50 60
0 10 20 30 40 50 60
nT/km
100 200 300 400 500 600 700 800 900 1000
0
50
100
150
0 50 100 150
0 200 400 600 800 1000 1200
Hình 2.5 - Kết quả tính tín hiệu giải tích mô hình 3 lăng trụ cùng độ sâu với I = 90 a) Dị thường từ 2D b) Dị thường từ 3D
c) Biên độ tín hiệu giải tích 2D d) Biên độ tín hiệu giải tích 3D
a) b)
c) C) d)
20 Góc nghiêng từ hóa I = 300
0 20 40 60 80 100 120 140
0 100 200 -600 -400 -200 0 200 400 600 800
x (km)
y (km)
0 10 20 30 40 50 60
0 10 20 30 40 50 60
nT
-400 -200 0 200 400 600
x (km)
y (km)
0 10 20 30 40 50 60
0 10 20 30 40 50 60
nT/km
100 200 300 400 500 600 700 800 900
0 20 40 60 80 100 120 140
0 50 100 150
0 200 400 600 800 1000 1200
a) b)
c) d)
Hình 2.6 - Kết quả tính tín hiệu giải tích mô hình 3 lăng trụ cùng độ sâu với I = 30 a) Dị thường từ 2D b) Dị thường từ 3D
c) Biên độ tín hiệu giải tích 2D d) Biên độ tín hiệu giải tích 3D
x (km)
y (km)
0 10 20 30 40 50 60
0 10 20 30 40 50 60
nT
-400 -200 0 200 400 600
a) b)
c) d)
21 Góc nghiêng từ hóa I = 50
x (km)
y (km)
0 10 20 30 40 50 60
0 10 20 30 40 50 60
nT
-400 -300 -200 -100 0 100 200 300
0 20 40 60 80 100 120 140
0 50 100 150 -600 -400 -200 0 200 400
x (km)
y (km)
0 10 20 30 40 50 60
0 10 20 30 40 50 60
nT/km
100 200 300 400 500 600
0 20 40 60 80 100 120 140
0 50 100 150
0 200 400 600 800
Hình 2.7 - Kết quả tính tín hiệu giải tích mô hình 3 lăng trụ cùng độ sâu với I = 5 a) Dị thường từ 2D b) Dị thường từ 3D
c) Biên độ tín hiệu giải tích 2D d) Biên độ tín hiệu giải tích 3D
a) b)
c) d)
22
II.1.2.2. Trường hợp 3 lăng trụ có độ sâu khác nhau a) Các thông số mô hình
Mô hình tiếp theo bao gồm ba lăng trụ được phân bố gần mặt quan sát với các thông số như sau:
- Độ sâu tới mặt trên của cả 3 lăng trụ lần lượt là: 0,2 km; 0,5 km và 1km - Kích thước theo trục x và trục y của lăng trụ là (15, 15); (10, 40); (15, 15) .
- Số điểm quan sát trên tuyến: 128 x 128.
- Khoảng cách giữa các điểm quan sát 0.5/63.5 km
- Tọa độ tâm trên mặt phẳng là (31.5, 12); (31.5, 31.5); (31.5, 51.5) km - Độ từ thiên của Trái đất và vecto từ hóa: -1
- Cường độ từ hóa: 5, 5, 5 (A/m) - Góc quay: 45 , 0 , 45 .
b) Kết quả tính toán
Với các thông số đã đưa ra, theo công thức (2.6), chúng tôi xác định được dị thường từ của 3 lăng trụ này lần lượt với các góc nghiêng từ hóa : I = 900 (hình 2.8a,b), I = 300 (hình 2.9a,b) , I = 50 (hình 2.10a,b) sau đó thực hiện việc tính tín hiệu giải tích. Kết quả tính toán tương ứng được đưa ra trên các hình vẽ: hình 2.8c,d; hình 2.9c,d và hình 2.10c,d.
23 Góc nghiêng từ hóa I = 900
x (km)
y (km)
0 10 20 30 40 50 60
0 10 20 30 40 50 60
nT
-600 -400 -200 0 200 400 600 800 1000 1200 1400
0
50
100
150
0 50 100 150 -1000 -500 0 500 1000 1500 2000
x (km)
y (km)
0 10 20 30 40 50 60
0 10 20 30 40 50 60
nT/km
500 1000 1500 2000 2500 3000
0 20 40 60 80 100 120 140 0
50 100
150
0 1000 2000 3000 4000
Hình 2.8 - Kết quả tính tín hiệu giải tích mô hình 3 lăng trụ khác độ sâu với I = 90 a) Dị thường từ 2D b) Dị thường từ 3D
c) Biên độ tín hiệu giải tích 2D d) Biên độ tín hiệu giải tích 3D
a) b)
c) d)
24 Góc nghiêng từ hóa I = 300
x (km)
y (km)
0 10 20 30 40 50 60
0 10 20 30 40 50 60
nT
-1500 -1000 -500 0 500 1000 1500 2000
0 20 40 60 80 100 120 140
0 100 200 -2000 -1500 -1000 -500 0 500 1000 1500 2000 2500
x (km)
y (km)
0 10 20 30 40 50 60
0 10 20 30 40 50 60
nT/km
500 1000 1500 2000 2500 3000 3500
0 20 40 60 80 100 120 140 0
50 100
150
0 1000 2000 3000 4000 5000
Hình 2.9 - Kết quả tính tín hiệu giải tích mô hình 3 lăng trụ khác độ sâu với I = 30 a) Dị thường từ 2D b) Dị thường từ 3D
c) Biên độ tín hiệu giải tích 2D d) Biên độ tín hiệu giải tích 3D
a) b)
c) d)
25 Góc nghiêng từ hóa I = 50
x (km)
y (km)
0 10 20 30 40 50 60
0 10 20 30 40 50 60
nT
-1000 -500 0 500 1000
0 20 40 60 80 100 120 140
0 100 200 -2000 -1500 -1000 -500 0 500 1000 1500
x (km)
y (km)
0 10 20 30 40 50 60
0 10 20 30 40 50 60
nT/km
500 1000 1500 2000 2500
0 20 40 60 80 100 120 140 0
50 100
150
0 500 1000 1500 2000 2500 3000
Hình 2.10 - Kết quả tính tín hiệu giải tích mô hình 3 lăng trụ khác độ sâu với I = 5 a) Dị thường từ 2D b) Dị thường từ 3D
c) Biên độ tín hiệu giải tích 2D d) Biên độ tín hiệu giải tích 3D
a) b)
c) d)
26