Lưu lượng ép và áp suất miệng giếng: tách tổn thất tuyến và vỉa

1. Ba vị trí đo không có cùng áp suất

Một xu hướng áp suất tại trạm tăng lên có thể đến từ đường ống bẩn, van đổi vị trí, nước thay đổi mật độ hoặc khả năng nhận ép của vỉa giảm. Kết luận “vỉa bị nghẽn” chỉ từ đồng hồ trạm là vượt quá thông tin quan sát. Tài liệu Injection inflow của Harmony xây dựng quan hệ tuyến tính từ áp suất vỉa và lưu lượng thử tại áp suất mặt cát, với chỉ số nhận ép không đổi [1]. Từ cơ sở đó, bài này tự lập một mạng thủy lực đơn giản để thấy vì sao vị trí đo quyết định ý nghĩa của hệ số khớp.

Ký hiệu pᵣ là áp suất vỉa quy về độ sâu tiếp nhận, pᵦ là áp suất nước tại mặt cát, pₕ là áp suất miệng giếng và pₛ là áp suất đầu tuyến ở trạm. Tất cả dùng MPa, cùng quy chiếu áp suất. Lưu lượng q dùng m³/ngày trên cùng cơ sở thể tích; nước giả định không nén được. Chỉ số J có đơn vị (m³/ngày)/MPa, và q = J(pᵦ−pᵣ) khi dòng ép đi vào vỉa trong miền tuyến tính đã chọn.

Đây là quan hệ tại một trạng thái vỉa xác định. Nếu áp suất vỉa thay đổi theo thời gian, cùng q có thể cần pᵦ khác dù J chưa đổi. Nếu nước và đá thay đổi tính chất hoặc xuất hiện cơ chế tiếp nhận khác, việc giữ J hằng cũng cần xem lại. Những biến đó phải được tách khỏi đường đặc tính của ống trước khi gán nguyên nhân.

2. Viết cân bằng áp suất theo chiều dòng đi xuống

Với độ chênh cao thẳng đứng H tính từ miệng đến mặt cát, mật độ ρ và gia tốc trọng trường g, áp suất thủy tĩnh tăng xuống dưới bằng ρgH. Tổn thất ma sát lại làm giảm áp suất dọc chiều chảy. Vì thế pᵦ = pₕ + ρgH/10⁶ − Δpₜ, trong đó Δpₜ là tổn thất ống ép tính bằng MPa. Dấu cộng thủy tĩnh và dấu trừ ma sát cũng được nêu trong phần Hall Plot của IHS [3]. H ở đây là chênh cao đứng, còn chiều dài dùng để tính ma sát là chiều dài đường chảy; giếng nghiêng không cho phép thay hai đại lượng này bằng nhau.

Nếu đầu tuyến trạm và miệng giếng cùng cao độ, pₛ = pₕ + Δpₗ. Giả sử Δpₜ = Kₜq² và Δpₗ = Kₗq² trên một khoảng lưu lượng hẹp, ta được pₛ = pᵣ − ρgH/10⁶ + q/J + (Kₜ+Kₗ)q². Các hệ số K có đơn vị MPa/(m³/ngày)², phụ thuộc đúng đơn vị q đang dùng. Đổi q sang m³/s mà giữ nguyên số K sẽ làm sai kết quả rất lớn.

Dạng bình phương là mô hình cục bộ giả định, không khẳng định hệ số ma sát luôn không đổi. Với dòng tầng, thay đổi độ nhớt lớn, van điều tiết hoặc đường kính biến đổi, cần xây dựng lại quan hệ tổn thất. Nếu trạm thấp hơn miệng giếng, phương trình pₛ còn phải cộng phần cột áp để đưa nước lên cao. Một sơ đồ ghi rõ cao độ và chiều dòng giúp kiểm dấu trước khi tính.

3. Ví dụ giả định từ mặt cát về trạm

Chọn pᵣ = 26 MPa, J = 500 (m³/ngày)/MPa, H = 2.400 m, ρ = 1.020 kg/m³ và g = 9,80665 m/s². Tại q = 2.000 m³/ngày, giả sử tổn thất ống ép bằng 1,10 MPa và tuyến mặt đất bằng 0,70 MPa. Phần thủy tĩnh là 24,006679 MPa, còn chênh áp để nước đi vào vỉa là q/J = 4 MPa.

Do đó pᵦ = 30 MPa, pₕ = 30−24,006679+1,10 = 7,093321 MPa và pₛ = 7,793321 MPa. Áp suất trạm thấp hơn nhiều so với áp suất vỉa vì cột nước hỗ trợ ép xuống sâu. Kết quả không có nghĩa dòng nước tự có năng lượng miễn phí: nó chuyển thế năng trọng trường và năng lượng áp suất trong ranh giới đã chọn. Muốn tính công suất bơm phải biết cả áp suất hút và trạng thái dòng vào trạm.

Tại q = 0, mô hình cân bằng cho pₕ = pᵣ−ρgH/10⁶ = 1,993321 MPa nếu cột nước liên tục và đã cân bằng với vỉa. Áp suất đọc ngay sau đóng van chưa chắc bằng giá trị này vì còn quá trình phục hồi, đàn hồi của hệ và phân bố nhiệt. Không nên dùng một số tức thời làm áp suất vỉa tĩnh mà bỏ lịch sử thời gian.

4. Tìm điểm làm việc dưới giới hạn áp suất giả định

