Giới thiệu dự án

Khu vực Bắc Bể Sông Hồng (BBSH), đặc biệt là diện tích Đông Bắc đứt gãy Sông Lô thuộc Lô 102 và 106 thềm lục địa Vịnh Bắc Bộ, là một trong những bồn trũng trầm tích Kainozoi có tiềm năng dầu khí quan trọng của Việt Nam. Dù đã có lịch sử tìm kiếm thăm dò từ đầu thập niên 1980 qua các chiến dịch thu nổ địa chấn 2D/3D và khoan thăm dò của Tổng cục Dầu khí Việt Nam, Total, Idemitsu, PCOSB và PVEP (với các phát hiện thương mại tại cấu tạo Thái Bình, Yên Tử, Hàm Rồng, Bạch Long), công tác dự báo các bẫy phi cấu tạo và bẫy địa tầng vẫn đối mặt với rủi ro địa chất cao.

Trong các lát cắt địa tầng Kainozoi, bên cạnh đối tượng móng đá vôi nứt nẻ trước Đệ Tam và cát kết Oligocen chôn vùi sâu bị suy giảm độ rỗng do nén ép, cát kết tuổi Miocen giữa (Hệ tầng Phù Cừ - $N_1^2 pc$) nổi lên như một tầng chứa có chất lượng tối ưu với độ rỗng đạt từ 19,3% đến 29,6%. Tuy nhiên, sự phân bố tướng thạch học và môi trường lắng đọng của Miocen giữa có tính dị hướng cao do chịu ảnh hưởng bởi các chu kỳ biển tiến - biển lùi và hoạt động kiến tạo trượt bằng nghịch đảo cục bộ.

Vấn đề cốt lõi của công tác thăm dò tại khu vực là: Làm thế nào để xác định chính xác ranh giới không gian của các thể cát chứa (reservoir sand bodies), bóc tách các pha phân tập trầm tích và khoanh vùng tướng đá có triển vọng dầu khí cao từ tập dữ liệu địa chấn và địa vật lý giếng khoan (ĐVLGK) không đồng nhất?

Đồ á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ô" được thực hiện nhằm giải quyết bài toán trên thông qua các mục tiêu cụ thể:

  1. Chuẩn hóa và tích hợp cơ sở dữ liệu địa chấn 2D/3D cùng tài liệu ĐVLGK của 11 giếng khoan trọng điểm (102-TB-1X, 102-CQ-1X, 106-HRN-1X, 107-BAL-1X,...).
  2. Thiết lập khung địa tầng phân tập chi tiết cho tầng Miocen giữa, phân chia thành 4 tập trầm tích thứ cấp: Tập I ($U_{260}-U_{240}$), Tập II ($U_{240}-U_{220}$), Tập III ($U_{220}-U_{210}$) và Tập IV ($U_{210}-U_{200}$).
  3. Phân tích tướng địa chấn (Seismic Facies) kết hợp giải đoán dạng đường cong Gamma Ray (GR Electrofacies) để khôi phục cổ địa lý và môi trường lắng đọng trầm tích.
  4. Xây dựng hệ thống bản đồ đẳng thời (Time map), bản đồ đẳng sâu (Depth map), bản đồ bề dày (Isochore map) và bản đồ phân bố tướng trầm tích nhằm dự báo các diện tích bẫy cát kết triển vọng.

Phạm vi nghiên cứu tập trung vào khoảng địa tầng Miocen giữa tại khu vực Đông Bắc đứt gãy Sông Lô, giới hạn giữa đứt gãy Sông Lô (phía Tây Nam) và đứt gãy Hải Dương / đới nâng Quảng Ninh (phía Đông Bắc), ở độ sâu mực nước biển dao động từ 20 đến 35 mét.

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 địa chất tại Bắc Bể Sông Hồng chủ yếu dựa vào phương pháp địa tầng thạch học truyền thống (Lithostratigraphy) và minh giải cấu trúc thuần túy, dẫn đến việc bỏ sót các bẫy địa tầng tinh vi như doi cát ngầm biển nông (offshore sand bars), nêm lấn tiền châu thổ (prodelta clinoforms) hay bẫy vát nhọn kề áp (onlap/toplap traps).

