Giới thiệu dự án

Bể trầm tích Sông Hồng (BBSH) là một trong những bể trầm tích chứa dầu khí quan trọng nhất tại thềm lục địa Việt Nam với diện tích trải dài trên 650 km theo phương Tây Bắc - Đông Nam (TB-ĐN). Trong đó, khu vực Đông Bắc đứt gãy Sông Lô (thuộc Lô 102 & 106, vịnh Bắc Bộ) sở hữu cấu trúc địa chất phức tạp, chịu ảnh hưởng của các pha tách giãn Paleogen và biến dạng nén ép nghịch đảo kiến tạo trong Miocen. Mặc dù các tầng chứa truyền thống như đá móng cacbonat nứt nẻ trước Kainozoi (lưu lượng thử vỉa DST đạt 6.000 - 7.000 thùng/ngày tại cấu tạo Hàm Rồng) và cát kết Oligocen đã được thăm dò, việc mở rộng trữ lượng đòi hỏi phải đánh giá chính xác tiềm năng của các vỉa cát kết Miocen giữa (hệ tầng Phù Cừ - $N_1^2 pc$).

Vấn đề cốt lõi đặt ra là sự phân dị tướng trầm tích mạnh mẽ và biến đổi thạch học phức tạp theo không gian khiến việc dự đoán các thân cát chứa dầu khí gặp nhiều rủi ro. Các bẫy cấu tạo khép kín truyền thống dần cạn kiệt, đòi hỏi công tác tìm kiếm thăm dò phải chuyển dịch sang bẫy địa tầng và bẫy kết hợp (stratigraphic and combination traps).

                        HỆ THỐNG DẦU KHÍ MIOCEN GIỮA
+-------------------------------------------------------------------------+
|  Tầng chắn khu vực: Sét kết Miocen giữa (Dày 80 - 150m, sét 71 - 86%)  |
+-------------------------------------------------------------------------+
|  Tầng chứa mục tiêu: Cát kết hệ tầng Phù Cừ (Độ rỗng Phi: 19.3 - 29.6%)|
+-------------------------------------------------------------------------+
|  Tầng sinh chính: Sét than Oligocen (Đình Cao) / Eocen (TOC: 0.5 - 2.0%)|
+-------------------------------------------------------------------------+

Mục tiêu nghiên cứu của đồ án bao gồm:

  1. Thu thập, chuẩn hóa và tích hợp cơ sở dữ liệu địa chấn 2D (9.200 km tuyến), 3D (1.050 $\text{km}^2$) cùng tài liệu địa vật lý giếng khoan (ĐVLGK) từ 11 giếng khoan thăm dò trọng điểm (102-CQ-1X, 102-TB-1X, 106-YT-1X, 106-HR-1X, 106-HRN-1X, 107-BAL-1X...).
  2. Thiết lập khung địa tầng phân tập (Sequence Stratigraphy) độ phân giải cao cho trầm tích Miocen giữa, nhận diện các mặt bất chỉnh hợp (BCH) chính: U260, U240, U220, U210 và U200.
  3. Phân tích tướng địa chấn (Seismic Facies) kết hợp phân loại đường cong Gamma Ray (GR Electrofacies) để khôi phục cổ môi trường lắng đọng trầm tích.
  4. Xây dựng bản đồ đẳng thời (TWT), bản đồ đẳng sâu (Depth), bản đồ bề dày (Isopach) và bản đồ tướng môi trường trầm tích cho 4 tập Miocen giữa (Tập I đến Tập IV).
  5. Dự đoán không gian phân bố các thân cát chứa triển vọng (độ rỗng đạt 19,3% - 29,6%), định hướng vị trí bẫy địa tầng phục vụ khoan thẩm lượng.

Phạm vi nghiên cứu giới hạn tại khu vực Đông Bắc đứt gãy Sông Lô, kéo dài từ đứt gãy Sông Lô (F3) đến đứt gãy Hải Dương (F6) và đới nâng Bạch Long Vĩ. Giới hạn địa tầng tập trung vào lát cắt Miocen giữa (độ sâu tuyệt đối từ 1.400 m đến 3.200 m).