Từ ví dụ, K = Kₜ+Kₗ = 1,80/2.000² = 4,5×10⁻⁷ MPa/(m³/ngày)². Nếu chỉ xét một giới hạn nghiên cứu pₛ = 8 MPa, nghiệm dương của Kq² + q/J + pᵣ−ρgH/10⁶−pₛ = 0 là q = 2.054,043392 m³/ngày. Áp suất mặt cát tương ứng bằng 30,108087 MPa. Đây là giao điểm của mô hình, chưa phải một lưu lượng vận hành được phép.

Giả sử riêng bài toán có ràng buộc pᵦ ≤ 31 MPa, thì nhánh vỉa cho q ≤ 2.500 m³/ngày. Để đạt 2.500 với đường ống trên cần pₛ = 9,805821 MPa, vượt giới hạn trạm giả định 8 MPa. Vì vậy trong hai điều kiện đã nêu, trạm giới hạn lưu lượng trước. Giá trị 31 MPa được đặt để minh họa cách kết hợp ràng buộc; nó không phải áp suất nứt vỉa hay định mức thiết bị lấy từ tiêu chuẩn.

SLB mô tả thử nhận ép là xác định quan hệ lưu lượng–áp suất trong phạm vi không gây nứt vỉa [2]. Trong một hồ sơ thực, giới hạn địa cơ học, áp suất làm việc cho phép của từng đoạn và đường đặc tính bơm phải đến từ dữ liệu riêng. Không thể lấy một giao điểm đại số rồi gọi nó là lưu lượng an toàn cho mọi giếng.

5. Tách nghẽn ống khỏi thay đổi khả năng nhận ép

Tại q = 1.000 m³/ngày, tổn thất ống ép theo mô hình bình phương chỉ còn 0,275 MPa và pₕ = 4,268321 MPa. Dùng hai điểm miệng giếng 1.000 và 2.000 để tính Δq/Δpₕ cho 353,982301 (m³/ngày)/MPa, thấp hơn J thật là 500. Chênh lệch không chứng minh vỉa kém hơn; hệ số bề mặt đã gộp ma sát ống.

Xét hai thay đổi độc lập tại cùng q = 2.000. Nếu J giảm từ 500 xuống 400 mà đường ống giữ nguyên, pᵦ tăng từ 30 lên 31 MPa và pₛ tăng lên 8,793321 MPa. Nếu J vẫn bằng 500 nhưng tổn thất ống ép tăng gấp đôi từ 1,10 lên 2,20 MPa, pᵦ vẫn bằng 30 MPa còn pₛ tăng lên 8,893321 MPa. Hai số trạm khá gần nhau, trong khi cơ chế khác hẳn. Đo áp suất đáy hoặc kiểm chứng mô hình ống bằng dữ liệu độc lập mới giúp phân biệt.

Một đường Hall thay đổi độ dốc cũng cần đọc cùng lịch sử áp suất vỉa, nước ép và điều kiện dòng, như phần giả thiết của phương pháp đã lưu ý [3]. Không nên coi một dấu hiệu hình học trên đồ thị là phép xác định nguyên nhân duy nhất. Chất lượng dữ liệu ngày ép, lưu lượng và vị trí áp suất quan trọng không kém việc vẽ đường hồi quy.

6. Độ nhạy và giới hạn suy rộng

Đạo hàm dpₛ/dq = 1/J + 2Kq bằng 0,0038 MPa/(m³/ngày) tại q = 2.000. Nghịch đảo cho dq/dpₛ = 263,157895 (m³/ngày)/MPa. Một thay đổi nhỏ áp suất trạm sẽ không tạo mức tăng lưu lượng bằng J×Δp, vì một phần chênh áp tăng thêm bị tiêu hao trong tuyến. Giữ pₛ = 8 MPa, thay J bằng 400, 500 và 600 cho q lần lượt 1.811,800131; 2.054,043392; 2.244,186634 m³/ngày. Lợi ích tăng J giảm dần khi tổn thất tuyến chiếm tỷ trọng lớn.

Sai khác mật độ 20 kg/m³ trên cột 2.400 m tương ứng 0,470719 MPa thủy tĩnh. Nếu vẫn dùng mật độ cũ để suy pᵦ từ pₕ, sai khác đó có thể bị gán nhầm cho vỉa. Mật độ thực có thể thay đổi theo nhiệt độ, áp suất và độ mặn; mô hình chính xác hơn cần tích phân ρ(z)g dz và tính ma sát theo điều kiện cục bộ.

Bảng tính cuối cùng nên giữ riêng áp suất tại bốn vị trí, độ sâu đứng, chiều dài ống, cơ sở thể tích, nhiệt độ nước và khoảng dữ liệu dùng khớp J. Mô hình ở đây bỏ qua dòng quá độ, đàn hồi, nhiều lớp tiếp nhận và biến đổi độ mở khe nứt. Nó hữu ích để kiểm ranh giới và loại những kết luận sai từ đồng hồ bề mặt; để đánh giá giếng thực cần bổ sung dữ liệu phản ánh các cơ chế đang hoạt động.

Nguồn tham khảo
[1] Injection inflow — Harmony Enterprise 2026.1 — S&P Global. Truy cập ngày 09/11/2026.
[2] injectivity test — Energy Glossary — SLB. Truy cập ngày 09/11/2026.
[3] Surveillance Analysis Theory — Hall Plot Analysis — IHS Markit. Truy cập ngày 09/11/2026.

Rà soát nội dung: 09/11/2026.
 
Top