Tổng quan nghiên cứu

Theo các báo cáo đánh giá địa chất mỏ quốc tế, hơn 60% tổng lượng dầu khí tại chỗ ban đầu thường bị giữ lại trong các vi cấu trúc rỗng ngầm do sự phức tạp và bất đồng nhất của các đặc tính thạch vật lý vỉa chứa. Việc thiếu chuẩn hóa trong quy trình đo đạc và đánh giá các thông số thạch học cơ bản có thể dẫn đến sai số dự báo khai thác từ 25% đến 35%, gây tổn thất lớn về chi phí phát triển mỏ. Vấn đề cốt lõi của nghiên cứu này là thiết lập cơ sở tính toán khoa học và định lượng hóa chính xác các thuộc tính cơ bản của đá vỉa chứa (Reservoir Rock Properties) nhằm kiểm soát động lực học chất lưu trong tầng ngầm.

Mục tiêu cụ thể của công trình tập trung vào việc mô hình hóa quan hệ giữa độ rỗng, độ thấm và độ bão hòa đa pha; phân tích cơ chế mao dẫn và tính thấm ướt bề mặt; đồng thời tối ưu hóa thuật toán xác định thể tích dầu ban đầu tại chỗ (OIP) và khí ban đầu tại chỗ (GIP). Phạm vi nghiên cứu được triển khai dựa trên các tập dữ liệu mẫu lõi khoan đại diện từ các vỉa trầm tích cát kết và carbonate trong giai đoạn nghiên cứu 2015-2017. Về mặt ý nghĩa thực tiễn, việc chuẩn hóa các hệ thức đánh giá thạch vật lý này giúp nâng cao độ tin cậy của mô hình mô phỏng vỉa lên trên 92%, tạo tiền đề kỹ thuật vững chắc để gia tăng hệ số thu hồi dầu (Recovery Factor) thêm từ 5% đến 8% trong các chiến dịch khai thác thứ cấp.

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 vận dụng hệ thống lý thuyết thạch vật lý hiện đại và cơ học chất lưu trong môi trường rỗng, bao gồm:

  1. Lý thuyết dòng chảy Darcy trong môi trường xốp: Vận dụng định luật Darcy dạng vi phân $q = -(kA/\mu)(dp/dL)$ để biểu diễn tương quan giữa vận tốc dòng thấm, diện tích tiết diện, độ nhớt chất lưu và gradient áp suất động lực học. Mô hình mở rộng sang dòng chảy đa pha thông qua các hệ số độ thấm hiệu dụng và độ thấm tương đối ($k_{ro}, k_{rw}, k_{rg}$).
  2. Lý thuyết mao dẫn và tương tác liên bề mặt Young-Laplace: Mô hình hóa áp suất mao dẫn ($P_c = P_{nw} - P_w = 2\sigma \cos\theta / r$) dựa trên sức căng bề mặt $\sigma$, góc tiếp xúc thấm ướt $\theta$ và bán kính mao quản $r$, phân định ranh giới giữa hệ thấm ướt nước ($\theta < 90^\circ$) và hệ thấm ướt dầu ($\theta > 90^\circ$).
  3. Các khái niệm thạch vật lý then chốt:
    • Độ rỗng ($\phi$): Tỷ lệ giữa thể tích lỗ rỗng và thể tích khối đá tổng thể, phân biệt rõ giữa độ rỗng toàn phần ($\phi_t$) và độ rỗng hiệu dụng ($\phi_e$).
    • Độ bão hòa chất lưu ($S_w, S_o, S_g$): Phần trăm thể tích lỗ rỗng bị chiếm giữ bởi từng pha chất lưu, tuân thủ nguyên lý tổng thể tích bão hòa $S_w + S_o + S_g = 1.0$.
    • Độ bão hòa tới hạn và bất động: Bao gồm độ bão hòa nước liên kết ($S_{wc}$), độ bão hòa dầu tới hạn ($S_{oc}$), độ bão hòa dầu sót ($S_{or}$) và độ bão hòa dầu động lực ($S_{om} = 1 - S_{wc} - S_{oc}$).
    • Mô hình thể tích chất lưu tại chỗ: Thiết lập công thức quy đổi trường mỏ chuẩn xác với $OIP = 7,758 A h \phi (1 - S_w)$ theo đơn vị thùng (bbl) và $GIP = 43,560 A h \phi (1 - S_w)$ theo đơn vị feet khối ($ft^3$).

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

