Tổng quan nghiên cứu

Thăm dò địa chấn phản xạ là phương pháp địa vật lý giữ vai trò then chốt trong công tác nghiên cứu cấu trúc sâu của vỏ Trái Đất và tìm kiếm thăm dò dầu khí với độ sâu khảo sát hiệu quả từ 300 mét đến hơn 10.000 mét (10 km). Trong quá trình thu nổ thực địa trên các tuyến địa chấn 2D và mạng lưới 3D, sự khác biệt về khoảng cách giữa nguồn phát và máy thu (offset) tạo ra độ trễ thời gian truyền sóng theo quy luật hypebol, gây biến dạng trục đồng pha và làm suy giảm nghiêm trọng tỷ số tín hiệu trên nhiễu khi thực hiện cộng sóng.

Luận văn thạc sĩ khoa học chuyên ngành Vật lý địa cầu (mã số 8440130.06) của học viên Nguyễn Quang Sơn, dưới sự hướng dẫn của Tiến sĩ Nguyễn Đức Vinh tại Trường Đại học Khoa học Tự nhiên – Đại học Quốc gia Hà Nội năm 2019, đã tập trung giải quyết bài toán cốt lõi: Tính toán hiệu chỉnh động trong thăm dò địa chấn. Mục tiêu trọng tâm của nghiên cứu là xây dựng quy trình tính toán số gia pháp tuyến, nắn thẳng các trục đồng pha về mức thời gian pháp tuyến chuẩn và ứng dụng phương pháp quét phổ tốc độ không đổi để tối ưu hóa chất lượng lát cắt địa chấn.

Phạm vi nghiên cứu bao gồm việc phân tích lý thuyết trường sóng phản xạ điểm sâu chung và thử nghiệm mô phỏng trên 3 mô hình địa chất phân lớp điển hình có bổ sung phông nhiễu ngẫu nhiên. Ý nghĩa thực tiễn của công trình thể hiện qua việc nâng cao biên độ mạch cộng tổng hợp đạt mức cực đại 100%, tăng tỷ số tín hiệu trên nhiễu từ 2,5 đến 3,8 lần và triệt tiêu hơn 90% các sai số pha, mang lại công cụ tính toán tin cậy cho công tác minh giải địa chất dầu khí tại thềm lục địa Việt Nam.

Cơ sở lý thuyết và phương pháp nghiên cứu

Khung lý thuyết áp dụng

Nghiên cứu được phát triển dựa trên lý thuyết truyền sóng đàn hồi trong môi trường địa chất phân lớp và nguyên lý phương pháp địa chấn điểm sâu chung (Common Depth Point - CDP). Khi sóng địa chấn gặp các mặt ranh giới có sự biến đổi về độ kháng âm (tích số giữa mật độ đất đá và tốc độ truyền sóng), sóng phản xạ thứ cấp hình thành và truyền ngược về bề mặt quan sát.

Các khái niệm và mô hình toán học trọng tâm bao gồm:

  1. Biểu đồ thời khoảng tuyến dọc: Biểu diễn mối quan hệ giữa thời gian truyền sóng và tọa độ điểm quan sát theo phương trình hypebol: $t = \sqrt{t_0^2 + \frac{x^2}{v^2}}$ Trong đó, $t_0 = \frac{2h}{v}$ là thời gian pháp tuyến tại điểm thu nổ trùng nhau ($x = 0$), $h$ là chiều sâu ranh giới và $v$ là tốc độ truyền sóng trung bình của môi trường phủ.

  2. Số gia pháp tuyến và hiệu chỉnh động (NMO): Số gia thời gian $\Delta t = t(x) - t_0$ phản ánh độ trễ do khoảng cách thu nổ $x$ tạo ra. Khai triển chuỗi Taylor cho thấy lượng hiệu chỉnh động biến thiên theo quy luật hàm mũ giảm dần theo thời gian $t_0$ và tăng dần theo bình phương khoảng cách $x^2$.

  3. Vận tốc điểm sâu chung ($v_{\text{ĐSC}}$): Khi mặt phản xạ nằm nghiêng một góc $\varphi$, vận tốc hiệu dụng được chuẩn hóa qua biểu thức $v_{\text{ĐSC}} = \frac{v}{\cos\varphi}$. Trục đồng pha chỉ được nắn thẳng tuyệt đối khi vận tốc lựa chọn trùng khớp hoàn toàn với vận tốc điểm sâu chung.

  4. Tương quan liên kết trong tinh chỉnh sai số: Sử dụng hàm tương quan chéo giữa các cặp đường ghi địa chấn trong cửa sổ thời gian từ 0,5 giây đến 1,0 giây để xác định sai số thời gian dư $\tau_n$, loại bỏ độ lệch pha cục bộ có độ phân tán dưới 30% chu kỳ biểu kiến.

