1, Vẽ bản đồ dị thường Fai (Gravity) và địa hình đáy biển(Topo) trong KVNC từ KĐ 105111,VĐ 623 từ nguồn số liệu dowload trên trang topex.ucsd.edu.2, Bản đồ dị thường Fai xây dựng với phiên bản bao nhiêu? Độ chính xác? Mật độ điểm đo ntn?3, Bản đồ dị thường Fai xây dựng với phiên bản bao nhiêu? Độ chính xác? Mật độ điểm đo ntn?4, Tính dị thường trọng lực Bughe đơn giản từ nguồn số liệu trên?
Thuyết trình thăm dò trọng lực CHỦ ĐỀ BÀI TẬP Xây dựng đồ trường dị thường trọng lực bughe đơn giản khu vực bồn trũng Sông Hồng lân cận *Giáo viên hướng dẫn: Th.s Phan Thị Hồng Địa vật Lý K60 Đề 1, Vẽ đồ dị thường Fai (Gravity) địa hình đáy biển(Topo) KVNC từ KĐ 105-111,VĐ 6-23 từ nguồn số liệu dowload trang topex.ucsd.edu 2, Bản đồ dị thường Fai xây dựng với phiên bao nhiêu? Độ xác? Mật độ điểm đo ntn? 3, Bản đồ dị thường Fai xây dựng với phiên bao nhiêu? Độ xác? Mật độ điểm đo ntn? 4, Tính dị thường trọng lực Bughe đơn giản từ nguồn số liệu trên? I.Tổng quan khu vực nghiên cứu -Bồn trũng Sông Hồng nằm phía Bắc Việt Nam trải dài từ 105030 – 110030 kinh độ Đông, 16030 – 21000 vĩ độ Bắc Bồn trũng có phần nhỏ diện tích nằm đất liền thuộc đồng Sơng Hồng, lại phần lớn thuộc vùng biển vịnh Bắc Bộ biển miền Trung -Bồn trầm tích có dạng hình thoi kéo dài từ miền võng Hà Nội đến Vịnh Bắc Bộ phần ven biển miền Trung với lớp phủ trầm tích Đệ Tam dày 14 km I.Tổng quan khu vực nghiên cứu -Phía Tây bể trồi lộ đá móng Paleozoi-Mesozoi -Phía Đơng Bắc tiếp giáp bể Tây Lơi Châu -Phía Đơng Nam bể Đơng Nam Hải Nam bể Hồng Sa, -Phía Nam giáp bể trầm tích Phú Khánh -Bồn trũng Sơng Hồng có cấu trúc địa chất phức tạp thay đổi từ đất liền biển theo hướng Tây Bắc - Đông Nam Nam, bao gồm vùng địa chất khác nhau, đối tượng tìm kiếm thăm dò khác 1,Cách lấy số liệu:http://topex.ucsd.edu/cgi-bin/get_data.cgi II, Bản đồ dị thường Fai xây dựng với phiên bao nhiêu? Độ xác? Mật độ điểm đo nào? - Bản đồ dị thường Fai vẽ từ số liệu dị thường Fai với phiên 24.1 - Độ xác số liệu: sai số khoảng mGal Mật độ điểm đo: phút * phút cho ta điểm ( điểm cách hải lý ~1.852km) ( trích Wikipedia) III, Bản đồ độ sâu đáy biển xây dựng với phiên bao nhiêu? Độ xác? Mật độ điểm đo? +Bản đồ độ sâu đáy biển vẽ từ số liệu địa hình với phiên 18.1 +Độ xác số liệu: khoảng 0.03m +Mật độ điểm đo : tương tự với giá trị dị thường Fai phút * phút ghi điểm điểm cách hải lý ~ 1.852km IV, Vẽ đồ dị thường Fai (Gravity) địa hình đáy biển(Topo) KVNC từ KĐ 105-111,VĐ 16-23 từ nguồn số liệu dowload trang topex.ucsd.edu 1,Code GMT vẽ đồ dị thường Fai khu vực bồn trũng sông Hồng surface baitaptl1.xyz -R105/111/16/23 -I0.5m -C0.1 -Gbaitaptl1.grd grd2cpt baitaptl1.grd -Crainbow –Z > mycolor.cpt grdgradient baitaptl1.grd -Nt1 -A45 -Gbaitaptl1.grd grdimage baitaptl1.grd -Ibaitaptl1i.grd -R -JM6 -Cmycolor.cpt -X1i -Y3 -P -K > baitaptl1.ps psscale -D3/-0.5/4/0.2h -Cmycolor.cpt -Ba25f5:“BAN DO DI THUONG FAI (mGal)": -O -P –K >> baitapl1.ps grdcontour baitaptl1.grd -R -JM6 -A25f10 –C25 -S10 -W0.7p -Q200 -O -P -K >> baitaptl1.ps pscoast -R -JM6 -B3g3 -W0.5p/120/60/30 -G255 -O -P -K >>baitaptl1.ps pause Hình 2: Sơ đồ dị thường Fai(mGal)) khu vực bồn trũng Sông Hồng lân cận 2, Code GMT vẽ địa hình đáy biển(Topo) khu vực bồn trũng sông Hồng surface btttl2.xyz -R105/111/16/23 -I0.5m -C0.1 -Gbtttl2.grd grd2cpt btttl2.grd -Crainbow –Z > mycolor.cpt grdgradient btttl2.grd -Nt1 -A45 -Gbtttl2.grd grdimage btttl2.grd -Ibtttl2i.grd -R -JM6 -Cmycolor.cpt -X1i -Y3 -P -K > btttl2.ps psscale -D3/-0.5/4/0.2h -Cmycolor.cpt –Ba400f50:“BAN DO DO SAU DIA HINH DAY BIEN(m)": -O -P –K >> btttl2.ps grdcontour btttl2.grd -R -JM6 –A55f9 –C55 -S10 -W0.1p -Q200 -O -P -K >> btttl2.ps pscoast -R -JM6 -B3g3 -W0.5p/120/60/30 -G255 -O -P -K >>btttl2.ps pause Hình 3: Sơ đồ độ sâu địa hình đáy biển(m) khu vực bồn trũng Sơng Hồng lân cận V,Tính dị thường trọng lực Bughe đơn giản từ nguồn số liệu trên? A, Code mathlab để tính tốn dị thường Buge đơn giản xuất liệu để vẽ GMT x=load('topo.xyz'); y=load('baitaptl1.xyz'); dithuongfai=y(:,3); dithuongbugedg=dithuongfai+0.0419.*(2.67-1.03).*x; z=y(:,1); l=y(:,2); Muay = [z(:),l(:),dithuongbugedg(:)]; save dithuongtl.xyz Muay -ascii B,Giải thích code Mathlab lệnh nhập liệu x=load(‘topo.xyz’) baitaptl1 bttl2 y=load(‘baitaptl1.xyz’) y(:,3); lệnh giữ nguyên cột ma trận y dithuongbugedg=dithuongfai+(0.0419.*(2.67-1.03).*x): CT tính dị thường buge đơn giản z=y(:,1); l=y(:,2); Muay = [z(:),l(:),dithuongbugedg(:) ]; lệnh tạo ma trận cách ghép hàng cột Save dithuongtl.xyz Muay –ascii:lệnh xuất liệu C,Code vẽ dị thường buge đơn giản GMT surface dithuongbugedg.xyz -R105/111/16/23 -I0.5m -T0.2 -C0.1 -Gdithuongbugedg.grd grd2cpt dithuongbugedg.grd -Crainbow -Z > mycolor.cpt grdgradient dithuongbugedg.grd -Nt1 -A45 -Gdithuongbugedgi.grd grdimage dithuongbugedg.grd -Idithuongbugedgi.grd -R -JM6 -Cmycolor.cpt -X1i -Y3 -P -K > dithuongbugedg.ps psscale -D3/-0.5/5/0.2h -Cmycolor.cpt -Ba50f25:"DI THUONG BUGE DON GIAN(mGal)": -O -P -K >> dithuongbugedg.ps grdcontour dithuongbugedg.grd -R -JM6 -A10f9 -C10 -S10 -W0.5p -Q200 -O -P -K >> dithuongbugedg.ps pscoast -R -JM6 -B1g1 -W0.5p/120/60/30 -G255 -O -P -K >>dithuongbugedg.ps pause Hình 4: Sơ đồ dị thường Buge đơn giản(mGal) khu vực bồn trũng sông Hồng lân cận 3,Nâng trường surface dithuongbugedg.xyz -R105/111/16/23 -I0.5m -T0.2 -C0.1 -Gdithuongbugedg.grd grdfft dithuongbugedg.grd -C1000 -M -Gnangtrg.grd grd2cpt nangtrg.grd -Crainbow -Z > mycolor.cpt grdgradient nangtrg.grd -Nt1 -A45 -Gnangtrgi.grd grdimage nangtrg.grd -Inangtrgi.grd -R -JM6 -Cmycolor.cpt -X1i -Y3 -P -K > nangtrg.ps psscale -D3/-0.5/5/0.2h -Cmycolor.cpt -Ba25f5:"Nang truong": -O -P -K >> nangtrg.ps grdcontour nangtrg.grd -R -JM6 -A10f9 -C10 -S10 -W0.5p -Q200 -O -P -K >> nangtrg.ps pscoast -R -JM6 -B1g1 -W0.5p/120/60/30 -G255 -O -P -K >>nangtrg.ps pause Sơ đồ dị thường trọng lực buge mức nâng trường lên 1km Sơ đồ dị thường trọng lực buge mức nâng trường lên 10km 4,Hạ trường surface dithuongbugedg.xyz -R105/111/16/23 -I0.5m -T0.2 -C0.1 -Gdithuongbugedg.grd grdfft dithuongbugedg.grd -C-500 -M -Ghatrg.grd grd2cpt hatrg.grd -Crainbow -Z > mycolor.cpt grdgradient hatrg.grd -Nt1 -A45 -Ghatrgi.grd grdimage hatrg.grd -Ihatrgi.grd -R -JM6 -Cmycolor.cpt -X1i -Y3 -P -K > hatrg.ps psscale -D3/-0.5/5/0.2h -Cmycolor.cpt -Ba50f10:"Ha truong": -O -P -K >> hatrg.ps grdcontour hatrg.grd -R -JM6 -A25f9 -C25 -S10 -W0.2p -Q200 -O -P -K >> hatrg.ps pscoast -R -JM6 -B1g1 -W0.5p/120/60/30 -G255 -O -P -K >>hatrg.ps pause Sơ đồ dị thường trọng lực buge mức hạ trường 500m Sơ đồ dị thường trọng lực buge mức hạ trường 1000m 5,Gradient ngang cực đại surface dithuongbugedg.xyz -R105/111/16/23 -I0.5m -T0.2 -C0.1 -Gbaitaptl.grd grdmath baitaptl.grd DDX = A.grd grdmath baitaptl.grd DDY = B.grd grdmath A.grd B.grd HYPOT 111 DIV = E.grd grd2cpt E.grd -Crainbow -Z > mycolor.cpt grdgradient E.grd -Nt1 -A45 -GEi.grd grdimage E.grd -IEi.grd -R -JM6 -Cmycolor.cpt -X1i -Y3 -P -K > E.ps psscale -D3/-0.5/5/0.2h -Cmycolor.cpt -Ba5f1:"Gradient ngang cuc dai(mGal/km)": -O -P -K >> E.ps grdcontour E.grd -R -JM6 -A25f7 -C25 -S10 -W0.5p -Q200 -O -P -K >> E.ps Sơ đồ dị thường trọng lực buge đơn giản thực ... 3: Sơ đồ độ sâu địa hình đáy biển(m) khu vực bồn trũng Sơng Hồng lân cận V,Tính dị thường trọng lực Bughe đơn giản từ nguồn số liệu trên? A, Code mathlab để tính tốn dị thường Buge đơn giản xuất... -K >>baitaptl1.ps pause Hình 2: Sơ đồ dị thường Fai(mGal)) khu vực bồn trũng Sông Hồng lân cận 2, Code GMT vẽ địa hình đáy biển(Topo) khu vực bồn trũng sông Hồng surface btttl2.xyz -R105/111/16/23... ntn? 3, Bản đồ dị thường Fai xây dựng với phiên bao nhiêu? Độ xác? Mật độ điểm đo ntn? 4, Tính dị thường trọng lực Bughe đơn giản từ nguồn số liệu trên? I.Tổng quan khu vực nghiên cứu -Bồn trũng