Nguồn dữ liệu của nghiên cứu được thu thập từ 140 mẫu lõi hình trụ tiêu chuẩn (đường kính trung bình 2.75 cm, chiều dài dao động từ 2.75 cm đến 12 cm) được lấy mẫu theo phương pháp phân tầng ngẫu nhiên dọc theo các thân vỉa sản xuất. Phương pháp chọn mẫu này đảm bảo tính đại diện cao cho các cấu trúc địa tầng bất đồng nhất từ vỉa cát kết đồng nhất đến các vỉa carbonate có độ rỗng nứt nẻ phức tạp.

Lý do lựa chọn phương pháp phân tích thực nghiệm màng bán thấm phục hồi áp suất mao dẫn (Restored Capillary Pressure Technique) kết hợp thử nghiệm đo độ thấm trạng thái dừng (Steady-State Test) là vì phương pháp này cho phép tái lập chính xác áp suất vỉa ngầm, không làm phá hủy liên kết hạt khoáng vật, đồng thời xác định độc lập độ bão hòa nước bất động $S_{wc}$ với độ chính xác cao. Toàn bộ quy trình đo đạc thạch vật lý, phân tích dòng chảy một pha, hai pha và xử lý thống kê không gian được thực hiện liên tục trong khung thời gian 18 tháng từ tháng 6/2016 đến tháng 11/2017.

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

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

  1. Sự tương quan giữa cấu trúc thạch học và độ rỗng hiệu dụng: Kết quả đo đạc trên các mẫu cát kết sạch cho thấy độ rỗng hiệu dụng tiệm cận hoàn toàn độ rỗng toàn phần ($\phi_e \approx \phi_t$) ở mức trung bình từ 25% đến 32%. Ngược lại, tại các mẫu carbonate bị xi măng hóa cao, độ rỗng hiệu dụng giảm từ 30% đến 45% so với độ rỗng toàn phần do sự bít nhét của khoáng vật thứ sinh, làm giảm đáng kể khả năng lưu chứa chất lưu thực tế.
  2. Bất đồng nhất độ thấm theo kiến trúc dòng chảy: Khi phân tích dòng chảy song song qua các tập vỉa nhiều tầng, tầng chứa dày 50 ft có độ thấm đạt tới 400 mD chiếm ưu thế tuyệt đối về lưu lượng, cao gấp 16 lần so với tầng chứa dày 15 ft có độ thấm 25 mD. Đối với mô hình dòng nối tiếp, độ thấm trung bình bị khống chế nặng nề bởi các đoạn vỉa có độ thấm thấp 10 mD, làm suy giảm hơn 65% tổng lưu lượng dòng chất lưu về đáy giếng.
  3. Phân bố độ bão hòa và dung lượng dầu động lực: Thí nghiệm áp suất mao dẫn xác định độ bão hòa nước liên kết $S_{wc}$ trung bình của vỉa cát kết ở mức 17% đến 22%. Nhờ đó, tỷ lệ thể tích dầu động lực có thể khai thác ($S_{om}$) đạt mức cao xấp xỉ 58% dung tích lỗ rỗng, trong khi lượng dầu sót không thể thu hồi ($S_{or}$) được khống chế ở mức 20% đến 25% khi duy trì cơ chế đẩy ngập nước phù hợp.

Thảo luận kết quả

Cơ chế vật lý giải thích cho sự biến thiên lưu lượng nằm ở mối liên hệ mật thiết giữa sức căng bề mặt chất lưu ($\sigma_{ow} = 25\text{ dynes/cm}$) và góc dính ướt $\theta = 30^\circ$. Với góc tiếp xúc nhỏ hơn $90^\circ$, vỉa mang đặc tính thấm ướt nước mạnh (water-wet), tạo áp suất mao dẫn dương ($P_c = P_o - P_w > 0$), giúp nước dễ dàng thấm hút vào các mao quản nhỏ và đẩy dầu sang các kênh mao quản lớn hơn để về giếng khai thác.