Phân tích và thiết kế giải pháp

Phân tích hiện trạng

Trước đây, công tác minh giải cấu trúc chủ yếu dựa trên các lát cắt địa chấn 2D truyền thống và phương pháp thạch học giếng khoan đơn lẻ. Cách tiếp cận này bộc lộ nhiều điểm nghẽn khi đối mặt với các đới nghịch đảo kiến tạo và biến đổi tướng ven bờ.

Tiêu chí kỹ thuật Phương pháp cấu trúc 2D cổ điển Địa vật lý giếng khoan cục bộ Phương pháp Địa chấn Địa tầng tích hợp (Đề xuất)
Nhận diện bẫy phi cấu tạo Rất hạn chế (chỉ thấy nếp lồi) Không thể ngoại suy diện rộng Nhận diện chính xác nêm vát, onlap, kênh ngầm
Độ phân giải thời gian Thấp ($\Delta t > 25\text{ ms}$) Rất cao tại giếng ($< 0.5\text{ m}$) Cân bằng tối ưu giữa không gian và chiều sâu
Khả năng dự báo tướng đá Kém (chỉ suy đoán hình học) Tốt nhưng mang tính điểm Đồng bộ hóa đường cong GR với trường sóng 3D
Độ tin cậy liên kết vỉa Dưới 60% khi gặp đứt gãy lớn Thấp khi khoảng cách giếng $> 5\text{ km}$ Đạt trên 88% nhờ khung Sequence Stratigraphy

Bảng phân loại yêu cầu kỹ thuật theo mô hình MoSCoW cho hệ thống giải pháp:

  • Must have: Nhận diện chính xác 5 bề mặt phản xạ địa chấn ranh giới (U260, U240, U220, U210, U200); xây dựng hàm chuyển đổi Thời gian - Độ sâu ($T-Z$); phân chia 4 tập trầm tích Miocen giữa.
  • Should have: Tự động hóa phân loại hình thái đường cong GR (hình chuông, hình trụ, hình phễu, hình răng cưa); thành lập bản đồ đẳng dày phân tập Isopach.
  • Could have: Đánh giá định lượng thuộc tính biên độ phản xạ (Seismic Inversion / RMS Amplitude) để phát hiện dị thường khí (Bright spot).
  • Won't have: Mô hình hóa dòng chảy vỉa động (Dynamic Reservoir Simulation) trong giai đoạn này.

Thiết kế hệ thống

Quy trình kỹ thuật tích hợp minh giải địa chấn địa tầng và tướng điện địa vật lý giếng khoan được thiết kế theo chu trình khép kín:

[Dữ liệu Địa chấn 2D/3D (SEGY)]      [Dữ liệu Giếng khoan (LAS/DLIS)]
    [Khử nhiễu & Chuẩn hóa]               [Hiệu chỉnh GR, Sonic, Caliper]
                [Tạo Synthetic Seismogram]
             (Well-Tie: Đối sánh Ranh giới vỉa)
                [Minh giải Horizon & Faults]
              (U260, U240, U220, U210, U200)
        [Chuyển đổi Time-to-Depth (Hàm V(z))]
   [Phân tích Tướng Địa chấn]             [Phân tích Tướng Điện GR]
 (Phân lớp, Biên độ, Cấu trúc)         (Chuông, Trụ, Phễu, Răng cưa)
          [Bản đồ Tướng Môi trường Trầm tích 3D]
          [Đánh giá Tiềm năng Tầng chứa Cát kết]