Phương pháp nghiên cứu

Nghiên cứu sử dụng nguồn dữ liệu mô phỏng tổng hợp gồm 48 mạch ghi địa chấn cho mỗi tập hợp điểm giữa chung (CMP gather), phân bố trên 3 mô hình cấu trúc địa tầng nhân tạo có dải độ sâu từ 500 mét đến 3.000 mét. Phương pháp chọn mẫu áp dụng hệ thống quan sát đối xứng điểm nổ và điểm thu dọc theo tuyến khảo sát 2D, khống chế tỷ lệ khoảng cách thu nổ cực đại trên chiều sâu ranh giới ($x_{\max}/h$) trong khoảng 0,5 đến 1,2 nhằm hạn chế hiện tượng méo dạng hình học khi góc nghiêng địa tầng vượt quá 10 độ.

Quy trình phân tích số liệu được tiến hành qua 3 giai đoạn chính:

  • Tiền xử lý và hiệu chỉnh tĩnh: Đưa các điểm phát và thu sóng về cùng một mặt phẳng chuẩn quy ước (với độ sâu nguồn phát trên biển khoảng 6 mét và cáp thu khoảng 8 mét), triệt tiêu ảnh hưởng của đới tốc độ nhỏ có chiều dày từ 5 mét đến 25 mét.
  • Phân tích phổ tốc độ quét không đổi (CVS): Chia băng địa chấn thành các đoạn thời gian với các đầu mút $t_{0i}$, quét đồng thời $m$ giá trị vận tốc tiên nghiệm với bước quét từ 10 m/s đến 20 m/s.
  • Tính toán hàm năng lượng cực đại: Đánh giá biên độ cộng sóng trên từng cửa sổ thời gian $\Delta t$ nhằm xác định chính xác vận tốc hiệu dụng $v_{\text{hd}}$. Phương pháp này được lựa chọn vì tính ổn định cao, khả năng tự động hóa vượt trội và triệt tiêu hiệu quả các xung nhiễu ngẫu nhiên trong môi trường tính toán số lượng lớn.

Kết quả nghiên cứu và thảo luận

Những phát hiện chính

Quá trình thử nghiệm tính toán trên các mô hình số đã mang lại 4 phát hiện khoa học quan trọng:

  1. Độ nhạy của biên độ cộng sóng theo vận tốc hiệu chỉnh: Trên Mô hình 1 với vận tốc thực $v = 700\text{ m/s}$, khi áp dụng chính xác vận tốc 700 m/s, biểu đồ thời khoảng được nắn phẳng hoàn toàn, giúp biên độ mạch tổng sau khi cộng sóng đạt cực đại 100%. Ngược lại, khi vận tốc chọn sai lệch ở mức $v = 680\text{ m/s}$ (nắn chưa tới) hoặc $v = 730\text{ m/s}$ (nắn quá mức), biên độ mạch tổng bị suy giảm từ 32% đến 46% do các dao động lệch pha triệt tiêu lẫn nhau.

  2. Tác động của góc nghiêng ranh giới phản xạ: Khi góc nghiêng địa tầng $\varphi$ tăng từ 5 độ lên 15 độ, biểu đồ thời khoảng điểm sâu chung bị doãng rộng rõ rệt. Độ dịch chuyển của yếu tố phản xạ $MA$ mở rộng vượt mức 0,15 lần khoảng cách thu nổ lớn nhất ($x_{\max}$), đòi hỏi vận tốc tính toán phải được nâng lên theo tỷ lệ nghịch với $\cos\varphi$ để tránh hiện tượng uốn cong dư dạng parabol $t = t_0 + ax^2 + bx^4$.

  3. Khả năng phân tách tín hiệu trong môi trường nhiễu cao: Trên Mô hình 2 và Mô hình 3 có bổ sung 20% đến 30% phông nhiễu trắng ngẫu nhiên, phương pháp quét phổ tốc độ dạng đồ thị năng lượng vẫn phân lập chính xác các cực đại phổ, khôi phục trục đồng pha với sai số thời gian dưới 2 mili-giây.

  4. Hiệu quả đồng bộ giữa hiệu chỉnh tĩnh và hiệu chỉnh động: Sự kết hợp giữa hiệu chỉnh tĩnh (xử lý sai số đới tốc độ nhỏ dày 10 đến 25 mét) và hiệu chỉnh động chuẩn xác giúp loại bỏ từ 92% đến 96% hiện tượng méo dạng hypebol, nâng cao độ phân giải thẳng đứng của lát cắt địa chấn thêm 35%.

Thảo luận kết quả

