1. Khí chiếm phần lưu lượng lớn chưa chắc chiếm cùng phần thể tích ống
Trong dòng khí–lỏng, khí có thể đi nhanh hơn lỏng. Vì vậy, một tỷ lệ khí cao trong lưu lượng thể tích không tự cho biết trong đoạn ống đang có bao nhiêu chất lỏng. SLB định nghĩa holdup là phần pha hiện diện trong một đoạn ống và phân biệt nó với cut, tức phần đóng góp vào lưu lượng.[1]
Phần trọng lực trong phương trình áp suất phụ thuộc khối lượng đang hiện diện trên một đơn vị thể tích ống. Nếu thay holdup lỏng bằng tỷ phần lưu lượng lỏng, kết quả có thể sai đáng kể ngay cả khi tổng lưu lượng khối lượng được tính đúng.
2. Vận tốc bề ngoài và vận tốc pha
Xét tiết diện A của ống có phần khí α và phần lỏng H_L=1−α. Gọi q_g,q_L là lưu lượng thể tích địa phương, m³/s, đã đưa về cùng áp suất và nhiệt độ tại tiết diện. Vận tốc bề ngoài là j_g=q_g/A, j_L=q_L/A; vận tốc trung bình trong phần tiết diện của từng pha là v_g=j_g/α và v_L=j_L/H_L.
Đặt S=v_g/v_L, với hai pha cùng chảy theo chiều dương. Tự biến đổi bảo toàn lưu lượng cho:
α = j_g/(j_g+Sj_L)
H_L = Sj_L/(j_g+Sj_L)
β = q_g/(q_g+q_L) = j_g/(j_g+j_L)
β là tỷ phần lưu lượng khí không trượt, không mặc nhiên bằng α. Chỉ khi S=1 thì α=β. Bài khảo sát áp suất trong học liệu NTNU yêu cầu tách riêng các đại lượng này, vận tốc pha và gradient thủy tĩnh.[2] Các số bên dưới được tự dựng để kiểm quan hệ giữa chúng.
3. Ví dụ tại một tiết diện
Chọn đường kính trong D=0,10 m, A=0,007853981634 m², q_L=0,003 m³/s và q_g=0,009 m³/s. Khi đó j_L=0,381971863 m/s, j_g=1,145915590 m/s và β=0,75. Giả định S=2 ở trạng thái đang xét, ta có α=0,60 và H_L=0,40.
Vận tốc pha lỏng là0,954929659 m/s, khí là1,909859317 m/s, đúng tỷ số2. Nếu dùng H_L=1−β=0,25, các vận tốc suy ra sẽ bằng nhau, trái với giả thiết S=2 vừa đặt. Kiểm lại lưu lượng qua H_LAv_L và αAv_g giúp phát hiện sự pha trộn hai định nghĩa.
Với khối lượng riêng địa phương ρ_L=900 vàρ_g=60 kg/m³, mật độ hỗn hợp theo thể tích hiện diện là:
ρ_h = H_Lρ_L+(1−H_L)ρ_g = 396 kg/m³
Trong khi mật độ tính từ tỷ số tổng khối lượng qua tiết diện trên tổng lưu lượng thể tích là(ρ_Lq_L+ρ_gq_g)/(q_L+q_g)=270 kg/m³. Số270 vẫn có ý nghĩa cho phép chia hai thông lượng, nhưng không phải mật độ khối chất lưu đang chiếm ống khi hai pha trượt.
4. Phần trọng lực cần tích phân theo chiều cao
Trong mô hình trung bình tiết diện, phần tăng áp suất khi đi xuống một đoạn thẳng đứng là d p_g=ρ_h g dz, lấy z dương xuống dưới. Với đoạn nghiêng, dz=cosθ ds khi θ đo từ phương thẳng đứng. g=9,81 m/s², s và z bằng m, ρ_h bằng kg/m³ thì áp suất thu được bằng Pa.
Nếu giả định trạng thái tiết diện của ví dụ được giữ trên một đoạn cao1.500 m, phần trọng lực là396×9,81×1.500=5,827140 MPa. Dùng mật độ không trượt270 cho3,973050 MPa, thấp hơn1,854090 MPa. Giả định giữ mật độ trên toàn đoạn chỉ phục vụ so sánh; trong giếng khí nén được, tính chất thường phải đổi theo áp suất và nhiệt độ.
Whitson mô tả việc tích phân gradient theo quỹ đạo giếng và PVT; phương trình áp suất đầy đủ còn các thành phần ma sát và gia tốc.[3] Bởi vậy, phần trọng lực đúng vẫn chưa phải toàn bộ chênh áp đáy–miệng của một giếng đang chảy.
5. Ba đoạn minh họa phải giữ lưu lượng khối lượng
Xét ba đoạn từ sâu lên nông, mỗi đoạn có chênh chiều sâu đứng500 m. Giả định không trao đổi khối lượng giữa pha, không dòng nhánh, ρ_L=900 kg/m³ và q_L=0,003 m³/s ở cả ba đoạn. Lưu lượng khối lượng khí cố định0,54 kg/s; khối lượng riêng khí đại diện lần lượt120;60;30 kg/m³. Vì vậy q_g lần lượt0,0045;0,009;0,018 m³/s, không được giữ cùng một lưu lượng thể tích khí.
Đặt các tỷ số trượt đại diện S lần lượt1,5;2;3. Công thức trên cho H_L lần lượt0,50;0,40;1/3 vàρ_h lần lượt510;396;320 kg/m³. Đây là bộ trạng thái giả định để tính tích phân từng đoạn, chưa phải nghiệm tự nhất quán của một tương quan PVT và holdup theo áp suất.
Δp_g = 9,81×500×(510+396+320)/10⁶
Δp_g = 6,013530 MPa
Nếu chiều dài đo dọc lỗ của ba đoạn là800;700;600 m, dùng các chiều dài này thay500 m trong phần trọng lực sẽ cho8,605332 MPa, sai vì bỏ góc nghiêng. Ma sát lại phụ thuộc chiều dài dọc đường đi, nên không thể dùng một loại chiều sâu cho cả hai thành phần mà không xét hình học.
6. Trung bình riêng các đầu vào có thể phá cân bằng
Trung bình H_L của ba đoạn là0,411111111; trung bìnhρ_g là70 kg/m³. Nhân các trung bình này choρ_h=411,222222 kg/m³, trong khi trung bình đúng của ba mật độ đoạn là408,666667 kg/m³. Phần chênh đến từ mối liên hệ giữa holdup và mật độ khí theo vị trí.
Giữ bảng các đoạn với trọng số chênh cao sẽ bảo toàn tích phân trọng lực. Nếu đoạn có độ dài khác nhau, cần trọng số đúng thay vì trung bình ba giá trị như nhau. Với quỹ đạo lên xuống, tích phân phải giữ dấu dz; dùng tổng chiều sâu đứng tuyệt đối có thể cộng cả các phần đáng ra bù nhau.
Trong một mô hình dự báo đầy đủ, áp suất quyết định mật độ và phần khí; các đại lượng đó lại quyết định gradient áp suất. Cần giải đồng thời hoặc lặp theo đoạn, rồi kiểm độ nhạy khi chia lưới nhỏ hơn. Chỉ tăng số chữ số của holdup nhập tay không giải quyết vòng phụ thuộc này.
7. Suy holdup từ chênh áp phải tách các thành phần khác
Xét một phép suy riêng với đoạn thẳng đứng1.500 m, mật độ pha cố định900 và60 kg/m³. Giả định đo tổng chênh áp7,20 MPa, mô hình độc lập ước tính ma sát1,10 MPa và gia tốc0,05 MPa. Phần trọng lực còn6,05 MPa, tương ứngρ_h=411,145090 kg/m³ và:
H_L = (ρ_h−ρ_g)/(ρ_L−ρ_g) = 0,418029869
Nếu gán toàn bộ7,20 MPa cho trọng lực, H_L suy ra thành0,511067424. Phép tính sau không sai đại số, nhưng dùng sai ranh giới của thành phần áp suất. Kết quả suy holdup chỉ đáng tin trong mức độ các thành phần đã trừ và mật độ pha đáng tin.
Ở mật độ pha cố định, dΔp_g/dH_L=(ρ_L−ρ_g)gH. Tăng H_L thêm0,05 làm phần trọng lực của đoạn tăng0,618030 MPa. Vì thế khác biệt holdup vài điểm phần trăm có thể đáng kể trong suy áp suất đáy.
8. Điều kiện giới hạn
Quan hệ theo S ở trên giả định hai pha cùng chiều và các vận tốc trung bình có nghĩa trên khoảng xét. Dòng slug dao động, dòng ngược chiều hoặc lỏng gần đứng yên có thể làm cách biểu diễn bằng một S cố định không phù hợp. Holdup theo diện tích một tiết diện cũng cần được lấy trung bình đúng khi dùng cho một đoạn có cấu trúc thay đổi theo thời gian.
Bảng tính nên lưu riêng lưu lượng quy chiếu, lưu lượng địa phương, tỷ phần thông lượng, holdup và vận tốc pha. Kiểm bảo toàn khối lượng và mốc chiều sâu trước sẽ giúp phân biệt sai khác mô hình nhiều pha với lỗi nhập liệu đơn giản, đồng thời giữ rõ giới hạn của việc dùng phần trọng lực để suy toàn bộ áp suất giếng.
Nguồn tham khảo
[1] holdup — SLB Energy Glossary. Truy cập ngày 09/11/2026.
[2] TPG4245 Regular Exam2023 — Problem5: pressure survey and phase velocities — Prof. Milan Stanko, Norwegian University of Science and Technology. Truy cập ngày 09/11/2026.
[3] Bottomhole Pressure Calculations — Governing Equations and well deviation survey — Whitson user manual. Truy cập ngày 09/11/2026.
Rà soát nội dung: 09/11/2026.
Trong dòng khí–lỏng, khí có thể đi nhanh hơn lỏng. Vì vậy, một tỷ lệ khí cao trong lưu lượng thể tích không tự cho biết trong đoạn ống đang có bao nhiêu chất lỏng. SLB định nghĩa holdup là phần pha hiện diện trong một đoạn ống và phân biệt nó với cut, tức phần đóng góp vào lưu lượng.[1]
Phần trọng lực trong phương trình áp suất phụ thuộc khối lượng đang hiện diện trên một đơn vị thể tích ống. Nếu thay holdup lỏng bằng tỷ phần lưu lượng lỏng, kết quả có thể sai đáng kể ngay cả khi tổng lưu lượng khối lượng được tính đúng.
2. Vận tốc bề ngoài và vận tốc pha
Xét tiết diện A của ống có phần khí α và phần lỏng H_L=1−α. Gọi q_g,q_L là lưu lượng thể tích địa phương, m³/s, đã đưa về cùng áp suất và nhiệt độ tại tiết diện. Vận tốc bề ngoài là j_g=q_g/A, j_L=q_L/A; vận tốc trung bình trong phần tiết diện của từng pha là v_g=j_g/α và v_L=j_L/H_L.
Đặt S=v_g/v_L, với hai pha cùng chảy theo chiều dương. Tự biến đổi bảo toàn lưu lượng cho:
α = j_g/(j_g+Sj_L)
H_L = Sj_L/(j_g+Sj_L)
β = q_g/(q_g+q_L) = j_g/(j_g+j_L)
β là tỷ phần lưu lượng khí không trượt, không mặc nhiên bằng α. Chỉ khi S=1 thì α=β. Bài khảo sát áp suất trong học liệu NTNU yêu cầu tách riêng các đại lượng này, vận tốc pha và gradient thủy tĩnh.[2] Các số bên dưới được tự dựng để kiểm quan hệ giữa chúng.
3. Ví dụ tại một tiết diện
Chọn đường kính trong D=0,10 m, A=0,007853981634 m², q_L=0,003 m³/s và q_g=0,009 m³/s. Khi đó j_L=0,381971863 m/s, j_g=1,145915590 m/s và β=0,75. Giả định S=2 ở trạng thái đang xét, ta có α=0,60 và H_L=0,40.
Vận tốc pha lỏng là0,954929659 m/s, khí là1,909859317 m/s, đúng tỷ số2. Nếu dùng H_L=1−β=0,25, các vận tốc suy ra sẽ bằng nhau, trái với giả thiết S=2 vừa đặt. Kiểm lại lưu lượng qua H_LAv_L và αAv_g giúp phát hiện sự pha trộn hai định nghĩa.
Với khối lượng riêng địa phương ρ_L=900 vàρ_g=60 kg/m³, mật độ hỗn hợp theo thể tích hiện diện là:
ρ_h = H_Lρ_L+(1−H_L)ρ_g = 396 kg/m³
Trong khi mật độ tính từ tỷ số tổng khối lượng qua tiết diện trên tổng lưu lượng thể tích là(ρ_Lq_L+ρ_gq_g)/(q_L+q_g)=270 kg/m³. Số270 vẫn có ý nghĩa cho phép chia hai thông lượng, nhưng không phải mật độ khối chất lưu đang chiếm ống khi hai pha trượt.
4. Phần trọng lực cần tích phân theo chiều cao
Trong mô hình trung bình tiết diện, phần tăng áp suất khi đi xuống một đoạn thẳng đứng là d p_g=ρ_h g dz, lấy z dương xuống dưới. Với đoạn nghiêng, dz=cosθ ds khi θ đo từ phương thẳng đứng. g=9,81 m/s², s và z bằng m, ρ_h bằng kg/m³ thì áp suất thu được bằng Pa.
Nếu giả định trạng thái tiết diện của ví dụ được giữ trên một đoạn cao1.500 m, phần trọng lực là396×9,81×1.500=5,827140 MPa. Dùng mật độ không trượt270 cho3,973050 MPa, thấp hơn1,854090 MPa. Giả định giữ mật độ trên toàn đoạn chỉ phục vụ so sánh; trong giếng khí nén được, tính chất thường phải đổi theo áp suất và nhiệt độ.
Whitson mô tả việc tích phân gradient theo quỹ đạo giếng và PVT; phương trình áp suất đầy đủ còn các thành phần ma sát và gia tốc.[3] Bởi vậy, phần trọng lực đúng vẫn chưa phải toàn bộ chênh áp đáy–miệng của một giếng đang chảy.
5. Ba đoạn minh họa phải giữ lưu lượng khối lượng
Xét ba đoạn từ sâu lên nông, mỗi đoạn có chênh chiều sâu đứng500 m. Giả định không trao đổi khối lượng giữa pha, không dòng nhánh, ρ_L=900 kg/m³ và q_L=0,003 m³/s ở cả ba đoạn. Lưu lượng khối lượng khí cố định0,54 kg/s; khối lượng riêng khí đại diện lần lượt120;60;30 kg/m³. Vì vậy q_g lần lượt0,0045;0,009;0,018 m³/s, không được giữ cùng một lưu lượng thể tích khí.
Đặt các tỷ số trượt đại diện S lần lượt1,5;2;3. Công thức trên cho H_L lần lượt0,50;0,40;1/3 vàρ_h lần lượt510;396;320 kg/m³. Đây là bộ trạng thái giả định để tính tích phân từng đoạn, chưa phải nghiệm tự nhất quán của một tương quan PVT và holdup theo áp suất.
Δp_g = 9,81×500×(510+396+320)/10⁶
Δp_g = 6,013530 MPa
Nếu chiều dài đo dọc lỗ của ba đoạn là800;700;600 m, dùng các chiều dài này thay500 m trong phần trọng lực sẽ cho8,605332 MPa, sai vì bỏ góc nghiêng. Ma sát lại phụ thuộc chiều dài dọc đường đi, nên không thể dùng một loại chiều sâu cho cả hai thành phần mà không xét hình học.
6. Trung bình riêng các đầu vào có thể phá cân bằng
Trung bình H_L của ba đoạn là0,411111111; trung bìnhρ_g là70 kg/m³. Nhân các trung bình này choρ_h=411,222222 kg/m³, trong khi trung bình đúng của ba mật độ đoạn là408,666667 kg/m³. Phần chênh đến từ mối liên hệ giữa holdup và mật độ khí theo vị trí.
Giữ bảng các đoạn với trọng số chênh cao sẽ bảo toàn tích phân trọng lực. Nếu đoạn có độ dài khác nhau, cần trọng số đúng thay vì trung bình ba giá trị như nhau. Với quỹ đạo lên xuống, tích phân phải giữ dấu dz; dùng tổng chiều sâu đứng tuyệt đối có thể cộng cả các phần đáng ra bù nhau.
Trong một mô hình dự báo đầy đủ, áp suất quyết định mật độ và phần khí; các đại lượng đó lại quyết định gradient áp suất. Cần giải đồng thời hoặc lặp theo đoạn, rồi kiểm độ nhạy khi chia lưới nhỏ hơn. Chỉ tăng số chữ số của holdup nhập tay không giải quyết vòng phụ thuộc này.
7. Suy holdup từ chênh áp phải tách các thành phần khác
Xét một phép suy riêng với đoạn thẳng đứng1.500 m, mật độ pha cố định900 và60 kg/m³. Giả định đo tổng chênh áp7,20 MPa, mô hình độc lập ước tính ma sát1,10 MPa và gia tốc0,05 MPa. Phần trọng lực còn6,05 MPa, tương ứngρ_h=411,145090 kg/m³ và:
H_L = (ρ_h−ρ_g)/(ρ_L−ρ_g) = 0,418029869
Nếu gán toàn bộ7,20 MPa cho trọng lực, H_L suy ra thành0,511067424. Phép tính sau không sai đại số, nhưng dùng sai ranh giới của thành phần áp suất. Kết quả suy holdup chỉ đáng tin trong mức độ các thành phần đã trừ và mật độ pha đáng tin.
Ở mật độ pha cố định, dΔp_g/dH_L=(ρ_L−ρ_g)gH. Tăng H_L thêm0,05 làm phần trọng lực của đoạn tăng0,618030 MPa. Vì thế khác biệt holdup vài điểm phần trăm có thể đáng kể trong suy áp suất đáy.
8. Điều kiện giới hạn
Quan hệ theo S ở trên giả định hai pha cùng chiều và các vận tốc trung bình có nghĩa trên khoảng xét. Dòng slug dao động, dòng ngược chiều hoặc lỏng gần đứng yên có thể làm cách biểu diễn bằng một S cố định không phù hợp. Holdup theo diện tích một tiết diện cũng cần được lấy trung bình đúng khi dùng cho một đoạn có cấu trúc thay đổi theo thời gian.
Bảng tính nên lưu riêng lưu lượng quy chiếu, lưu lượng địa phương, tỷ phần thông lượng, holdup và vận tốc pha. Kiểm bảo toàn khối lượng và mốc chiều sâu trước sẽ giúp phân biệt sai khác mô hình nhiều pha với lỗi nhập liệu đơn giản, đồng thời giữ rõ giới hạn của việc dùng phần trọng lực để suy toàn bộ áp suất giếng.
Nguồn tham khảo
[1] holdup — SLB Energy Glossary. Truy cập ngày 09/11/2026.
[2] TPG4245 Regular Exam2023 — Problem5: pressure survey and phase velocities — Prof. Milan Stanko, Norwegian University of Science and Technology. Truy cập ngày 09/11/2026.
[3] Bottomhole Pressure Calculations — Governing Equations and well deviation survey — Whitson user manual. Truy cập ngày 09/11/2026.
Rà soát nội dung: 09/11/2026.