Tiêu chí so sánh Phương pháp Địa tầng Thạch học truyền thống Phương pháp Địa chấn Địa tầng & Phân tập (Nghiên cứu áp dụng) Ứng dụng Học máy/AI phân loại tướng địa chấn
Cơ sở ranh giới Thay đổi thành phần đất đá tại giếng khoan Các mặt bất chỉnh hợp (Unconformities) & chu kỳ mực nước biển Thuật toán phân cụm biên độ sóng đa chiều
Độ tin cậy liên kết Thấp ở vùng khoảng cách giếng khoan xa Rất cao nhờ liên kết pha sóng địa chấn liên tục Trung bình (phụ thuộc vào dữ liệu gán nhãn)
Khả năng dự báo tướng Cục bộ tại thân giếng khoan Rộng khắp không gian 2D/3D toàn bồn trũng Tự động hóa cao nhưng thiếu kiểm soát địa chất
Phát hiện bẫy phi cấu tạo Hạn chế, dễ nhầm lẫn ranh giới giả Tối ưu cho bẫy ngầm biển thấp (LST), bẫy nêm vát Khá tốt đối với các dị thường biên độ

Yêu cầu kỹ thuật của quy trình minh giải được chuẩn hóa theo mô hình MoSCoW:

  • Must have (Bắt buộc): Buộc phải liên kết chính xác giếng khoan với địa chấn (Well-to-seismic tie) qua băng địa chấn tổng hợp (Synthetic Seismogram); xác định 5 mặt ranh giới phản xạ chính ($U_{260}, U_{240}, U_{220}, U_{210}, U_{200}$); phân loại dạng đường cong Gamma Ray theo chuẩn Emery & Kenneth.
  • Should have (Nên có): Thiết lập hàm chuyển đổi thời gian - độ sâu (Time-to-Depth conversion) từ tài liệu đo vận tốc Checkshot giếng khoan 106-HRN-1X; tính toán bản đồ chiều dày thực tế cho 4 tập trầm tích.
  • Could have (Có thể có): Tích hợp phân tích thạch học từ mẫu vụn (ditch cuttings) và mẫu lõi sườn (sidewall cores) để định lượng tỷ lệ khoáng vật sét (Illite, Kaolinite).
  • Won't have (Chưa thực hiện): Nghịch đảo địa chấn đàn hồi trước cộng (Pre-stack seismic inversion) do hạn chế về góc tới của dữ liệu địa chấn 2D cũ.

Thiết kế hệ thống minh giải tích hợp

Hệ thống quy trình làm việc (Interpretation Workflow) kết hợp giữa địa chấn địa tầng học (Mitchum & Vail) và phân tích log giếng khoan được thiết kế theo cấu trúc modular:

[Dữ liệu Địa chấn 2D/3D (SEG-Y)] + [Tài liệu Well Logs (LAS/DLIS)]

Technology Stack sử dụng trong nghiên cứu:

  • Nền tảng minh giải địa chấn chuyên dụng: Schlumberger Petrel E&P Software Platform (v2018.2) / Landmark DecisionSpace 5000.10.
  • Xử lý tài liệu địa vật lý giếng khoan: Interactive Petrophysics (IP v4.5) và Schlumberger Techlog (v2017.1).
  • Công cụ phân tích dữ liệu & tự động hóa: Python 3.9 kết hợp các thư viện chuyên ngành (lasio v0.30, obspy v1.3.0, numpy v1.23, matplotlib v3.6, scipy v1.9).
  • Tiêu chuẩn dữ liệu: Định dạng SEG-Y chuẩn SEG Rev-1 cho dữ liệu trường sóng; LAS 2.0/3.0 cho dữ liệu đường cong carota giếng khoan.

Methodology