Nguyên nhân cốt lõi dẫn đến sự suy giảm biên độ mạch tổng khi chọn sai vận tốc là do sự biến dạng phi tuyến của trường sóng ở khoảng cách thu nổ lớn. Một sai số vận tốc nhỏ khoảng 4,28% (tương đương chênh lệch 30 m/s) đã đủ tạo ra độ trễ pha vượt quá 1/4 chu kỳ sóng biểu kiến, dẫn đến hiện tượng giao thoa triệt tiêu năng lượng khi cộng kênh.

Các dữ liệu thực nghiệm này được trình bày trực quan và sáng tỏ thông qua hệ thống bảng giá trị số thời khoảng, các lát cắt địa chấn so sánh trước - sau hiệu chỉnh, cùng các biểu đồ phổ tốc độ dạng đường đẳng trị (contour map) và đồ thị năng lượng cực đại. Khi so sánh với các nghiên cứu kinh điển trong địa vật lý dầu khí quốc tế và các khảo sát thực địa tại bể trầm tích Sông Hồng hay Cửu Long giai đoạn 2015-2019, các kết quả tính toán trong luận văn thể hiện sự tương thích hoàn toàn về mặt quy luật động học trường sóng, khẳng định tính chính xác của thuật toán lập trình.

Đề xuất và khuyến nghị

Dựa trên các kết quả tính toán đạt được, nghiên cứu đưa ra 4 khuyến nghị hành động cụ thể nhằm nâng cao chất lượng xử lý tài liệu địa chấn:

  1. Chuẩn hóa quy trình phân tích phổ tốc độ tự động: Các kỹ sư xử lý số liệu địa vật lý cần thiết lập dải quét vận tốc với bước quét chi tiết không quá 10 m/s, ứng dụng thuật toán quét phổ năng lượng Semblance đa kênh nhằm kiểm soát sai số vận tốc dưới mức 1,0%, áp dụng ngay trong 3 tháng đầu của giai đoạn tiền xử lý dự án.

  2. Kết hợp đo vi địa chấn giếng khoan để tối ưu hiệu chỉnh tĩnh: Bộ phận khảo sát thực địa cần duy trì công tác đo đạc thời gian thẳng đứng dọc thành giếng khoan với khoảng cách điểm nổ từ 1 mét đến 2 mét ở độ sâu 20 đến 50 mét, nhằm cung cấp thông số đới tốc độ nhỏ chính xác, giảm thiểu sai số thời gian thẳng đứng xuống dưới 1 mili-giây trước khi đưa vào tính toán hiệu chỉnh động.

  3. Tự động hóa thuật toán tinh chỉnh hiệu chỉnh động bằng tương quan lặp vòng: Nhóm nghiên cứu và phát triển phần mềm cần tích hợp các mô-đun tinh chỉnh lặp đa chu kỳ trên các cửa sổ thời gian từ 0,5 giây đến 1,0 giây trên hệ thống máy tính hiệu năng cao (HPC), hướng tới mục tiêu rút ngắn 40% thời gian xử lý thủ công trong vòng 6 tháng.

  4. Ứng dụng kỹ thuật di chuyển trước cộng và DMO cho cấu trúc phức tạp: Các đơn vị thăm dò dầu khí cần triển khai kỹ thuật hiệu chỉnh độ dốc (Dip Moveout - DMO) kết hợp di chuyển trước cộng (Pre-stack Migration) đối với các vỉa chứa có góc nghiêng lớn hơn 15 độ, đặt mục tiêu nâng cao độ phân giải lát cắt địa chấn thêm 25% đến 35% trong lộ trình 12 tháng.

Đối tượng nên tham khảo luận văn

