Một hệ số tổng cần cùng cơ sở thể tích
Trong mô hình áp suất vỉa, độ nén tổng thường nằm trong một ô nhỏ nhưng quyết định lượng tích trữ thay đổi khi áp suất giảm. Việc cộng số của dầu, nước, khí và đá chỉ có ý nghĩa khi các hệ số cùng đơn vị, cùng điều kiện và cùng định nghĩa thể tích. Độ nén của vật liệu hạt rắn không tự bằng độ nén thể tích lỗ rỗng.
SLB mô tả compressibility là thay đổi tương đối thể tích theo một đơn vị thay đổi áp suất. [1] Bài này dùng định nghĩa vi phân để tự dựng mô hình lưu trữ cục bộ, với các số liệu giả định. Không lấy kết quả làm hệ số đã hiệu chuẩn cho một vỉa thực hoặc suy trực tiếp dự báo thu hồi.
Độ nén pha và dấu của phần lỗ rỗng
Với một lượng lưu chất cố định, đẳng nhiệt và không đổi thành phần, định nghĩa c_i=−(1/V_i)(∂V_i/∂p). Hệ số dương nghĩa là tăng áp suất làm thể tích lưu chất giảm. Tương đương, c_i=(1/ρ_i)(∂ρ_i/∂p) vì khối lượng không đổi.
Đối với không gian chứa, chọn c_p=(1/V_p)(∂V_p/∂p) theo một đường ứng suất ngoài đã xác định. Trong ví dụ ta giữ ứng suất bao ngoài cố định và lấy c_p dương: tăng áp suất lỗ rỗng làm thể tích lỗ rỗng tăng; giảm áp suất làm nó giảm. Tài liệu iTOUGH2 xác định rõ tham số là độ nén không gian lỗ rỗng. [3]
Hai dấu khác nhau không mâu thuẫn. Khi áp suất giảm, lưu chất muốn nở nhưng chỗ chứa co lại; cả hai góp vào lượng có thể được đẩy ra. Nếu thí nghiệm báo hệ số theo ứng suất hữu hiệu, cần xét phép đổi biến và đường tải thay vì chép nguyên dấu. Nếu thể tích khối đá cũng thay đổi, không được thay V_p bằng độ rỗng φ mà bỏ sự thay đổi của thể tích khối.
Cộng có trọng số tại một trạng thái tham chiếu
Gọi S_o, S_w, S_g là tỷ phần thể tích dầu, nước và khí trong tổng thể tích lỗ rỗng, với S_o+S_w+S_g=1. Trong mô hình cục bộ mà các pha chịu cùng gia số áp suất, không có chuyển pha hay trao đổi thành phần, độ nén tổng là:
c_t = c_p + S_o c_o + S_w c_w + S_g c_g.
WellTest trình bày cấu trúc tổng gồm phần thành hệ và phần lưu chất có trọng số bão hòa. [2] c_t và từng c có đơn vị Pa⁻¹ nếu áp suất dùng Pa. Các S là tỷ phần từ 0 đến 1, không phải số phần trăm từ 0 đến 100. c_p đã quy về V_p nên được cộng một lần, không nhân thêm độ rỗng.
Xét giảm áp nhỏ Δp_dương=p_đầu−p_cuối. Lượng thể tích giải phóng xấp xỉ là ΔV≈V_p c_t Δp_dương. Phần nở lưu chất là V_pΣS_i c_i Δp_dương; phần giảm chỗ chứa là V_p c_p Δp_dương. Đây là khai triển bậc nhất quanh trạng thái đã chọn, chưa mô tả những thay đổi lớn của bão hòa, thành phần và hệ số trong một quá trình dài.
Ví dụ dầu–nước chưa có khí tự do
Cho V_khối=1.000.000 m³, φ=0,18, nên V_p=180.000 m³. Chọn S_w=0,25 và S_o=0,75; c_w=4,4×10⁻¹⁰ Pa⁻¹, c_o=1,0×10⁻⁹ Pa⁻¹, c_p=5,0×10⁻¹⁰ Pa⁻¹. Mọi giá trị là giả định của bài.
Khi đó c_t=1,36×10⁻⁹ Pa⁻¹. Với giảm áp 1 MPa, ΔV≈244,8 m³. Có thể kiểm riêng ba khoản: dầu nở 135 m³, nước nở 19,8 m³ và chỗ chứa giảm 90 m³. Tổng đúng 244,8 m³.
Nếu quên phần lỗ rỗng, kết quả còn 154,8 m³. Nếu nhân c_p thêm φ lần nữa, phần đá chỉ còn 16,2 m³, tổng thành 171 m³. Cả hai sai khác đều xuất phát từ định nghĩa cơ sở thể tích, không phải sai số làm tròn. Độ rỗng đã xuất hiện trong V_p=φV_khối.
Một ít khí có thể thay đổi tổng đáng kể
Ở trạng thái giả định thứ hai, giữ S_w=0,25, đưa S_g=0,02 và giảm S_o còn 0,73. Không giữ S_o=0,75 rồi thêm khí vì như vậy tổng bão hòa vượt một. Giả định khí lý tưởng tại p=20 MPa tuyệt đối, c_g=1/p=5,0×10⁻⁸ Pa⁻¹.
Tổng mới c_t=2,34×10⁻⁹ Pa⁻¹, lớn hơn trường hợp không khí 72,0588%. Lượng giải phóng bậc nhất cho cùng 1 MPa là 421,2 m³. Riêng phần S_g c_g chiếm 42,7350% c_t: nó rất đáng kể nhưng chưa chiếm quá một nửa. Vì thế cần tính trọng số, không thay phép tính bằng phát biểu rằng cứ có khí là mọi phần khác đều bỏ được.
Với S_w cố định, viết c_t=c_t,0+(c_g−c_o)S_g. Đạo hàm theo S_g bằng c_g−c_o=4,9×10⁻⁸ Pa⁻¹ trên một đơn vị bão hòa. Phần khí vượt tổng các phần còn lại khi S_g>c_t,0/(c_g+c_o)=0,026667 trong chính bộ hệ số này. Ngưỡng 2,6667% là kết quả của ví dụ, không phải ranh giới dùng chung cho vỉa.
Độ nén khí không phải hệ số Z
Z là hệ số không thứ nguyên trong pV=ZnRT. Còn c_g mang đơn vị nghịch đảo áp suất. Giữ n và T, lấy log rồi đạo hàm cho c_g=1/p−(1/Z)(dZ/dp). Do đó chỉ khi Z không đổi theo p mới được dùng chính xác 1/p trong mô hình đó.
Chọn riêng Z=0,90 và dZ/dp=0,003 MPa⁻¹ tại 20 MPa. Đổi đạo hàm thành 3×10⁻⁹ Pa⁻¹, ta được c_g=4,666667×10⁻⁸ Pa⁻¹. Thay vào trạng thái S_g=0,02 cho c_t=2,273333×10⁻⁹ Pa⁻¹. Chỉ biết một giá trị Z tại một điểm chưa đủ suy đạo hàm của nó; cần quan hệ nhiệt động hoặc các điểm lân cận đáng tin.
Đối chiếu với một nghiệm hữu hạn đơn pha
Xét một bài toán riêng, tất cả lỗ rỗng chứa nước, c_p và c_w giữ không đổi; không có dòng vào và giữ giả thiết đường tải. Khi giảm áp một lượng D, tích phân các định nghĩa cho V_p,1=V_p,0 exp(−c_pD) và ρ_1=ρ_0 exp(−c_wD). Khối lượng còn lại bằng M_1=M_0 exp[−(c_p+c_w)D].
Với D=10 MPa và V_p,0=180.000 m³, khối lượng lấy ra quy về thể tích ở mật độ ban đầu là V_p,0{1−exp[−(c_p+c_w)D]}=1.684,072459 m³. Nếu quy về mật độ cuối, cùng khối lượng đó tương ứng V_p,0[exp(c_wD)−exp(−c_pD)]=1.691,498704 m³.
Xấp xỉ tuyến tính cho 1.692 m³. Hai thể tích hữu hạn khác nhau vì dùng hai mật độ tham chiếu, dù đều diễn tả một khối lượng. Ví dụ này chỉ kiểm sự nhất quán toán học của hệ số hằng; không khẳng định hệ số thực giữ nguyên trên 10 MPa hay dùng cùng nghiệm cho hệ ba pha có khí thoát ra.
Nối với hệ số chứa và kiểm đơn vị
Với bão hòa nước hoàn toàn, iTOUGH2 nêu S_s=φρg(c_w+c_p). [3] Ở φ=0,18, ρ=1.050 kg/m³ và các hệ số trên, S_s=1,7428446×10⁻⁶ m⁻¹. Đơn vị này khác Pa⁻¹ của c_t vì ρg đã chuyển biến áp suất sang cột áp. Không thay hai hệ số cho nhau trong mô hình mà bỏ hệ số chuyển đổi.
Chẳng hạn c_p=5×10⁻¹⁰ Pa⁻¹ tương đương 3,447379×10⁻⁶ psi⁻¹. Khi đổi từ psi⁻¹ về Pa⁻¹ phải chia 6.894,757293, không nhân. Một bảng nhập liệu tốt cần lưu định nghĩa V_p, đường ứng suất, đơn vị, nhiệt độ, thành phần, trạng thái pha và phạm vi áp suất. Khi có chuyển pha hoặc nén không tuyến tính, cập nhật mô hình tích trữ theo trạng thái mới có ý nghĩa hơn việc giữ một c_t cố định chỉ vì nó khớp tại điểm đầu.
Nguồn tham khảo
[1] compressibility — Energy Glossary — SLB. Truy cập ngày 09/11/2026.
[2] Formation Properties — WellTest 2019.1, Total Compressibility — IHS Markit Ltd., revised November 9, 2021. Truy cập ngày 09/11/2026.
[3] Manual Page for Command >> COMPRESSIBILITY — Lawrence Berkeley National Laboratory, iTOUGH2. Truy cập ngày 09/11/2026.
Rà soát nội dung: 09/11/2026.
Trong mô hình áp suất vỉa, độ nén tổng thường nằm trong một ô nhỏ nhưng quyết định lượng tích trữ thay đổi khi áp suất giảm. Việc cộng số của dầu, nước, khí và đá chỉ có ý nghĩa khi các hệ số cùng đơn vị, cùng điều kiện và cùng định nghĩa thể tích. Độ nén của vật liệu hạt rắn không tự bằng độ nén thể tích lỗ rỗng.
SLB mô tả compressibility là thay đổi tương đối thể tích theo một đơn vị thay đổi áp suất. [1] Bài này dùng định nghĩa vi phân để tự dựng mô hình lưu trữ cục bộ, với các số liệu giả định. Không lấy kết quả làm hệ số đã hiệu chuẩn cho một vỉa thực hoặc suy trực tiếp dự báo thu hồi.
Độ nén pha và dấu của phần lỗ rỗng
Với một lượng lưu chất cố định, đẳng nhiệt và không đổi thành phần, định nghĩa c_i=−(1/V_i)(∂V_i/∂p). Hệ số dương nghĩa là tăng áp suất làm thể tích lưu chất giảm. Tương đương, c_i=(1/ρ_i)(∂ρ_i/∂p) vì khối lượng không đổi.
Đối với không gian chứa, chọn c_p=(1/V_p)(∂V_p/∂p) theo một đường ứng suất ngoài đã xác định. Trong ví dụ ta giữ ứng suất bao ngoài cố định và lấy c_p dương: tăng áp suất lỗ rỗng làm thể tích lỗ rỗng tăng; giảm áp suất làm nó giảm. Tài liệu iTOUGH2 xác định rõ tham số là độ nén không gian lỗ rỗng. [3]
Hai dấu khác nhau không mâu thuẫn. Khi áp suất giảm, lưu chất muốn nở nhưng chỗ chứa co lại; cả hai góp vào lượng có thể được đẩy ra. Nếu thí nghiệm báo hệ số theo ứng suất hữu hiệu, cần xét phép đổi biến và đường tải thay vì chép nguyên dấu. Nếu thể tích khối đá cũng thay đổi, không được thay V_p bằng độ rỗng φ mà bỏ sự thay đổi của thể tích khối.
Cộng có trọng số tại một trạng thái tham chiếu
Gọi S_o, S_w, S_g là tỷ phần thể tích dầu, nước và khí trong tổng thể tích lỗ rỗng, với S_o+S_w+S_g=1. Trong mô hình cục bộ mà các pha chịu cùng gia số áp suất, không có chuyển pha hay trao đổi thành phần, độ nén tổng là:
c_t = c_p + S_o c_o + S_w c_w + S_g c_g.
WellTest trình bày cấu trúc tổng gồm phần thành hệ và phần lưu chất có trọng số bão hòa. [2] c_t và từng c có đơn vị Pa⁻¹ nếu áp suất dùng Pa. Các S là tỷ phần từ 0 đến 1, không phải số phần trăm từ 0 đến 100. c_p đã quy về V_p nên được cộng một lần, không nhân thêm độ rỗng.
Xét giảm áp nhỏ Δp_dương=p_đầu−p_cuối. Lượng thể tích giải phóng xấp xỉ là ΔV≈V_p c_t Δp_dương. Phần nở lưu chất là V_pΣS_i c_i Δp_dương; phần giảm chỗ chứa là V_p c_p Δp_dương. Đây là khai triển bậc nhất quanh trạng thái đã chọn, chưa mô tả những thay đổi lớn của bão hòa, thành phần và hệ số trong một quá trình dài.
Ví dụ dầu–nước chưa có khí tự do
Cho V_khối=1.000.000 m³, φ=0,18, nên V_p=180.000 m³. Chọn S_w=0,25 và S_o=0,75; c_w=4,4×10⁻¹⁰ Pa⁻¹, c_o=1,0×10⁻⁹ Pa⁻¹, c_p=5,0×10⁻¹⁰ Pa⁻¹. Mọi giá trị là giả định của bài.
Khi đó c_t=1,36×10⁻⁹ Pa⁻¹. Với giảm áp 1 MPa, ΔV≈244,8 m³. Có thể kiểm riêng ba khoản: dầu nở 135 m³, nước nở 19,8 m³ và chỗ chứa giảm 90 m³. Tổng đúng 244,8 m³.
Nếu quên phần lỗ rỗng, kết quả còn 154,8 m³. Nếu nhân c_p thêm φ lần nữa, phần đá chỉ còn 16,2 m³, tổng thành 171 m³. Cả hai sai khác đều xuất phát từ định nghĩa cơ sở thể tích, không phải sai số làm tròn. Độ rỗng đã xuất hiện trong V_p=φV_khối.
Một ít khí có thể thay đổi tổng đáng kể
Ở trạng thái giả định thứ hai, giữ S_w=0,25, đưa S_g=0,02 và giảm S_o còn 0,73. Không giữ S_o=0,75 rồi thêm khí vì như vậy tổng bão hòa vượt một. Giả định khí lý tưởng tại p=20 MPa tuyệt đối, c_g=1/p=5,0×10⁻⁸ Pa⁻¹.
Tổng mới c_t=2,34×10⁻⁹ Pa⁻¹, lớn hơn trường hợp không khí 72,0588%. Lượng giải phóng bậc nhất cho cùng 1 MPa là 421,2 m³. Riêng phần S_g c_g chiếm 42,7350% c_t: nó rất đáng kể nhưng chưa chiếm quá một nửa. Vì thế cần tính trọng số, không thay phép tính bằng phát biểu rằng cứ có khí là mọi phần khác đều bỏ được.
Với S_w cố định, viết c_t=c_t,0+(c_g−c_o)S_g. Đạo hàm theo S_g bằng c_g−c_o=4,9×10⁻⁸ Pa⁻¹ trên một đơn vị bão hòa. Phần khí vượt tổng các phần còn lại khi S_g>c_t,0/(c_g+c_o)=0,026667 trong chính bộ hệ số này. Ngưỡng 2,6667% là kết quả của ví dụ, không phải ranh giới dùng chung cho vỉa.
Độ nén khí không phải hệ số Z
Z là hệ số không thứ nguyên trong pV=ZnRT. Còn c_g mang đơn vị nghịch đảo áp suất. Giữ n và T, lấy log rồi đạo hàm cho c_g=1/p−(1/Z)(dZ/dp). Do đó chỉ khi Z không đổi theo p mới được dùng chính xác 1/p trong mô hình đó.
Chọn riêng Z=0,90 và dZ/dp=0,003 MPa⁻¹ tại 20 MPa. Đổi đạo hàm thành 3×10⁻⁹ Pa⁻¹, ta được c_g=4,666667×10⁻⁸ Pa⁻¹. Thay vào trạng thái S_g=0,02 cho c_t=2,273333×10⁻⁹ Pa⁻¹. Chỉ biết một giá trị Z tại một điểm chưa đủ suy đạo hàm của nó; cần quan hệ nhiệt động hoặc các điểm lân cận đáng tin.
Đối chiếu với một nghiệm hữu hạn đơn pha
Xét một bài toán riêng, tất cả lỗ rỗng chứa nước, c_p và c_w giữ không đổi; không có dòng vào và giữ giả thiết đường tải. Khi giảm áp một lượng D, tích phân các định nghĩa cho V_p,1=V_p,0 exp(−c_pD) và ρ_1=ρ_0 exp(−c_wD). Khối lượng còn lại bằng M_1=M_0 exp[−(c_p+c_w)D].
Với D=10 MPa và V_p,0=180.000 m³, khối lượng lấy ra quy về thể tích ở mật độ ban đầu là V_p,0{1−exp[−(c_p+c_w)D]}=1.684,072459 m³. Nếu quy về mật độ cuối, cùng khối lượng đó tương ứng V_p,0[exp(c_wD)−exp(−c_pD)]=1.691,498704 m³.
Xấp xỉ tuyến tính cho 1.692 m³. Hai thể tích hữu hạn khác nhau vì dùng hai mật độ tham chiếu, dù đều diễn tả một khối lượng. Ví dụ này chỉ kiểm sự nhất quán toán học của hệ số hằng; không khẳng định hệ số thực giữ nguyên trên 10 MPa hay dùng cùng nghiệm cho hệ ba pha có khí thoát ra.
Nối với hệ số chứa và kiểm đơn vị
Với bão hòa nước hoàn toàn, iTOUGH2 nêu S_s=φρg(c_w+c_p). [3] Ở φ=0,18, ρ=1.050 kg/m³ và các hệ số trên, S_s=1,7428446×10⁻⁶ m⁻¹. Đơn vị này khác Pa⁻¹ của c_t vì ρg đã chuyển biến áp suất sang cột áp. Không thay hai hệ số cho nhau trong mô hình mà bỏ hệ số chuyển đổi.
Chẳng hạn c_p=5×10⁻¹⁰ Pa⁻¹ tương đương 3,447379×10⁻⁶ psi⁻¹. Khi đổi từ psi⁻¹ về Pa⁻¹ phải chia 6.894,757293, không nhân. Một bảng nhập liệu tốt cần lưu định nghĩa V_p, đường ứng suất, đơn vị, nhiệt độ, thành phần, trạng thái pha và phạm vi áp suất. Khi có chuyển pha hoặc nén không tuyến tính, cập nhật mô hình tích trữ theo trạng thái mới có ý nghĩa hơn việc giữ một c_t cố định chỉ vì nó khớp tại điểm đầu.
Nguồn tham khảo
[1] compressibility — Energy Glossary — SLB. Truy cập ngày 09/11/2026.
[2] Formation Properties — WellTest 2019.1, Total Compressibility — IHS Markit Ltd., revised November 9, 2021. Truy cập ngày 09/11/2026.
[3] Manual Page for Command >> COMPRESSIBILITY — Lawrence Berkeley National Laboratory, iTOUGH2. Truy cập ngày 09/11/2026.
Rà soát nội dung: 09/11/2026.