Quy trình địa tầng phân tập được triển khai theo các pha kiểm soát chất lượng nghiêm ngặt:

  1. Pha 1 - Đánh giá ranh giới tập (Sequence Boundaries - SB): Nhận diện các mặt bất chỉnh hợp đáy (Baselap: onlap, downlap) và bất chỉnh hợp nóc (Toplap: erosional truncation) trên các tuyến địa chấn ngang và dọc bồn trũng.
  2. Pha 2 - Phân tích tướng địa chấn (Seismic Facies Mapping): Phân loại cấu trúc phản xạ nội bộ thành: dạng phân lớp song song (parallel), dạng nêm lấn hình sigma / chữ S (clinoform progradation), dạng xiên chéo (oblique/shingled), và dạng hỗn độn (chaotic).
  3. Pha 3 - Phân tích điện sai & Tướng ĐVLGK (Electrofacies Analysis): Giải đoán cơ chế bồi lắng thông qua xu thế hạt thô dần lên (coarsening-upward / funnel shape) hay hạt mịn dần lên (fining-upward / bell shape) trên đường cong Gamma Ray.

Implementation và kết quả

Development process & Core Algorithms

Để minh giải chính xác mối liên hệ giữa thời gian phản xạ sóng địa chấn hai chiều (Two-Way Time - TWT, ms) và chiều sâu địa tầng thực tế (Depth, m), thuật toán tính toán hệ số phản xạ và tích chập sóng âm địa chấn tổng hợp được triển khai kết hợp với hàm chuyển đổi vận tốc từ giếng khoan kiểm chứng 106-HRN-1X.

Hệ số phản xạ sóng âm $R_i$ tại ranh giới giữa lớp $i$ và $i+1$ được xác định theo trở sóng âm thanh ($Z = \rho \cdot V$): $$R_i = \frac{\rho_{i+1}V_{i+1} - \rho_i V_i}{\rho_{i+1}V_{i+1} + \rho_i V_i} = \frac{Z_{i+1} - Z_i}{Z_{i+1} + Z_i}$$

Băng địa chấn tổng hợp $S(t)$ được tạo ra bằng phép tích chập giữa chuỗi xung phản xạ $R(t)$ với sóng con Ricker $W(t)$: $$S(t) = R(t) * W(t) = \int_{-\infty}^{\infty} R(\tau) W(t - \tau) d\tau$$

Hàm chuyển đổi vận tốc Time-to-Depth được xây dựng dựa trên kết quả đo vận tốc giếng khoan 106-HRN-1X thông qua phương trình hồi quy đa thức bậc hai: $$Depth(m) = 0.00042 \cdot TWT^2 + 0.8125 \cdot TWT + 125.4$$

Dưới đây là đoạn mã Python chuẩn hóa xử lý tín hiệu đường cong Gamma Ray, tự động bóc tách chỉ số sét $V_{shale}$ theo mô hình Larionov cho đá trầm tích Đệ Tam và phân loại tướng điện sai phục vụ minh giải:

import numpy as np
import lasio
import pandas as pd

def process_gamma_ray_electrofacies(las_file_path: str, gr_clean: float = 25.0, gr_shale: float = 120.0):
    """
    Xử lý đường cong Gamma Ray (GR), tính toán thể tích sét (Vshale)
    theo mô hình Larionov cho đá Đệ Tam và nhận diện tướng trầm tích.
    """
    # Đọc dữ liệu giếng khoan chuẩn LAS 2.0
    las = lasio.read(las_file_path)
    df = las.df().reset_index()
    
    # Kiểm tra sự tồn tại của kênh Gamma Ray
    if 'GR' not in df.columns:
        raise ValueError("Kênh Gamma Ray (GR) không tồn tại trong tệp LAS.")
        
    # Tính toán chỉ số Gamma Ray phân đoạn (IGR)
    df['IGR'] = (df['GR'] - gr_clean) / (gr_shale - gr_clean)
    df['IGR'] = df['IGR'].clip(lower=0.0, upper=1.0)
    
    # Tính toán thể tích sét Vshale (Larionov Tertiary Formula)
    # Vsh = 0.083 * (2^(3.7 * IGR) - 1)
    df['VSHALE_LARIONOV'] = 0.083 * (np.power(2.0, 3.7 * df['IGR']) - 1.0)
    df['VSHALE_LARIONOV'] = df['VSHALE_LARIONOV'].clip(lower=0.0, upper=1.0)
    
    # Phân loại tướng điện sai sơ bộ dựa trên đặc tính Gamma Ray
    conditions = [
        (df['GR'] < 45.0),                                      # Cát sạch ven bờ / Lòng sông cổ
        (df['GR'] >= 45.0) & (df['GR'] < 75.0),                 # Cát bột kết xen kẽ / Cồn cát cửa sông
        (df['GR'] >= 75.0) & (df['GR'] < 105.0),                # Bột kết / Sét tiền châu thổ
        (df['GR'] >= 105.0)                                     # Sét dày biển nông / Đầm hồ
    ]
    facies_labels = [
        'Clean_Sand_Facies',
        'Interbedded_Sand_Silt_Facies',
        'Prodelta_Mud_Facies',
        'Marine_Shale_Seal_Facies'
    ]
    df['ELECTRO_FACIES'] = np.select(conditions, facies_labels, default='Undetermined')
    
    return df[['DEPT', 'GR', 'IGR', 'VSHALE_LARIONOV', 'ELECTRO_FACIES']]

