TÓM TẮT LUẬN VĂN THẠC SĨ PHƯƠNG PHÁP BIẾN ĐỔI TRƯỜNG TRỌNG LỰC VÀ VIỆC ÁP DỤNG CHÚNG CHO KHU VỰC X THUỘC THỀM LỤC ĐỊA VIỆT NAM

CHƯƠNG CÁC PHÉP BIẾN ĐỔI TRƯỜNG TRONG MIỀN KHÔNG GIAN Các dị thường trọng lực quan sát phản ánh toàn hiệu ứng trọng lực yếu tố địa chất gây Trong trường tổng cộng yếu tố địa chất có đóng góp phần định Vì vậy, giải nhiệm vụ địa chất cụ thể, từ trường tổng phải tách thành phần trường riêng biệt có liên hệ trực tiếp đến đối tượng cần nghiên cứu Muốn vậy, người ta phải tiến hành biến đổi trường quan sát nhằm nhấn mạnh thành phần trường cần thiết (được coi phần hữu ích) làm yếu thành phần khác (được coi nhiễu) Như vậy, phép biến đổi trường dị thường trọng lực có điểm chung phép lọc nhiễu, phân tách tín hiệu lý thuyết truyền tin Mục đích phép biến đổi trường trọng lực (hoặc từ) tách trường quan sát thành thành phần tương ứng với đối tượng địa chất nằm độ sâu khác Dưới dạng toán học, tất phép biến đổi biểu diễn công thức sau [3]: Vbđ(x0,y0,z0)= ∫∫ V xp (ξ ,η ,0) K ( x0 − ξ , y − η , z 0) dξdη (1.1) Trong trường hợp toán ba chiều, và: Vbđ(x0,z0)= ∫∫ V xp (ξ ,0) K ( x0 − ξ , z 0)dξ Trong trường hợp tốn hai chiều, đó: (1.2) + Vbd ( x0, y 0, z 0) Vbd ( x0, z 0) hàm số biến đổi + V xp (ξ ,η ,0) V xp (ξ ,0) hàm số xuất phát (trường tổng) + K ( x0 − ξ , y0 − η , z 0) K ( x0 − ξ , z 0) nhân biến đổi (đơi cịn gọi hàm trọng số) Vì K ( x0 − ξ , y0 − η , z 0) K ( x0 − ξ , z 0) thường toán tử tuyến tính nên tất biến đổi tương ứng gọi biến đổi tuyến tính Phép biến đổi trường trọng lực từ miền không gian chia làm ba nhóm chính: + Trung bình hố + Tiếp tục giải tích dị thường trọng lực (xem hàm điều hồ) + Tính đạo hàm bậc cao trọng lực Chúng ta xét đến nhóm phương pháp 1.1 Phương pháp trung bình hóa Việc phân chia dị thường trọng lực thành thành phần khu vực địa phương nhờ phương pháp trung bình hố sử dụng rộng rãi thực tế Bản chất phương pháp trung bình hố sau: Xem trường trọng lực quan sát gồm hai thành phần, thành phần khu vực Vr thành phần địa phương Vl V = V r + Vl (1.3) Lấy trung bình trường quan sát phạm vi đường trịn bán kính R Giá trị trung bình biểu diễn tích phân sau [1,2,3]: V (0,0,0)= πR 2π ∞ ∫ ∫V (r,α ,0)rdrdα (1.4) 0 Bán kính R chọn cho lớn nhiều so với kích thước dị thường địa phương nhỏ nhiều so với kích thước dị thường khu vực Khi thoả mãn điều kiện thành phần khu vực tách riêng từ trường quan sát Do dị thường địa phương âm dương bù trừ lẫn thành phần khu vực bị thay đổi Do V ≈ Vr , trường hợp đặc biệt trường khu vực thay đổi theo quy luật tuyến tính hồn tồn khơng bị thay đổi lấy trung bình, tức: V (0,0,0) = Vr(0,0,0) (1.5) Sau tính trường khu vực Vr, trường dị thường địa phương tính theo cơng thức: Vl = V - V (1.6) Để làm sáng tỏ ý nghĩa vật lý phương pháp trung bình hố, người ta đưa vào khái niệm mức độ trung bình hố, tỷ số trường trung bình hố trường xuất phát ε= V V (1.7) Mức độ trung bình hố đồng thời đặc trưng cho mức độ xác việc tách trường địa phương Trong phương pháp trung bình hố, ngồi cách lấy trung bình theo vịng trịn người ta cịn lấy trung bình theo hình khác Một hình hay dùng hình vng , nhờ có Pa-lét vng mà khối lượng phép tính giảm nhiều Phương pháp trung bình hố phép biến đổi trường miền không gian thường thực Pa-lét Người ta đưa tiêu chuẩn để đánh giá khả lọc Pa-lét qua việc đánh giá độ sâu Đó q trình theo dõi biến đổi dị thường theo chiều sâu đơn vị nguồn điểm nằm độ sâu Z chứa toàn nguồn dị thường cần tách Trong phương pháp trung bình hố, để đánh giá độ sâu người ta thường đưa vào đại lượng gọi đại lượng đặc trưng tương đối ký hiệu N(z) Đại lượng định nghĩa sau: Đặc trưng độ sâu tương đối N(z) tỷ số dị thường trọng lực biến đổi chưa biến đổi M ( z) N(z) = Mbt ( z ) Nhờ biểu thức ta biết vật thể độ sâu Z sau phép biến đổi dị thường biến đổi M(z) - Dị thường trọng lực sau biến đổi Mbt(z) - Dị thường trọng lực chưa biến đổi Các công thức Andrêep Klusin xây dựng cơng thức tính sau: ∞ ∞ ∞ 0 −αz Mbt(z) = ∫ αdα ∫ ∆ g ( ρ ) J (αρ ) ρdρ = ∫ e αdα = z2 (1.8) + ∆g ( ρ ) giá trị trung bình trường dị thường quan sát đường trịn bán kính ρ nhận sau tính tốn với giá trị đọc điểm nút Palét + J (α , ρ ) hàm Bessel loại cấp 0, khác với hàm J 0(x) J1(x) giống khác biệt COS (αx) SIN (αx) với COS(x) SIN(x) + α đóng vai trị tần số vịng trường hợp hàm điều hồ ∞ Còn M(z) = J (αR ) −αR ∫ αR e = Z + R (Z + Z + R ) (1.9) Thay M(z) Mbt(z) vào N(z) ta có: N(z) = 2Z Z + R (Z + Z + R ) (1.10) với : Z- Độ sâu đến vật thể gây dị thường R- Bán kính trung bình hố Theo cơng thức ta thấy N(z) hàm phụ thuộc vào độ sâu nằm Z Khảo sát hàm N(z) thấy Z → N(z) → Z → ∞ N(z) → Nhìn đồ thị ta thấy dị thường gây vật thể nằm độ sâu khoảng lấy trung bình Z = 2R độ sâu thay đổi N( z r ) bắt đầu tiệm cận với N( z r ) = từ z r = Có thể chọn Z=2R làm độ sâu nghiên cứu R bán kính trung bình hố tối ưu Có thể xác định R theo cách trình bày trên, biết R tìm Z 1.2 Phương pháp tiếp tục giải tích trường Cơ sở phương pháp tiếp tục giải tích trường dị thường trọng lực từ là: Hàm xem hàm điều hoà Phương pháp tiếp tục giải tích dị thường trọng lực từ sử dụng rộng rãi để tách dị thường mà đơi cịn sử dụng để xác định thông số vật thể gây nên dị thường Các dị thường vật thể có kích thước khác nằm độ sâu khác bị biến đổi khác q trình tiếp tục giải tích Để thấy rõ ý nghĩa việc tiếp tục giải tích dị thường trọng lực ta xét ví dụ sau: + Hai cầu, nằm độ sâu h có khối lượng M, nằm độ sâu nh có khối lượng n M Các dị thường trọng lực cầu gây mặt đất điểm tâm cầu tương ứng là: Vz1(0,0,0) = kM h2 (1.11) Vz2(0,0,0) = kMn kn M = 2 h n h (1.12) Tức là: Vz =n V z1 (1.13) tiếp tục giải tích dị thường lên độ cao H = h, thì: Vz1(0,0,-H) = kn M kM Vz2(0,0,-H) = h (n + 1) 4h Nên Vz (0,0, − H ) 4n3 = f n Vz1 (0,0, − H ) (n + 1) (1.14) Khi n  Như vậy, tiếp tục giải tích dị thường trọng lực lên nửa khơng gian dị thường khối vật chất nằm nông giảm nhiều so với dị thường có nguồn gốc sâu Bây ta lại tiếp tục giải tích dị thường xuống nửa khơng gian bên đến độ sâu H=0.5h Tương ứng ta có: kMh 4kM Vz1(0,0,H) = Vz2(0,0,H) = h (2n − 1) h V z (0,0, H ) n3 =  n V z1 (0,0, H ) (2n − 1) (1.15) n>1 Điều chứng tỏ rằng, tiếp tục giải tích dị thường trọng lực xuống nửa không gian bên dị thường khu vực V z2 tăng lên chậm so với dị thường địa phương Vz1 Dị thường địa phương làm rõ qua phép biến đổi Sau số toán tiếp tục giải tích cụ thể 1.2.1 Bài tốn tiếp tục giải tích trường lên nửa khơng gian Nếu hàm điều hồ cho trước hình cầu hay mặt phẳng để xác định hàm khơng gian ngồi người ta sử dụng cơng thức Poisson Trong hệ toạ độ vng góc, trục Z hướng xuống dưới, tích phân Poisson có dạng: 2π ∞ z V(x,y,-z) = 2π Trong đó: V (ξ ,η ,0) dηdξ ∫ ∫ [(ξ - x) − ∞− ∞ + (η − y ) + z ] 32 (1.16) V(x,y,-z) giá trị hàm điều hoà điểm (x,y,-z) V( ξ ,η ,0 ) giá trị hàm điều hoà điểm mặt phẳng x0y Trong hệ toạ độ trụ thẳng đứng (r, α , z) có gốc toạ độ nằm hình chiếu điểm cần tính hàm mặt phẳng x0y tích phân Poisson có dạng: V(x,y,-z) = z 2π 2π ∞ V (r , α ,0) rdrdα + z )3 ∫ ∫ (r 0 (1.17) Nếu biến đổi tích phân Poisson hệ toạ độ vng góc cách lấy tích phân theo biến η từ − ∞ → ∞ ta thu tích phân Poisson trường hợp hai chiều: z V(x,-z) = π ∞ V (ξ ,0) dξ + z2 ∫ (ξ − x) −∞ (1.18) Đối với điểm có toạ độ (0,-h) nằm độ cao h thì: ∞ h V (ξ ,0) V(0,-h) = ∫ 2 dξ π − ∞ξ + h (1.19) 1.2.2 Bài tốn tiếp tục giải tích trường xuống nửa khơng gian Việc tiếp tục giải tích hàm điều hồ xuống nửa khơng gian bên phức tạp nhiều so với việc tiếp tục giải tích xuống nửa khơng gian bên Bài tốn tiếp tục giải tích xuống nửa khơng gian bên tốn khơng ổn định với biến đổi nhỏ hàm xuất phát cho giá trị hàm tiếp tục giải tích xuống bị sai lệch nhiều Hiện có nhiều phương pháp để tính chuyển trường xuống nửa khơng gian bên ta xét đến vài phương pháp thường sử dụng * Phương pháp thứ : Để tiếp tục giải tích trường ta sử dụng trực tiếp công thức Poisson tương tự (1.18) [2]: ∞ h V (ξ , h) V(0,0) = ∫ 2 dξ π − ∞ξ + h Do hàm số chưa biết V( ξ , h ) nằm dấu tích phân nên để giải tốn ta phải giải phương trình tích phân Có thể giải phương trình phương pháp gần liên tiếp Bản chất phương pháp sau: Trường độ cao h so với mặt quan sát có độ lớn nhỏ so với độ lớn trường mặt quan sát, ta có độ chênh lệch đó, tạm gọi ∆ 1V ( giả sử trục z hướng xuống dưới): ∆ 1V = V(0,0) – V(0,-h) Tương tự, trường mặt quan sát trường độ sâu h: ∆ 1V = V(0,h) – V(0,0) ta giả sử hiệu số khơng đổi h thay đổi [3,5], lúc với mức gần bậc thì: V(0,h) = V(0,0) + ∆ 1V (1.25) V(0,h) = 2V(0,0) – V(0,-h) (1.26) hay: Sau tiếp tục tính chuyển trường lên mức –2h ta có hiệu số giới nội khoảng (-h,-2h) ký hiệu hiệu số giới nội ∆ V, ta có cơng thức gần đúng: V(0,h) = V(0,0) + ∆ 1V + ∆ V (1.27) hay : V(0,h)=3V(0,0) – 3V(0,-h) + V(0,-2h) (1.28) Tiếp tục làm đến hiệu số thứ n thu công thức gần đúng: n V(0,h) = ∑ (−1) C k =0 k k +1 n +1 V (0,−kh) V(0,h) = (n+1) V(0,0) + (1.29) ∞ ∑V K i = −∞ i i Với n k k +1 K i = ∑ (−1) C n +1 k =1 n ξ ξ kh ξi +1 dξ k +1 = ∑ ( −1) k C n +1 (arctg i +1 − arctg i ) ∫ξi ξ + k h k =1 π π kh kh Dựa vào cơng thức thành lập Pa-let để tính chuyển trường xuống nửa khơng gian * Phương pháp thứ Phương pháp sử dụng định lý trung bình Gauss Trong trường hợp tốn hai chiều, ơng cho giá trị trung bình hàm vịng trịn giá trị hàm tâm vịng trịn đó: V(0,0) = V ( r , α )dα 2πr ∫ Thay phép lấy tích phân phép lấy tổng giới hạn số điểm lấy tổng 4, hình 1.2: V(0,0) = [V(-h,0) + V(h,0) + V(0,-h) + V(0,h)]/4 V(0,h) = 4V(0,0) – V(-h,0) - V(h,0) - V(0,-h) 1.3 Tính đạo hàm bậc cao trọng lực: Khi phân tích số liệu trọng lực, việc tính tốn đạo hàm bậc cao trọng lực đóng vai trị quan trọng Trong nhiều trường hợp đạo hàm bậc cao cho phép đơn giản nhiều thông số vật thể Việc tính đạo hàm bậc cao so với thành phần đo cho phép người ta phân chia trường thành thành phần khu vực địa phương riêng biệt Hàm trọng lực hàm thoả mãn điều kiện Trong hệ toạ độ Đề các, nghiệm toán Neuman ngồi xác định cơng thức: V(x,y,z) = 2π ∞ ∞ V z (ξ ,η ,0)dξdη ∫ ∫ [(ξ - x) − ∞− ∞ + (η − y ) + z ] Như xác định V z mặt phẳng z=0 (tức xác định ∆g mặt phẳng quan sát) người ta dựa vào cơng thức để tính hàm V Sau xác định hàm V, cách lấy đạo hàm theo toạ độ tương ứng người ta xác định đạo hàm với bậc khác trọng lực Trên sở công thức tổng quát để tính đạo hàm bậc cao trọng lực có dạng: ∂ m + n + pV ( x, y, z ) = m n p 2π ∂x ∂y ∂z ∞ ∞ ∫ ∫ Vz (ξ ,η ,0) − ∞− ∞  ∂ m+n+ p  dξdη  m n p  ∂x ∂y ∂z  [ (ξ − x) + (η − y ) + z ]  Để tính gần tích phân trường hợp chiều, người ta chia tồn diện tích lấy tích phân vòng tròn đồng tâm tia xuyên tâm Với hướng tính vậy, nhiều palet xây dựng giúp cho việc tính tốn đạo hàm tiện lợi palet Malovisco, palet Vexelop, palet Chepkin để tính đạo hàm thẳng đứng palet để tính đạo hàm ngang [5] Hình 1.3 Palet Malovisko (a), Vexelop (b), Chepkin (c) Nhìn chung, sách giáo khoa kinh điển thăm dò từ trọng lực [1,2,3,4,5,6] khơng đưa cơng thức giải tích để tính tốn đạo hàm bậc cao trọng lực Trong trường hợp toán hai chiều, số tác giả đề xuất cách tính đơn giản theo định nghĩa đạo hàm [6,8]: Cơng thức để tính đạo hàm theo phương nằm ngang thẳng đứng 1.4 Tính đạo hàm ngang cực đại Việc tính đạo hàm ngang khơng có lạ, tốn nhằm xác định điểm uốn đường cong quan sát Điểm uốn đường cong trường trọng lực quan sát thường nơi chuyển tiếp hai khối có mật độ khác Với số liệu quan sát diện, ta tính đào hàm theo chiều x y sau: Hướng véc tơ gradient tổng xác định theo qui tắc hình bình hành: Giá trị cực đại H[G(x,y,z)] điểm có cực đại X max xác định đa thức bậc dạng: a X2max + b Xmax + c Đây đường cong hồi qui qua điểm xem xét hai điểm lân cận CHƯƠNG CÁC PHÉP BIẾN ĐỔI TRƯỜNG TRONG MIỀN TẦN SỐ Ta chuyển sang việc nghiên cứu phép biến đổi trường phương pháp phổ Để thực việc phải sử dụng phép biến đổi Fourier ứng dụng miền tần số ta xem xét qua phép biến đổi Fourier định lý phổ [3,4] 2.1 Phép biến đổi Fourier Định nghĩa Một hàm F(x) không tuần hồn có giới hạn biểu diễn dạng tích phân Fourier: ∞ F(x)= S (ω )e jωx dω π −∫ ∞ (2.1) Trong đó: S( ω ) = ∞ ∫ F ( x )e − jωx dx (2.2) −∞ Theo (II-1) (II-2) với j= − , ω tham số đo Radian/m 2.2 Biến đổi trường miền tần số Lý thuyết chung biến đổi trường miền tần số Ta dùng tính chất ứng dụng phép biến đổi Fourier để chuyển sang xây dựng công thức biến đổi trường miền tần số Trong chương trước, xét phép biến đổi trường miền không gian ghi nhận: Dưới dạng toán học, tất phép biến đổi biểu diễn dạng: Vbd ( x0, y 0, z 0) = ∫∫ Vxp (ξ ,η ,0) K ( x0 − ξ , y − η , z 0)d ξ dη (2.43) trường hợp toán chiều, và: Vbd ( x0, z 0) = ∫ V xp (ξ ,0) K ( x0 − ξ , z 0)dξ (2.44) trường hợp toán chiều Chúng ta chứng tỏ công thức tổng quát phép biến đổi trường miền tần số có dạng trên, tức là: Vbd ( x0, z 0) = ∫ S xp (ω )F (ω )dω (2.45) đây, để đơn giản xét trường hợp chiều Theo lý thuyết, mối liên hệ tích chập đặc trưng hàm lối lối vào dạng tần số (2.33),(2.41) mô tả giống trình lọc tần số Và vậy, phép biến đổi khác sử dụng phân chia trường giống trình lọc tần số có đặc trưng khác Trong cơng thức tích chập (2.30) hàm F1 ( x − ξ ) coi hàm xuất phát (tín hiệu vào), cịn hàm F(x) hàm biến đổi (tín hiệu ra), theo (2.30), ta có Fbd ( x) = ∞ ∫ −∞ Fxp ( x − ξ ) P(ξ )dξ (2.46) Với toán biến đổi trường mà ta xét, hàm xuất phát hàm biến đổi trọng lực chúng phụ thuộc vào biến x z P( ξ ) đặc trưng phép biến đổi Giả sử hàm biến đổi theo x, ta có: ∞ Vbd ( x ) = ∫ V xf ( x − ξ ) P (ξ ) dξ (2.47) −∞ Giống (2.33) hay (2.41), miền tần số liên hệ với hệ thức: SVbd (ω ) = SVxf (ω ) ∗ F (ω ) (2.48) Với SVbd (ω ), SVxf (ω ), F (ω ) tương ứng đặc trưng phổ hàm biến đổi, hàm xuát phát nhân biến đổi Biểu thức (2.47) mô tả mối liên hệ Vbd V xp miền không gian, (2.48) mô tả quan hệ chúng miền tần số Để có cơng thức liên hệ miền không gian miền tần số ta biến đổi Fourier ngược hai vế (2.48) ta có: 2π ∞ jωx ∫ SVbd e dω = −∞ 2π ∞ ∫S Vbd (ω ) F (ω )e jωx dω (2.49) −∞ hay Vbd = 2π ∞ ∫S Vxf (ω ) F (ω )e jωx dω (2.50) −∞ Đây cơng thức mơ tả phép biến đổi trường miền tần số Như vậy, thực toán biến đổi trường miền tần số ta phải thực theo bước sau: Biến đổi Fourier thuận để tính phổ hàm xuất phát theo giá trị cho trước Tính đặc trưng tần số phép biến đổi (dạng biểu thức đặc trưng tần số phụ thuộc vào phép biến đổi cụ thể) Nhân đặc trưng phổ hàm xuất phát với đặc trưng tần số phép biến đổi Biến đổi Fourier ngược tích hàm xuất phát với đặc trưng tần số ta tìm giá trị tiếp tục giải tích mức CHƯƠNG MỘT SỐ KẾT QUẢ THỬ NGHIỆM TRÊN MƠ HÌNH VÀ KHU VỰC X THỀM LỤC ĐỊA VIỆT NAM 3.1 Mơ hình kết thử nghiệm Chúng tơi chọn mơ hình khu vực quan sát giả định có chiều dài theo x chiều rộng theo y, chiều 64 km Giữa khu vực, độ sâu 12 km có vật thể dạng cầu, có bán kính km, có mật độ dư 0.2 g/cm3 Nằm tuyến xuyên tâm theo trục X (cách tâm 15 km) phía, vị trí có cầu thể, có bán kính km, độ sâu km, mật độ dư hai cho 0.2 g/cm3 Mơ hình khu vực giả định nhằm mô khu vực quan sát có trường khu vực (do cầu thể lớn trung tâm) dị thường địa phương (do cầu thể nhỏ nông) tạo nên 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 Hình 3.1 Trường trọng lực cầu thể trung tâm 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 Hình 3.2 Trường trọng lực mơ hình có ba cầu thể 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 Hình 3.3 Tính trung bình trường (bán kính km) 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 Hình 3.4 Tính đạo hàm thẳng đứng (bậc 1) 0 0 0 0 0 0 0 0 0 0 0 0 0 0 Hình 3.5 Tính hạ trường xuống km 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 Hình 3.6 Tính nâng trường lên km Trên sở số liệu mơ hình khu vực giả định trình bày trên, chúng tơi thử nghiệm chương trình đánh giá khả tính tốn, hiệu số phép biến đổi trường trình bày chương 2., hiệu phép tính trung bình hố nâng trường Tất nhiên, phụ thuộc vào bán kính trung bình hố mức nâng mà độ trơn (mức loại nhiễu địa phương) khác (hình 3.4, 3.5, 3.8, 3.9, 3.10) Hiệu phép hạ trường tính đạo hàm theo phương thẳng đứng giống Thành phần địa phương (hai cầu thể nhỏ) rõ sau phép tính (hình 3.6, 3.7) Hiệu việc tính đạo hàm ngang cực đại lại khác, cho ta thấy rõ đường biên, phần tiếp xúc khối vật chất chênh lệch mật độ Với kết thử nghiệm mơ trên, n tâm sử dụng qui trình nói việc phân tích xử lý tài liệu thực tế 3.2 Kết thử nghiệm cho vùng X thuộc thềm lục địa Việt nam 9 0 8 0 7 0 6 0 0 0 0 Hình 3.7 Bản đồ trọng lực Bughe khu vực X Khu vực nghiên cứu nằm thềm lục địa phía Đơng nam Việt nam Nó giới hạn từ kinh độ 106 o 30’ đến kinh độ 1090 20’ từ vĩ độ 6o 0’ đến 9o 30’ Về mặt địa chất khoáng sản, khu vực có tiềm lớn lĩnh vực dầu khí [4] Phía Bắc Tây bắc khu vực tiếp giáp với khối nâng Côn sơn, phía Tây Nam giáp khối nâng Khorat-Natuna Phía Đơng Đơng bắc với bể Tư chính-Vũng mây, vùng sụt lún sâu [4] Bản đồ trọng lực Bughe khu vực trình bày hình 3.12 Mục đích tìm hiểu đặc điểm cấu trúc sâu khu vực nên phép nâng trường tiến hành 9 0 8 0 7 0 6 0 0 0 0 Hình 3.18 Bản đồ nâng trường lên 30 km Trên hình từ 3.14 đến 3.18 đồ trường nâng lên độ cao khác Ta thấy đến mức nâng 20 km 30 km, giá trị trường thay đổi, việc tính nâng cao khơng cần thiết Kết hợp nâng với việc tính đạo hàm ngang cực đại ta có tranh 0 8 0 7 0 0 0 0 Hình 3.21 Tính đạo hàm ngang cực đại mức nâng 20 30 km Qua hình từ 3.19 đến 3.21 ta thấy “băm nát” hình dần bị Như vậy, nâng trường lên độ cao lớn yếu tố địa phương dần Ở mức nâng 20 đến 30 km tồn ảnh hưởng cấu trúc sâu sâu, địa chất coi đứt gãy cấp Rõ ràng, với công cụ biến đổi trường (kết hợp nâng đạo hàm ngang cực đại) phát đặc điểm cấu trúc tầng nơng sâu khác Trên hình vẽ đạo hàm ngang cực đại hướng mũi tên hướng véc tơ đạo hàm ngang từ cao xuống thấp Có thể dựa vào để nói nâng lên hay hạ xuống đối tượng địa chất KẾT LUẬN Thơng qua việc tìm hiểu tốn biến đổi trường cụ thể phép tính giải tích trường đến nhận xét sau: - Có thể thực phép biến đổi trường miền tần số cách nhanh chóng nhờ máy tính điện tử - Kết biến đổi trường thử nghiệm số liệu mơ hình cho thấy độ tin cậy ổn định thuật toán thử nghiệm - Các phép biến đổi trường miền tần số yêu cầu số lượng điểm quan sát phải phù hợp sử dụng biến đổi Furier nhanh - Kết hợp phép nâng hạ trường với tính đạo hàm ngang cực đại cơng cụ hỗ trợ tốt cho chuyên gia phân tích tài liệu địa chất - Cần phải khảo sát nhiều kỹ lưỡng vấn đề kết hợp với nguồn thơng tin khác, ví dụ địa chấn thăm dò  ... thường trọng lực (xem hàm điều hoà) + Tính đạo hàm bậc cao trọng lực Chúng ta x? ?t đến nhóm phương pháp 1.1 Phương pháp trung bình hóa Việc phân chia dị thường trọng lực thành thành phần khu vực địa. .. sử dụng qui trình nói việc phân tích x? ?? lý tài liệu thực tế 3.2 Kết thử nghiệm cho vùng X thuộc thềm lục địa Việt nam 9 0 8 0 7 0 6 0 0 0 0 Hình 3.7 Bản đồ trọng lực Bughe khu vực X Khu vực. .. mà ta x? ?t, hàm xuất phát hàm biến đổi trọng lực chúng phụ thuộc vào biến x z P( ξ ) đặc trưng phép biến đổi Giả sử hàm biến đổi theo x, ta có: ∞ Vbd ( x ) = ∫ V xf ( x − ξ ) P (ξ ) dξ (2.47) −∞