1. Một bán kính có thể đang đại diện cho những việc khác nhau
Trong bảng tính năng suất, r_e thường được gọi là bán kính tháo khô. Nhưng ký hiệu ấy chỉ có ý nghĩa sau khi xác định điều kiện biên: áp suất được giữ ở đó, biên kín hay khoảng cách tới một vùng chịu ảnh hưởng giếng khác. Bán kính hình tròn cùng diện tích với ô khai thác là một quy đổi hình học; bán kính khảo sát từ thử áp suất lại gắn với thời gian và độ nhạy tín hiệu.
Bài này dùng miền tròn có áp suất biên không đổi để tách các ý nghĩa. Học liệu Groundwater Project xây dựng dòng ổn định từ Darcy và bảo toàn lưu lượng qua mặt trụ.[1] Toán học đó dùng được cho chất lỏng đơn pha có tính chất không đổi khi thay đúng độ thấm và độ nhớt. Không sử dụng hệ số thấm thủy lực của nước như thể nó là độ thấm riêng của đá dầu khí.
2. Vì sao kích thước miền đi vào logarithm
Giả định giếng thẳng đứng xuyên hết h, đá đồng nhất đẳng hướng, bỏ skin và chênh cao. Áp suất đáy p_wf và áp suất p_e ở r_e được giữ không đổi. Từ q_r=2πrhk(dp/dr)/μ, tích phân cho:
q_s = 2πkh(p_e−p_wf)/[μB_o ln(r_e/r_w)]
p(r) = p_wf + (p_e−p_wf)ln(r/r_w)/ln(r_e/r_w)
Ở đây q_s tính bằng m³/s bồn, k bằng m², h và bán kính bằng m, μ bằng Pa·s, áp suất bằng Pa. B_o là thể tích dầu vỉa trên thể tích dầu bồn. Logarithm tự nhiên xuất hiện vì diện tích dẫn dòng tăng tuyến tính theo r; mỗi khoảng bán kính tăng theo cùng một tỷ số tạo cùng phần chênh áp trong miền đồng nhất.
Chọn k=1,1×10⁻¹³ m², h=20 m, μ=0,002 Pa·s, B_o=1,20, r_w=0,10 m, r_e=500 m, p_e=25 MPa và p_wf=20 MPa. Kết quả q_s=292,131613 m³/ngày. Áp suất ở r=1;10;100 m lần lượt là 21,351728;22,703455;24,055183 MPa. Mỗi thập kỷ bán kính tạo khoảng 1,351728 MPa, dù bề rộng vật lý của các vành rất khác nhau.
3. Tăng đôi bán kính khác tăng đôi diện tích
Giữ nguyên áp suất hai đầu và mọi thuộc tính, tăng r_e lên 1.000 m cho q_s=270,146518 m³/ngày, giảm 7,525750%. Bán kính tăng đôi làm đường đi xa hơn; tại cùng chênh áp, trở lực tăng và lưu lượng giảm. Đây không phải kết luận rằng có thêm tài nguyên làm giếng kém hơn trong mọi lịch sử khai thác.
Thể tích lỗ rỗng hình học tỷ lệ với r_e²−r_w². Trong hai miền vừa xét, tỷ số thể tích gần 4,00000012, trong khi lưu lượng ổn định chỉ giảm khoảng 7,5%. Lượng chứa và khả năng truyền dòng là hai đại lượng khác nhau. Dùng độ nhạy năng suất nhỏ để kết luận kích thước miền không quan trọng cho trữ lượng là một suy luận sai.
Nếu chỉ tăng diện tích gấp đôi, bán kính mới là 500√2 m và q_s=280,709256 m³/ngày. Đạo hàm cục bộ, với r_w nhỏ so với r_e, cho:
d ln q_s / d ln r_e = −1/ln(r_e/r_w)
d ln q_s / d ln A ≈ −1/[2ln(r_e/r_w)]
Hai hệ số tại điểm cơ sở là −0,117409571 và −0,058704786. Đây là độ nhạy vi phân; khi thay đầu vào nhiều, cần tính lại biểu thức đầy đủ thay vì nhân tuyến tính.
4. Áp suất trung bình không bằng áp suất biên
Với độ rỗng và chiều dày không đổi, áp suất trung bình theo thể tích của vành đá là tích phân p(r)2r dr chia r_e²−r_w². Đặt C=(p_e−p_wf)/ln(r_e/r_w), tích phân trực tiếp được:
p̄ = p_e − C[1/2 − r_w²ln(r_e/r_w)/(r_e²−r_w²)]
Trong ví dụ, p̄=24,706476272 MPa, thấp hơn p_e. Nếu thay p_e bằng p̄ mà giữ nguyên mẫu số ln(r_e/r_w), sẽ giảm drawdown từ 5 xuống 4,706476272 MPa và làm sai lưu lượng. Để dùng p̄ trong đúng mô hình ổn định này, mẫu số cần đổi thành ln(r_e/r_w)−1/2+r_w²ln(r_e/r_w)/(r_e²−r_w²).
Số 1/2 ở đây được suy từ phép lấy trung bình một trường áp suất ổn định. Nó không phải hằng số được tự do đổi sang 3/4, hoặc gắn vào mọi giếng có dữ liệu áp suất trung bình. Một miền kín đang suy giảm áp suất có phương trình và điều kiện khác; cần thiết lập lại bài toán trước khi chọn hằng số hình học.
5. Cùng diện tích chưa chắc cùng năng suất
Tạo một mô hình khác gồm hai hình quạt độc lập, ngăn bằng các vách xuyên tâm không thấm. Hình quạt thứ nhất chiếm 1/4 góc tròn, bán kính ngoài 200 m; phần còn lại chiếm 3/4, bán kính 800 m. Cả hai cùng k,h,μ,B_o,p_e,p_wf và r_w như trên. Vách ngăn bảo đảm không có trao đổi theo phương góc, nên có thể giải chính xác từng phần như các đường dẫn song song.
Gọi K=2πkh(p_e−p_wf)×86.400/(μB_o), đơn vị m³/ngày. Tổng lưu lượng bằng:
q_s = K[0,25/ln(200/0,1) + 0,75/ln(800/0,1)]
q_s = 289,477538 m³/ngày
Bán kính tròn cùng diện tích là √(0,25×200²+0,75×800²)=700 m. Thay miền bằng hình tròn 700 m cho 281,029524 m³/ngày, thấp hơn mô hình hình quạt 2,918366%. Bán kính tái tạo đúng trở lực của mô hình song song lại là 540,609922 m, không phải 700 m.
Ví dụ này không đề xuất dùng các vách tưởng tượng cho vỉa thực. Nó chứng minh diện tích và trở lực không xác định cùng một bán kính tương đương. Trong miền thực không có vách, dòng có thành phần theo phương góc; khi đó cần giải trường áp suất của chính hình học đó hoặc dùng mô hình đã kiểm chứng.
6. Suy ngược bán kính có thể rất nhạy
Ở mô hình tròn không skin, từ một lưu lượng cố định có thể viết r_e=r_w exp[K/q_s]. Vì vậy, độ nhạy thuận logarithm nhỏ không đồng nghĩa bài toán ngược ổn định. Giữ q_s đo được của trường hợp cơ sở và tăng k giả định 10%, bán kính suy ra trở thành 1.171,836456 m, hơn gấp đôi 500 m.
Nếu đồng thời cho phép skin thay đổi, một điểm q–drawdown chỉ giữ cố định tổng ln(r_e/r_w)+s. Tăng r_e gấp đôi có thể được bù bằng giảm s đi ln2=0,693147181. Không thể xác định độc lập hai đầu vào từ cùng một phương trình.
IHS yêu cầu nhận diện giai đoạn xuyên tâm qua dạng đường áp suất và đạo hàm trước khi diễn giải kh và skin.[2] Điều này nhấn mạnh vai trò của hình dạng đáp ứng theo thời gian, thay vì cố suy mọi tham số từ một điểm năng suất.
7. Bán kính khảo sát không phải một mặt biên được đo trực tiếp
TestWells trình bày nhiều cách định nghĩa bán kính khảo sát và lưu ý nghiệm khuếch tán đầy đủ không có một bán kính hữu hạn nơi thay đổi áp suất bằng đúng không.[3] Vì vậy, r_i từ một quy ước tín hiệu không thể tự thay cho r_e của bài toán giữ áp suất biên.
Một hồ sơ tính cần ghi rõ bán kính dùng cho diện tích chứa, cho điều kiện biên hay cho mức nhạy của thử áp suất. Khi hình học, thời gian và áp suất được giữ nhất quán, độ nhạy logarithm giúp đánh giá độ vững của dự báo năng suất; nó không loại bỏ nhu cầu xác định miền liên thông và lịch sử suy giảm áp suất của vỉa.
Nguồn tham khảo
[1] 7.1 Steady-State Conditions in a Confined Groundwater System — The Groundwater Project, An Introduction to Hydraulic Testing in Hydrogeology. Truy cập ngày 09/11/2026.
[2] Radial Flow Analysis — S&P Global, WellTest2024.1; revised6June2024. Truy cập ngày 09/11/2026.
[3] Radius of Investigation — TestWells. Truy cập ngày 09/11/2026.
Rà soát nội dung: 09/11/2026.
Trong bảng tính năng suất, r_e thường được gọi là bán kính tháo khô. Nhưng ký hiệu ấy chỉ có ý nghĩa sau khi xác định điều kiện biên: áp suất được giữ ở đó, biên kín hay khoảng cách tới một vùng chịu ảnh hưởng giếng khác. Bán kính hình tròn cùng diện tích với ô khai thác là một quy đổi hình học; bán kính khảo sát từ thử áp suất lại gắn với thời gian và độ nhạy tín hiệu.
Bài này dùng miền tròn có áp suất biên không đổi để tách các ý nghĩa. Học liệu Groundwater Project xây dựng dòng ổn định từ Darcy và bảo toàn lưu lượng qua mặt trụ.[1] Toán học đó dùng được cho chất lỏng đơn pha có tính chất không đổi khi thay đúng độ thấm và độ nhớt. Không sử dụng hệ số thấm thủy lực của nước như thể nó là độ thấm riêng của đá dầu khí.
2. Vì sao kích thước miền đi vào logarithm
Giả định giếng thẳng đứng xuyên hết h, đá đồng nhất đẳng hướng, bỏ skin và chênh cao. Áp suất đáy p_wf và áp suất p_e ở r_e được giữ không đổi. Từ q_r=2πrhk(dp/dr)/μ, tích phân cho:
q_s = 2πkh(p_e−p_wf)/[μB_o ln(r_e/r_w)]
p(r) = p_wf + (p_e−p_wf)ln(r/r_w)/ln(r_e/r_w)
Ở đây q_s tính bằng m³/s bồn, k bằng m², h và bán kính bằng m, μ bằng Pa·s, áp suất bằng Pa. B_o là thể tích dầu vỉa trên thể tích dầu bồn. Logarithm tự nhiên xuất hiện vì diện tích dẫn dòng tăng tuyến tính theo r; mỗi khoảng bán kính tăng theo cùng một tỷ số tạo cùng phần chênh áp trong miền đồng nhất.
Chọn k=1,1×10⁻¹³ m², h=20 m, μ=0,002 Pa·s, B_o=1,20, r_w=0,10 m, r_e=500 m, p_e=25 MPa và p_wf=20 MPa. Kết quả q_s=292,131613 m³/ngày. Áp suất ở r=1;10;100 m lần lượt là 21,351728;22,703455;24,055183 MPa. Mỗi thập kỷ bán kính tạo khoảng 1,351728 MPa, dù bề rộng vật lý của các vành rất khác nhau.
3. Tăng đôi bán kính khác tăng đôi diện tích
Giữ nguyên áp suất hai đầu và mọi thuộc tính, tăng r_e lên 1.000 m cho q_s=270,146518 m³/ngày, giảm 7,525750%. Bán kính tăng đôi làm đường đi xa hơn; tại cùng chênh áp, trở lực tăng và lưu lượng giảm. Đây không phải kết luận rằng có thêm tài nguyên làm giếng kém hơn trong mọi lịch sử khai thác.
Thể tích lỗ rỗng hình học tỷ lệ với r_e²−r_w². Trong hai miền vừa xét, tỷ số thể tích gần 4,00000012, trong khi lưu lượng ổn định chỉ giảm khoảng 7,5%. Lượng chứa và khả năng truyền dòng là hai đại lượng khác nhau. Dùng độ nhạy năng suất nhỏ để kết luận kích thước miền không quan trọng cho trữ lượng là một suy luận sai.
Nếu chỉ tăng diện tích gấp đôi, bán kính mới là 500√2 m và q_s=280,709256 m³/ngày. Đạo hàm cục bộ, với r_w nhỏ so với r_e, cho:
d ln q_s / d ln r_e = −1/ln(r_e/r_w)
d ln q_s / d ln A ≈ −1/[2ln(r_e/r_w)]
Hai hệ số tại điểm cơ sở là −0,117409571 và −0,058704786. Đây là độ nhạy vi phân; khi thay đầu vào nhiều, cần tính lại biểu thức đầy đủ thay vì nhân tuyến tính.
4. Áp suất trung bình không bằng áp suất biên
Với độ rỗng và chiều dày không đổi, áp suất trung bình theo thể tích của vành đá là tích phân p(r)2r dr chia r_e²−r_w². Đặt C=(p_e−p_wf)/ln(r_e/r_w), tích phân trực tiếp được:
p̄ = p_e − C[1/2 − r_w²ln(r_e/r_w)/(r_e²−r_w²)]
Trong ví dụ, p̄=24,706476272 MPa, thấp hơn p_e. Nếu thay p_e bằng p̄ mà giữ nguyên mẫu số ln(r_e/r_w), sẽ giảm drawdown từ 5 xuống 4,706476272 MPa và làm sai lưu lượng. Để dùng p̄ trong đúng mô hình ổn định này, mẫu số cần đổi thành ln(r_e/r_w)−1/2+r_w²ln(r_e/r_w)/(r_e²−r_w²).
Số 1/2 ở đây được suy từ phép lấy trung bình một trường áp suất ổn định. Nó không phải hằng số được tự do đổi sang 3/4, hoặc gắn vào mọi giếng có dữ liệu áp suất trung bình. Một miền kín đang suy giảm áp suất có phương trình và điều kiện khác; cần thiết lập lại bài toán trước khi chọn hằng số hình học.
5. Cùng diện tích chưa chắc cùng năng suất
Tạo một mô hình khác gồm hai hình quạt độc lập, ngăn bằng các vách xuyên tâm không thấm. Hình quạt thứ nhất chiếm 1/4 góc tròn, bán kính ngoài 200 m; phần còn lại chiếm 3/4, bán kính 800 m. Cả hai cùng k,h,μ,B_o,p_e,p_wf và r_w như trên. Vách ngăn bảo đảm không có trao đổi theo phương góc, nên có thể giải chính xác từng phần như các đường dẫn song song.
Gọi K=2πkh(p_e−p_wf)×86.400/(μB_o), đơn vị m³/ngày. Tổng lưu lượng bằng:
q_s = K[0,25/ln(200/0,1) + 0,75/ln(800/0,1)]
q_s = 289,477538 m³/ngày
Bán kính tròn cùng diện tích là √(0,25×200²+0,75×800²)=700 m. Thay miền bằng hình tròn 700 m cho 281,029524 m³/ngày, thấp hơn mô hình hình quạt 2,918366%. Bán kính tái tạo đúng trở lực của mô hình song song lại là 540,609922 m, không phải 700 m.
Ví dụ này không đề xuất dùng các vách tưởng tượng cho vỉa thực. Nó chứng minh diện tích và trở lực không xác định cùng một bán kính tương đương. Trong miền thực không có vách, dòng có thành phần theo phương góc; khi đó cần giải trường áp suất của chính hình học đó hoặc dùng mô hình đã kiểm chứng.
6. Suy ngược bán kính có thể rất nhạy
Ở mô hình tròn không skin, từ một lưu lượng cố định có thể viết r_e=r_w exp[K/q_s]. Vì vậy, độ nhạy thuận logarithm nhỏ không đồng nghĩa bài toán ngược ổn định. Giữ q_s đo được của trường hợp cơ sở và tăng k giả định 10%, bán kính suy ra trở thành 1.171,836456 m, hơn gấp đôi 500 m.
Nếu đồng thời cho phép skin thay đổi, một điểm q–drawdown chỉ giữ cố định tổng ln(r_e/r_w)+s. Tăng r_e gấp đôi có thể được bù bằng giảm s đi ln2=0,693147181. Không thể xác định độc lập hai đầu vào từ cùng một phương trình.
IHS yêu cầu nhận diện giai đoạn xuyên tâm qua dạng đường áp suất và đạo hàm trước khi diễn giải kh và skin.[2] Điều này nhấn mạnh vai trò của hình dạng đáp ứng theo thời gian, thay vì cố suy mọi tham số từ một điểm năng suất.
7. Bán kính khảo sát không phải một mặt biên được đo trực tiếp
TestWells trình bày nhiều cách định nghĩa bán kính khảo sát và lưu ý nghiệm khuếch tán đầy đủ không có một bán kính hữu hạn nơi thay đổi áp suất bằng đúng không.[3] Vì vậy, r_i từ một quy ước tín hiệu không thể tự thay cho r_e của bài toán giữ áp suất biên.
Một hồ sơ tính cần ghi rõ bán kính dùng cho diện tích chứa, cho điều kiện biên hay cho mức nhạy của thử áp suất. Khi hình học, thời gian và áp suất được giữ nhất quán, độ nhạy logarithm giúp đánh giá độ vững của dự báo năng suất; nó không loại bỏ nhu cầu xác định miền liên thông và lịch sử suy giảm áp suất của vỉa.
Nguồn tham khảo
[1] 7.1 Steady-State Conditions in a Confined Groundwater System — The Groundwater Project, An Introduction to Hydraulic Testing in Hydrogeology. Truy cập ngày 09/11/2026.
[2] Radial Flow Analysis — S&P Global, WellTest2024.1; revised6June2024. Truy cập ngày 09/11/2026.
[3] Radius of Investigation — TestWells. Truy cập ngày 09/11/2026.
Rà soát nội dung: 09/11/2026.