# Khởi chạy hàm mẫu với giếng khoan kiểm tra
# well_data = process_gamma_ray_electrofacies("106_HRN_1X_Logs.las")

Testing và validation

Quy trình minh giải địa tầng phân tập được kiểm chứng chéo (Cross-Validation) qua 11 giếng khoan trong toàn khu vực nghiên cứu:

  • Độ chính xác Well-Tie: Hệ số tương quan pha sóng (Wavelet Cross-Correlation Coefficient) giữa băng địa chấn tổng hợp và dữ liệu thực tế tại giếng 106-HRN-1X và 102-TB-1X đạt $r = 0.892$.
  • Giảm sai số đóng vòng tuyến địa chấn (Mistie Resolution): Sai số đóng vòng tại các điểm giao cắt giữa mạng lưới tuyến địa chấn 2D (khảo sát Total 1989-1990 và PCOSB 2007) được hiệu chỉnh từ $\pm 18\text{ ms}$ xuống $< 3.5\text{ ms}$.
  • Độ phân giải tầng chứa: Độ phân giải thẳng đứng địa chấn (Vertical Resolution) tại khoảng độ sâu Miocen giữa ($1600 - 2800\text{ m}$) đạt xấp xỉ $\lambda/4 \approx 18 - 22\text{ m}$, cho phép nhận diện rõ nét các thân cát bồi tụ dày từ 15m trở lên.

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

Nghiên cứu đã làm sáng tỏ cấu trúc không gian và môi trường trầm tích của 4 tập địa chấn Miocen giữa ($N_1^2 pc$) khu vực Đông Bắc đứt gãy Sông Lô:

  1. Tập Miocen giữa I ($U_{260} - U_{240}$):
    • Bề dày dao động từ 150m đến 450m, dày nhất tại trũng Đông Quan và mỏng dần về phía nâng Tiên Lãng.
    • Trường sóng đặc trưng bởi phản xạ xiên chéo, biên độ từ trung bình đến cao, ranh giới đáy gá đáy (onlap) lên mặt bất chỉnh hợp $U_{260}$.
    • Môi trường lắng đọng: Đồng bằng châu thổ (Delta plain) chuyển tiếp sang đồng bằng ngập lũ; cát kết dạng thấu kính có độ rỗng trung bình 22 - 25%.
  2. Tập Miocen giữa II ($U_{240} - U_{220}$):
    • Bề dày phát triển từ 180m đến 520m; phản ánh thời kỳ phát triển cực thịnh của châu thổ sông.
    • Thể hiện các phản xạ nêm lấn dạng chữ S (Sigmoid clinoforms), góc dốc nghiêng thoái lui dần về phía Đông Nam. Tướng GR dạng phễu (funnel) hạt thô dần lên.
    • Môi trường: Sườn châu thổ (Delta front) và cồn cát cửa sông (Mouth bars). Đây là tầng chứa cát kết có tiềm năng tốt nhất với diện phân bố rộng khắp lô 102 và Tây Lô 106.
  3. Tập Miocen giữa III ($U_{220} - U_{210}$):
    • Bề dày từ 100m đến 320m. Trường sóng thể hiện dạng gò đồi xen kẽ phân lớp song song gián đoạn.
    • Dạng đường cong GR có dạng hình trụ (cylindrical) kết hợp dạng răng cưa, minh chứng cho các đê cát ngập nước và bãi cát biển nông ven bờ.
  4. Tập Miocen giữa IV ($U_{210} - U_{200}$):
    • Bề dày từ 80m đến 250m. Phản xạ song song, liên tục tốt, biên độ trung bình - thấp, đại diện cho pha biển tiến mở rộng (Transgressive phase) trước khi kết thúc Miocen giữa. Phía nóc bị bào mòn cục bộ bởi mặt bất chỉnh hợp $U_{200}$.
    • Tướng sét bột kết chiếm ưu thế, đóng vai trò tầng chắn khu vực trực tiếp cho các thân cát phía dưới.

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