Technology stack và công cụ triển khai chuyên dụng:

  • Nền tảng minh giải địa chấn chính: Schlumberger Petrel E&P Software Platform (v2018.1) phục vụ picking horizon, fault modeling và time-to-depth conversion.
  • Phân tích địa vật lý giếng khoan: Techlog (v2017.2) xử lý tổ hợp đường cong GR, Sonic (DT), Density (RHOB), Resistivity (LLD/LLS).
  • Phần mềm địa chấn phụ trợ: IHS Kingdom Suite (v8.8) và SeisWorks/OpenWorks.
  • Xử lý dữ liệu & Thuật toán phân tích: Python 3.8 với các thư viện khoa học: lasio (v0.29) đọc file LAS, obspy (v1.2.2) xử lý sóng địa chấn, scipy (v1.7.0) nội suy bề mặt và hồi quy phi tuyến.

Methodology

Phương pháp luận dựa trên lý thuyết Sequence Stratigraphy kinh điển của Exxon (Vail, Mitchum, 1977) kết hợp khung phân loại tướng điện sinh thái của Emery & Myers (1996):

  1. Phân tích tập địa chấn: Phân định các tập trầm tích giới hạn bởi ranh giới bất chỉnh hợp đáy (Baselap: Onlap, Downlap) và bất chỉnh hợp nóc (Toplap: Bào mòn cắt xén - Erosional truncation, Chống nóc).
  2. Khảo sát thuộc tính trường sóng: Xác định kiểu cấu trúc phân lớp (song song, phân kỳ, dạng sigma, chữ S, xếp ngói, hỗn loạn) để xác định năng lượng thủy động lực học môi trường.
  3. Phân tích tướng điện GR: Liên kết biến thiên độ phóng xạ tự nhiên của đất đá với chu kỳ biển tiến (TST), biển lùi (RST) và mực nước biển thấp (LST).

Implementation và kết quả

Development process

Quá trình minh giải được thực hiện qua 4 pha liên hoàn:

  • Pha 1: Chuẩn hóa dữ liệu và liên kết Giếng - Địa chấn (Well-to-Seismic Tie): Sử dụng băng địa chấn tổng hợp (Synthetic Seismogram) tính từ đường cong đo carota âm thanh ($DT$) và mật độ ($\rho$) theo công thức hệ số phản xạ: $$R_i = \frac{\rho_{i+1}V_{i+1} - \rho_i V_i}{\rho_{i+1}V_{i+1} + \rho_i V_i}$$ Nhân chập xung sóng Ricker tần số trội 28 Hz với dãy hệ số phản xạ để tạo trace địa chấn nhân tạo, đạt độ tương quan pha $> 85%$ tại giếng 106-HRN-1X.

  • Pha 2: Minh giải bề mặt địa chấn và chuyển đổi $T-Z$: Xây dựng mô hình vận tốc chuyển đổi từ miền thời gian (Two-Way Time - TWT) sang miền độ sâu thực tế (Depth) dựa trên số liệu đo vận tốc giếng khoan: $$Z(T) = aT^2 + bT + c$$ (với hệ số thực nghiệm hiệu chỉnh cho khu vực Lô 106: $a = 0.00042$, $b = 1.125$, $c = -15.4$).

  • Pha 3: Tự động hóa phân loại Electrofacies bằng Python:

import numpy as np
import pandas as pd
import lasio

def classify_gr_facies(las_file_path, gr_mnemonic='GR', depth_step=0.5):
    """
    Thuật toán nhận diện xu hướng đường cong Gamma Ray (GR Electrofacies)
    để xác định môi trường lắng đọng trầm tích Miocen giữa.
    """
    las = lasio.read(las_file_path)
    df = las.df().dropna(subset=[gr_mnemonic])
    
    # Tính gradient GR trên cửa sổ trượt 5m (10 điểm đo)
    window_size = int(5.0 / depth_step)
    df['GR_smooth'] = df[gr_mnemonic].rolling(window=window_size, center=True).mean()
    df['GR_slope'] = np.gradient(df['GR_smooth'])
    
    facies_labels = []
    for i, row in df.iterrows():
        gr = row['GR_smooth']
        slope = row['GR_slope']
        
        if gr < 45.0:
            facies_labels.append("Clean_Sand_Cylindrical (Dải cát chắn / Lòng sông)")
        elif gr > 90.0:
            facies_labels.append("Shale_Marine (Bùn sét biển nông / Tầng chắn)")
        elif slope > 0.8:
            facies_labels.append("Bell_Shape (Mịn dần lên trên - Phù sa lòng sông / Biển tiến)")
        elif slope < -0.8:
            facies_labels.append("Funnel_Shape (Thô dần lên trên - Tiền châu thổ / Cồn cát cửa sông)")
        else:
            facies_labels.append("Serrated (Răng cưa - Đồng bằng châu thổ xen kẹp sét than)")
            
    df['Facies'] = facies_labels
    return df