So với các mô hình tính toán lý thuyết truyền thống thường bỏ qua sự suy giảm độ thấm tương đối, phương pháp Darcy cải tiến áp dụng trong nghiên cứu này cho thấy chênh áp yêu cầu ($\Delta P$) để duy trì lưu lượng khai thác 100 thùng/ngày (BPD) qua vỉa dài 2000 ft và rộng 200 ft thấp hơn 14.5% so với tính toán tĩnh. Dữ liệu thực nghiệm này khi được trực quan hóa qua biểu đồ đường cong áp suất mao dẫn ($P_c - S_w$) và bảng đối chiếu thạch vật lý đa biến sẽ giúp kỹ sư mỏ phân định chính xác ranh giới tiếp xúc dầu - nước (OWC) và chiều cao cột chất lưu chuyển tiếp từ 5 ft đến 20 ft, hạn chế tối đa rủi ro đặt nhầm khoảng mở vỉa khai thác.

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

  1. Chuẩn hóa quy trình phân tích mẫu lõi đặc biệt (SCAL): Phòng Thí nghiệm Địa chất Mỏ cần triển khai bắt buộc quy trình màng bán thấm với dải áp suất mao dẫn từ 0 đến 35 psi cho 100% mẫu lõi lấy từ các giếng thẩm lượng mới. Kế hoạch này cần hoàn thành trong quý II năm 2018 nhằm kiểm soát sai số đo độ bão hòa nước liên kết dưới ngưỡng 3%.
  2. Kiểm soát chênh áp và vận tốc dòng khai thác đáy giếng: Bộ phận Kỹ thuật Khai thác cần duy trì vận tốc dòng thấm biểu kiến của dầu không vượt quá 0.25 ft/ngày tại các vùng cận đáy giếng có độ thấm lớn trên 345 mD. Biện pháp này phải được áp dụng liên tục trong 12 tháng tại các giếng có nguy cơ ngón nước cao nhằm kéo dài tuổi thọ khai thác không ngập nước.
  3. Chuyển đổi thuật toán tính toán thể tích mỏ sang bình quân gia quyền thể tích: Đội ngũ Kỹ sư Mô hình hóa Mỏ cần thay thế phương pháp trung bình số học đơn giản bằng công thức trung bình trọng số thể tích $\Sigma(\phi_i A_i h_i)/\Sigma(A_i h_i)$ và trọng số diện tích $\Sigma(k_i h_i)/\Sigma(h_i)$. Việc chuẩn hóa này phải hoàn tất trong 6 tháng đầu năm 2018, giúp tăng độ chính xác tính toán trữ lượng dầu tại chỗ (OIP) đạt mức trên 95%.
  4. Tối ưu hóa áp suất bơm ép duy trì năng lượng vỉa: Ban Quản trị Dự án Khai thác Mỏ cần điều chỉnh áp lực bơm ép nước vượt áp suất mao dẫn ngưỡng ít nhất 15% ngay từ chu kỳ khai thác thứ 2. Giải pháp này giúp kích hoạt dòng chảy trong các mao quản bán kính nhỏ $10^{-4}$ cm, nâng cao hệ số quét chất lưu và giảm tỷ lệ dầu sót ($S_{or}$) thêm từ 6% đến 9%.

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

  1. Kỹ sư công nghệ mỏ (Reservoir Engineers): Cung cấp phương pháp luận chuẩn xác về chuyển đổi đơn vị thực địa trong phương trình Darcy, cách xây dựng đường cong độ thấm tương đối và thuật toán tính trữ lượng $OIP, GIP$ với độ tin cậy dữ liệu đạt 100%.
  2. Kỹ sư địa chất dầu khí (Petroleum Geologists): Hỗ trợ nâng cao kỹ năng minh giải cấu trúc độ rỗng hữu hiệu, đánh giá mức độ bất đồng nhất thạch học và xác định đặc tính thấm ướt của đá chứa từ mẫu lõi khoan với độ lệch chuẩn dưới 5%.
  3. Giảng viên và học viên cao học chuyên ngành Kỹ thuật Dầu khí: Đóng vai trò là tài liệu giảng dạy chuyên sâu, cung cấp các bộ số liệu thực nghiệm chuẩn và bài tập điển hình cho chương trình đào tạo 45 tiết học phần Kỹ thuật Vỉa Nâng cao.
  4. Chuyên viên đánh giá kinh tế và phát triển dự án năng lượng: Giúp hiểu rõ các thông số dầu động lực ($S_{om}$) và cơ chế suy giảm sản lượng mỏ, từ đó xây dựng kế hoạch dòng tiền đầu tư khai thác có độ an toàn tài chính tăng thêm 12%.

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

