Chỉ số năng suất trả lời một câu hỏi có điều kiện
Một giếng đạt lưu lượng cao chưa đủ chứng minh vỉa có khả năng cung cấp tốt hơn giếng khác. Cần biết lưu lượng được tạo ra với mức chênh áp nào, loại lưu chất nào, tại thời điểm nào và theo điều kiện quy chiếu nào. Chỉ số năng suất PI hoặc J thường được dùng để liên hệ lưu lượng với drawdown tại vùng tiếp xúc giếng–vỉa. SLB mô tả PI theo khả năng cung cấp lưu chất của vỉa trên một đơn vị chênh áp. [1]
Trong bài này, mục tiêu là phân tích một mô hình ổn định đơn pha để hiểu vai trò độ thấm, chiều dày và skin. Ví dụ dùng số liệu giả định. Không dự báo sản lượng của mỏ cụ thể và không suy ra chương trình can thiệp giếng. Trọng tâm là làm rõ mỗi tham số xuất hiện ở đâu trong phương trình và những kết luận nào không thể rút ra chỉ từ một điểm đo.
Mô hình dòng xuyên tâm được lựa chọn
Xét một giếng thẳng đứng xuyên hết chiều dày hữu hiệu h của một vỉa đồng nhất, đẳng hướng. Dòng đi ngang, hội tụ theo hướng xuyên tâm. Độ thấm k, độ nhớt μ và tính chất lưu chất được coi không đổi. Áp suất p_e được giữ tại bán kính ngoài r_e; áp suất đáy giếng đang chảy p_wf được xác định tại bán kính r_w. Không có ảnh hưởng ranh giới không đều, dòng nhiều pha hoặc dòng quán tính đáng kể.
Darcy liên hệ lưu lượng riêng qua môi trường rỗng với gradient áp suất và tỷ số k/μ. Học liệu Penn State trình bày dạng quan hệ này và phân biệt vận tốc Darcy với vận tốc trong lỗ rỗng. [2] Trong hình học xuyên tâm, tiết diện dòng chảy qua một mặt trụ có bán kính r là 2πrh. Như vậy diện tích không cố định theo vị trí; càng đến gần giếng, cùng lưu lượng đi qua diện tích càng nhỏ.
Từ Darcy đến công thức lưu lượng
Lấy q_r là lưu lượng thể tích ở điều kiện vỉa, chiều dương hướng vào giếng. Với áp suất tăng theo r, có thể viết q_r = (2πrhk/μ)(dp/dr). Tách biến và tích phân từ r_w đến r_e cho quan hệ q_r μ ln(r_e/r_w) = 2πkh(p_e − p_wf). Logarithm là log tự nhiên, không phải log cơ số mười.
Đưa skin s vào như một phần trở lực không thứ nguyên bổ sung của mô hình, ta có:
q_r = 2πkhΔp / {μ[ln(r_e/r_w) + s]}
q_s = q_r/B; J_s = q_s/Δp
Δp = p_e − p_wf, Pa; k tính bằng m²; h, r_e và r_w bằng m; μ bằng Pa·s; q_r và q_s bằng m³/s theo điều kiện tương ứng. B là tỷ số thể tích vỉa trên thể tích bề mặt của cùng lượng lưu chất, không có đơn vị. J_s có đơn vị m³/(s·Pa), hoặc có thể quy đổi nhất quán sang m³/ngày/MPa.
Skin dương biểu thị tác động làm suy giảm năng suất so với trạng thái lý tưởng được chọn, còn skin âm thường phản ánh khả năng cung cấp được cải thiện. Đây là cách diễn giải của SLB. [3] Tuy nhiên, giá trị skin không tự cho biết duy nhất cơ chế hư hại hoặc giải pháp xử lý. Ý nghĩa của nó phụ thuộc mô hình đối chiếu và dữ liệu được dùng để suy ra.
Ví dụ giả định ở hệ SI
Chọn k = 1,0 × 10⁻¹³ m²; h = 18 m; μ = 0,002 Pa·s; r_e = 500 m; r_w = 0,10 m; s = 3; Δp = 5 MPa và B = 1,20. Khi thay số phải đổi Δp thành 5.000.000 Pa. Không sử dụng một hệ số chuyển đổi có sẵn cho mD, cP, ft và psi trong cùng phép tính SI này.
Tỷ số r_e/r_w = 5.000; ln(5.000) xấp xỉ 8,5172. Tổng trở lực không thứ nguyên là 11,5172. Tử số 2πkhΔp bằng khoảng 5,6549 × 10⁻⁵ theo đơn vị kết hợp của phương trình. Chia cho μ và phần logarithm cho q_r xấp xỉ 0,00245497 m³/s tại điều kiện vỉa.
Đổi sang điều kiện bề mặt: q_s = 0,00245497/1,20, rồi nhân 86.400 s/ngày, được khoảng 176,758 m³/ngày. Chỉ số năng suất theo cơ sở này bằng 176,758/5 = 35,352 m³/ngày/MPa. Nếu quên chia B, lưu lượng bề mặt bị ghi cao hơn 20% trong ví dụ, dù phần giải dòng trong vỉa hoàn toàn đúng.
So sánh skin mà không thay các tham số khác
Giữ toàn bộ dữ liệu và cho s = 0, lưu lượng bề mặt của mô hình tăng lên khoảng 239,017 m³/ngày. Với s = 6, lưu lượng giảm còn khoảng 140,230 m³/ngày. Như vậy việc chuyển từ s = 3 về s = 0 trong bài toán giả định làm lưu lượng tăng khoảng 35,2%, không phải tăng theo một tỷ lệ tuyến tính trực tiếp với số skin.
Lý do nằm ở mẫu số: s được cộng với ln(r_e/r_w), chứ không thay thế toàn bộ trở lực. Nếu phần logarithm rất lớn, cùng một thay đổi skin tạo tỷ lệ tăng lưu lượng khác. Do đó không nên dùng một quy tắc như “giảm một đơn vị skin sẽ tăng một tỷ lệ sản lượng cố định” cho mọi giếng.
Cũng không thể xem kết quả 239,017 m³/ngày là mục tiêu can thiệp đã được chứng minh. Trong thực tế, thay đổi vùng gần giếng có thể đồng thời làm thay đổi dòng nhiều pha, điều kiện nâng, áp suất đáy và giới hạn thiết bị. Bài so sánh chỉ cô lập một tham số để giải thích cơ chế của phương trình.
Độ nhạy độ thấm và bán kính ngoài
Ở cùng các giả thiết, q tỷ lệ thuận với k và h, nghịch với μ. Tăng k 10% cho lưu lượng khoảng 194,433 m³/ngày, đúng bằng 1,10 lần trường hợp cơ sở. Nhưng tăng r_e từ 500 lên 1.000 m chỉ làm lưu lượng giảm xuống khoảng 166,724 m³/ngày vì r_e xuất hiện bên trong logarithm. Độ nhạy của các tham số vì vậy không giống nhau dù chúng cùng xuất hiện trong danh sách đầu vào.
Một phép suy luận ngược cần đặc biệt thận trọng. Nếu chỉ biết một cặp q và Δp, có thể tìm nhiều tổ hợp kh và s cho cùng J. Tăng kh đồng thời tăng skin có thể giữ nguyên lưu lượng. Bởi vậy, một điểm vận hành không đủ tách riêng độ thấm và hư hại gần giếng; cần thêm bằng chứng về địa chất, chiều dày, tính chất lưu chất và thử giếng phù hợp.
Không trộn lẫn các chế độ dòng
Áp suất p_e trong mô hình là áp suất biên tại r_e, không tự động bằng áp suất vỉa trung bình của mọi chế độ khai thác. Bài toán giả ổn định trong vùng kín hoặc bài toán chuyển tiếp có dạng quan hệ và tham số thời gian khác. Không nên thêm hay bỏ một hằng số hình học trong mẫu số mà không xác định mình đang dùng mô hình nào.
Tương tự, khi áp suất giảm qua điểm bọt hoặc có nước và khí cùng chuyển động, độ thấm tương đối và tính chất pha có thể làm quan hệ lưu lượng–drawdown không còn tuyến tính như ở đây. Dùng một PI duy nhất để ngoại suy quá xa dữ liệu sẽ che giấu các thay đổi ấy. Áp suất đo tại miệng giếng cũng không thể thay thẳng cho p_wf nếu chưa mô tả cột lưu chất và tổn thất trong giếng.
Một hồ sơ có thể kiểm tra được nên lưu cả điều kiện quy chiếu lưu lượng, áp suất biên, nhiệt độ, B, μ, cơ sở kh và cách suy ra skin. Khi các đầu vào chưa chắc chắn, nên báo một tập kịch bản có giải thích thay vì một lưu lượng chính xác tới nhiều chữ số. Công thức xuyên tâm hữu ích nhất khi giúp người đọc thấy rõ dữ liệu nào chi phối kết quả và câu hỏi nào vẫn còn mở.
Nguồn tham khảo
[1] productivity index (PI) — SLB Energy Glossary. Truy cập ngày 09/11/2026.
[2] Lesson 6: Flow and Solute Transport Processes in 1D Systems — 6.2 Advection — Pennsylvania State University, College of Earth and Mineral Sciences. Truy cập ngày 09/11/2026.
[3] skin — SLB Energy Glossary. Truy cập ngày 09/11/2026.
Rà soát nội dung: 09/11/2026.
Một giếng đạt lưu lượng cao chưa đủ chứng minh vỉa có khả năng cung cấp tốt hơn giếng khác. Cần biết lưu lượng được tạo ra với mức chênh áp nào, loại lưu chất nào, tại thời điểm nào và theo điều kiện quy chiếu nào. Chỉ số năng suất PI hoặc J thường được dùng để liên hệ lưu lượng với drawdown tại vùng tiếp xúc giếng–vỉa. SLB mô tả PI theo khả năng cung cấp lưu chất của vỉa trên một đơn vị chênh áp. [1]
Trong bài này, mục tiêu là phân tích một mô hình ổn định đơn pha để hiểu vai trò độ thấm, chiều dày và skin. Ví dụ dùng số liệu giả định. Không dự báo sản lượng của mỏ cụ thể và không suy ra chương trình can thiệp giếng. Trọng tâm là làm rõ mỗi tham số xuất hiện ở đâu trong phương trình và những kết luận nào không thể rút ra chỉ từ một điểm đo.
Mô hình dòng xuyên tâm được lựa chọn
Xét một giếng thẳng đứng xuyên hết chiều dày hữu hiệu h của một vỉa đồng nhất, đẳng hướng. Dòng đi ngang, hội tụ theo hướng xuyên tâm. Độ thấm k, độ nhớt μ và tính chất lưu chất được coi không đổi. Áp suất p_e được giữ tại bán kính ngoài r_e; áp suất đáy giếng đang chảy p_wf được xác định tại bán kính r_w. Không có ảnh hưởng ranh giới không đều, dòng nhiều pha hoặc dòng quán tính đáng kể.
Darcy liên hệ lưu lượng riêng qua môi trường rỗng với gradient áp suất và tỷ số k/μ. Học liệu Penn State trình bày dạng quan hệ này và phân biệt vận tốc Darcy với vận tốc trong lỗ rỗng. [2] Trong hình học xuyên tâm, tiết diện dòng chảy qua một mặt trụ có bán kính r là 2πrh. Như vậy diện tích không cố định theo vị trí; càng đến gần giếng, cùng lưu lượng đi qua diện tích càng nhỏ.
Từ Darcy đến công thức lưu lượng
Lấy q_r là lưu lượng thể tích ở điều kiện vỉa, chiều dương hướng vào giếng. Với áp suất tăng theo r, có thể viết q_r = (2πrhk/μ)(dp/dr). Tách biến và tích phân từ r_w đến r_e cho quan hệ q_r μ ln(r_e/r_w) = 2πkh(p_e − p_wf). Logarithm là log tự nhiên, không phải log cơ số mười.
Đưa skin s vào như một phần trở lực không thứ nguyên bổ sung của mô hình, ta có:
q_r = 2πkhΔp / {μ[ln(r_e/r_w) + s]}
q_s = q_r/B; J_s = q_s/Δp
Δp = p_e − p_wf, Pa; k tính bằng m²; h, r_e và r_w bằng m; μ bằng Pa·s; q_r và q_s bằng m³/s theo điều kiện tương ứng. B là tỷ số thể tích vỉa trên thể tích bề mặt của cùng lượng lưu chất, không có đơn vị. J_s có đơn vị m³/(s·Pa), hoặc có thể quy đổi nhất quán sang m³/ngày/MPa.
Skin dương biểu thị tác động làm suy giảm năng suất so với trạng thái lý tưởng được chọn, còn skin âm thường phản ánh khả năng cung cấp được cải thiện. Đây là cách diễn giải của SLB. [3] Tuy nhiên, giá trị skin không tự cho biết duy nhất cơ chế hư hại hoặc giải pháp xử lý. Ý nghĩa của nó phụ thuộc mô hình đối chiếu và dữ liệu được dùng để suy ra.
Ví dụ giả định ở hệ SI
Chọn k = 1,0 × 10⁻¹³ m²; h = 18 m; μ = 0,002 Pa·s; r_e = 500 m; r_w = 0,10 m; s = 3; Δp = 5 MPa và B = 1,20. Khi thay số phải đổi Δp thành 5.000.000 Pa. Không sử dụng một hệ số chuyển đổi có sẵn cho mD, cP, ft và psi trong cùng phép tính SI này.
Tỷ số r_e/r_w = 5.000; ln(5.000) xấp xỉ 8,5172. Tổng trở lực không thứ nguyên là 11,5172. Tử số 2πkhΔp bằng khoảng 5,6549 × 10⁻⁵ theo đơn vị kết hợp của phương trình. Chia cho μ và phần logarithm cho q_r xấp xỉ 0,00245497 m³/s tại điều kiện vỉa.
Đổi sang điều kiện bề mặt: q_s = 0,00245497/1,20, rồi nhân 86.400 s/ngày, được khoảng 176,758 m³/ngày. Chỉ số năng suất theo cơ sở này bằng 176,758/5 = 35,352 m³/ngày/MPa. Nếu quên chia B, lưu lượng bề mặt bị ghi cao hơn 20% trong ví dụ, dù phần giải dòng trong vỉa hoàn toàn đúng.
So sánh skin mà không thay các tham số khác
Giữ toàn bộ dữ liệu và cho s = 0, lưu lượng bề mặt của mô hình tăng lên khoảng 239,017 m³/ngày. Với s = 6, lưu lượng giảm còn khoảng 140,230 m³/ngày. Như vậy việc chuyển từ s = 3 về s = 0 trong bài toán giả định làm lưu lượng tăng khoảng 35,2%, không phải tăng theo một tỷ lệ tuyến tính trực tiếp với số skin.
Lý do nằm ở mẫu số: s được cộng với ln(r_e/r_w), chứ không thay thế toàn bộ trở lực. Nếu phần logarithm rất lớn, cùng một thay đổi skin tạo tỷ lệ tăng lưu lượng khác. Do đó không nên dùng một quy tắc như “giảm một đơn vị skin sẽ tăng một tỷ lệ sản lượng cố định” cho mọi giếng.
Cũng không thể xem kết quả 239,017 m³/ngày là mục tiêu can thiệp đã được chứng minh. Trong thực tế, thay đổi vùng gần giếng có thể đồng thời làm thay đổi dòng nhiều pha, điều kiện nâng, áp suất đáy và giới hạn thiết bị. Bài so sánh chỉ cô lập một tham số để giải thích cơ chế của phương trình.
Độ nhạy độ thấm và bán kính ngoài
Ở cùng các giả thiết, q tỷ lệ thuận với k và h, nghịch với μ. Tăng k 10% cho lưu lượng khoảng 194,433 m³/ngày, đúng bằng 1,10 lần trường hợp cơ sở. Nhưng tăng r_e từ 500 lên 1.000 m chỉ làm lưu lượng giảm xuống khoảng 166,724 m³/ngày vì r_e xuất hiện bên trong logarithm. Độ nhạy của các tham số vì vậy không giống nhau dù chúng cùng xuất hiện trong danh sách đầu vào.
Một phép suy luận ngược cần đặc biệt thận trọng. Nếu chỉ biết một cặp q và Δp, có thể tìm nhiều tổ hợp kh và s cho cùng J. Tăng kh đồng thời tăng skin có thể giữ nguyên lưu lượng. Bởi vậy, một điểm vận hành không đủ tách riêng độ thấm và hư hại gần giếng; cần thêm bằng chứng về địa chất, chiều dày, tính chất lưu chất và thử giếng phù hợp.
Không trộn lẫn các chế độ dòng
Áp suất p_e trong mô hình là áp suất biên tại r_e, không tự động bằng áp suất vỉa trung bình của mọi chế độ khai thác. Bài toán giả ổn định trong vùng kín hoặc bài toán chuyển tiếp có dạng quan hệ và tham số thời gian khác. Không nên thêm hay bỏ một hằng số hình học trong mẫu số mà không xác định mình đang dùng mô hình nào.
Tương tự, khi áp suất giảm qua điểm bọt hoặc có nước và khí cùng chuyển động, độ thấm tương đối và tính chất pha có thể làm quan hệ lưu lượng–drawdown không còn tuyến tính như ở đây. Dùng một PI duy nhất để ngoại suy quá xa dữ liệu sẽ che giấu các thay đổi ấy. Áp suất đo tại miệng giếng cũng không thể thay thẳng cho p_wf nếu chưa mô tả cột lưu chất và tổn thất trong giếng.
Một hồ sơ có thể kiểm tra được nên lưu cả điều kiện quy chiếu lưu lượng, áp suất biên, nhiệt độ, B, μ, cơ sở kh và cách suy ra skin. Khi các đầu vào chưa chắc chắn, nên báo một tập kịch bản có giải thích thay vì một lưu lượng chính xác tới nhiều chữ số. Công thức xuyên tâm hữu ích nhất khi giúp người đọc thấy rõ dữ liệu nào chi phối kết quả và câu hỏi nào vẫn còn mở.
Nguồn tham khảo
[1] productivity index (PI) — SLB Energy Glossary. Truy cập ngày 09/11/2026.
[2] Lesson 6: Flow and Solute Transport Processes in 1D Systems — 6.2 Advection — Pennsylvania State University, College of Earth and Mineral Sciences. Truy cập ngày 09/11/2026.
[3] skin — SLB Energy Glossary. Truy cập ngày 09/11/2026.
Rà soát nội dung: 09/11/2026.