Testing và validation

Mô hình minh giải và phân chia ranh giới tập địa tầng được kiểm chứng chéo (Cross-validation) trên toàn bộ mạng lưới 11 giếng khoan trong khu vực nghiên cứu:

Tiêu chuẩn kiểm thử Giá trị mục tiêu Kết quả thực nghiệm đạt được Trạng thái
Độ tương quan Well-Tie ($R$) $\ge 0.80$ $0.88$ (Giếng 106-HRN-1X) Đạt chuẩn xuất sắc
Sai số khép góc lưới tuyến (Closure Error) $< 8\text{ ms}$ $3.5\text{ ms}$ (Lưới 2D Total & PCOSB) Đạt chuẩn
Độ lệch độ sâu tuyệt đối ($Z_{\text{calc}} - Z_{\text{well}}$) $< 10\text{ m}$ $4.2\text{ m}$ tại chiều sâu 2.450 m Chính xác cao
Độ phân giải nhận diện vỉa cát mỏng $\le 12\text{ m}$ $8.5\text{ m}$ trên địa chấn 3D Vượt mục tiêu

Kết quả đạt được

Đồ án đã phân chia và lập bản đồ chi tiết 4 tập địa chấn - trầm tích Miocen giữa (hệ tầng Phù Cừ):

  1. Tập Miocen giữa I (Ranh giới U260 - U240): Bề dày dao động 80 - 220 m. Đặc trưng bởi các pha phản xạ biên độ trung bình - cao, dạng nêm lấn xếp ngói (shingled clinoforms). Môi trường thành tạo: Đồng bằng châu thổ chuyển tiếp sang tiền châu thổ (Prodelta).
  2. Tập Miocen giữa II (Ranh giới U240 - U220): Bề dày 120 - 310 m, phát triển mạnh về phía trũng trung tâm lô 102. Xuất hiện các phản xạ dạng lòng máng xói lở (channel-fill) đặc trưng cho hệ thống kênh ngầm biển thấp (Lowstand Channels).
  3. Tập Miocen giữa III (Ranh giới U220 - U210): Bề dày 60 - 180 m. Pha phản xạ song song, liên tục tốt, biên độ cao. Tướng trầm tích chủ đạo: Cồn cát ven bờ (barrier bars) và bãi cát biển nông, có độ chọn lọc cơ học cao.
  4. Tập Miocen giữa IV (Ranh giới U210 - U200): Bề dày 100 - 250 m. Đặc trưng bởi trường sóng phân lớp song song mỏng đến phản xạ tự do (sét đồng nhất), thể hiện một pha biển tiến cực đại mở rộng toàn bể.
                           LÁT CẮT ĐỊA TẦNG MIOCEN GIỮA
  Độ sâu (m)

Đặc tính vỉa chứa cát kết Miocen giữa đạt chất lượng hồ chứa lý tưởng:

  • Độ rỗng hiệu dụng ($\Phi_e$): Đạt từ 19,3% đến 29,6% (theo số liệu đo điện trở và mật độ giếng khoan).
  • Hàm lượng sét ($V_{\text{shale}}$): Dưới 15% trong các thân cát kênh ngầm và dải cát chắn.
  • Tầng chắn khu vực (Seal): Lớp sét kết nằm trên ranh giới U200 dày 80 - 150 m với hàm lượng khoáng vật sét đạt 71% - 86% (chủ yếu là Illit 55% và Kaolinit 28%), đảm bảo độ kín tuyệt đối cho các bẫy dầu khí.

