1. Đường IPR mô tả một trạng thái của giếng và vỉa
IPR liên hệ lưu lượng đi vào giếng với áp suất đáy đang chảy. Khi dầu và khí cùng chuyển động trong vùng vỉa gần giếng, một đường thẳng với chỉ số năng suất cố định có thể không mô tả được toàn dải áp suất. Whitson trình bày dạng Vogel chuẩn hóa cho dầu bão hòa, cùng dạng nối với đường thẳng khi áp suất vỉa cao hơn điểm bọt.[1]
Bài này xem Vogel như mô hình quan hệ lưu lượng tại một trạng thái đã xác định. Các hệ số không tạo thành một mô hình lịch sử vỉa đầy đủ. Lưu lượng giả định là dầu bồn, không phải tổng lỏng gồm nước, và mọi áp suất trong tỷ số đều là áp suất tuyệt đối ở cùng mốc chiều sâu.
2. Chuẩn hóa trước khi thay số
Gọi p_r là áp suất vỉa đại diện, p_wf là áp suất đáy đang chảy, x=p_wf/p_r và q_max là tham số giới hạn của mô hình khi p_wf=0. Với 0≤x≤1:
q/q_max = f(x) = 1−0,2x−0,8x²
q_max = q_test/f(p_test/p_r)
q và q_max phải có cùng đơn vị, chẳng hạn m³ dầu bồn/ngày. x và f không có đơn vị. q_max được gọi là AOF theo mô hình; điều kiện áp suất đáy tuyệt đối bằng không không phải một chế độ vận hành khả thi được đề xuất.[1]
S&P Global nêu các đầu vào cần có gồm áp suất vỉa, cặp lưu lượng–áp suất đáy từ thử giếng và áp suất bọt.[2] Ghi một lưu lượng mà không biết nó đi cùng áp suất nào sẽ không đủ dựng đường. Cũng không ghép lưu lượng một ngày với áp suất vỉa của thời kỳ đã suy giảm khác.
3. Ví dụ từ một điểm thử giả định
Chọn p_r=24 MPa, p_b=27 MPa, p_test=16 MPa và q_test=400 m³/ngày. Vỉa đang ở dưới điểm bọt trong ví dụ. Tỷ số x_test=2/3 và f_test=0,511111111, nên:
q_max = 400/0,511111111 = 782,608696 m³/ngày
Tại p_wf=12 MPa, x=0,5, f=0,7 và q=547,826087 m³/ngày. Tại 8 MPa, q=660,869565 m³/ngày. Giảm thêm 4 MPa ở vùng áp suất thấp tạo mức tăng lưu lượng khác với giảm 4 MPa ở gần áp suất vỉa; độ cong đã thay đổi khả năng đáp ứng với drawdown.
Muốn tìm áp suất cho một tỷ lệ y=q/q_max, giải phương trình bậc hai và lấy nghiệm trong miền vật lý:
x = [−0,2 + √(0,04+3,2(1−y))]/1,6
Với y=0,5, x=0,675390530. Nghiệm còn lại âm nên không thuộc miền áp suất của bài. Một chương trình tính phải kiểm 0≤y≤1 trước khi dùng căn; không tự trả một áp suất âm khi người dùng nhập lưu lượng vượt giới hạn đường đang xét.
4. Chỉ số qua hai điểm khác độ dốc tiếp tuyến
Từ điểm thử, chỉ số secant J_sec=q/(p_r−p_wf)=50 m³/ngày/MPa. Trong khi đó, độ nhạy lưu lượng theo drawdown ở một điểm là:
J_tan = −dq/dp_wf = (q_max/p_r)(0,2+1,6x)
Tại điểm thử, J_tan=41,304348 m³/ngày/MPa. Ngay gần đóng giếng, x=1, giá trị tiếp tuyến là 58,695652 m³/ngày/MPa. Cả ba số đều được tính đúng nhưng trả lời ba câu hỏi khác nhau.
Nếu dùng J_sec=50 để ngoại suy tới p_wf=12 MPa, đường thẳng cho 600 m³/ngày, cao hơn Vogel khoảng 52,17 m³/ngày. Nếu dùng tiếp tuyến tại đóng giếng cho cả dải, sai khác còn khác nữa. Vì thế, nhãn PI trong dữ liệu cần nói rõ là hệ số của nhánh tuyến tính, tỷ số tại điểm thử hay đạo hàm của đường cong.
5. Một cách nhìn qua hàm độ linh động
Học liệu NTNU chỉ ra rằng dạng Vogel có thể xuất hiện khi hàm M(p)=k_ro/(B_oμ_o) biến thiên tuyến tính theo áp suất và M(0)=M(p_r)/9.[3] Đây là một giả thiết để hiểu hình dạng, không phải quan hệ PVT chung cho mọi dầu.
Đặt M_r=M(p_r), khi đó M(p)=M_r[1+8p/p_r]/9. Nếu q tỷ lệ với tích phân M(p)dp từ p_wf đến p_r, chia tích phân đó cho giá trị khi cận dưới bằng0 thu được:
q/q_max = [(1−x)+4(1−x²)]/5
Biểu thức trở lại 1−0,2x−0,8x². Cách suy này cho thấy đường cong đã gói một giả thiết về độ linh động vào hai hệ số cố định. Nếu độ thấm tương đối, thành phần hoặc cơ chế dòng thay đổi mạnh, giữ nguyên hình dạng chuẩn hóa cần được kiểm bằng dữ liệu.
6. Điểm thử gần đóng giếng làm phép suy ngược nhạy hơn
Giữ q_test=400 và p_r=24, đổi p_test từ16 thành15,8 hoặc16,2 MPa cho q_max lần lượt 766,854830 và799,200799 m³/ngày. Sai lệch áp suất nhỏ đã ảnh hưởng tham số ngoại suy, dù lưu lượng đo không đổi.
Trong một phép thử riêng rất gần đóng giếng, mô hình gốc cho q=57,608696 m³/ngày ở23 MPa. Nếu cùng lưu lượng đó bị gán áp suất22,8 hoặc23,2 MPa, q_max suy ra thành654,644269 hoặc974,583197 m³/ngày. Vì f(x) gần0 khi x gần1, chia cho f khuếch đại sai lệch. Đây là tính chất toán học của phép hiệu chỉnh, không phải bằng chứng vỉa thật đột ngột đổi năng suất.
Sai số p_r cũng cần xét. Giữ điểm thử400 tại16 MPa nhưng dùng p_r=23,5 hoặc24,5 MPa cho q_max=811,386593 hoặc757,293802 m³/ngày. Áp suất vỉa đại diện vì thế không nên được coi là một số chính xác tuyệt đối chỉ vì bảng tính nhập được nhiều chữ số.
7. Kiểm nhiều điểm và nối qua điểm bọt
Với p_r biết trước và các điểm cùng trạng thái, có thể khớp một tham số A=q_max bằng cách tối thiểu tổng(q_i−Af_i)². Nghiệm là A=tổng(f_iq_i)/tổng(f_i²). Ba điểm giả định tại21;16;8 MPa với lưu lượng180;400;650 m³/ngày cho A=776,460037. Phần dư theo thứ tự là15,002242;3,142648;−5,677364 m³/ngày. Mẫu phần dư có ý nghĩa hơn việc chỉ báo một đường đã đi gần các điểm; trọng số cần thay đổi nếu độ tin cậy từng phép đo khác nhau.
Với một vỉa khác có p_r=30 MPa, p_b=20 MPa và J=35 m³/ngày/MPa, dùng nhánh tuyến tính q=J(p_r−p_wf) phía trên p_b. Phía dưới, một dạng ghép liên tục là:
q = J(p_r−p_b)+(Jp_b/1,8)f(p_wf/p_b)
Ở20 MPa, cả hai nhánh cho350 m³/ngày và cùng đạo hàm−J. Ở15;10;0 MPa, nhánh dưới cho505,555556;622,222222;738,888889 m³/ngày. Việc ghép giữ cả giá trị và độ dốc ở điểm nối, nhưng vẫn cần dữ liệu kiểm chứng ngoài điểm nối.
8. Phạm vi sử dụng
Nước có thể cần một nhánh riêng trước khi cộng thành tổng lỏng; không tự dùng lưu lượng dầu cộng nước để hiệu chỉnh đường chỉ dành cho dầu. Áp suất đáy phải cùng mốc và đã được xử lý nhất quán với gradient trong giếng. Khi p_r, nước hoặc trạng thái hoàn thiện thay đổi, cần đánh giá lại đường IPR thay vì chỉ dịch một điểm trên đồ thị cũ.
Một hồ sơ tốt lưu dữ liệu thử, áp suất bọt, giả thiết nhánh, phần dư và các kịch bản sai số. Đường Vogel khi đó hỗ trợ so sánh khả năng cung cấp vào giếng; lưu lượng thực còn cần đồng thời đáp ứng đường nâng và điều kiện hệ thống phía trên.
Nguồn tham khảo
[1] Nodal Analysis - IPR / VLP — Reservoir IPR — Whitson user manual. Truy cập ngày 09/11/2026.
[2] Vogel inflow — S&P Global, Harmony Enterprise2026.1. Truy cập ngày 09/11/2026.
[3] TPG4245 Regular Exam2023 — Problems1 and6 — Prof. Milan Stanko, Norwegian University of Science and Technology. Truy cập ngày 09/11/2026.
Rà soát nội dung: 09/11/2026.
IPR liên hệ lưu lượng đi vào giếng với áp suất đáy đang chảy. Khi dầu và khí cùng chuyển động trong vùng vỉa gần giếng, một đường thẳng với chỉ số năng suất cố định có thể không mô tả được toàn dải áp suất. Whitson trình bày dạng Vogel chuẩn hóa cho dầu bão hòa, cùng dạng nối với đường thẳng khi áp suất vỉa cao hơn điểm bọt.[1]
Bài này xem Vogel như mô hình quan hệ lưu lượng tại một trạng thái đã xác định. Các hệ số không tạo thành một mô hình lịch sử vỉa đầy đủ. Lưu lượng giả định là dầu bồn, không phải tổng lỏng gồm nước, và mọi áp suất trong tỷ số đều là áp suất tuyệt đối ở cùng mốc chiều sâu.
2. Chuẩn hóa trước khi thay số
Gọi p_r là áp suất vỉa đại diện, p_wf là áp suất đáy đang chảy, x=p_wf/p_r và q_max là tham số giới hạn của mô hình khi p_wf=0. Với 0≤x≤1:
q/q_max = f(x) = 1−0,2x−0,8x²
q_max = q_test/f(p_test/p_r)
q và q_max phải có cùng đơn vị, chẳng hạn m³ dầu bồn/ngày. x và f không có đơn vị. q_max được gọi là AOF theo mô hình; điều kiện áp suất đáy tuyệt đối bằng không không phải một chế độ vận hành khả thi được đề xuất.[1]
S&P Global nêu các đầu vào cần có gồm áp suất vỉa, cặp lưu lượng–áp suất đáy từ thử giếng và áp suất bọt.[2] Ghi một lưu lượng mà không biết nó đi cùng áp suất nào sẽ không đủ dựng đường. Cũng không ghép lưu lượng một ngày với áp suất vỉa của thời kỳ đã suy giảm khác.
3. Ví dụ từ một điểm thử giả định
Chọn p_r=24 MPa, p_b=27 MPa, p_test=16 MPa và q_test=400 m³/ngày. Vỉa đang ở dưới điểm bọt trong ví dụ. Tỷ số x_test=2/3 và f_test=0,511111111, nên:
q_max = 400/0,511111111 = 782,608696 m³/ngày
Tại p_wf=12 MPa, x=0,5, f=0,7 và q=547,826087 m³/ngày. Tại 8 MPa, q=660,869565 m³/ngày. Giảm thêm 4 MPa ở vùng áp suất thấp tạo mức tăng lưu lượng khác với giảm 4 MPa ở gần áp suất vỉa; độ cong đã thay đổi khả năng đáp ứng với drawdown.
Muốn tìm áp suất cho một tỷ lệ y=q/q_max, giải phương trình bậc hai và lấy nghiệm trong miền vật lý:
x = [−0,2 + √(0,04+3,2(1−y))]/1,6
Với y=0,5, x=0,675390530. Nghiệm còn lại âm nên không thuộc miền áp suất của bài. Một chương trình tính phải kiểm 0≤y≤1 trước khi dùng căn; không tự trả một áp suất âm khi người dùng nhập lưu lượng vượt giới hạn đường đang xét.
4. Chỉ số qua hai điểm khác độ dốc tiếp tuyến
Từ điểm thử, chỉ số secant J_sec=q/(p_r−p_wf)=50 m³/ngày/MPa. Trong khi đó, độ nhạy lưu lượng theo drawdown ở một điểm là:
J_tan = −dq/dp_wf = (q_max/p_r)(0,2+1,6x)
Tại điểm thử, J_tan=41,304348 m³/ngày/MPa. Ngay gần đóng giếng, x=1, giá trị tiếp tuyến là 58,695652 m³/ngày/MPa. Cả ba số đều được tính đúng nhưng trả lời ba câu hỏi khác nhau.
Nếu dùng J_sec=50 để ngoại suy tới p_wf=12 MPa, đường thẳng cho 600 m³/ngày, cao hơn Vogel khoảng 52,17 m³/ngày. Nếu dùng tiếp tuyến tại đóng giếng cho cả dải, sai khác còn khác nữa. Vì thế, nhãn PI trong dữ liệu cần nói rõ là hệ số của nhánh tuyến tính, tỷ số tại điểm thử hay đạo hàm của đường cong.
5. Một cách nhìn qua hàm độ linh động
Học liệu NTNU chỉ ra rằng dạng Vogel có thể xuất hiện khi hàm M(p)=k_ro/(B_oμ_o) biến thiên tuyến tính theo áp suất và M(0)=M(p_r)/9.[3] Đây là một giả thiết để hiểu hình dạng, không phải quan hệ PVT chung cho mọi dầu.
Đặt M_r=M(p_r), khi đó M(p)=M_r[1+8p/p_r]/9. Nếu q tỷ lệ với tích phân M(p)dp từ p_wf đến p_r, chia tích phân đó cho giá trị khi cận dưới bằng0 thu được:
q/q_max = [(1−x)+4(1−x²)]/5
Biểu thức trở lại 1−0,2x−0,8x². Cách suy này cho thấy đường cong đã gói một giả thiết về độ linh động vào hai hệ số cố định. Nếu độ thấm tương đối, thành phần hoặc cơ chế dòng thay đổi mạnh, giữ nguyên hình dạng chuẩn hóa cần được kiểm bằng dữ liệu.
6. Điểm thử gần đóng giếng làm phép suy ngược nhạy hơn
Giữ q_test=400 và p_r=24, đổi p_test từ16 thành15,8 hoặc16,2 MPa cho q_max lần lượt 766,854830 và799,200799 m³/ngày. Sai lệch áp suất nhỏ đã ảnh hưởng tham số ngoại suy, dù lưu lượng đo không đổi.
Trong một phép thử riêng rất gần đóng giếng, mô hình gốc cho q=57,608696 m³/ngày ở23 MPa. Nếu cùng lưu lượng đó bị gán áp suất22,8 hoặc23,2 MPa, q_max suy ra thành654,644269 hoặc974,583197 m³/ngày. Vì f(x) gần0 khi x gần1, chia cho f khuếch đại sai lệch. Đây là tính chất toán học của phép hiệu chỉnh, không phải bằng chứng vỉa thật đột ngột đổi năng suất.
Sai số p_r cũng cần xét. Giữ điểm thử400 tại16 MPa nhưng dùng p_r=23,5 hoặc24,5 MPa cho q_max=811,386593 hoặc757,293802 m³/ngày. Áp suất vỉa đại diện vì thế không nên được coi là một số chính xác tuyệt đối chỉ vì bảng tính nhập được nhiều chữ số.
7. Kiểm nhiều điểm và nối qua điểm bọt
Với p_r biết trước và các điểm cùng trạng thái, có thể khớp một tham số A=q_max bằng cách tối thiểu tổng(q_i−Af_i)². Nghiệm là A=tổng(f_iq_i)/tổng(f_i²). Ba điểm giả định tại21;16;8 MPa với lưu lượng180;400;650 m³/ngày cho A=776,460037. Phần dư theo thứ tự là15,002242;3,142648;−5,677364 m³/ngày. Mẫu phần dư có ý nghĩa hơn việc chỉ báo một đường đã đi gần các điểm; trọng số cần thay đổi nếu độ tin cậy từng phép đo khác nhau.
Với một vỉa khác có p_r=30 MPa, p_b=20 MPa và J=35 m³/ngày/MPa, dùng nhánh tuyến tính q=J(p_r−p_wf) phía trên p_b. Phía dưới, một dạng ghép liên tục là:
q = J(p_r−p_b)+(Jp_b/1,8)f(p_wf/p_b)
Ở20 MPa, cả hai nhánh cho350 m³/ngày và cùng đạo hàm−J. Ở15;10;0 MPa, nhánh dưới cho505,555556;622,222222;738,888889 m³/ngày. Việc ghép giữ cả giá trị và độ dốc ở điểm nối, nhưng vẫn cần dữ liệu kiểm chứng ngoài điểm nối.
8. Phạm vi sử dụng
Nước có thể cần một nhánh riêng trước khi cộng thành tổng lỏng; không tự dùng lưu lượng dầu cộng nước để hiệu chỉnh đường chỉ dành cho dầu. Áp suất đáy phải cùng mốc và đã được xử lý nhất quán với gradient trong giếng. Khi p_r, nước hoặc trạng thái hoàn thiện thay đổi, cần đánh giá lại đường IPR thay vì chỉ dịch một điểm trên đồ thị cũ.
Một hồ sơ tốt lưu dữ liệu thử, áp suất bọt, giả thiết nhánh, phần dư và các kịch bản sai số. Đường Vogel khi đó hỗ trợ so sánh khả năng cung cấp vào giếng; lưu lượng thực còn cần đồng thời đáp ứng đường nâng và điều kiện hệ thống phía trên.
Nguồn tham khảo
[1] Nodal Analysis - IPR / VLP — Reservoir IPR — Whitson user manual. Truy cập ngày 09/11/2026.
[2] Vogel inflow — S&P Global, Harmony Enterprise2026.1. Truy cập ngày 09/11/2026.
[3] TPG4245 Regular Exam2023 — Problems1 and6 — Prof. Milan Stanko, Norwegian University of Science and Technology. Truy cập ngày 09/11/2026.
Rà soát nội dung: 09/11/2026.
