1. Tỷ số Horner không phải một số giờ
Sau khi đóng giếng, áp suất có thể phục hồi theo một đường cong; đổi trục thời gian đúng giúp nhận ra một đoạn xuyên tâm có ý nghĩa. Tài liệu IHS định nghĩa Horner là trường hợp riêng của chồng chất theo thời gian cho một giai đoạn lưu lượng không đổi rồi đóng giếng [2]. Vì vậy thao tác vẽ một đường thẳng không thể thay cho việc xác nhận lịch sử lưu lượng và chế độ dòng.
Gọi tₚ là thời gian khai thác trước đóng, Δt là thời gian kể từ đóng giếng. Tỷ số H = (tₚ+Δt)/Δt không có đơn vị. Nếu tₚ nhập theo giờ thì Δt cũng phải theo giờ; có thể đổi cả hai sang giây mà H giữ nguyên. H giảm từ rất lớn về 1 khi thời gian đóng tăng. Mốc H = 1 tương ứng giới hạn Δt tiến tới vô hạn, không phải một giờ sau đóng.
Thời gian tương đương Agarwal là một đại lượng khác: tₑ = tₚΔt/(tₚ+Δt) = tₚ/H, có cùng đơn vị với tₚ. Whitson trình bày phép biến đổi này để liên hệ dữ liệu phục hồi với dạng phân tích giảm áp [1]. Việc gọi chung cả H và tₑ là “thời gian tương đương” mà bỏ công thức dễ khiến bảng dữ liệu dùng nhầm trục.
2. Chồng chất làm xuất hiện logarithm của tỷ số
Xét chất lỏng một pha hơi nén được, vỉa đồng nhất, dòng xuyên tâm và chưa chịu ảnh hưởng biên. Dùng xấp xỉ logarithm trong miền thời gian thích hợp, đáp ứng giảm áp do lưu lượng tại vỉa qᵣ bắt đầu lúc 0 có dạng A[ln(t/t₀)+C], với A = qᵣμ/(4πkh). Ở đây qᵣ dùng m³/s, μ dùng Pa·s, k dùng m², h dùng m; A có đơn vị Pa. t₀ chỉ là mốc để đối số logarithm không thứ nguyên, còn C gom các tham số cố định trong mô hình.
Đóng giếng ở tₚ được biểu diễn bằng bước lưu lượng −qᵣ chồng lên bước +qᵣ ban đầu [1]. Tại t = tₚ+Δt, hai thành phần hằng triệt tiêu và pᵢ−p_ws = A ln[(tₚ+Δt)/Δt]. Vì vậy p_ws = pᵢ−M log₁₀H, trong đó M = A ln10 là độ lớn độ dốc theo một thập kỷ logarithm. Độ dốc có dấu của đường p theo log₁₀H là −M, còn p tăng theo thời gian đóng.
Sự triệt tiêu trên thuộc mô hình lý tưởng, không nói mọi hiệu ứng gần giếng biến mất ngay sau đóng van. Tích trữ trong giếng và dòng tiếp tục ở đáy có thể làm đoạn sớm lệch khỏi biểu thức. Việc dùng đúng trục Horner cũng không biến một giai đoạn dòng tuyến tính hoặc song tuyến tính thành dòng xuyên tâm; IHS nêu rõ cần hàm thời gian phù hợp với chế độ đó [2].
3. Ví dụ giả định có thể tính ngược
Chọn lưu lượng dầu ở điều kiện báo cáo qₛ = 500 m³/ngày, hệ số thể tích B = 1,20, μ = 1,5×10⁻³ Pa·s, h = 25 m và k = 8×10⁻¹⁴ m². Khi đó qᵣ = qₛB/86.400 = 0,006944444 m³/s. Giả sử pᵢ = 28 MPa và khai thác hằng tₚ = 72 giờ trước đóng. Tính được A = 414.465,997635 Pa và M = 0,954343228 MPa trên một thập kỷ.
Sau 1 giờ đóng, H = 73 và p_ws = 26,221750 MPa. Sau 4 giờ, H = 19 và p_ws = 26,779630 MPa. Sau 12 giờ, H = 7 và p_ws = 27,193486 MPa. Sau 24 giờ, H = 4 và p_ws = 27,425428 MPa. Các điểm này được tạo từ mô hình, không phải dữ liệu của một giếng thực tế.
Từ độ dốc đo trên đồ thị đúng, k = qᵣμ ln10/(4πhM), với M phải đổi về Pa. Thay số trả lại 8×10⁻¹⁴ m². Nếu lấy lưu lượng bề mặt mà quên B, độ thấm suy ra giảm theo cùng tỷ lệ; nếu dùng M tính bằng MPa trong một công thức SI đòi Pa, sai số lên tới hệ số một triệu. Không cần ghi nhớ một hệ số chuyển đổi gộp khi từng đại lượng và đơn vị đã rõ.
Giao điểm ngoại suy tại log₁₀H = 0 cho pᵢ = 28 MPa trong hệ vô hạn đã giả định. Ở vỉa hữu hạn đã suy giảm, một giao điểm Horner không tự động bằng áp suất trung bình hiện thời. Tên và ý nghĩa của áp suất ngoại suy cần gắn với mô hình biên, lịch sử sản xuất và thời điểm khảo sát.
4. Vì sao không thay H bằng thời gian đóng?
Với tₚ = 72 giờ, tₑ ở Δt bằng 1; 4; 12; 24; 72 giờ lần lượt là 0,986301; 3,789474; 10,285714; 18; 36 giờ. Khi Δt rất nhỏ so với tₚ, tₑ gần Δt. Khi đóng rất lâu, tₑ tiến tới 72 giờ và không tăng vô hạn. Giới hạn này xuất phát từ lịch sử khai thác hữu hạn, không phải đồng hồ vật lý ngừng chạy.
Do H = tₚ/tₑ, phương trình có thể viết p_ws = pᵢ−M log₁₀tₚ+M log₁₀tₑ, nếu hai logarithm được hiểu theo cùng đơn vị tham chiếu. Độ dốc theo log₁₀tₑ là +M. Ngược lại, độ dốc cục bộ theo log₁₀Δt bằng M·tₚ/(tₚ+Δt), giảm dần theo Δt. Trong ví dụ, nó bằng 0,904115 tại 4 giờ và 0,715757 MPa mỗi thập kỷ tại 24 giờ.
Nếu hồi quy ba điểm 4, 12 và 24 giờ theo log₁₀Δt thay cho trục đúng, độ dốc chỉ còn 0,833323 MPa mỗi thập kỷ, thấp hơn 12,681019%. Đường vẫn có thể nhìn khá thẳng trên một cửa sổ hẹp. Một hệ số khớp tốt vì vậy chưa đủ xác minh phép biến đổi thời gian; cần kiểm phương trình đã dùng và các đoạn chế độ dòng.
5. Lịch sử nhiều mức lưu lượng cần được giữ
Hướng dẫn Saphir mô tả cách lấy cộng dồn chia lưu lượng cuối để tạo thời gian sản xuất quy đổi, đồng thời nhắc giới hạn của Horner chuẩn [3]. Phép quy đổi này giữ một đại lượng thể tích, nhưng không bảo toàn toàn bộ đáp ứng logarithm của mọi lịch sử. Có thể thấy điều đó bằng một ví dụ mới: sản xuất 300 m³/ngày trong 48 giờ, rồi 600 m³/ngày trong 24 giờ và đóng giếng.
Tổng thể tích báo cáo là 1.200 m³; chia cho 600 m³/ngày cho t_c = 2 ngày = 48 giờ. Tuy nhiên, với cùng tính chất của ví dụ trước, chồng chất chính xác trong xấp xỉ logarithm cho pᵢ−p_ws = A₆₀₀S, trong đó A₆₀₀ = 0,497359197 MPa và S = 0,5ln[(72+Δt)/Δt] + 0,5ln[(24+Δt)/Δt]. Hai hạng 0,5 đến từ bước tăng 300 lúc bắt đầu và bước tăng thêm 300 sau 48 giờ, chuẩn hóa bằng lưu lượng cuối 600.
Tại Δt = 4 giờ, S = 2,445174564, tương đương H_eff = exp(S) = 11,532563. Dùng thời gian quy đổi 48 giờ lại cho H = 13. Áp suất từ chồng chất bằng 26,783870 MPa; phép thay thế cho 26,724299 MPa, thấp hơn 0,059571 MPa. Ở 24 giờ, sai khác còn 0,029290 MPa. Không có một thời gian quy đổi duy nhất tạo đúng cả đường cong trong ví dụ này.
Các logarithm nên dùng tỷ số như trên để đơn vị tự triệt tiêu. Với lịch sử nhiều bước hơn, giữ thời điểm và độ lớn từng thay đổi lưu lượng rồi chồng chất theo mô hình đã chọn. Không nên tùy ý gộp những khoảng ngừng, thay choke và khởi động lại chỉ để đồ thị trông thẳng hơn.
6. Độ nhạy và các dấu kiểm cuối
Vì k tỷ lệ nghịch với M, tăng độ dốc 5% trong khi giữ các đầu vào khác làm k thành 7,619048×10⁻¹⁴ m²; giảm độ dốc 5% cho 8,421053×10⁻¹⁴ m². Khoảng không đối xứng là hệ quả của phép nghịch đảo. Sai lệch áp suất cộng cố định 0,20 MPa làm giao điểm tăng 0,20 MPa nhưng không đổi độ dốc; một sai lệch trôi theo thời gian thì có thể tác động cả hai.
Từ dp_ws/dΔt = A tₚ/[Δt(tₚ+Δt)], lệch mốc thời gian 30 giây quanh Δt = 4 giờ tạo thay đổi bậc nhất khoảng 0,000818 MPa trong ví dụ. Ở sát thời điểm đóng, độ nhạy thời gian lớn hơn, đồng thời ảnh hưởng tích trữ thường cần được xem xét riêng. Ghi chính xác thời điểm sự kiện và đồng bộ đồng hồ là một phần của dữ liệu, không chỉ là định dạng trình bày.
Một kết quả có thể rà soát cần lưu lịch sử lưu lượng, độ sâu quy chiếu áp suất, cơ sở thể tích, đoạn chọn để khớp và lý do nhận dạng dòng xuyên tâm. Tỷ số Horner và thời gian Agarwal là phép biến đổi phục vụ mô hình; chúng không tự xác nhận chất lượng phép thử, và không thay cho việc kiểm những cơ chế khiến mô hình không còn phù hợp.
Nguồn tham khảo
[1] PTA — Pressure Transient Analysis, sections 6.1–6.3 — Whitson. Truy cập ngày 09/11/2026.
[2] Horner Time — WellTest 2019.1, revised 9 November 2021 — IHS Markit. Truy cập ngày 09/11/2026.
[3] Saphir Guided Session #1 — G01 Horner plot, Ecrin v4.20, Doc v4.20.01 — KAPPA. Truy cập ngày 09/11/2026.
Rà soát nội dung: 09/11/2026.
Sau khi đóng giếng, áp suất có thể phục hồi theo một đường cong; đổi trục thời gian đúng giúp nhận ra một đoạn xuyên tâm có ý nghĩa. Tài liệu IHS định nghĩa Horner là trường hợp riêng của chồng chất theo thời gian cho một giai đoạn lưu lượng không đổi rồi đóng giếng [2]. Vì vậy thao tác vẽ một đường thẳng không thể thay cho việc xác nhận lịch sử lưu lượng và chế độ dòng.
Gọi tₚ là thời gian khai thác trước đóng, Δt là thời gian kể từ đóng giếng. Tỷ số H = (tₚ+Δt)/Δt không có đơn vị. Nếu tₚ nhập theo giờ thì Δt cũng phải theo giờ; có thể đổi cả hai sang giây mà H giữ nguyên. H giảm từ rất lớn về 1 khi thời gian đóng tăng. Mốc H = 1 tương ứng giới hạn Δt tiến tới vô hạn, không phải một giờ sau đóng.
Thời gian tương đương Agarwal là một đại lượng khác: tₑ = tₚΔt/(tₚ+Δt) = tₚ/H, có cùng đơn vị với tₚ. Whitson trình bày phép biến đổi này để liên hệ dữ liệu phục hồi với dạng phân tích giảm áp [1]. Việc gọi chung cả H và tₑ là “thời gian tương đương” mà bỏ công thức dễ khiến bảng dữ liệu dùng nhầm trục.
2. Chồng chất làm xuất hiện logarithm của tỷ số
Xét chất lỏng một pha hơi nén được, vỉa đồng nhất, dòng xuyên tâm và chưa chịu ảnh hưởng biên. Dùng xấp xỉ logarithm trong miền thời gian thích hợp, đáp ứng giảm áp do lưu lượng tại vỉa qᵣ bắt đầu lúc 0 có dạng A[ln(t/t₀)+C], với A = qᵣμ/(4πkh). Ở đây qᵣ dùng m³/s, μ dùng Pa·s, k dùng m², h dùng m; A có đơn vị Pa. t₀ chỉ là mốc để đối số logarithm không thứ nguyên, còn C gom các tham số cố định trong mô hình.
Đóng giếng ở tₚ được biểu diễn bằng bước lưu lượng −qᵣ chồng lên bước +qᵣ ban đầu [1]. Tại t = tₚ+Δt, hai thành phần hằng triệt tiêu và pᵢ−p_ws = A ln[(tₚ+Δt)/Δt]. Vì vậy p_ws = pᵢ−M log₁₀H, trong đó M = A ln10 là độ lớn độ dốc theo một thập kỷ logarithm. Độ dốc có dấu của đường p theo log₁₀H là −M, còn p tăng theo thời gian đóng.
Sự triệt tiêu trên thuộc mô hình lý tưởng, không nói mọi hiệu ứng gần giếng biến mất ngay sau đóng van. Tích trữ trong giếng và dòng tiếp tục ở đáy có thể làm đoạn sớm lệch khỏi biểu thức. Việc dùng đúng trục Horner cũng không biến một giai đoạn dòng tuyến tính hoặc song tuyến tính thành dòng xuyên tâm; IHS nêu rõ cần hàm thời gian phù hợp với chế độ đó [2].
3. Ví dụ giả định có thể tính ngược
Chọn lưu lượng dầu ở điều kiện báo cáo qₛ = 500 m³/ngày, hệ số thể tích B = 1,20, μ = 1,5×10⁻³ Pa·s, h = 25 m và k = 8×10⁻¹⁴ m². Khi đó qᵣ = qₛB/86.400 = 0,006944444 m³/s. Giả sử pᵢ = 28 MPa và khai thác hằng tₚ = 72 giờ trước đóng. Tính được A = 414.465,997635 Pa và M = 0,954343228 MPa trên một thập kỷ.
Sau 1 giờ đóng, H = 73 và p_ws = 26,221750 MPa. Sau 4 giờ, H = 19 và p_ws = 26,779630 MPa. Sau 12 giờ, H = 7 và p_ws = 27,193486 MPa. Sau 24 giờ, H = 4 và p_ws = 27,425428 MPa. Các điểm này được tạo từ mô hình, không phải dữ liệu của một giếng thực tế.
Từ độ dốc đo trên đồ thị đúng, k = qᵣμ ln10/(4πhM), với M phải đổi về Pa. Thay số trả lại 8×10⁻¹⁴ m². Nếu lấy lưu lượng bề mặt mà quên B, độ thấm suy ra giảm theo cùng tỷ lệ; nếu dùng M tính bằng MPa trong một công thức SI đòi Pa, sai số lên tới hệ số một triệu. Không cần ghi nhớ một hệ số chuyển đổi gộp khi từng đại lượng và đơn vị đã rõ.
Giao điểm ngoại suy tại log₁₀H = 0 cho pᵢ = 28 MPa trong hệ vô hạn đã giả định. Ở vỉa hữu hạn đã suy giảm, một giao điểm Horner không tự động bằng áp suất trung bình hiện thời. Tên và ý nghĩa của áp suất ngoại suy cần gắn với mô hình biên, lịch sử sản xuất và thời điểm khảo sát.
4. Vì sao không thay H bằng thời gian đóng?
Với tₚ = 72 giờ, tₑ ở Δt bằng 1; 4; 12; 24; 72 giờ lần lượt là 0,986301; 3,789474; 10,285714; 18; 36 giờ. Khi Δt rất nhỏ so với tₚ, tₑ gần Δt. Khi đóng rất lâu, tₑ tiến tới 72 giờ và không tăng vô hạn. Giới hạn này xuất phát từ lịch sử khai thác hữu hạn, không phải đồng hồ vật lý ngừng chạy.
Do H = tₚ/tₑ, phương trình có thể viết p_ws = pᵢ−M log₁₀tₚ+M log₁₀tₑ, nếu hai logarithm được hiểu theo cùng đơn vị tham chiếu. Độ dốc theo log₁₀tₑ là +M. Ngược lại, độ dốc cục bộ theo log₁₀Δt bằng M·tₚ/(tₚ+Δt), giảm dần theo Δt. Trong ví dụ, nó bằng 0,904115 tại 4 giờ và 0,715757 MPa mỗi thập kỷ tại 24 giờ.
Nếu hồi quy ba điểm 4, 12 và 24 giờ theo log₁₀Δt thay cho trục đúng, độ dốc chỉ còn 0,833323 MPa mỗi thập kỷ, thấp hơn 12,681019%. Đường vẫn có thể nhìn khá thẳng trên một cửa sổ hẹp. Một hệ số khớp tốt vì vậy chưa đủ xác minh phép biến đổi thời gian; cần kiểm phương trình đã dùng và các đoạn chế độ dòng.
5. Lịch sử nhiều mức lưu lượng cần được giữ
Hướng dẫn Saphir mô tả cách lấy cộng dồn chia lưu lượng cuối để tạo thời gian sản xuất quy đổi, đồng thời nhắc giới hạn của Horner chuẩn [3]. Phép quy đổi này giữ một đại lượng thể tích, nhưng không bảo toàn toàn bộ đáp ứng logarithm của mọi lịch sử. Có thể thấy điều đó bằng một ví dụ mới: sản xuất 300 m³/ngày trong 48 giờ, rồi 600 m³/ngày trong 24 giờ và đóng giếng.
Tổng thể tích báo cáo là 1.200 m³; chia cho 600 m³/ngày cho t_c = 2 ngày = 48 giờ. Tuy nhiên, với cùng tính chất của ví dụ trước, chồng chất chính xác trong xấp xỉ logarithm cho pᵢ−p_ws = A₆₀₀S, trong đó A₆₀₀ = 0,497359197 MPa và S = 0,5ln[(72+Δt)/Δt] + 0,5ln[(24+Δt)/Δt]. Hai hạng 0,5 đến từ bước tăng 300 lúc bắt đầu và bước tăng thêm 300 sau 48 giờ, chuẩn hóa bằng lưu lượng cuối 600.
Tại Δt = 4 giờ, S = 2,445174564, tương đương H_eff = exp(S) = 11,532563. Dùng thời gian quy đổi 48 giờ lại cho H = 13. Áp suất từ chồng chất bằng 26,783870 MPa; phép thay thế cho 26,724299 MPa, thấp hơn 0,059571 MPa. Ở 24 giờ, sai khác còn 0,029290 MPa. Không có một thời gian quy đổi duy nhất tạo đúng cả đường cong trong ví dụ này.
Các logarithm nên dùng tỷ số như trên để đơn vị tự triệt tiêu. Với lịch sử nhiều bước hơn, giữ thời điểm và độ lớn từng thay đổi lưu lượng rồi chồng chất theo mô hình đã chọn. Không nên tùy ý gộp những khoảng ngừng, thay choke và khởi động lại chỉ để đồ thị trông thẳng hơn.
6. Độ nhạy và các dấu kiểm cuối
Vì k tỷ lệ nghịch với M, tăng độ dốc 5% trong khi giữ các đầu vào khác làm k thành 7,619048×10⁻¹⁴ m²; giảm độ dốc 5% cho 8,421053×10⁻¹⁴ m². Khoảng không đối xứng là hệ quả của phép nghịch đảo. Sai lệch áp suất cộng cố định 0,20 MPa làm giao điểm tăng 0,20 MPa nhưng không đổi độ dốc; một sai lệch trôi theo thời gian thì có thể tác động cả hai.
Từ dp_ws/dΔt = A tₚ/[Δt(tₚ+Δt)], lệch mốc thời gian 30 giây quanh Δt = 4 giờ tạo thay đổi bậc nhất khoảng 0,000818 MPa trong ví dụ. Ở sát thời điểm đóng, độ nhạy thời gian lớn hơn, đồng thời ảnh hưởng tích trữ thường cần được xem xét riêng. Ghi chính xác thời điểm sự kiện và đồng bộ đồng hồ là một phần của dữ liệu, không chỉ là định dạng trình bày.
Một kết quả có thể rà soát cần lưu lịch sử lưu lượng, độ sâu quy chiếu áp suất, cơ sở thể tích, đoạn chọn để khớp và lý do nhận dạng dòng xuyên tâm. Tỷ số Horner và thời gian Agarwal là phép biến đổi phục vụ mô hình; chúng không tự xác nhận chất lượng phép thử, và không thay cho việc kiểm những cơ chế khiến mô hình không còn phù hợp.
Nguồn tham khảo
[1] PTA — Pressure Transient Analysis, sections 6.1–6.3 — Whitson. Truy cập ngày 09/11/2026.
[2] Horner Time — WellTest 2019.1, revised 9 November 2021 — IHS Markit. Truy cập ngày 09/11/2026.
[3] Saphir Guided Session #1 — G01 Horner plot, Ecrin v4.20, Doc v4.20.01 — KAPPA. Truy cập ngày 09/11/2026.
Rà soát nội dung: 09/11/2026.