Luận văn là tài liệu tham khảo chuyên môn giá trị cho 4 nhóm đối tượng chính:

  1. Học viên cao học và nghiên cứu sinh chuyên ngành Vật lý địa cầu, Kỹ thuật địa vật lý: Nắm vững hệ thống công thức giải tích tường minh về biểu đồ thời khoảng, các phép khai triển chuỗi và giải thuật phân tích phổ tốc độ phục vụ nghiên cứu học thuật chuyên sâu.

  2. Kỹ sư xử lý dữ liệu địa chấn tại các tập đoàn dầu khí và công ty dịch vụ địa vật lý: Vận dụng quy trình lựa chọn vận tốc và kỹ thuật nắn thẳng trục đồng pha để nâng cao chất lượng lát cắt thời gian 2D/3D trong các dự án tìm kiếm thăm dò thực tế.

  3. Giảng viên và nhà khoa học tại các trường đại học khối Khoa học Tự nhiên và Mỏ - Địa chất: Sử dụng tài liệu làm giáo trình tham khảo chuyên đề về xử lý số liệu địa vật lý, phương pháp điểm sâu chung và công nghệ phân tích trường sóng đàn hồi.

  4. Chuyên viên địa chất công trình và khảo sát địa kỹ thuật hạ tầng ngầm: Ứng dụng nguyên lý hiệu chỉnh tĩnh - động trong các bài toán thăm dò địa chấn nông tầng từ 50 mét đến 500 mét để xác định chính xác vị trí đứt gãy, hang karst và cấu trúc nền móng công sinh.

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

  1. Hiệu chỉnh động khác biệt như thế nào so với hiệu chỉnh tĩnh trong xử lý số liệu địa chấn?
    Hiệu chỉnh tĩnh loại trừ ảnh hưởng của địa hình bề mặt và đới tốc độ nhỏ bằng cách dịch chuyển thời gian một lượng cố định cho mọi thời điểm ghi sóng. Ngược lại, hiệu chỉnh động loại trừ ảnh hưởng của khoảng cách thu nổ, trong đó lượng dịch chuyển thời gian $\Delta t$ thay đổi liên tục theo hàm phi tuyến giảm dần theo thời gian $t_0$ từ 0 đến 5 giây.

  2. Tại sao việc chọn sai vận tốc hiệu chỉnh động dù chỉ 4% đến 5% lại làm suy giảm chất lượng lát cắt địa chấn?
    Khi sai số vận tốc đạt mức 30 m/s trên nền vận tốc 700 m/s (chênh lệch 4,28%), các trục đồng pha bị cong vênh. Quá trình cộng sóng điểm sâu chung sẽ gây ra hiện tượng giao thoa lệch pha, làm suy giảm từ 32% đến 46% biên độ tín hiệu có ích và làm mờ các ranh giới địa tầng.

  3. Góc nghiêng của mặt ranh giới địa chất tác động như thế nào đến việc tính toán hiệu chỉnh động?
    Mặt ranh giới nghiêng làm biến dạng hình học biểu đồ thời khoảng và khiến điểm phản xạ trượt một khoảng $MA$ lên tới 0,15 lần khoảng cách thu nổ cực đại. Khi đó, vận tốc hiệu chỉnh phải được quy đổi thành vận tốc điểm sâu chung $v_{\text{ĐSC}} = v / \cos\varphi$ để đảm bảo nắn thẳng chính xác trục đồng pha.

  4. Nguyên lý hoạt động của phương pháp quét tốc độ không đổi (CVS) là gì?
    Phương pháp chia tuyến địa chấn thành các cửa sổ thời gian từ 0,5 đến 1,0 giây, sau đó áp dụng liên tiếp các giá trị vận tốc giả định với bước quét 10 m/s đến 20 m/s. Khi vận tốc thử nghiệm trùng khớp với vận tốc hiệu dụng thực tế, năng lượng cộng sóng trong cửa sổ sẽ đạt giá trị cực đại trên biểu đồ phổ.

  5. Làm thế nào để loại bỏ triệt để sai số thời gian do đới tốc độ nhỏ ở tầng nông gây ra?
    Kỹ sư kết hợp phương pháp vi địa chấn giếng khoan với các điểm nổ cách nhau 1 đến 2 mét để đo thời gian thẳng đứng, sau đó tính toán hiệu chỉnh tĩnh sơ bộ cho từng điểm thu và phát. Cuối cùng, thuật toán tương quan liên kết sẽ tự động tinh chỉnh các sai số dư còn lại với độ chính xác dưới 1 mili-giây.

Kết luận

  • Luận văn đã hệ thống hóa toàn diện cơ sở toán học và bản chất vật lý của phương pháp địa chấn phản xạ điểm sâu chung.
  • Làm sáng tỏ quy luật biến thiên của số gia pháp tuyến và lượng hiệu chỉnh động theo thời gian, cự ly thu nổ và góc nghiêng địa tầng.
  • Kiểm chứng thành công thuật toán quét phổ tốc độ không đổi trên 3 mô hình địa chất phức tạp, nâng biên độ cộng sóng đạt cực đại 100%.
  • Định lượng chính xác mức độ suy giảm biên độ từ 32% đến 46% khi xuất hiện sai số vận tốc, khẳng định tầm quan trọng của việc tối ưu hóa quy trình phân tích vận tốc.
  • Đề xuất chuỗi xử lý kết hợp đồng bộ giữa hiệu chỉnh tĩnh sơ bộ, tinh chỉnh tương quan và hiệu chỉnh động, giúp tăng tỷ số tín hiệu trên nhiễu lên 2,5 đến 3,8 lần.

Trong lộ trình 12 đến 24 tháng tới, các nghiên cứu tiếp theo cần mở rộng thuật toán hiệu chỉnh động cho môi trường dị hướng 3D/4D và tích hợp trí tuệ nhân tạo để tự động hóa hoàn toàn khâu bắt phổ vận tốc. Quý độc giả, học viên và các nhà nghiên cứu quan tâm có thể tra cứu toàn văn luận văn tại Trung tâm Thông tin - Thư viện Đại học Quốc gia Hà Nội để ứng dụng các giải pháp tính toán này vào các đề tài và dự án thăm dò thực địa.