Đổi mới và đóng góp

  1. Đổi mới phương pháp phân chia địa tầng: Thay vì chia gộp hệ tầng Phù Cừ thành một khối đồng nhất như các nghiên cứu trước năm 2009, đồ án đã bóc tách chi tiết thành 4 phân tập độc lập dựa trên 5 mặt bất chỉnh hợp chuẩn quốc tế, giúp nâng cao độ chính xác minh giải thêm 35%.
  2. Khám phá mô hình tướng kênh ngầm biển thấp (Lowstand Submarine Channels): Chứng minh sự tồn tại của các bẫy địa tầng dạng thấu kính cát lòng sông cổ Miocen giữa tại khu vực trũng Đông Quan và sườn nâng Tiên Lãng - Yên Tử, mở ra đối tượng tìm kiếm dầu khí mới ngoài bẫy nếp lồi truyền thống.
  3. Tích hợp liên ngành độ phân giải cao: Kết hợp thành công dữ liệu địa chấn 3D băng thông rộng với phân tích vi cổ sinh và hình thái điện vi sai GR, giúp giảm thời gian khoanh vùng bẫy chứa từ vài tuần xuống còn vài ngày làm việc.

Ứng dụng thực tế và triển khai

  • Triển khai tại các cấu tạo trọng điểm: Kết quả nghiên cứu áp dụng trực tiếp cho công tác thẩm lượng cấu tạo Thái Bình (Lô 102) và mở rộng diện tích chứa dầu khí tại cấu tạo Yên Tử - Hàm Rồng (Lô 106).
  • Tối ưu hóa vị trí giếng khoan thăm dò: Định vị chính xác tọa độ các đới cát sạch ven biển, giúp giảm thiểu nguy cơ khoan vào đới sét tiền châu thổ bất lợi, tiết kiệm ước tính 8 - 12 triệu USD cho mỗi giếng khoan ngoài khơi.
  • Lộ trình áp dụng:
    • Giai đoạn 1 (Hiện tại): Bàn giao bộ bản đồ đẳng dày, bản đồ tướng trầm tích cho tổ vận hành thăm dò Lô 102 & 106.
    • Giai đoạn 2 (Mở rộng): Áp dụng quy trình chuẩn Sequence Stratigraphy cho các lô lân cận (Lô 103, 107) thuộc miền võng Hà Nội và trũng Bạch Long Vĩ.

Hạn chế và hướng phát triển

  • Hạn chế kỹ thuật: Mật độ mạng lưới địa chấn 3D chưa bao phủ đồng đều toàn bộ phần phía Đông Lô 106; một số tuyến 2D thu nổ trước năm 1990 có tỷ số tín hiệu trên tạp âm ($S/N$) thấp gây khó khăn khi bám sát ranh giới U240 tại đới đứt gãy sâu.
  • Hướng phát triển tiếp theo:
    1. Ứng dụng kỹ thuật biến đổi thuộc tính địa chấn nâng cao (Pre-stack Inversion và AVO - Amplitude Versus Offset) để phân biệt trực tiếp vỉa chứa bão hòa khí và bão hòa nước.
    2. Ứng dụng mô hình học máy (Machine Learning / Deep Learning) để tự động hóa quá trình nhận diện bề mặt nêm lấn và phân vùng tướng địa chấn 3D (Seismic Facies Segmentation).

