Suy giảm mũ: dự báo lưu lượng với giả thiết tốc độ suy giảm cố định

1. Tốc độ suy giảm cố định có hai cách diễn đạt

“Giảm 18% mỗi năm” thường được hiểu là sau một năm còn 82% lưu lượng đầu năm. Đưa trực tiếp số 0,18 vào số mũ lại tạo một mô hình khác. Tài liệu Whitson phân biệt tốc độ suy giảm danh nghĩa với mức suy giảm hữu hiệu trên một khoảng thời gian [1]. Sự phân biệt này cần được giải quyết trước khi khớp dữ liệu, vì sai khác sẽ tích lũy trong dự báo nhiều năm.

Định nghĩa D = −(1/q)dq/dt, với q là lưu lượng và D có đơn vị nghịch đảo thời gian. Nếu D không đổi thì dq/dt = −Dq, suy ra q(t) = qᵢ exp(−Dt). Đây là dạng mũ được trình bày trong tài liệu ResInsight [2]. Từ nghiệm có thể thấy cùng một khoảng Δt sẽ có cùng tỷ số q(t+Δt)/q(t), dù lượng giảm tuyệt đối nhỏ dần. Mô hình không nói giếng mất cùng một số m³/ngày sau mỗi năm.

Mức suy giảm hữu hiệu trên khoảng Δt là d = 1−exp(−DΔt), nên D = −ln(1−d)/Δt. Cần ghi rõ d là số thập phân hay phần trăm và Δt tính bằng ngày, tháng hay năm. Số mũ Dt phải không thứ nguyên. Việc ghi một cột “decline” mà bỏ loại tốc độ và đơn vị làm cho bảng tính khó kiểm toán dù các ô vẫn cho kết quả trơn đẹp.

2. Ví dụ giả định với 18% hữu hiệu mỗi năm

Chọn qᵢ = 600 m³ dầu ở điều kiện báo cáo mỗi ngày, tại t = 0 của đoạn dự báo. Lấy một năm mô hình bằng đúng 365 ngày và d = 0,18 trên một năm. Khi đó D = −ln(0,82) = 0,198450939 năm⁻¹, tương đương 0,000543701 ngày⁻¹. Sau một năm q = 492 m³/ngày; sau ba năm q = 330,820800 m³/ngày; sau năm năm q = 222,443906 m³/ngày.

Nếu nhầm D = 0,18 năm⁻¹ thì lưu lượng năm thứ ba thành 349,648951 m³/ngày. Sai khác này không đến từ mô hình vỉa phức tạp mà chỉ do đổi nghĩa tham số. Mức hữu hiệu cho một phần mười hai năm là 1−0,82^(1/12) = 0,016401583, khoảng 1,640158%. Lấy 18% chia 12 cho 1,5% sẽ không ghép lại đúng mức giảm năm ban đầu.

Trong hồ sơ dữ liệu, qᵢ phải gắn với thời điểm bắt đầu đoạn đang khớp. Nó có thể là lưu lượng sau một đợt sửa giếng hoặc sau khi kết thúc giai đoạn tăng tải, không bắt buộc là đỉnh lịch sử. Sản lượng đã khai thác trước t = 0 được lưu riêng để cộng vào tổng lịch sử, tránh tính lặp khi nối nhiều đoạn.

3. Tích phân phải giữ hệ số đổi ngày và năm

Quan hệ giữa lưu lượng, thời gian và cộng dồn đã là trọng tâm của công trình gốc Arps [3]. Với cách chọn đơn vị trong ví dụ, sản lượng mới trên 0 đến T năm là N(T) = 365∫q(t)dt = 365qᵢ[1−exp(−DT)]/D, đơn vị m³. Không có hệ số 365 thì kết quả vẫn còn đơn vị m³·năm/ngày, chưa phải một thể tích.

Sau ba năm, N = 495.086,637694 m³. Năm đầu cho 198.638,516167 m³, tương ứng lưu lượng trung bình 544,215113 m³/ngày. Lấy lưu lượng cuối năm nhân 365 cho 179.580 m³ sẽ đánh giá thấp phần đã sản xuất khi lưu lượng còn cao. Lấy lưu lượng đầu năm nhân 365 lại cho 219.000 m³ và đánh giá cao. Với mô hình liên tục đã chọn, tích phân mới là phép tính nhất quán.

Nếu dữ liệu nhập là lưu lượng trung bình tháng, cần so với trung bình mô hình trên đúng khoảng tháng. Công thức trung bình từ t₁ đến t₂ là q̄ = qᵢ[exp(−Dt₁)−exp(−Dt₂)]/[D(t₂−t₁)]. Không nên tự xem mọi số trung bình là lưu lượng tức thời tại cuối tháng. Có thể xấp xỉ bằng giá trị giữa tháng khi biến động nhỏ, nhưng ngày đầu, ngày cuối và số ngày của từng kỳ vẫn phải được giữ trong dữ liệu gốc.