Nghiên cứu mang lại những đóng góp khoa học và thực tiễn mang tính đột phá cho công tác địa chất dầu khí tại bể Sông Hồng:

                         CẢI TIẾN ĐỘ CHÍNH XÁC
                     CỦA QUY TRÌNH MINH GIẢI MỚI
                     
   Độ chính xác liên kết địa tầng:  [====================] +42.5%
   Giảm sai số độ sâu ranh giới:    [==================  ] -38.0%
   Hiệu quả khoanh vùng bẫy cát:    [====================] +55.0%
   Độ tin cậy mô hình cổ địa lý:    [=================   ] +35.0%
  • Đổi mới phương pháp luận: Thay vì xem Miocen giữa như một đối tượng chứa đồng nhất, đề tài đã thành công trong việc bóc tách lát cắt thành 4 đơn vị phân tập chu kỳ bậc 3 (3rd-order sequences), làm sáng tỏ quy luật dịch chuyển của đường bờ cổ từ Tây Bắc xuống Đông Nam.
  • Minh chứng bẫy phi cấu tạo: Khẳng định sự tồn tại của hệ thống bẫy địa tầng dạng cát nêm lấn tiền châu thổ (Prodelta sand pinches) và cát lòng sông cổ ngầm biển thấp (Lowstand channel sands) tại Lô 102 & 106 – tương tự mô hình đã mang lại phát hiện khí tại giếng 107-BAL-1X và Yên Tử-1X.
  • Tối ưu hóa tham số thạch học: Cung cấp bộ thông số định lượng chính xác về độ rỗng ($\Phi = 19,3 - 29,6%$) và thể tích sét ($V_{sh} = 15 - 35%$) của các thân cát Miocen giữa, phục vụ trực tiếp cho công tác mô phỏng vỉa (Reservoir Simulation).

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

Tình huống ứng dụng thực tế trong E&P

  1. Định vị giếng khoan thăm dò/thẩm lượng: Hỗ trợ các nhà thầu dầu khí (PVEP, Vietsovpetro, Eni,...) xác định tọa độ thả giếng tối ưu nhằm xuyên qua tâm của các cồn cát cửa sông (mouth bar sands) có bề dày hiệu dụng cao nhất trong Tập II ($U_{240}-U_{220}$).
  2. Đánh giá rủi ro địa chất (Geological Risk Assessment): Giảm thiểu xác suất khoan vào tầng cát chứa nước hoặc sét tiền châu thổ dày thông qua bản đồ phân bố tướng trầm tích chi tiết.

Phân tích Hiệu quả Đầu tư & Chi phí (Cost-Benefit Analysis)

Chi phí khoan một giếng thăm dò ngoài khơi tại khu vực BBSH dao động từ 15 đến 25 triệu USD. Việc áp dụng mô hình phân tích tướng địa chấn địa tầng phân tập giúp nâng tỷ lệ thành công của các giếng khoan thăm dò bẫy địa tầng từ 40% lên 65-70%, mang lại giá trị thặng dư kinh tế hàng triệu USD cho các chiến dịch thăm dò.

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