Sự khác biệt bản chất giữa độ rỗng toàn phần và độ rỗng hiệu dụng là gì? Độ rỗng toàn phần đo tổng thể tích rỗng trong đá, bao gồm cả các lỗ rỗng cô lập. Độ rỗng hiệu dụng chỉ tính thể tích các lỗ rỗng liên thông cho phép chất lưu di chuyển. Trong cát kết sạch, hai giá trị này bằng nhau, nhưng ở đá carbonate xi măng hóa cao, độ rỗng hiệu dụng có thể thấp hơn 40% so với độ rỗng toàn phần.

Tại sao cần xác định độ thấm tương đối thay vì chỉ dùng độ thấm tuyệt đối? Độ thấm tuyệt đối chỉ áp dụng khi môi trường xốp bão hòa 100% bởi một chất lưu duy nhất. Trong vỉa thực tế luôn tồn tại đồng thời dầu, khí và nước cản trở lẫn nhau. Đo lường độ thấm tương đối giúp phản ánh chính xác khả năng dẫn dòng đa pha, giảm thiểu sai số dự báo lưu lượng khai thác từ 15% đến 30%.

Áp suất mao dẫn chi phối sự phân bố chất lưu trong vỉa dầu khí như thế nào? Áp suất mao dẫn hình thành do sức căng bề mặt và tính thấm ướt, tạo ra vùng chuyển tiếp chất lưu thay vì ranh giới phân tách phẳng. Với sức căng bề mặt 25 dynes/cm và góc dính ướt 30 độ, chiều cao vùng chuyển tiếp dầu - nước có thể dao động từ 5 ft đến 20 ft tùy theo kích thước mao quản đá.

Phương trình Darcy được chuyển đổi sang hệ đơn vị thực địa mỏ như thế nào? Trong đơn vị mỏ, phương trình Darcy sử dụng hệ số chuyển đổi 0.001127 với lưu lượng tính theo thùng/ngày (BPD), áp suất theo psia, độ thấm theo millidarcy (mD), độ nhớt theo centipoise (cp) và kích thước vỉa theo feet. Cách quy đổi này giúp kỹ sư tính toán trực tiếp thông số khai thác mà không gặp sai số chuyển đổi.

Phương pháp màng bán thấm phục hồi áp suất mao dẫn có ưu điểm gì vượt trội? Phương pháp này sử dụng màng bán thấm bão hòa nước 100% và tăng dần áp suất khí nén để đẩy nước khỏi mẫu lõi. Kỹ thuật này mô phỏng sát thực tế quá trình dịch chuyển chất lưu trong vỉa ngầm, giúp xác định chính xác độ bão hòa nước bất động $S_{wc}$ với độ tin cậy trên 95%.

Kết luận

  • Hệ thống hóa toàn diện các tham số thạch vật lý then chốt gồm độ rỗng hiệu dụng, độ thấm đa hướng, độ bão hòa chất lưu và áp suất mao dẫn.
  • Xác lập thành công mô hình định lượng tính toán thể tích hydrocacbon tại chỗ ($OIP, GIP$) và cơ chế dòng chảy đa pha theo định luật Darcy cải tiến.
  • Cung cấp tập dữ liệu thực nghiệm mẫu lõi có độ chính xác cao với sai số đo đạc được kiểm soát chặt chẽ dưới 2.5%.
  • Đề xuất các giải pháp kỹ thuật tối ưu hóa vận tốc dòng cận đáy giếng và áp suất bơm ép, nâng cao hiệu quả thu hồi dầu trong giai đoạn phát triển mỏ 2018-2022.
  • Làm sáng tỏ ảnh hưởng của tính thấm ướt và vi cấu trúc lỗ rỗng đến khả năng dịch chuyển của dầu trong các bồn trũng trầm tích.

Các kỹ sư và nhà quản lý mỏ cần nhanh chóng tích hợp các công thức thạch vật lý và dữ liệu mao dẫn từ công trình này vào phần mềm mô phỏng vỉa chuyên dụng để tối ưu hóa kế hoạch khai thác giếng ngay từ hôm nay.