Xây dựng mô hình lịch sử dòng nhiệt để xác định sự phân bố dầu khí ở vùng nước sâu bể phú khánh

ĐẠI HỌC QUỐC GIA TP HCM TRƯỜNG ĐẠI HỌC BÁCH KHOA NGUYỄN THỊ THU TRANG XÂY DỰNG MƠ HÌNH LỊCH SỬ DÒNG NHIỆT ĐỂ XÁC ĐỊNH SỰ PHÂN BỐ DẦU KHÍ Ở VÙNG NƯỚC SÂU BỂ PHÚ KHÁNH Chuyên ngành: Kỹ Thuật Dầu Khí Mã số: 8520604 LUẬN VĂN THẠC SĨ TP HỒ CHÍ MINH, tháng 10 năm 2020 Cơng trình hồn thành tại: Trường Đại học Bách Khoa – ĐHQG-HCM Cán hướng dẫn khoa học: TS NGUYỄN XUÂN HUY (Ghi rõ họ, tên, học hàm, học vị chữ ký) Cán chấm nhận xét 1: TS TRẦN NHƯ HUY (Ghi rõ họ, tên, học hàm, học vị chữ ký) Cán chấm nhận xét 2: TS BÙI THANH VÂN (Ghi rõ họ, tên, học hàm, học vị chữ ký) Luận văn thạc sĩ bảo vệ Trường Đại học Bách Khoa, ĐHQG Tp HCM ngày 15 tháng 10 năm 2020 Thành phần Hội đồng đánh giá luận văn thạc sĩ gồm: (Ghi rõ họ, tên, học hàm, học vị Hội đồng chấm bảo vệ luận văn thạc sĩ) PGS.TS TRẦN VĂN XUÂN – Chủ tịch Hội đồng TS TẠ QUỐC DŨNG – Thư ký Hội đồng TS TRẦN NHƯ HUY – Cán Phản biện TS BÙI THANH VÂN – Cán Phản biện TS TRẦN LÊ PHƯƠNG – Ủy viên Hội đồng Xác nhận Chủ tịch Hội đồng đánh giá LV Trưởng Khoa quản lý chuyên ngành sau luận văn sửa chữa (nếu có) CHỦ TỊCH HỘI ĐỒNG TRƯỞNG KHOA CỘNG HÒA XÃ HỘI CHỦ NGHĨA VIỆT NAM ĐẠI HỌC QUỐC GIA TP.HCM TRƯỜNG ĐẠI HỌC BÁCH KHOA Độc lập - Tự - Hạnh phúc NHIỆM VỤ LUẬN VĂN THẠC SĨ Họ tên học viên: NGUYỄN THỊ THU TRANG MSHV: 1970038 Ngày, tháng, năm sinh: 22/06/1996 Nơi sinh: Tiền Giang Chuyên ngành: KỸ THUẬT DẦU KHÍ Mã số: 8520604 I TÊN ĐỀ TÀI: XÂY DỰNG MƠ HÌNH LỊCH SỬ DÒNG NHIỆT ĐỂ XÁC ĐỊNH SỰ PHÂN BỐ DẦU KHÍ Ở VÙNG NƯỚC SÂU BỂ PHÚ KHÁNH NHIỆM VỤ VÀ NỘI DUNG: ▪ Nghiên cứu sở lý thuyết dòng nhiệt dòng nhiệt cổ ▪ Xây dựng sở quy trình tính tốn dịng nhiệt dịng nhiệt cổ từ liệu giếng, từ đó, biện luận dòng nhiệt khu vực nước sâu bể Phú Khánh ▪ Xây dựng mơ hình bồn trầm tích để đánh giá độ trưởng thành đá mẹ lịch sử hình thành hydrocarbon theo mặt cắt Tây Bắc-Đơng Nam, từ xác định phân bố dầu khí vùng nước sâu bể Phú Khánh II NGÀY GIAO NHIỆM VỤ: 24/02/2020 III NGÀY HOÀN THÀNH NHIỆM VỤ: 30/09/2020 IV CÁN BỘ HƯỚNG DẪN: TS NGUYỄN XUÂN HUY Tp HCM, ngày… tháng năm 2020 CÁN BỘ HƯỚNG DẪN CHỦ NHIỆM BỘ MÔN ĐÀO TẠO (Họ tên chữ ký) (Họ tên chữ ký) TRƯỞNG KHOA (Họ tên chữ ký) LỜI C M ƠN Luận văn hoàn thành sau thời gian học tập, nghiên cứu làm việc lâu dài, nghiêm túc Để có thành ngày hôm nay, tác giả xin gửi lời cảm ơn chân thành đến trường Đại học Bách Khoa, ĐHQG-HCM hỗ trợ thời gian phương tiện vật chất cho nghiên cứu Cám ơn giúp đỡ tận tình, kiến thức, kinh nghiệm vơ quý báu mà Quý Thầy Cô Khoa Kỹ thuật Địa chất & Dầu khí truyền đạt cho tác giả suốt khoảng thời gian học đại học cao học Với lòng biết ơn sâu sắc nhất, tác giả xin gửi lời cám ơn chân thành đến Thầy PGS.TS Trần Văn Xuân, Thầy TS Đỗ Quang Khánh dành thời gian có góp ý vơ q báu để tác giả chỉnh sửa hoàn thiện luận văn cách tốt Đặc biệt tác giả xin gửi đến Thầy TS Nguyễn Xuân Huy lòng biết ơn chân thành nhất, Thầy tận tình bảo, giảng dạy, chia sẻ kinh nghiệm cho tác giả bước trình học tập từ việc định hướng đề tài luận văn đến việc hoàn thành luận văn cách tốt Tác giả xin cảm ơn tài trợ Đại học Quốc gia Thành phố Hồ Chí Minh (ĐHQG-HCM) khuôn khổ đề tài mã số C2018-20-24 Tác giả xin gửi lời cám ơn đến Tổng công ty thăm dị Khai thác dầu khí (PVEP) cho phép tác giả sử dụng thông tin tài liệu quý báu Luận văn Mặc dù có nhiều nỗ lực kinh nghiệm, kiến thức thời gian nhiều hạn chế nên luận văn chắn khơng thể tránh khỏi sai sót Tác giả mong nhận ý kiến nhận xét quý báu Quý Thầy Cô bạn để luận văn hoàn thiện Thành phố Hồ Chí Minh, ngày 30 tháng 09 năm 2020 Học viên cao học NGUYỄN THỊ THU TRANG i TÓM TẮT LUẬN VĂN THẠC SĨ Bể Phú Khánh coi bồn trũng nước sâu tiềm dầu khí cịn lại thềm lục địa Việt Nam đến chưa nghiên cứu đầy đủ Lịch sử thăm dò bể Phú Khánh thưa thớt chưa đồng Với giếng khoan thăm dị vùng nơng từ 2011, đến chưa có giếng khoan thực sâu 1,000 m mực nước biển Bằng chứng lớp sapropel giàu rong tảo vết thấm dầu Đầm Thị Nại nghiên cứu Saurin (1944-1964) tạo quan tâm nguồn gốc dầu khu vực Vì vậy, dù sở liệu thăm dị khiêm tốn, bể Phú Khánh thu hút nhiều đầu tư nghiên cứu nhà khoa học nước lịch sử kiến tạo, chế hình thành bồn, hay minh giải địa chấn địa tầng mơ hình hệ thống dầu khí Luận văn hịa vào khơng khí chung thơng qua đóng góp nghiên cứu tính tốn dịng nhiệt vùng nước sâu bể Phú Khánh Dịng nhiệt yếu tố đóng vai trị quan trọng q trình sinh dầu khí Dịng nhiệt ảnh hưởng trực tiếp đến mức độ trưởng thành thể địa chất xem tầng đá mẹ sinh dầu tiềm Tuy nhiên, nghiên cứu dòng nhiệt ảnh hưởng chúng đến q trình sinh dầu khí bể Phú Khánh chưa quan tâm nghiên cứu nghiêm túc, thường quan tâm đến thông số nhiệt độ tập trầm tích nên chưa đủ sở thuyết phục Nghiên cứu làm sáng tỏ lịch sử dịng nhiệt theo thời gian địa chất, thơng qua: (1) Tổng hợp lại lịch sử kiến tạo với giai đoạn (Eocene-Miocene Sớm, Miocene Giữa-Muộn, Pliocene-nay), làm sở để ngoại suy giá trị dòng nhiệt vùng nước sâu (2) Xác định phân bố dòng nhiệt khu vực bể Phú khánh từ 42 điểm liệu liệu dịng nhiệt tồn cầu Goutorbe et al (2011) (3) Tính tốn dịng nhiệt khoảng độ sâu theo thành phần thạch học giá trị dịng nhiệt trung bình cho giếng TH-1X CMT-1X vùng nước nông bể ii (4) Giá trị dịng nhiệt cổ tính cho giếng CMT-1X từ quy trình bắt đầu việc giải nén, backstripping, tính sụt lún kiến tạo, hệ số tách giãn kết thúc với giá trị dòng nhiệt cổ theo mơ hình McKenzie (5) Từ đó, làm sở ngoại suy giá trị min-max điều kiện biên theo thời gian địa chất vùng nước sâu bể, gồm: dòng nhiệt, độ sâu mực nước biển cổ, nhiệt độ ranh giới bề mặt trầm tích (6) Xây dựng mơ hình bồn trầm tích 1D 2D để xác định thời điểm hình thành hydrocarbon, độ trưởng thành đá mẹ dự báo phân bố dầu khí vùng nước sâu bể Phú Khánh iii ABSTRACT Phu Khanh basin is considered as a potential deepwater basin for oil and gas remaining in Vietnam's continental shelf that has not been fully studied Exploration history in the Phu Khanh basin is still sparse and asynchronous With three exploration wells in the shallow areas since 2011, so far, no wells have been made more than 1,000 m below sea level The evidence of the algae-rich sapropel layers in oil seeps in Dam Thi Nai Lagoon in the Saurin study (1944-1964) has brought interests about the oil origin here Although the exploration database is still quite modest, the Phu Khanh basin still attracts a lot of research investment from domestic and foreign researchers on tectonic history, basin formation mechanism, or explanations of sequence stratigraphy and the model of the petroleum system This thesis also blends into that general atmosphere through research contributions in calculating heat flow in the Phu Khanh Basin’s deepwater Heat flow is one of the factors playing a very important role in oil and gas generation Heat flow directly affects a geological body's maturity that is considered a potential oil-generating source rock layer However, the study of heat and their effects on the oil and gas generation process in the Phu Khanh basin has not been seriously studied and usually only interested in the the sediments' current temperature, so there is not have enough convincing evidence This study will shed light on the history of heat flow over geological time, through: (1) Synthesizing the tectonic history with three main stages (Eocene – Early Miocene, Middle – Late Miocene, Pliocene – present), as a basis for extrapolating heat flow values in deepwater (2) Determining the distribution of heat flow in the Phu Khanh basin area from 42 data points in the Goutorbe global heat flow dataset (2011) (3) Calculating the heat flow in each depth range according to the lithological composition and average present heat flow for the wells TH-1X and CMT-1X in the shallow water of the Phu Khanh Basin iv (4) Paleo heat flow value is calculated for CMT-1X well following the process beginning with calculations of decompaction, backstripping, tectonic subsidence, stretching factor, and ending with heat flow value according to the McKenzie model (5) From there, as the basis for extrapolating the min-max value of boundary conditions according to geological time in deepwater Phu Khanh Basin, including heat flow, paleo water depth, sediment-water interface temperature (6) Building basin models of 1D and 2D to determine the hydrocarbon generation's time, the maturity of source rock, and predict oil and gas distribution in deepwater of the Phu Khanh Basin v LỜI CAM ĐOAN Tác giả xin cam đoan cơng trình nghiên cứu thân tác giả Các kết nghiên cứu kết luận luận văn trung thực không chép từ nguồn hình thức Việc tham khảo nguồn tài liệu (nếu có) thực trích dẫn ghi nguồn tài liệu tham khảo quy định TPHCM, ngày 30 tháng 09 năm 2020 Học viên cao học NGUYỄN THỊ THU TRANG vi MỤC LỤC LỜI C M ƠN i TÓM TẮT LUẬN VĂN THẠC SĨ ii ABSTRACT iv LỜI CAM ĐOAN vi MỤC LỤC vii DANH MỤC HÌNH NH xi DANH MỤC B NG BIỂU xiv DANH MỤC THUẬT NGỮ VÀ CHỮ VIẾT TẮT xvi MỞ ĐẦU xvii CHƯƠNG I TỔNG QUAN VỀ KHU VỰC NGHIÊN CỨU 1.1 Vị trí địa lí bể Phú Khánh 1.2 Lịch sử tìm kiếm thăm dị dầu khí bể Phú Khánh 1.3 Lịch sử kiến tạo 1.4 Hệ thống dầu khí bể Phú Khánh 1.4.1 Đá mẹ 1.4.2 Đá chứa 1.4.3 Đá chắn 1.4.4 Di cư, nạp bẫy 1.5 Hệ thống sở liệu sử dụng nghiên cứu CHƯƠNG CƠ SỞ LÝ THUYẾT VỀ DÒNG NHIỆT VÀ QUY TRÌNH XÂY DỰNG MƠ HÌNH DỊNG NHIỆT 11 2.1 Tổng quan dòng nhiệt 11 2.2 Tính tốn dịng nhiệt từ liệu giếng 12 vii − cy1 ' − cy2 ' 0 ( e − e ) = c y2' − y1' (A.2.8) A.3 Backstripping Độ sâu lớp địa tầng theo thời gian phần kết q trình tải trọng trầm tích bên Backstripping việc loại bỏ ảnh hưởng tải trọng trầm tích để tiết lộ động lực gây sụt lún kiến tạo, với giả định diễn bồn trầm tích chứa đầy nước Backstripping cho phép giếng khoan mặt cắt địa tầng khác so sánh độ lún kiến tạo tải trọng nước chúng (Hình A.3) Hình A Cân đẳng tĩnh minh họa kỹ thuật backstripping tải trọng trầm tích (được chỉnh sửa từ Allen Allen, 2013) Nói cách đơn giản nhất, theo Sclater Christie (1980), ảnh hưởng tải trọng trầm tích đánh sau Độ rỗng lớp trầm tích độ sâu là: 78 − cy1 ' − cy2 ' 0 ( e − e ) = c y2' − y1' (A.3.1) Vì mật độ khối lớp trầm tích (ρ b) phụ thuộc vào độ rỗng mật độ hạt trầm tích (ρg): b =  w + (1 −  )  g (A.3.2) Mật độ khối tồn cột trầm tích ( b ) tính đến lớp thứ i là: ( )  i  w + − i  gi   yi' b =   S  i    (A.2.3) ϕi độ rỗng trung bình lớp thứ i, ρgi mật độ hạt trầm tích lớp, yi bề dày lớp trầm tích thứ i, S tổng bề dày cột hiệu chỉnh cho nén ép nh hưởng tải trọng trầm tích sau coi vấn đề cân đẳng tĩnh cục (Airy) Nơi trầm tích thay cột nước:   −  Y =S m b   m −  w  (A.3.4) Y độ sâu tầng móng hiệu chỉnh theo tải trọng trầm tích, ρm, b ρw mật độ lớp manti, mật độ cột trầm tích trung bình mật độ nước Việc kết hợp tác động khác phép đo độ sâu theo chiều biển, thay đổi mực nước biển tuyệt đối tải trọng trầm tích tạo sụt lún kiến tạo cân đẳng tĩnh Airy:   −   w  Y = S  m b  −  SL   + (Wd −  SL )  −   −   m w  m w 79 (A.3.5) Với ΔSL mực nước biển cổ so với nay, Wd độ sâu mực nước biển cổ Đây độ sụt lún so với mặt cố định (mực nước biển ngày nay) xảy bồn trầm tích chứa đầy nước 80 PHỤ LỤC B: DỮ LIỆU DÒNG NHIỆT TỪ 108o-115o ĐÔNG ĐẾN 10 o-15o BẮC TỔNG HỢP TỪ CƠ SỞ DỮ LIỆU GOUTORBE ET AL (2011) STT 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 Kinh độ Vĩ độ Dòng nhiệt từ Dòng phương pháp nhiệt đo tương tự đạc (mW/m ) (mW/m2) 108.5 15.5 75 109.5 15.5 74 110.5 15.5 73.4 111.5 15.5 68.3 112.5 15.5 76.5 113.5 15.5 61.3 60.0901 114.5 15.5 67.8 108.5 14.5 77.1 109.5 14.5 77.9 110.5 14.5 78.9 111.5 14.5 64.7 64.2656 112.5 14.5 74.5 113.5 14.5 56.2 56.1389 114.5 14.5 60.7 108.5 13.5 74.4 109.5 13.5 85 110.5 13.5 85.7 111.5 13.5 71.7 112.5 13.5 84.9 113.5 13.5 53.2 54.284 114.5 13.5 56.2 50.7556 108.5 12.5 78.6 109.5 12.5 84.2 110.5 12.5 85.8 111.5 12.5 76.7 112.5 12.5 64.4 113.5 12.5 60.9 114.5 12.5 64.6 108.5 11.5 77.6 109.5 11.5 66.8 110.5 11.5 78.1 111.5 11.5 77.1 112.5 11.5 54.6 113.5 11.5 67.8 114.5 11.5 86.1 108.5 10.5 76.9 109.5 10.5 70.7 81 38 110.5 10.5 39 111.5 10.5 40 112.5 10.5 41 113.5 10.5 42 114.5 10.5 Giá trị trung bình (mW/m2) 82 74.5 85.9 89.1 82.5 77.6 73.2 57.1 PHỤ LỤC C: CODE VBA TRONG EXCEL MACRO PHỤC VỤ CHO TÍNH TỐN SỤT LÚN KIẾN TẠO Phần mềm có sẵn để tạo đường cong sụt lún bị hạn chế khả nhập liệu gốc từ bảng tính tệp văn thiếu kiểm sốt người dùng tham số đầu vào, truy cập vào giá trị tính tốn hiển thị đường cong sụt lún cuối Việc tính tốn sụt lún thực thơng qua Excel để giúp kiểm sốt giai đoạn tính tốn Tuy nhiên, có số lượng lớn giếng khoan với nhiều lớp địa tầng khác nhau, dẫn đến số lượng tính tốn lặp lặp lại nhiều lần việc tính tốn sụt lún Excel bảng tính tốn khơng đơn giản Việc tự động hóa đơn giản hóa q trình nhiều tốt làm giảm khả xảy sai số q trình tính tốn Và VBA (Visual Basic for Applications) giải pháp cho vấn đề VBA ngơn ngữ lập trình Excel nói riêng ứng dụng khác Microsoft Office nói chung Với VBA, tất thao tác cố định lặp lặp lại Excel tự động thực vòng vài giây Điều giúp tiết kiệm thời gian tăng hiệu suất tính tốn, đồng thời dễ dàng cập nhật lại liệu đầu vào tạo kết cách nhanh chóng Đây động lực ban đầu để tác giả thực ý tưởng Tuy nhiên, hạn chế mặt thời gian, kiến thức kinh nghiệm nên đoạn code áp dụng cho giếng để tự động lặp lại thao tác giải nén, backstripping tính tốn sụt lún kiến tạo Về phần dự tính ban đầu, tác giả dự định cần nhập liệu tên giếng, tên horizon độ sâu nóc, phần liệu khác độ rỗng bề mặt, hệ số nén ép, mật độ hạt, tuổi,… tự động điền vào để thực q trình tính tốn tốn lặp chưa thực Và đoạn code bên chưa xét đến kiện địa chất bào mịn, ngừng trầm tích Những hạn chế cần khắc phục tương lai để áp dụng cho việc tính tốn nhanh sụt lún kiến bồn trầm tích Tham chiếu R1C1 (Hình C.1) chọn sử dụng đoạn code tham chiếu hữu ích cho việc tính tốn vị trí dịng cột VBA Đối với kiểu R1C1, Excel ký hiệu “R” để dòng theo sau số thứ tự dòng (Row) ký hiệu "C" để cột (Column) theo sau số thứ tự cột Các số đại diện cho vị trí, khoảng cách tính 83 từ chọn tới tham chiếu khác dòng, cột Chỉ số sau R hay C để móc vng [ ] tham chiếu tương đối Chỉ số sau R hay C không để móc vng [ ] tham chiếu tuyệt đối R (chính R[0]) hay C (chính C[0]) khơng có số theo sau biểu thị dịng hay cột với hành Hình C Minh họa tham chiếu R1C1 VBA Các macro sử dụng để tính tốn độ sụt lún liệt kê bên dưới, với thích macro liệt kê bên phải Chú thích macro hiển thị màu xanh Các lệnh Visual Basic hiển thị màu xanh dương C.1 Định dạng biến Option Explicit Dim myCell Dim myNum Dim i As Integer Dim j As Integer Dim k As Integer Dim Rng As Range Lệnh Option Explicit có nghĩa tất biến sử dụng macro phải xác định trước chúng sử dụng Lệnh Dim đặt tên cho biến C.2 Start Sub Start() ‘Xac dinh so lop cua gieng, dien vao cot No., chep tieu de cot vao gieng Worksheets("Data").Activate Set myCell = Application.InputBox(prompt:="Chon o dau tien cho cot dau tien (ngoai cung ben trai) cho gieng", Type:=8) myNum = Application.InputBox("Gieng co bao nhieu horizons") myCell.Value = 84 Đặt ô tham chiếu (myCell) số lớp (myNum) cột trầm tích For i = To myNum - myCell.Offset(i, 0).Value = i + Next Cung cấp cho lớp số thứ tự cột bảng tính myCell.Offset(-1, 0).Value = "No." myCell.Offset(-1, 0).HorizontalAlignment = xlCenter myCell.Offset(-1, 0).VerticalAlignment = xlCenter Dán nhãn cột số lớp địa tầng myCell.Offset(-1, 1).Value = "Well ID" myCell.Offset(-1, 1).HorizontalAlignment = xlCenter myCell.Offset(-1, 1).VerticalAlignment = xlCenter Dán nhãn cột tên giếng myCell.Offset(-1, 2).Value = "Horizon" & Chr(10) & "Name" myCell.Offset(-1, 2).HorizontalAlignment = xlCenter myCell.Offset(-1, 2).VerticalAlignment = xlCenter Dán nhãn cột tên lớp myCell.Offset(-1, 3).Value = "Present Upper" & Chr(10) & "Burial Depth" & Chr(10) & "(m)" myCell.Offset(-1, 3).HorizontalAlignment = xlCenter myCell.Offset(-1, 3).VerticalAlignment = xlCenter Dán nhãn cột Độ sâu myCell.Offset(-1, 4).Value = "Present Lower" & Chr(10) & "Burial Depth" & Chr(10) & "(m)" myCell.Offset(-1, 4).HorizontalAlignment = xlCenter myCell.Offset(-1, 4).VerticalAlignment = xlCenter Dán nhãn cột Độ sâu đáy myCell.Offset(-1, 5).Value = "Lithology" & Chr(10) & "Type" myCell.Offset(-1, 5).HorizontalAlignment = xlCenter myCell.Offset(-1, 5).VerticalAlignment = xlCenter Dán nhãn cột thạch học myCell.Offset(-1, 6).Value = "Compaction" & Chr(10) & "Coefficient" myCell.Offset(-1, 6).HorizontalAlignment = xlCenter myCell.Offset(-1, 6).VerticalAlignment = xlCenter Dán nhãn cột hệ số nén ép myCell.Offset(-1, 7).Value = "Porosity" myCell.Offset(-1, 7).HorizontalAlignment = xlCenter myCell.Offset(-1, 7).VerticalAlignment = xlCenter Dán nhãn cột độ rỗng myCell.Offset(-1, 8).Value = "Density" & Chr(10) & "(kg/m^3)" myCell.Offset(-1, 8).HorizontalAlignment = xlCenter myCell.Offset(-1, 8).VerticalAlignment = xlCenter Dán nhãn cột mật độ myCell.Offset(-1, 9).Value = "Age" & Chr(10) & "(Ma)" myCell.Offset(-1, 9).HorizontalAlignment = xlCenter myCell.Offset(-1, 9).VerticalAlignment = xlCenter Dán nhãn cột tuổi myCell.Offset(-1, 10).Value = "Total" & Chr(10) & "Thickness" & Chr(10) & "(m)" myCell.Offset(-1, 10).HorizontalAlignment = xlCenter myCell.Offset(-1, 10).VerticalAlignment = xlCenter Dán nhãn cột bề d myCell.Offset(-1, 11).Value = "Total" & Chr(10) & "Subsidence" & Chr(10) & "(m)" myCell.Offset(-1, 11).HorizontalAlignment = xlCenter myCell.Offset(-1, 11).VerticalAlignment = xlCenter Dán nhãn cột Tổng sụt lún myCell.Offset(-1, 12).Value = "Tectonic" & Chr(10) & "Subsidence" & Chr(10) & "(m)" myCell.Offset(-1, 12).HorizontalAlignment = xlCenter myCell.Offset(-1, 12).VerticalAlignment = xlCenter Dán nhãn cột Sụt lún kiến 85 tạo myCell.Offset(-1, 13).Value = "Total" & Chr(10) & "Column" & Chr(10) & "Density" & Chr(10) & "(kg/m^3)" myCell.Offset(-1, 13).HorizontalAlignment = xlCenter myCell.Offset(-1, 13).VerticalAlignment = xlCenter Dán nhãn cột Tổng mật độ myCell.Offset(-1, 14).Value = "Present" & Chr(10) & "Thickness" & Chr(10) & "(m)" myCell.Offset(-1, 14).HorizontalAlignment = xlCenter myCell.Offset(-1, 14).VerticalAlignment = xlCenter Dán nhãn cột bề dày myCell.Offset(-1, 15).Value = "Porosity" myCell.Offset(-1, 15).HorizontalAlignment = xlCenter myCell.Offset(-1, 15).VerticalAlignment = xlCenter Dán nhãn cột độ rỗng tính tốn myCell.Offset(-1, 16).Value = "Density" & Chr(10) & "(kg/m^3)" myCell.Offset(-1, 16).HorizontalAlignment = xlCenter myCell.Offset(-1, 16).VerticalAlignment = xlCenter Dán nhãn cột mật độ tính tốn myCell.Offset(-1, 17).Value = "Proportional" & Chr(10) & "Density" & Chr(10) & "(kg/m^3)" myCell.Offset(-1, 17).HorizontalAlignment = xlCenter myCell.Offset(-1, 17).VerticalAlignment = xlCenter Dán nhãn cột mật độ tỉ lệ Range((myCell.Offset(-2, 15)), (myCell.Offset(-2, 17))).Select Selection.Merge myCell.Offset(-2, 15).Value = "T= 0Ma" myCell.Offset(-2, 15).HorizontalAlignment = xlCenter End Sub Nhóm độ rỗng, mật độ mật độ tỉ lệ T=0Ma lại với C.3 Độ sâu đ y Sub pLBD() 'Tinh toan sau day hien tai cho moi lop Điền độ sâu đáy cho lớp (nóc lớp bên dưới) For i = To myNum - myCell.Offset(i, 4).Value = "=R[1]C[-1]" Next End Sub C.4 Bề dày Sub pThick() ‘Tinh be day hien tai cua moi lop Tính tốn bề dày lớp: y2-y1 For i = To myNum - myCell.Offset(i, 14).Formula = "= y2_-y1_" Next End Sub C.5 Bề dày tổng Sub tThick() 'Hoan Cot be day tong Tính tốn bề dày cộng dồn cột For i = To myNum - Set Rng = Range(myCell.Offset(i, 14), myCell.Offset(myNum 2, 14)) 86 myCell.Offset(i, 10).Formula = "=SUM(" & Rng.Address() & ")" Next Set Rng = Range(myCell.Offset(0, 14), myCell.Offset(myNum - 2, 14)) myCell.Offset(0, 11).Formula = _ "=SUM(" & Rng.Address() & ")" End Sub Tổng sụt lún T=0Ma bề dày cộng dồn cột C.7 Độ rỗng mật độ Sub Time_0MA() 'Tinh Do rong tai T=0Ma Tính toán độ rỗng lớp: − cy − e − cy 0 ( e = c ) y − y1 For i = To myNum - myCell.Offset(i, 15).Formula = "=phi/cc)*((EXP(-cc*y1_)EXP(-cc*y2_))/(y2_-y1_)))" myCell.Offset(i, 15).NumberFormat = "0.000" Next 'Tinh mat tai T=0Ma For i = To myNum - myCell.Offset(i, 16).Formula = "=(RC[-1]*rho_h2o)+((1-RC[1])*rho)" myCell.Offset(i, 16).NumberFormat = "0.000" Tính tốn mật độ lớp T=0Ma sử dụng công thức: 'Tinh toan mat ti le tai T=0Ma For j = To myNum - myCell.Offset(i, 17).Formula = "=(RC[-1]*RC[-3])/" & myCell.Offset(0, 10) & "" myCell.Offset(i, 17).NumberFormat = "0.000" Next Tính toán mật độ tỉ lệ lớp T=0Ma:  S * =  S*y Set Rng = Range(myCell.Offset(0, 17), myCell.Offset(myNum 2, 17)) myCell.Offset(0, 13).Formula = "=SUM(" & Rng.Address() & ")" Next End Sub Tổng mật độ tỉ lệ điền vào  s =  w + (1 −  )  sg S cột Tổng mật độ T=0Ma C.8 Các tính tốn khoảng thời gian Vì số lượng khoảng thời gian phụ thuộc vào số lớp, số lượng tính tốn cần thiết khác cột trầm tích Số lượng khoảng thời gian xác định từ số lớp cột trầm tích Phần thiết lập cột trống để điền vào phép tính Sub D_TS_col() ' Thiet lap cac cot cho Cac tinh toan Giai nen va Sut lun kien tao tai moi buoc thoi gian For i = To myNum myCell.Offset(-1, 17 Chr(10) & "to Top" & myCell.Offset(-1, 17 xlCenter + ((6 * i) + 1)).Formula = "Depth" & Chr(10) & "(m)" + ((6 * i) + 1)).HorizontalAlignment = 87 Dán nhãn cột “Độ Sâu Tới Nóc” cho khoảng thời gian myCell.Offset(-1, 17 + ((6 * i) + 1)).VerticalAlignment = xlCenter myCell.Offset(-1, 17 + ((6 * i) Chr(10) & "to Bottom" & Chr(10) myCell.Offset(-1, 17 + ((6 * i) xlCenter myCell.Offset(-1, 17 + ((6 * i) xlCenter + 2)).Formula = "Depth" & & " (m)" + 2)).HorizontalAlignment = Dán nhãn cột “Độ Sâu Tới Đáy” cho khoảng thời gian + 2)).VerticalAlignment = myCell.Offset(-1, 17 + ((6 * i) + 3)).Formula = "Decompacted" & Chr(10) & "Thickness" & Chr(10) & "(m)" myCell.Offset(-1, 17 + ((6 * i) + 3)).HorizontalAlignment = xlCenter myCell.Offset(-1, 17 + ((6 * i) + 3)).VerticalAlignment = xlCenter myCell.Offset(-1, 17 + ((6 * i) + 3)).ColumnWidth = 13.29 Dán nhãn cột “Bề dày giải nén” khoảng thời gian myCell.Offset(-1, 17 + ((6 * i) + 4)).Formula = "Porosity" myCell.Offset(-1, 17 + ((6 * i) + 4)).HorizontalAlignment = xlCenter myCell.Offset(-1, 17 + ((6 * i) + 4)).VerticalAlignment = xlCenter Dán nhãn cột độ rỗng tính tốn khoảng thời gian myCell.Offset(-1, 17 + ((6 * i) + 5)).Formula = "Density" & Chr(10) & "(kg/m^3)" myCell.Offset(-1, 17 + ((6 * i) + 5)).HorizontalAlignment = xlCenter myCell.Offset(-1, 17 + ((6 * i) + 5)).VerticalAlignment = xlCenter Dán nhãn cột mật độ tính tốn khoảng thời gian myCell.Offset(-1, 17 + ((6 "Proportional" & Chr(10) & myCell.Offset(-1, 17 + ((6 xlCenter myCell.Offset(-1, 17 + ((6 xlCenter myCell.Offset(-1, 17 + ((6 Dán nhãn cột mật độ tỉ lệ bước thời gian * i) + 6)).Formula = "Density" & Chr(10) & "(kg/m^3)" * i) + 6)).HorizontalAlignment = * i) + 6)).VerticalAlignment = * i) + 6)).ColumnWidth = 12.29 Range((myCell.Offset(-2, 17 + ((6 * i) + 1))), (myCell.Offset(-2, 17 + ((6 * i) + 6)))).Select Selection.Merge myCell.Offset(-2, 17 + ((6 * i) + 1)).Value = "T=" & myCell.Offset(i + 1, 9) & "Ma" myCell.Offset(-2, 17 + ((6 * i) + 1)).HorizontalAlignment = xlCenter Nhóm Độ sâu tới Nóc, Độ sâu tới Đáy, bề dày giải nén, độ rỗng, mật độ mật độ tỉ lệ với khoảng thời gian Next End Sub C.8.1 C c độ sâu bề dày giải nén Sub NewDepths() ' Tinh toan cac sau va be day da giai nen Tính tốn độ sâu tới đáy lớp bề dày lớp khoảng thời gian For i = To myNum - For j = i + To myNum - myCell.Offset(j, 17 + ((6 * i) + 1)).Formula = "=R[-1]C+R[1]C[2]" myCell.Offset(j, 17 + ((6 * i) + 1)).NumberFormat = "0.00" Độ sâu tới cùa lớp bên độ sâu lớp bên + bề dày lớp 88 myCell.Offset(j, 17 + ((6 * i) + 2)).Formula = "=Tn+RC[1]+((phi/cc)*(EXP(-cc*RC[-1])-EXP(-cc*RC)-EXP(-cc*y1_)+EXP(cc*y2_))" myCell.Offset(j, 17 + ((6 * i) + 2)).NumberFormat = "0.00" Đối với tất lớp khác sử dụng phép tính lặp cơng thức: y2' − y1' = y2 − y1 0 − cy ( e − e −cy c  + e − cy − e −cy c − ( myCell.Offset(j, 17 + ((6 * i) + 3)).Formula = "=RC[-1]-RC[2]" myCell.Offset(j, 17 + ((6 * i) + 3)).NumberFormat = "0.00" Next j Next i End Sub ) ' ' ) Bề dày lớp khoảng thời gian tính sử dụng cơng thức: y ' = y2' − y1' C.8.2 Độ rỗng mật độ giải nén Sub NewDensity() 'Tinh toan rong va mat da giai nen For i = To myNum - For j = i + To myNum - myCell.Offset(j, 17 + ((6 * i) + 4)).Formula = _ "=phi/cc)*((EXP(-cc*RC[-3])-EXP(-cc*RC[-2]))/(RC[-2]-RC[3])))" myCell.Offset(j, 17 + ((6 * i) + 4)).NumberFormat = "0.000" myCell.Offset(j, 17 + ((6 * i) + 5)).Formula = "=(RC[1]*rho_h2o)+((1-RC[-1])*rho)" myCell.Offset(j, 17 + ((6 * i) + 5)).NumberFormat = "0.000" Next j Next i End Sub Tính tốn độ rỗng giải nén mật độ cho lớp khoảng thời gian định cách sử dụng phương trình: = − cy ' − e − cy 0 ( e c ' ) y −y ' ' Mật độ lớp khoảng thời gian cụ thể tính cơng thức:  s =  '  w + (1 −  ' ) sg C.8.3 Tổng sụt lún Sub ToSub() 'Hoan cot Tong sut lun For i = To myNum - Set Rng = Range(myCell.Offset(0, 11 + ((6 * i) + 3)), (myCell.Offset(myNum - 2, 11 + ((6 * i) + 3)))) myCell.Offset(i, 11).Formula = "=SUM(" & Rng.Address() & ")" myCell.Offset(i, 11).NumberFormat = "0.00" Next End Sub C.8.4 Mật độ tỉ lệ Sub PropDens() 'Tinh toan mat ti le va dien vao cot mat tong Dim total As Long 89 Tổng sụt lún khoảng thời gian định tổng bề dày giải nén lớp khoảng thời gian For i = To myNum - total = myCell.Offset(i + 1, 11) For j = i + To myNum - myCell.Offset(j, 17 + ((6 * i) + 3]/" & total & "" myCell.Offset(j, 17 + ((6 * i) + If IsError(myCell.Offset(j, 17 + myCell.Offset(j, 17 + ((6 * i) + Next j Next i 6)).Formula = "=RC[-1]*RC[- Tính tốn mật độ tỉ lệ lớp khoảng thời gian sử dụng công thức  y  S * = S* S 6)).NumberFormat = "0.000" ((6 * i) + 6)).Value) Then 6)).Value = Else For i = To myNum - Set Rng = Range(myCell.Offset(0, 12 + ((6 * (myCell.Offset(myNum - 2, 12 + ((6 * i) + 5)))) myCell.Offset(i, 13).Formula = _ "=SUM(" & Rng.Address() & ")" myCell.Offset(i, 13).NumberFormat = "0.000" Next End Sub i) + 5)), Tổng mật độ tỉ lệ trong khoảng thời gian C.8.5 Sụt lún kiến tạo Sub TeSub() 'Tinh sut lun kien tao For i = To myNum - myCell.Offset(i, 12).Formula = "=(ToS*(rho_mantleTCD))/(rho_mantle-rho_h2o)" myCell.Offset(i, 12).NumberFormat = "0.00" Next End Sub 90 Tính tốn sụt lún kiến tạo cột trầm tích khoảng thời gian cụ thể:   − s  Z = S  m   m −  w  LÝ LỊCH TRÍCH NGANG Họ tên: NGUYỄN THỊ THU TRANG Ngày, tháng, năm sinh: 22/06/1996 Nơi sinh: Tiền Giang Địa liên lạc: 3/51 Thành Thái, Phường 14, Quận 10, Tp Hồ Chí Minh Q TRÌNH ĐÀO TẠO ▪ 9/2014-11/2018: sinh viên Khoa Kỹ thuật Địa chất Dầu khí, Bộ mơn Địa chất Dầu khí, Đại học Bách Khoa TP.HCM ▪ 3/2019 – nay: học viên cao học chuyên ngành Kỹ thuật Dầu khí, Khoa Kỹ thuật Địa chất & Dầu khí, Đại học Bách Khoa TP.HCM Q TRÌNH CƠNG TÁC ▪ 11/2018 – nay: Nghiên Cứu viên Phịng Thí Nghiệm Mơ Phỏng, Khoa Kỹ thuật Địa chất & Dầu khí, Đại học Bách Khoa TP.HCM CÁC GI I THƯỞNG VÀ CHỨNG NHẬN 2019 Đạt Học bổng TOSHIBA dành cho Học viên cao học, Nghiên cứu sinh có thành tích học thuật nghiên cứu xuất sắc 2019 Giải Best Poster Award cho báo “Testing Neural Networks Assessment Based on Data-Driven Using Well Log Data: a Case Study in Cuu Long Basin, Offshore Vietnam" Hội Nghị Khoa Học lần thứ 16 trường ĐH Bách Khoa TpHCM 2019 Đạt Học bổng tham dự Trường hè Khoa Học Việt Nam (VSSS’07) Quy Nhơn, Bình Định 2019 Đội thi đạt Giải thi học thuật Science Hackathon Trường hè Khoa Học Việt Nam (VSSS’07) Quy Nhơn, Bình Định 2018 Tốt nghiệp Huy Chương Bạc, ngành Kỹ thuật Dầu khí 12/2017- Đội thi vào vịng B n kết thi Imperial Barrel Award (IBA) 2018 03/2018 (khu vực Châu Á-Thái Bình Dương) AAPG tổ chức 2017, 2018 Đạt Học bổng VALLET (2 lần) Tổ chức gặp gỡ Việt Nam (Rencontres du Vietnam) trao tặng 2015, 2016, 2017 Đạt Học bổng SPE (3 lần) cho năm học 2016-2017, 2015-2016, 20142015 2017 Giải Best Poster Award cho báo “Application of Seismic Attributes to Determine Fracture and Fault System in Fractured Basement Reservoir” Hội Nghị Khoa Học lần thứ 15 trường ĐH Bách Khoa TpHCM 91 2017 Đội thi đạt Giải nhì thi Học thuật ROCK THE SOIL 2017 Khoa Kỹ thuật Địa chất & Dầu khí, trường ĐH Bách Khoa Tp.HCM tổ chức 2016 Đội thi đạt Giải nhì thi Học thuật PETROSMART 2016 SPE ĐH Khoa Học Tự Nhiên Tp.HCM tổ chức 2016 Đội thi đạt Giải thi Học thuật CHẤN ĐỘNG PANGAEA 2016 ĐH Khoa Học Tự Nhiên Tp.HCM tổ chức 2016 Đạt Học bổng KrisEnergy 92 ... TÀI: XÂY DỰNG MƠ HÌNH LỊCH SỬ DÒNG NHIỆT ĐỂ XÁC ĐỊNH SỰ PHÂN BỐ DẦU KHÍ Ở VÙNG NƯỚC SÂU BỂ PHÚ KHÁNH NHIỆM VỤ VÀ NỘI DUNG: ▪ Nghiên cứu sở lý thuyết dòng nhiệt dòng nhiệt cổ ▪ Xây dựng sở quy... Cơ sở lý thuyết dòng nhiệt quy trình xây dựng mơ hình dịng nhiệt CHƯƠNG 3: Tính tốn dịng nhiệt vùng nước sâu bể Phú Khánh CHƯƠNG 4: Xác định phân bố dầu khí vùng nước sâu bể Phú Khánh KẾT LUẬN... “XÂY DỰNG MƠ HÌNH DỊNG NHIỆT ĐỂ XÁC ĐỊNH SỰ PHÂN BỐ DẦU KHÍ Ở VÙNG NƯỚC SÂU BỂ PHÚ KHÁNH” Nghiên cứu thay đổi dịng nhiệt theo thời gian địa chất cơng việc quan trọng, sở để xây dựng mơ hình bồn