Dù đạt được các kết quả mang tính khoa học và thực tiễn cao, đề tài vẫn ghi nhận một số hạn chế kỹ thuật:

  • Mật độ mạng lưới địa chấn không đồng đều: Vùng phía Đông Lô 106 chủ yếu dựa trên lưới địa chấn 2D thưa ($4 \times 6\text{ km}$ của Total từ năm 1989), dẫn đến độ chi tiết của bản đồ đẳng thời tại khu vực này thấp hơn khu vực trung tâm đã có địa chấn 3D.
  • Hạn chế về tài liệu mẫu lõi: Số lượng giếng khoan có mẫu lõi liên tục (conventional cores) qua tầng Miocen giữa còn ít, việc định danh thạch học chi tiết chủ yếu dựa trên mẫu vụn và log địa vật lý.

Hướng phát triển tiếp theo:

  1. Nghịch đảo địa chấn biên độ đàn hồi (AVO & Seismic Inversion): Ứng dụng công nghệ Simultaneous Inversion để dự báo trực tiếp tích số mật độ - vận tốc ($V_p/V_s$) và xác suất có mặt của chất lưu khí/condensate trong vỉa cát.
  2. Ứng dụng Trí tuệ Nhân tạo (AI & Deep Learning): Xây dựng mạng nơ-ron tích chập 3D (3D-CNN) để tự động nhận diện và phân vùng ranh giới các đứt gãy kiến tạo và tướng địa chấn phức tạp.

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

  • Sinh viên & Học viên cao học ngành Kỹ thuật Địa chất / Địa vật lý: Cung cấp tài liệu tham khảo bài bản, mẫu mực về phương pháp luận kết hợp giữa địa chấn địa tầng học cổ điển và kỹ thuật minh giải hiện đại.
  • Kỹ sư Địa vật lý & Địa chất mỏ (Geoscientists): Nắm bắt các quy tắc giải đoán đường cong Gamma Ray, quy trình tính toán chuyển đổi Time-Depth và kỹ năng nhận diện bất chỉnh hợp trên lát cắt sóng địa chấn thực tế.
  • Các Công ty Dầu khí & Cơ quan Quản lý (PVEP, PVN, VPI): Cơ sở dữ liệu và bản đồ chuyên đề đáng tin cậy phục vụ đánh giá tiềm năng dầu khí, xây dựng kế hoạch mở thầu Lô và tối ưu hóa vị trí các giếng khoan phát triển mỏ tại Vịnh Bắc Bộ.
  • Nhà nghiên cứu Khoa học Trái đất: Củng cố mô hình tiến hóa bồn trũng Kainozoi và cơ chế tác động của hoạt động nghịch đảo kiến tạo đứt gãy Sông Lô đối với hệ sinh thái trầm tích.

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

1. Yêu cầu cấu hình phần cứng và phần mềm tối thiểu để chạy quy trình minh giải này là gì?

Hệ thống yêu cầu máy trạm Workstation chuyên dụng chạy Red Hat Enterprise Linux 7/8 hoặc Windows 10/11 Pro 64-bit; CPU tối thiểu 16 Cores (Intel Xeon hoặc AMD Threadripper); RAM tối thiểu 64 GB (khuyến nghị 128 GB cho khối dữ liệu địa chấn 3D); Card đồ họa chuyên dụng NVIDIA Quadro RTX 4000 hoặc RTX A4000 (tối thiểu 8 GB VRAM hỗ trợ OpenGL 4.5); Bộ nhớ lưu trữ NVMe SSD tối thiểu 2 TB để đảm bảo tốc độ đọc/ghi dữ liệu SEG-Y khối lượng lớn. Các phần mềm chuyên dụng cần có: Petrel 2018+, Techlog, Interactive Petrophysics.

2. Làm thế nào để phân biệt chính xác giữa hiện tượng gá đáy (Onlap) và phủ đáy (Downlap) trong điều kiện bồn trũng bị biến dạng nén ép nghịch đảo?

Trong các khu vực chịu ảnh hưởng của pha nghịch đảo kiến tạo cuối Miocen như BBSH, thế nằm nguyên thủy của các lớp đá bị uốn nếp hoặc nghiêng dốc, dễ gây nhầm lẫn giữa Onlap (gá đáy ở đầu nguồn/đường bờ) và Downlap (phủ đáy ở chân sườn thềm/hướng ra biển). Giải pháp là phải thực hiện khôi phục mặt cắt địa chất cân bằng (Structural Restoration / Flattening) về mặt ranh giới đáy chuẩn trước biến dạng ($U_{260}$ hoặc nóc móng), từ đó xác định chính xác góc dốc lắng đọng ban đầu của các pha phản xạ.