Đối tượng hưởng lợi

  • Sinh viên ngành Địa chất Dầu khí / Địa vật lý: Nguồn tài liệu tham khảo chuẩn mực về quy trình minh giải địa chấn địa tầng và ứng dụng phần mềm chuyên ngành (Petrel, Techlog).
  • Kỹ sư địa chất & địa vật lý (Geoscientists): Bộ quy trình làm việc (workflow) chi tiết, công thức chuyển đổi $T-Z$ và mã nguồn phân tích Electrofacies có khả năng tái sử dụng cao.
  • Các công ty E&P (PVEP, PVN, Nhà thầu quốc tế): Cung cấp cơ sở khoa học tin cậy để thẩm định trữ lượng, tối ưu hóa danh mục giếng khoan và giảm thiểu rủi ro tài chính trong công tác thăm dò thềm lục địa phía Bắc.

Câu hỏi thường gặp

  1. Yêu cầu kỹ thuật tối thiểu về mặt dữ liệu để triển khai quy trình minh giải này là gì? Cần tối thiểu lưới địa chấn 2D mật độ từ $2 \times 2\text{ km}$ hoặc địa chấn 3D kèm số liệu đo carota (GR, Sonic, Density) từ ít nhất 3 giếng khoan đã hiệu chuẩn độ sâu để thực hiện liên kết Well-Tie chuẩn xác.

  2. Tại sao cát kết Miocen giữa lại có độ rỗng cao (19.3 - 29.6%) dù nằm ở độ sâu trên 2.000 m? Do môi trường lắng đọng năng lượng cao (đồng bằng ven biển, cồn cát chắn) tạo ra độ chọn lọc hạt rất tốt, đồng thời quá trình nén ép cơ học được hạn chế bởi sự xuất hiện sớm của lớp bao bọc xi măng gắn kết mềm và áp suất lỗ rỗng thuận lợi.

  3. Làm thế nào để phân biệt bất chỉnh hợp bào mòn (Erosional Truncation) và giả bất chỉnh hợp trên địa chấn? Cần kết hợp kiểm tra góc tiếp xúc của pha phản xạ cắt cụt trên nhiều tuyến vuông góc (lát cắt 3D) và đối chiếu bước nhảy đột ngột của đường cong vận tốc Sonic/Density tại ranh giới giếng khoan tương ứng.

  4. Tầng chắn sét Miocen giữa có đủ năng lực giữ các thân khí áp suất cao không? Hoàn toàn đủ năng lực. Lớp sét kết Miocen giữa dày 80 - 150 m với tỷ lệ khoáng vật sét mịn $> 71%$ tạo nên áp suất đẩy mao quản cực lớn, đã được chứng minh qua thực tế giữ thân dầu khí tại mỏ Thái Bình và Bạch Long.

  5. Quy trình này có thể áp dụng cho các bể trầm tích khác như Cửu Long hay Nam Côn Sơn không? Nguyên lý Sequence Stratigraphy và phân tích Electrofacies mang tính phổ quát toàn cầu, hoàn toàn áp dụng được cho bể Cửu Long và Nam Côn Sơn bằng cách hiệu chỉnh lại các mốc thời địa tầng và hàm vận tốc $V(z)$ cục bộ.

Kết luận

Đồ án tốt nghiệp "Phân tích đặc điểm môi trường trầm tích Miocen giữa khu vực Đông Bắc đứt gãy Sông Lô" đã giải quyết triệt để bài toán khôi phục cổ môi trường và dự báo không gian phân bố vỉa chứa dầu khí thông qua tích hợp địa chấn địa tầng và địa vật lý giếng khoan. Việc xác lập thành công 4 phân tập địa tầng chi tiết (Tập I đến IV) cùng hệ thống bản đồ tướng trầm tích 3D không chỉ khẳng định tiềm năng chứa vượt trội của cát kết hệ tầng Phù Cừ (độ rỗng 19,3% - 29,6%) mà còn đóng góp luận cứ khoa học then chốt cho chiến lược gia tăng trữ lượng dầu khí tại thềm lục địa vịnh Bắc Bộ. Quý độc giả và các đơn vị nghiên cứu quan tâm có thể ứng dụng trực tiếp khung phân tích và thuật toán này vào các dự án thăm dò khai thác có điều kiện địa chất tương đồng.