1. Một cột số chưa cho biết cách tích phân
Dữ liệu sản xuất có thể là lưu lượng tức thời, lưu lượng trung bình kỳ hoặc thể tích đã cộng dồn. Ba loại này đòi hỏi ba cách xử lý khác nhau. Với số tức thời cần mô tả diễn biến giữa các lần đo; với trung bình kỳ chỉ cần nhân đúng độ dài kỳ; với bộ đếm thể tích lấy chênh sau khi xử lý các lần reset. Áp dụng hình thang cho cả ba sẽ tạo một tổng có vẻ hợp lý nhưng không bảo toàn ý nghĩa đo lường.
Giả sử q(t) là lưu lượng dầu đã quy đổi về cùng điều kiện báo cáo, đơn vị m³/ngày, còn dấu thời gian tính bằng giờ. Sản lượng trên một ngày là N = (1/24)∫q(t)dt, đơn vị m³. Hệ số 1/24 không được bỏ khi dữ liệu gốc ghi giờ. Nếu cơ sở nhiệt độ, áp suất thay đổi giữa hai nguồn dữ liệu, cần quy đổi trước khi nối thành một chuỗi duy nhất.
2. Ví dụ hình thang trên lưới không đều
Quy tắc hình thang xấp xỉ bằng đoạn thẳng nối hai điểm; SciPy cho phép cung cấp trực tiếp các vị trí lấy mẫu và lưu ý chúng không được tự sắp xếp [1]. Với hai điểm (tᵢ,qᵢ), (tᵢ₊₁,qᵢ₊₁), phần thể tích là ΔNᵢ = (qᵢ+qᵢ₊₁)(tᵢ₊₁−tᵢ)/48 khi t dùng giờ và q dùng m³/ngày.
Xét dữ liệu giả định tại giờ 0, 2, 5, 12, 18, 24 với q lần lượt 600, 660, 600, 480, 540, 500. Nếu chấp nhận đoạn thẳng giữa mọi cặp, năm phần thể tích lần lượt là 52,50; 78,75; 157,50; 127,50; 130,00 m³. Tổng N = 546,25 m³, và lưu lượng bình quân của đúng 24 giờ này cũng có trị số 546,25 m³/ngày.
Không nên lấy trung bình số học sáu lưu lượng rồi nhân một ngày. Các số đo không đại diện cho những khoảng bằng nhau; điểm ở đầu và cuối còn chỉ đóng góp một phía trong hình thang. Tương tự, dùng khoảng cách mặc định bằng một cho tất cả cặp sẽ biến lưới giờ không đều thành một lưới khác.
3. Khoảng thiếu không tự có nghĩa là giếng đóng
Giả sử nhật ký chất lượng cho biết từ giờ 5 đến 12 không có dữ liệu liên tục; hai số ở biên vẫn hợp lệ. Khi ấy 157,50 m³ là phần ước tính theo nội suy, không phải phần đã được quan sát đủ. Bốn khoảng còn lại dài tổng cộng 17 giờ, đóng góp 388,75 m³ theo phương pháp đang dùng. Cần báo riêng thể tích có dữ liệu và thể tích phụ thuộc giả thiết lấp chỗ trống.
MathWorks cung cấp phương pháp điền theo tọa độ thời gian và trả cờ nhận diện điểm đã được điền [2]. Khả năng đó chỉ thực hiện lựa chọn toán học; nó không xác nhận giếng có chạy hay không. Với cùng hai điểm biên, ta có thể tạo nhiều đường nằm giữa, cho nhiều tích phân khác nhau. Xóa các hàng thiếu rồi nối thẳng hai điểm còn lại vẫn chính là chấp nhận một nội suy qua khoảng thiếu.
Giả sử thông tin vận hành độc lập chỉ ràng buộc q nằm từ 450 đến 650 m³/ngày suốt khoảng 7 giờ. Thể tích khoảng đó nằm từ 131,25 đến 189,583333 m³. Tổng ngày tương ứng từ 520,00 đến 578,333333 m³. Đây là khoảng theo ràng buộc đã giả định, không phải khoảng tin cậy thống kê. Nếu ràng buộc không có bằng chứng, chính khoảng này cũng không nên được trình bày như đã biết.
4. Nhật ký đóng giếng thay đổi cách dựng lại đường
Xét kịch bản thay thế: nhật ký xác nhận đóng từ giờ 8 đến 10. Để minh họa, giả sử lưu lượng giữ 600 từ giờ 5 đến sát giờ 8, chuyển ngay về 0, giữ 0 đến giờ 10 rồi tăng tuyến tính đến 480 ở giờ 12. Ba phần thể tích là 75, 0 và 20 m³, tổng khoảng thiếu 95 m³. Tổng ngày thành 483,75 m³, thấp hơn nội suy xuyên suốt 62,50 m³.
Bước nhảy tại đúng thời điểm đóng không có diện tích, nhưng việc kéo một đường dốc từ số trước đóng đến số sau đóng có thể tạo ra sản lượng trong khoảng đã dừng. Dữ liệu sự kiện nên chia đoạn tích phân và xác định giá trị phía trái, phía phải tại điểm chuyển. Không thể sửa sai này chỉ bằng cách làm lưới thời gian dày hơn trên cùng đường nội suy không phù hợp.
Nếu muốn ngoại suy trung bình 17 giờ có dữ liệu ra toàn ngày, lưu lượng suy ra là 388,75×24/17 = 548,823529 m³/ngày. Đây là một giả thiết khác, cho rằng thời gian thiếu có mức trung bình giống thời gian quan sát. Nó không trùng kết quả 546,25 và cũng không tận dụng thông tin đóng giếng. Mọi lựa chọn phải giữ dấu vết để khi có dữ liệu bổ sung có thể tính lại.
5. Trung bình kỳ và ngày hoạt động
Tài liệu dữ liệu tháng của Whitson phân biệt đơn vị lưu lượng và số ngày giếng hoạt động qua trường Days On [3]. Trong một ví dụ độc lập, báo cáo tháng có thể tích 12.000 m³, tháng dài 30 ngày, trong đó chạy 24 ngày. Bình quân lịch là 400 m³/ngày, còn bình quân khi chạy là 500 m³/ngày. Cả hai đúng nếu tên cột rõ ràng.
Nhân 500 với 30 sẽ tạo 15.000 m³, thừa 3.000. Ngược lại, nhân 400 với 24 chỉ còn 9.600, thiếu 2.400. Nếu hệ thống đã lưu thể tích tháng được kiểm đếm, chuyển qua lưu lượng rồi tích phân lại phải tái tạo đúng 12.000 trước khi ghép vào dự báo. Các số trung bình tháng cũng không phải các điểm tức thời cuối tháng để tự động nối bằng hình thang.
Với bộ đếm cộng dồn, chênh từ 9.800 lên 10.150 m³ cho 350 m³ trong kỳ, nhưng chênh âm sau reset chưa có nghĩa dòng chảy đảo chiều. Cần sự kiện reset hoặc giá trị chốt trước reset. Nếu thiếu thông tin đó, giữ khoảng chưa xác định sẽ trung thực hơn việc tự ép mọi chênh âm về không.
6. Độ cong thật và độ không đảm bảo khác với thiếu dữ liệu
Ngay cả khi hai đầu đều đo đúng, hình thang vẫn có sai số do đường thật cong. Dựng ví dụ độc lập q(τ) = 600exp(−0,1τ), với τ tính giờ, trong 7 giờ. Giá trị cuối là 297,951182 m³/ngày. Tích phân giải tích chia 24 cho 125,853674 m³; hình thang hai đầu cho 130,951214 m³, cao hơn 4,050370%. Dạng mũ lồi nằm dưới dây nối hai đầu, giải thích dấu sai số trong trường hợp cụ thể này.
Chèn thêm các điểm do chính nội suy tuyến tính sinh ra không làm đổi 130,951214 m³. Muốn giảm sai số phải có thêm đo lường hoặc một mô hình có căn cứ về độ cong. Chọn một công thức tích phân bậc cao hơn trên dữ liệu thiếu không tự tạo ra thông tin mới về các lần đổi chế độ.
Với bộ sáu điểm ban đầu, N là tổng có trọng số wᵢqᵢ. Các trọng số theo ngày lần lượt 1/24, 2,5/24, 5/24, 6,5/24, 6/24 và 3/24, tổng bằng 1 ngày. Nếu từng q có độ không đảm bảo chuẩn 10 m³/ngày độc lập, độ không đảm bảo N từ riêng các số đo là 10√Σwᵢ² = 4,554836 m³. Nếu đó là cùng một độ lệch hiệu chuẩn tác động lên mọi điểm, thành phần tương ứng là 10Σwᵢ = 10 m³. Không được coi các điểm nội suy thêm là các phép đo độc lập để làm giảm con số này.
7. Từ tổng số đến một kết quả có thể truy nguyên
Bảng kết quả nên giữ giờ bắt đầu, giờ kết thúc, đơn vị, loại số đo, cờ chất lượng, trạng thái chạy, phương pháp lấp khoảng và đóng góp thể tích từng đoạn. Cần xử lý bản ghi trùng theo nguồn ưu tiên đã công bố; sắp thời gian phải đi cùng lưu lượng tương ứng. Tích phân ngược thứ tự có thể tạo diện tích âm dù mọi q đều dương [1].
Sau cùng, so tổng với kiểm đếm độc lập hoặc cân bằng bồn, nhưng phải xét thay đổi tồn chứa và ranh giới đo. Sai lệch giữa lưu lượng giếng và xuất bán không nhất thiết nằm ở công thức tích phân: chúng có thể thuộc các dòng phụ, cơ sở quy đổi hoặc thời gian lưu khác nhau. Giữ các phần “đã đo”, “đã ước tính” và “chưa xác định” cho phép cập nhật sản lượng mà không làm mất lịch sử ra quyết định.
Nguồn tham khảo
[1] scipy.integrate.trapezoid — SciPy Manual — SciPy community. Truy cập ngày 09/11/2026.
[2] fillmissing — Fill missing entries — MathWorks. Truy cập ngày 09/11/2026.
[3] Production Data — Monthly Data — Whitson. Truy cập ngày 09/11/2026.
Rà soát nội dung: 09/11/2026.
Dữ liệu sản xuất có thể là lưu lượng tức thời, lưu lượng trung bình kỳ hoặc thể tích đã cộng dồn. Ba loại này đòi hỏi ba cách xử lý khác nhau. Với số tức thời cần mô tả diễn biến giữa các lần đo; với trung bình kỳ chỉ cần nhân đúng độ dài kỳ; với bộ đếm thể tích lấy chênh sau khi xử lý các lần reset. Áp dụng hình thang cho cả ba sẽ tạo một tổng có vẻ hợp lý nhưng không bảo toàn ý nghĩa đo lường.
Giả sử q(t) là lưu lượng dầu đã quy đổi về cùng điều kiện báo cáo, đơn vị m³/ngày, còn dấu thời gian tính bằng giờ. Sản lượng trên một ngày là N = (1/24)∫q(t)dt, đơn vị m³. Hệ số 1/24 không được bỏ khi dữ liệu gốc ghi giờ. Nếu cơ sở nhiệt độ, áp suất thay đổi giữa hai nguồn dữ liệu, cần quy đổi trước khi nối thành một chuỗi duy nhất.
2. Ví dụ hình thang trên lưới không đều
Quy tắc hình thang xấp xỉ bằng đoạn thẳng nối hai điểm; SciPy cho phép cung cấp trực tiếp các vị trí lấy mẫu và lưu ý chúng không được tự sắp xếp [1]. Với hai điểm (tᵢ,qᵢ), (tᵢ₊₁,qᵢ₊₁), phần thể tích là ΔNᵢ = (qᵢ+qᵢ₊₁)(tᵢ₊₁−tᵢ)/48 khi t dùng giờ và q dùng m³/ngày.
Xét dữ liệu giả định tại giờ 0, 2, 5, 12, 18, 24 với q lần lượt 600, 660, 600, 480, 540, 500. Nếu chấp nhận đoạn thẳng giữa mọi cặp, năm phần thể tích lần lượt là 52,50; 78,75; 157,50; 127,50; 130,00 m³. Tổng N = 546,25 m³, và lưu lượng bình quân của đúng 24 giờ này cũng có trị số 546,25 m³/ngày.
Không nên lấy trung bình số học sáu lưu lượng rồi nhân một ngày. Các số đo không đại diện cho những khoảng bằng nhau; điểm ở đầu và cuối còn chỉ đóng góp một phía trong hình thang. Tương tự, dùng khoảng cách mặc định bằng một cho tất cả cặp sẽ biến lưới giờ không đều thành một lưới khác.
3. Khoảng thiếu không tự có nghĩa là giếng đóng
Giả sử nhật ký chất lượng cho biết từ giờ 5 đến 12 không có dữ liệu liên tục; hai số ở biên vẫn hợp lệ. Khi ấy 157,50 m³ là phần ước tính theo nội suy, không phải phần đã được quan sát đủ. Bốn khoảng còn lại dài tổng cộng 17 giờ, đóng góp 388,75 m³ theo phương pháp đang dùng. Cần báo riêng thể tích có dữ liệu và thể tích phụ thuộc giả thiết lấp chỗ trống.
MathWorks cung cấp phương pháp điền theo tọa độ thời gian và trả cờ nhận diện điểm đã được điền [2]. Khả năng đó chỉ thực hiện lựa chọn toán học; nó không xác nhận giếng có chạy hay không. Với cùng hai điểm biên, ta có thể tạo nhiều đường nằm giữa, cho nhiều tích phân khác nhau. Xóa các hàng thiếu rồi nối thẳng hai điểm còn lại vẫn chính là chấp nhận một nội suy qua khoảng thiếu.
Giả sử thông tin vận hành độc lập chỉ ràng buộc q nằm từ 450 đến 650 m³/ngày suốt khoảng 7 giờ. Thể tích khoảng đó nằm từ 131,25 đến 189,583333 m³. Tổng ngày tương ứng từ 520,00 đến 578,333333 m³. Đây là khoảng theo ràng buộc đã giả định, không phải khoảng tin cậy thống kê. Nếu ràng buộc không có bằng chứng, chính khoảng này cũng không nên được trình bày như đã biết.
4. Nhật ký đóng giếng thay đổi cách dựng lại đường
Xét kịch bản thay thế: nhật ký xác nhận đóng từ giờ 8 đến 10. Để minh họa, giả sử lưu lượng giữ 600 từ giờ 5 đến sát giờ 8, chuyển ngay về 0, giữ 0 đến giờ 10 rồi tăng tuyến tính đến 480 ở giờ 12. Ba phần thể tích là 75, 0 và 20 m³, tổng khoảng thiếu 95 m³. Tổng ngày thành 483,75 m³, thấp hơn nội suy xuyên suốt 62,50 m³.
Bước nhảy tại đúng thời điểm đóng không có diện tích, nhưng việc kéo một đường dốc từ số trước đóng đến số sau đóng có thể tạo ra sản lượng trong khoảng đã dừng. Dữ liệu sự kiện nên chia đoạn tích phân và xác định giá trị phía trái, phía phải tại điểm chuyển. Không thể sửa sai này chỉ bằng cách làm lưới thời gian dày hơn trên cùng đường nội suy không phù hợp.
Nếu muốn ngoại suy trung bình 17 giờ có dữ liệu ra toàn ngày, lưu lượng suy ra là 388,75×24/17 = 548,823529 m³/ngày. Đây là một giả thiết khác, cho rằng thời gian thiếu có mức trung bình giống thời gian quan sát. Nó không trùng kết quả 546,25 và cũng không tận dụng thông tin đóng giếng. Mọi lựa chọn phải giữ dấu vết để khi có dữ liệu bổ sung có thể tính lại.
5. Trung bình kỳ và ngày hoạt động
Tài liệu dữ liệu tháng của Whitson phân biệt đơn vị lưu lượng và số ngày giếng hoạt động qua trường Days On [3]. Trong một ví dụ độc lập, báo cáo tháng có thể tích 12.000 m³, tháng dài 30 ngày, trong đó chạy 24 ngày. Bình quân lịch là 400 m³/ngày, còn bình quân khi chạy là 500 m³/ngày. Cả hai đúng nếu tên cột rõ ràng.
Nhân 500 với 30 sẽ tạo 15.000 m³, thừa 3.000. Ngược lại, nhân 400 với 24 chỉ còn 9.600, thiếu 2.400. Nếu hệ thống đã lưu thể tích tháng được kiểm đếm, chuyển qua lưu lượng rồi tích phân lại phải tái tạo đúng 12.000 trước khi ghép vào dự báo. Các số trung bình tháng cũng không phải các điểm tức thời cuối tháng để tự động nối bằng hình thang.
Với bộ đếm cộng dồn, chênh từ 9.800 lên 10.150 m³ cho 350 m³ trong kỳ, nhưng chênh âm sau reset chưa có nghĩa dòng chảy đảo chiều. Cần sự kiện reset hoặc giá trị chốt trước reset. Nếu thiếu thông tin đó, giữ khoảng chưa xác định sẽ trung thực hơn việc tự ép mọi chênh âm về không.
6. Độ cong thật và độ không đảm bảo khác với thiếu dữ liệu
Ngay cả khi hai đầu đều đo đúng, hình thang vẫn có sai số do đường thật cong. Dựng ví dụ độc lập q(τ) = 600exp(−0,1τ), với τ tính giờ, trong 7 giờ. Giá trị cuối là 297,951182 m³/ngày. Tích phân giải tích chia 24 cho 125,853674 m³; hình thang hai đầu cho 130,951214 m³, cao hơn 4,050370%. Dạng mũ lồi nằm dưới dây nối hai đầu, giải thích dấu sai số trong trường hợp cụ thể này.
Chèn thêm các điểm do chính nội suy tuyến tính sinh ra không làm đổi 130,951214 m³. Muốn giảm sai số phải có thêm đo lường hoặc một mô hình có căn cứ về độ cong. Chọn một công thức tích phân bậc cao hơn trên dữ liệu thiếu không tự tạo ra thông tin mới về các lần đổi chế độ.
Với bộ sáu điểm ban đầu, N là tổng có trọng số wᵢqᵢ. Các trọng số theo ngày lần lượt 1/24, 2,5/24, 5/24, 6,5/24, 6/24 và 3/24, tổng bằng 1 ngày. Nếu từng q có độ không đảm bảo chuẩn 10 m³/ngày độc lập, độ không đảm bảo N từ riêng các số đo là 10√Σwᵢ² = 4,554836 m³. Nếu đó là cùng một độ lệch hiệu chuẩn tác động lên mọi điểm, thành phần tương ứng là 10Σwᵢ = 10 m³. Không được coi các điểm nội suy thêm là các phép đo độc lập để làm giảm con số này.
7. Từ tổng số đến một kết quả có thể truy nguyên
Bảng kết quả nên giữ giờ bắt đầu, giờ kết thúc, đơn vị, loại số đo, cờ chất lượng, trạng thái chạy, phương pháp lấp khoảng và đóng góp thể tích từng đoạn. Cần xử lý bản ghi trùng theo nguồn ưu tiên đã công bố; sắp thời gian phải đi cùng lưu lượng tương ứng. Tích phân ngược thứ tự có thể tạo diện tích âm dù mọi q đều dương [1].
Sau cùng, so tổng với kiểm đếm độc lập hoặc cân bằng bồn, nhưng phải xét thay đổi tồn chứa và ranh giới đo. Sai lệch giữa lưu lượng giếng và xuất bán không nhất thiết nằm ở công thức tích phân: chúng có thể thuộc các dòng phụ, cơ sở quy đổi hoặc thời gian lưu khác nhau. Giữ các phần “đã đo”, “đã ước tính” và “chưa xác định” cho phép cập nhật sản lượng mà không làm mất lịch sử ra quyết định.
Nguồn tham khảo
[1] scipy.integrate.trapezoid — SciPy Manual — SciPy community. Truy cập ngày 09/11/2026.
[2] fillmissing — Fill missing entries — MathWorks. Truy cập ngày 09/11/2026.
[3] Production Data — Monthly Data — Whitson. Truy cập ngày 09/11/2026.
Rà soát nội dung: 09/11/2026.