3. Tầng chắn (Cap rock) cho các vỉa cát kết Miocen giữa tại Lô 102 & 106 có đảm bảo độ kín để giữ các tích tụ khí/dầu không?

Tầng chắn khu vực cho Miocen giữa là các tập sét và sét-bột kết phân bố đồng nhất với chiều dày hiệu dụng dao động từ 80 đến 150m (gặp tại các giếng 112-BT-1XR, 106-HR-1X). Phân tích khoáng vật cho thấy hàm lượng khoáng vật sét đạt từ 71% đến 86% (chủ yếu là Illite 55% và Kaolinite 28%), tạo nên áp suất bẻ vỡ thủy tĩnh cao, đảm bảo khả năng chắn kín tuyệt đối cho các bẫy cát kết bên dưới.

4. Chiều dày tối thiểu của một thân cát để có thể quan sát thấy trên băng địa chấn 2D/3D là bao nhiêu?

Theo tiêu chuẩn Rayleigh, giới hạn độ phân giải thẳng đứng có thể phát hiện được (limit of resolution) xấp xỉ bằng $1/4$ bước sóng ($\lambda/4$). Với tần số ưu thế của trường sóng địa chấn tại tầng Miocen giữa là $f = 30 - 35\text{ Hz}$ và vận tốc truyền sóng trung bình trong cát kết $V_p \approx 2600\text{ m/s}$, bước sóng $\lambda = V/f \approx 74 - 86\text{ m}$. Do đó, chiều dày tối thiểu của thân cát để phân biệt được nóc và đáy trên địa chấn là khoảng $18 - 22\text{ m}$. Các lớp mỏng hơn ($< 10\text{ m}$) chỉ có thể phát hiện gián tiếp thông qua dị thường biên độ sóng (tuning effect) kết hợp với minh giải chi tiết từ log Gamma Ray giếng khoan.

5. Dữ liệu vận tốc từ giếng khoan 106-HRN-1X có thể áp dụng chung cho toàn bộ Lô 102 và 106 không?

Hàm chuyển đổi vận tốc từ giếng khoan 106-HRN-1X có độ tin cậy cao nhất tại khu vực đới nâng và các địa hào lân cận thuộc Lô 106. Khi mở rộng sang trung tâm trũng Đông Quan (Lô 102) - nơi có bề dày trầm tích Kainozoi sâu hơn và tốc độ nén chặt nhanh hơn, cần phải hiệu chỉnh mô hình vận tốc cục bộ bằng cách nội suy kết hợp với số liệu Checkshot từ các giếng 102-TB-1X và 102-CQ-1X để tránh sai lệch độ sâu tuyệt đối (depth error) vượt quá sai số cho phép ($< 1.5%$).

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ô" của tác giả Hà Diệu Linh đã hoàn thành xuất sắc các mục tiêu nghiên cứu đặt ra. Bằng việc kết hợp nhuần nhuyễn giữa lý thuyết địa chấn địa tầng học, địa tầng phân tập của trường phái Exxon và kỹ thuật minh giải địa vật lý giếng khoan hiện đại, công trình đã:

  • Tái hiện chi tiết bức tranh cổ địa lý bồn trũng BBSH trong giai đoạn Miocen giữa với 4 tập trầm tích chu kỳ rõ rệt.
  • Khẳng định cát kết tướng sườn châu thổ (Delta front) và cồn cát cửa sông thuộc Tập II ($U_{240}-U_{220}$) là tầng chứa tối ưu với độ rỗng vỉa lên tới 29,6%.
  • Mở ra hướng tiếp cận mới trong việc săn tìm các bẫy dầu khí phi cấu tạo, góp phần gia tăng tài nguyên và định hướng hiệu quả cho các chiến dịch khoan thăm dò tương lai tại thềm lục địa Việt Nam.