4. Ngưỡng lưu lượng tạo một điểm dừng hữu hạn

Giả sử chỉ nghiên cứu đoạn cho đến q_c = 80 m³/ngày. Đây là ngưỡng giả định để minh họa, chưa được tính từ doanh thu hay chi phí. Giải q_c = qᵢexp(−Dt_c) được t_c = ln(qᵢ/q_c)/D = 10,153154 năm. Sản lượng đến ngưỡng là N_c = 365(qᵢ−q_c)/D = 956.407,670432 m³.

Tích phân tới vô hạn cho 365qᵢ/D, nhưng không thể gọi toàn bộ kết quả đó là trữ lượng thương mại. Lưu lượng dương rất nhỏ vẫn đóng góp toán học dù giếng có thể đã dừng vì nước, thiết bị, thị trường hoặc điều kiện hợp đồng. Một giới hạn thời gian dự án cũng có thể đến trước t_c; khi ấy phải tích phân đến thời điểm sớm hơn. Nếu q_c lớn hơn qᵢ, mô hình không có một giai đoạn tương lai nằm trên ngưỡng để ngoại suy theo cách này.

5. Độ nhạy tham số và chất lượng đoạn khớp

Giữ qᵢ = 600 và q_c = 80, thay mức giảm hữu hiệu năm bằng 16%, 18% và 20% cho sản lượng đến ngưỡng lần lượt 1.088.593,706771; 956.407,670432; 850.573,538344 m³. Lưu lượng năm thứ năm tương ứng 250,927165; 222,443906; 196,608000 m³/ngày. Các kịch bản này thể hiện ảnh hưởng của tham số, không phải các mức xác suất dự báo.

Từ hai điểm dương q₁ tại t₁ và q₂ tại t₂ có thể tính D = ln(q₁/q₂)/(t₂−t₁). Nhưng một cặp điểm không kiểm tra được giả thiết D hằng. Nên xem nhiều khoảng ln(qᵢ/qᵢ₊₁)/Δt, đọc nhật ký điều kiện khai thác và kiểm phần dư theo thời gian. Giá trị q bằng không do đóng giếng không thể đưa thẳng vào logarithm; xóa nó cũng không có nghĩa thời gian đóng đã biến mất khỏi sản lượng thực tế.

Khớp bình phương sai số của q và khớp bình phương sai số của ln q ưu tiên sai lệch khác nhau. Cách thứ nhất nhạy với số tuyệt đối ở giai đoạn cao; cách thứ hai gần với sai lệch tương đối ở giai đoạn thấp. Với dữ liệu nhiễu, hai bộ tham số có thể khác. Cần công bố hàm mục tiêu và cửa sổ khớp, đồng thời xem tổng sản lượng dự báo có tái tạo được tổng đã quan sát không.

6. Thời gian lịch và thời gian hoạt động

Giả sử q(t) ở trên là năng lực khi đang chạy nhưng suy giảm theo thời gian lịch, còn hệ số hoạt động đều U = 0,90. Khi không có hiệu ứng khởi động lại, sản lượng ba năm bằng U×495.086,637694 = 445.577,973925 m³. Đây là một giả thiết vận hành riêng, không phải phép sửa luôn đúng cho mọi bộ dữ liệu.

Trong mô hình khác, suy giảm chỉ tiến theo thời gian hoạt động τ = Ut. Lưu lượng bình quân lịch là Uqᵢexp(−DUt), và N(T) = 365qᵢ[1−exp(−DUT)]/D. Với cùng U và T = 3 năm, kết quả là 457.761,709069 m³. Khác biệt xuất hiện vì mô hình thứ hai giữ lại năng lực trong thời gian nghỉ. Không được thay đổi đồng hồ thời gian mà vẫn dùng nguyên tham số đã khớp theo đồng hồ cũ.

Nếu q đầu vào vốn là sản lượng kỳ chia cho toàn bộ ngày lịch, hệ số hoạt động đã nằm trong số đo; nhân U thêm lần nữa có thể đếm tổn thất hai lần. Vì vậy kết quả cần lưu kèm định nghĩa lưu lượng, lịch dừng, đoạn khớp, quy ước năm và ngưỡng kết thúc. Một đường mũ chỉ trở thành dự báo có thể rà soát khi các quyết định đó được viết rõ bên cạnh công thức.

Nguồn tham khảo
[1] Decline Curve Analysis — sections 2.4 and 2.5 — Whitson. Truy cập ngày 09/11/2026.
[2] Decline Curve Analysis — ResInsight User Documentation — Open Porous Media / ResInsight. Truy cập ngày 09/11/2026.
[3] Analysis of Decline Curves — Transactions of the AIME, 160, 228–247 (1945) — J. J. Arps; AIME, bản lưu tại Texas A&M University. Truy cập ngày 09/11/2026.

Rà soát nội dung: 09/11/2026.
 
Top