1. Đường thẳng bắt đầu từ một hệ kín
Đồ thị p/Z theo khí cộng dồn thường được dùng để ước lượng lượng khí ban đầu. Tuy nhiên, một đường hồi quy đẹp chỉ có ý nghĩa vật lý khi các số đo và giả thiết mô hình cùng phù hợp. Tài liệu Harmony trình bày trường hợp cơ bản bằng cân bằng khí khô, trong đó giãn nở khí là cơ chế cung cấp thể tích chính.[1]
Ở đây xét một khoang khí liên thông, không nhận khí từ ngoài, không bơm ép, không ngưng tụ lỏng và không có khí hấp phụ đáng kể. Thể tích khí chiếm trong lỗ rỗng được giữ không đổi ở phần đầu. Các con số là dữ liệu tổng hợp tự dựng, không đại diện một mỏ hoặc một phép xác nhận trữ lượng.
2. Giữ cơ sở mol khi chuyển trạng thái
Gọi G là lượng khí đầu và G_p là lượng khí đã khai thác, cùng quy về một điều kiện chuẩn p_s,T_s,Z_s. Chúng có đơn vị m³ quy chiếu. Nếu thành phần khí không đổi, tỷ số mol còn lại trên mol đầu bằng (G−G_p)/G.
IADC định nghĩa Z là hệ số dùng để mô tả độ lệch quan hệ khí so với khí lý tưởng.[2] Dùng phương trình pV=nZRT, với p là áp suất tuyệt đối, T theo kelvin, ta có cho cùng thể tích và nhiệt độ:
p/Z = (p_i/Z_i)(1−G_p/G)
Đặt y=p/Z, y_i=p_i/Z_i thì y=y_i−(y_i/G)G_p. Giao điểm trục hoành của đường lý tưởng là G. Đơn vị của y vẫn là đơn vị áp suất; Z không có đơn vị. Dùng áp suất dư sẽ làm sai tỷ số, đặc biệt khi áp suất cuối thấp.
G không phải thể tích lỗ rỗng. Muốn tính thể tích khí đầu ở vỉa phải dùng V_gi=G(p_s/Z_sT_s)(Z_iT_i/p_i), với áp suất hai trạng thái cùng đơn vị.
3. Dữ liệu giả định có Z thay đổi
Chọn G=8 tỷ m³ quy chiếu, p_i=36 MPa và Z_i=1,05. Khi G_p lần lượt bằng 0; 1; 2; 3 tỷ m³, giá trị y theo mô hình là 34,285714; 30; 25,714286; 21,428571 MPa.
Để biểu diễn số đo áp suất, chọn các Z tương ứng 1,05; 0,98; 0,92; 0,88. Đây là bảng giả định cho bài toán, không phải một tương quan khí thực. Nhân y với Z cho p=36; 29,4; 23,657143; 18,857143 MPa.
Nếu dùng x=G_p tính bằng tỷ m³, phương trình đúng là:
y = 34,285714286 − 4,285714286 x
G = −a/b = 8 tỷ m³, khi y=a+bx
Với p_s=0,101325 MPa, T_s=288,15 K, Z_s=1 và T_i=370 K, thể tích khí đầu ở vỉa là 30.358.233,558910 m³. Tỷ số hàng trăm lần giữa khí quy chiếu và khí vỉa là hệ quả chuyển trạng thái, không phải tăng hoặc mất số mol.
4. Bỏ Z tạo một phép ngoại suy khác
Hồi quy trực tiếp p theo x của bốn điểm trên bằng bình phương tối thiểu không trọng số cho a=35,554285714 MPa và b=−5,717142857 MPa trên tỷ m³. Giao điểm suy ra là 6,218890555 tỷ m³, khác rõ với tám tỷ.
Các áp suất đều đúng trong mô hình; sai lệch xuất phát từ chọn biến phụ thuộc. Một đường p có vẻ gần thẳng trên khoảng ngắn vẫn không thay được quan hệ p/Z khi Z biến đổi.
Ngược lại, nếu mọi y đều bị nhân cùng một hệ số không đổi và cả hệ số chặn lẫn độ dốc đều được hồi quy tự do, tỷ số −a/b không đổi. Điều này không làm sai số Z vô hại: sai lệch phụ thuộc áp suất làm đổi độ dốc, còn việc cố định một hệ số chặn từ bộ Z khác phá mất sự triệt tiêu đó.
5. Dữ liệu đầu đời có thể rất nhạy
Khi y_i đã biết, một điểm cho G=G_p/(1−y/y_i). Mẫu số là mức suy giảm tương đối của y. Nếu mới khai thác G_p=0,5 tỷ m³, mức suy giảm chỉ là 0,0625; một lệch nhỏ ở y dễ bị khuếch đại.
Tăng y của điểm này thêm 0,5 MPa làm G suy ra thành 10,434782609 tỷ m³. Cùng lệch y=0,5 MPa tại G_p=3 tỷ m³ chỉ làm G thành 8,323699422 tỷ m³. Đây là hai thử độ nhạy riêng, không phải hai lần cộng dồn sai số.
Lệch p và lệch y cũng khác nhau. Tại điểm cuối có Z=0,88, tăng p thêm 0,5 MPa làm y tăng 0,568181818 MPa và cho G=8,369881110 tỷ m³. Cần truyền sai số qua phép chia trước khi đánh giá ảnh hưởng đến lượng khí.
Đạo hàm cục bộ ∂G/∂y=G_p/[y_i(1−y/y_i)²] giải thích mức khuếch đại. Nó không tự bao gồm sai số y_i, G_p hoặc tương quan giữa các lần khảo sát.
6. Thể tích khí giảm vì nước vào
Từ bảo toàn mol nhưng cho phép thể tích khí đổi, quan hệ tổng quát trong mô hình đồng nhất là:
(p/Z)/(p_i/Z_i) = (1−G_p/G)(T/T_i)(V_gi/V_g)
Giữ nhiệt độ không đổi, tại G_p=3 tỷ m³ giả định nước đã chiếm chỗ làm V_g chỉ còn 95% ban đầu. Khi đó y=21,428571429/0,95=22,556390977 MPa. Bỏ thay đổi thể tích và dùng y_i cố định sẽ suy G biểu kiến 8,769230769 tỷ m³.
Ví dụ này không xây một mô hình tầng nước theo thời gian; nó chỉ cho biết một thể tích khí nhỏ hơn tạo áp suất lớn hơn với cùng mol còn lại. Muốn dự báo diễn biến cần quan hệ dòng nước, kích thước liên thông và lịch sử áp suất tương ứng.
7. Nhiệt độ không đổi cũng là một giả thiết
Xét riêng trường hợp thể tích khí không đổi nhưng T giảm từ 370 xuống 360 K tại cùng mức khai thác ba tỷ m³. Khi đó y=20,849420849 MPa và phương pháp đẳng nhiệt suy G=7,655172414 tỷ m³.
Nếu nhiệt độ đại diện được đo đúng, đồ thị p/(ZT) sẽ giữ quan hệ tuyến tính với G_p cho hệ có thể tích không đổi. Tuy nhiên, một nhiệt độ tại giếng chưa chắc bằng nhiệt độ trung bình của toàn bộ khí chứa. Không nên thêm một phép hiệu chỉnh hình thức bằng số đo không đại diện rồi cho rằng mô hình đã đầy đủ.
Thay đổi thành phần, khí hòa tan trong nước hoặc ngưng tụ hydrocarbon còn làm quan hệ giữa thể tích quy chiếu và mol cần được xem lại. Chúng không được sửa bằng một hệ số Z duy nhất nếu bản thân sổ cấu tử chưa đóng.
8. Áp suất khảo sát phải đại diện cho khoang
Mô-đun của AGI do Bureau of Economic Geology lưu trữ chỉ ra ảnh hưởng của thời gian đóng giếng, giao thoa giếng lân cận và giếng mở nhiều tầng đối với áp suất dùng trong phân tích.[3] Các vấn đề này làm sai lệch ý nghĩa của từng điểm trước cả bước hồi quy.
Nên ghép mỗi áp suất với lượng khí cộng dồn đúng thời điểm, cùng ranh giới các giếng và cùng điều kiện quy chiếu. Nếu thêm một tầng hoặc thay phân bổ sản lượng giữa các khoang, phải đánh giá lại chuỗi dữ liệu. Không ghép nối các điểm thuộc những hệ vật chất khác nhau chỉ vì tên giếng không đổi.
Một đường cong có thể do nước hỗ trợ, biến dạng đá, liên thông muộn hoặc dữ liệu chưa ổn định. Đồ thị đơn lẻ thường không phân biệt duy nhất các nguyên nhân đó; cần đối chiếu lịch sử hoàn thiện giếng và dữ liệu địa chất.
9. Giao điểm hình học và lượng khai thác có điều kiện
Với mô hình cơ sở tám tỷ, nếu đặt một trạng thái cuối giả định p=5 MPa, Z=0,95, lượng đã khai thác đến đó là 6.771.929.824,561403 m³ quy chiếu. Phần còn lại vẫn nằm trong khoang ở áp suất cuối.
Đây chỉ là tính toán theo trạng thái áp suất giả định. Nó chưa xét khả năng nâng, công suất nén, hợp đồng, chi phí hoặc giới hạn vận hành để xác định điểm dừng kinh tế. Vì vậy không gọi giao điểm p/Z=0 là lượng chắc chắn sẽ sản xuất, và không chuyển trực tiếp kết quả mô hình thành phân loại trữ lượng thương mại.
Nguồn tham khảo
[1] Material balance analysis theory — Gas material balance — S&P Global / IHS, Harmony Enterprise documentation. Truy cập ngày 09/11/2026.
[2] Z factor — International Association of Drilling Contractors, IADC UBO / MPD Glossary, December 2011. Truy cập ngày 09/11/2026.
[3] AGI Training Module, Production History Analysis — Analysis of Pressure History — AGI, hosted by Bureau of Economic Geology, University of Texas at Austin, 2000. Truy cập ngày 09/11/2026.
Rà soát nội dung: 09/11/2026.
Đồ thị p/Z theo khí cộng dồn thường được dùng để ước lượng lượng khí ban đầu. Tuy nhiên, một đường hồi quy đẹp chỉ có ý nghĩa vật lý khi các số đo và giả thiết mô hình cùng phù hợp. Tài liệu Harmony trình bày trường hợp cơ bản bằng cân bằng khí khô, trong đó giãn nở khí là cơ chế cung cấp thể tích chính.[1]
Ở đây xét một khoang khí liên thông, không nhận khí từ ngoài, không bơm ép, không ngưng tụ lỏng và không có khí hấp phụ đáng kể. Thể tích khí chiếm trong lỗ rỗng được giữ không đổi ở phần đầu. Các con số là dữ liệu tổng hợp tự dựng, không đại diện một mỏ hoặc một phép xác nhận trữ lượng.
2. Giữ cơ sở mol khi chuyển trạng thái
Gọi G là lượng khí đầu và G_p là lượng khí đã khai thác, cùng quy về một điều kiện chuẩn p_s,T_s,Z_s. Chúng có đơn vị m³ quy chiếu. Nếu thành phần khí không đổi, tỷ số mol còn lại trên mol đầu bằng (G−G_p)/G.
IADC định nghĩa Z là hệ số dùng để mô tả độ lệch quan hệ khí so với khí lý tưởng.[2] Dùng phương trình pV=nZRT, với p là áp suất tuyệt đối, T theo kelvin, ta có cho cùng thể tích và nhiệt độ:
p/Z = (p_i/Z_i)(1−G_p/G)
Đặt y=p/Z, y_i=p_i/Z_i thì y=y_i−(y_i/G)G_p. Giao điểm trục hoành của đường lý tưởng là G. Đơn vị của y vẫn là đơn vị áp suất; Z không có đơn vị. Dùng áp suất dư sẽ làm sai tỷ số, đặc biệt khi áp suất cuối thấp.
G không phải thể tích lỗ rỗng. Muốn tính thể tích khí đầu ở vỉa phải dùng V_gi=G(p_s/Z_sT_s)(Z_iT_i/p_i), với áp suất hai trạng thái cùng đơn vị.
3. Dữ liệu giả định có Z thay đổi
Chọn G=8 tỷ m³ quy chiếu, p_i=36 MPa và Z_i=1,05. Khi G_p lần lượt bằng 0; 1; 2; 3 tỷ m³, giá trị y theo mô hình là 34,285714; 30; 25,714286; 21,428571 MPa.
Để biểu diễn số đo áp suất, chọn các Z tương ứng 1,05; 0,98; 0,92; 0,88. Đây là bảng giả định cho bài toán, không phải một tương quan khí thực. Nhân y với Z cho p=36; 29,4; 23,657143; 18,857143 MPa.
Nếu dùng x=G_p tính bằng tỷ m³, phương trình đúng là:
y = 34,285714286 − 4,285714286 x
G = −a/b = 8 tỷ m³, khi y=a+bx
Với p_s=0,101325 MPa, T_s=288,15 K, Z_s=1 và T_i=370 K, thể tích khí đầu ở vỉa là 30.358.233,558910 m³. Tỷ số hàng trăm lần giữa khí quy chiếu và khí vỉa là hệ quả chuyển trạng thái, không phải tăng hoặc mất số mol.
4. Bỏ Z tạo một phép ngoại suy khác
Hồi quy trực tiếp p theo x của bốn điểm trên bằng bình phương tối thiểu không trọng số cho a=35,554285714 MPa và b=−5,717142857 MPa trên tỷ m³. Giao điểm suy ra là 6,218890555 tỷ m³, khác rõ với tám tỷ.
Các áp suất đều đúng trong mô hình; sai lệch xuất phát từ chọn biến phụ thuộc. Một đường p có vẻ gần thẳng trên khoảng ngắn vẫn không thay được quan hệ p/Z khi Z biến đổi.
Ngược lại, nếu mọi y đều bị nhân cùng một hệ số không đổi và cả hệ số chặn lẫn độ dốc đều được hồi quy tự do, tỷ số −a/b không đổi. Điều này không làm sai số Z vô hại: sai lệch phụ thuộc áp suất làm đổi độ dốc, còn việc cố định một hệ số chặn từ bộ Z khác phá mất sự triệt tiêu đó.
5. Dữ liệu đầu đời có thể rất nhạy
Khi y_i đã biết, một điểm cho G=G_p/(1−y/y_i). Mẫu số là mức suy giảm tương đối của y. Nếu mới khai thác G_p=0,5 tỷ m³, mức suy giảm chỉ là 0,0625; một lệch nhỏ ở y dễ bị khuếch đại.
Tăng y của điểm này thêm 0,5 MPa làm G suy ra thành 10,434782609 tỷ m³. Cùng lệch y=0,5 MPa tại G_p=3 tỷ m³ chỉ làm G thành 8,323699422 tỷ m³. Đây là hai thử độ nhạy riêng, không phải hai lần cộng dồn sai số.
Lệch p và lệch y cũng khác nhau. Tại điểm cuối có Z=0,88, tăng p thêm 0,5 MPa làm y tăng 0,568181818 MPa và cho G=8,369881110 tỷ m³. Cần truyền sai số qua phép chia trước khi đánh giá ảnh hưởng đến lượng khí.
Đạo hàm cục bộ ∂G/∂y=G_p/[y_i(1−y/y_i)²] giải thích mức khuếch đại. Nó không tự bao gồm sai số y_i, G_p hoặc tương quan giữa các lần khảo sát.
6. Thể tích khí giảm vì nước vào
Từ bảo toàn mol nhưng cho phép thể tích khí đổi, quan hệ tổng quát trong mô hình đồng nhất là:
(p/Z)/(p_i/Z_i) = (1−G_p/G)(T/T_i)(V_gi/V_g)
Giữ nhiệt độ không đổi, tại G_p=3 tỷ m³ giả định nước đã chiếm chỗ làm V_g chỉ còn 95% ban đầu. Khi đó y=21,428571429/0,95=22,556390977 MPa. Bỏ thay đổi thể tích và dùng y_i cố định sẽ suy G biểu kiến 8,769230769 tỷ m³.
Ví dụ này không xây một mô hình tầng nước theo thời gian; nó chỉ cho biết một thể tích khí nhỏ hơn tạo áp suất lớn hơn với cùng mol còn lại. Muốn dự báo diễn biến cần quan hệ dòng nước, kích thước liên thông và lịch sử áp suất tương ứng.
7. Nhiệt độ không đổi cũng là một giả thiết
Xét riêng trường hợp thể tích khí không đổi nhưng T giảm từ 370 xuống 360 K tại cùng mức khai thác ba tỷ m³. Khi đó y=20,849420849 MPa và phương pháp đẳng nhiệt suy G=7,655172414 tỷ m³.
Nếu nhiệt độ đại diện được đo đúng, đồ thị p/(ZT) sẽ giữ quan hệ tuyến tính với G_p cho hệ có thể tích không đổi. Tuy nhiên, một nhiệt độ tại giếng chưa chắc bằng nhiệt độ trung bình của toàn bộ khí chứa. Không nên thêm một phép hiệu chỉnh hình thức bằng số đo không đại diện rồi cho rằng mô hình đã đầy đủ.
Thay đổi thành phần, khí hòa tan trong nước hoặc ngưng tụ hydrocarbon còn làm quan hệ giữa thể tích quy chiếu và mol cần được xem lại. Chúng không được sửa bằng một hệ số Z duy nhất nếu bản thân sổ cấu tử chưa đóng.
8. Áp suất khảo sát phải đại diện cho khoang
Mô-đun của AGI do Bureau of Economic Geology lưu trữ chỉ ra ảnh hưởng của thời gian đóng giếng, giao thoa giếng lân cận và giếng mở nhiều tầng đối với áp suất dùng trong phân tích.[3] Các vấn đề này làm sai lệch ý nghĩa của từng điểm trước cả bước hồi quy.
Nên ghép mỗi áp suất với lượng khí cộng dồn đúng thời điểm, cùng ranh giới các giếng và cùng điều kiện quy chiếu. Nếu thêm một tầng hoặc thay phân bổ sản lượng giữa các khoang, phải đánh giá lại chuỗi dữ liệu. Không ghép nối các điểm thuộc những hệ vật chất khác nhau chỉ vì tên giếng không đổi.
Một đường cong có thể do nước hỗ trợ, biến dạng đá, liên thông muộn hoặc dữ liệu chưa ổn định. Đồ thị đơn lẻ thường không phân biệt duy nhất các nguyên nhân đó; cần đối chiếu lịch sử hoàn thiện giếng và dữ liệu địa chất.
9. Giao điểm hình học và lượng khai thác có điều kiện
Với mô hình cơ sở tám tỷ, nếu đặt một trạng thái cuối giả định p=5 MPa, Z=0,95, lượng đã khai thác đến đó là 6.771.929.824,561403 m³ quy chiếu. Phần còn lại vẫn nằm trong khoang ở áp suất cuối.
Đây chỉ là tính toán theo trạng thái áp suất giả định. Nó chưa xét khả năng nâng, công suất nén, hợp đồng, chi phí hoặc giới hạn vận hành để xác định điểm dừng kinh tế. Vì vậy không gọi giao điểm p/Z=0 là lượng chắc chắn sẽ sản xuất, và không chuyển trực tiếp kết quả mô hình thành phân loại trữ lượng thương mại.
Nguồn tham khảo
[1] Material balance analysis theory — Gas material balance — S&P Global / IHS, Harmony Enterprise documentation. Truy cập ngày 09/11/2026.
[2] Z factor — International Association of Drilling Contractors, IADC UBO / MPD Glossary, December 2011. Truy cập ngày 09/11/2026.
[3] AGI Training Module, Production History Analysis — Analysis of Pressure History — AGI, hosted by Bureau of Economic Geology, University of Texas at Austin, 2000. Truy cập ngày 09/11/2026.
Rà soát nội dung: 09/11/2026.
