1. Một độ trễ quan sát cần một định nghĩa
Khi giếng quan sát đổi áp suất sau một tác động ở giếng khác, thời gian ghi nhận có thể được định nghĩa theo lúc vượt ngưỡng, lúc đạo hàm đổi dấu hoặc lúc áp suất đạt cực đại. Ba cách đó không nhất thiết cho cùng kết quả. Giáo trình NTNU dùng biến không thứ nguyên để đưa thời gian và khoảng cách về cùng thang khuếch tán, đồng thời chỉ rõ bán kính điều tra là một quy ước mô hình [1]. Bài này xây dựng một xung thể tích giả định để làm rõ ý nghĩa của từng thời điểm.
Xét chất lỏng một pha hơi nén được trong lớp đồng nhất, đẳng hướng, chiều dày không đổi, mở rộng vô hạn theo phương ngang. Bỏ qua thay đổi tính chất theo áp suất, trọng lực trong mặt phẳng và ảnh hưởng biên. Phương trình tuyến tính là ∂p/∂t = D[(1/r)∂/∂r(r∂p/∂r)], với D = k/(φμcₜ). k tính bằng m², φ là độ rỗng, μ dùng Pa·s và cₜ dùng Pa⁻¹; do đó D có đơn vị m²/s.
Cấu trúc tỷ số giữa khả năng dẫn và khả năng tích trữ cũng xuất hiện trong phương trình nước ngầm [2]. Nó cho thấy tăng k làm áp lực phân bố nhanh hơn, còn tăng tích trữ làm cần nhiều chất lưu hơn để tạo cùng biến đổi áp suất. D không phải vận tốc của một phần tử nước; đơn vị m²/s đã khác m/s.
2. Chọn chiều dài đặc trưng trước khi tính thời gian
Nếu lấy bán kính giếng rᵥ làm chiều dài đặc trưng, đặt rD = r/rᵥ và tD = Dt/rᵥ². Đây là một cách quy đổi dùng trong phân tích thử giếng [1]. Nếu quan tâm một giếng quan sát cách R, có thể dùng θ = Dt/R². Hai biến liên hệ bởi θ = tD/rD² tại r = R. Một giá trị tD lớn tại bán kính giếng không đồng nghĩa đáp ứng ở một khoảng cách rất xa đã đủ rõ.
Ví dụ giả định k = 5×10⁻¹⁴ m², φ = 0,18, μ = 1,2×10⁻³ Pa·s và cₜ = 1,5×10⁻⁹ Pa⁻¹ cho D = 0,154320988 m²/s. Với rᵥ = 0,10 m, sau một giờ tD = 55.555,555556. Nhưng tại R = 100 m, θ chỉ bằng 0,055556. Hai con số mô tả cùng một trạng thái, được chuẩn hóa bằng hai chiều dài khác nhau.
Thang chiều dài √(4Dt) ở một giờ bằng 47,140452 m. Đây là đường mức của nhóm r²/(4Dt), chưa phải một mặt phân cách giữa vùng “đã tới” và “chưa tới”. Nghiệm khuếch tán toán học có đuôi rất nhỏ ngoài đường mức đó. Mô hình tuyến tính này không dùng để giải quá trình truyền sóng áp suất âm học ở thời gian cực ngắn.
3. Xung thể tích cho một ví dụ có cực đại rõ ràng
Giả sử tại t = 0 có thể tích V được bổ sung tức thời vào nguồn đường xuyên suốt lớp dày h. Phân bố áp suất tăng tại r > 0, t > 0 là Δp(r,t) = V exp[−r²/(4Dt)]/(4πDthφcₜ). Ta có thể kiểm hệ số bằng tích phân tích trữ: ∫₀∞hφcₜΔp·2πr dr = V. Như vậy tổng thể tích bổ sung được bảo toàn trong mô hình, dù áp suất tại từng điểm tăng rồi giảm.
Dạng hàm mũ chia thời gian này cũng thu được khi lấy đạo hàm thời gian của đáp ứng lưu lượng bậc thang. Tài liệu ORION trình bày nghiệm Theis, nhóm r²/(4Dt) và phép chồng chất các bước lưu lượng [3]. Ở đây xung tức thời chỉ là mô hình giới hạn; một đợt ép có độ dài hữu hạn cần tích phân các đáp ứng theo lịch sử lưu lượng thực.
Lấy V = 20 m³ và h = 20 m, giữ D cùng các tính chất ở trên. Tại R = 100 m, đạo hàm của lnΔp theo t là −1/t + R²/(4Dt²). Cho bằng không được tₚ = R²/(4D) = 16.200 s = 4,5 giờ. Đạo hàm dương trước tₚ và âm sau tₚ, nên đây là cực đại duy nhất của đáp ứng tại điểm đó. Không cần giả định một vận tốc lan truyền cố định để tìm thời gian này.
Thay tₚ vào nghiệm được Δpₚ = V/(πR²hφcₜe) = 43.370,245574 Pa, tức 43,370246 kPa. Trong cùng mô hình, giếng ở 200 m có tₚ = 18 giờ và đỉnh chỉ còn 10,842561 kPa. Gấp đôi khoảng cách làm thời gian cực đại tăng bốn lần và biên độ đỉnh giảm bốn lần. Đây là dấu hiệu khác với chuyển động đối lưu đều, vốn tạo thời gian tỷ lệ bậc nhất với khoảng cách.
4. Thời điểm phát hiện phụ thuộc ngưỡng
Đặt x = t/tₚ thì Δp/Δpₚ = exp(1−1/x)/x. Tại x bằng 0,10; 0,25; 0,50; 1; 2; 4, áp suất ở 100 m lần lượt bằng 0,053523; 8,637110; 31,910043; 43,370246; 35,752723; 22,953703 kPa. Chuỗi này cho thấy tại cùng một vị trí, đáp ứng không bắt đầu bằng một bước nhảy từ không lên cực đại.
Giả sử chỉ để phân tích toán học, ngưỡng phát hiện là 10 kPa trên nền ổn định. Giải exp(1−1/x)/x = 10/43,370246 trên nhánh tăng được x = 0,262952295 và thời gian 1,183285 giờ. Nếu ngưỡng là 1 kPa, thời điểm thành 0,674973 giờ; nếu là 20 kPa, thành 1,603883 giờ. D vẫn không đổi. Vì vậy lấy “lúc đầu nhìn thấy tín hiệu” làm thời gian truyền mà không ghi ngưỡng có thể tạo ra các ước lượng khác nhau từ cùng một mô hình.
Phương trình ngưỡng còn có nghiệm trên nhánh giảm: với 10 kPa, tín hiệu xuống dưới ngưỡng sau khoảng 48,335520 giờ. Một thuật toán tìm nghiệm phải xác định đang tìm lần vượt lên hay lần đi xuống. Nếu ngưỡng cao hơn đỉnh, hoàn toàn không có nghiệm dù hai điểm vẫn liên thông trong mô hình. Không phát hiện trong cửa sổ đo vì thế chưa đủ chứng minh không có liên thông.
5. Độ nhạy giúp nhận ra tham số bị lẫn
Từ tₚ = R²φμcₜ/(4k), biến thiên nhỏ thỏa Δtₚ/tₚ ≈ 2ΔR/R + Δφ/φ + Δμ/μ + Δcₜ/cₜ − Δk/k. Nếu khoảng cách thực là 105 m thay vì 100 m, tₚ bằng 4,961250 giờ, tăng 10,25%; xấp xỉ tuyến tính cho 10%. Độ nhạy bình phương theo khoảng cách cần được giữ khi tọa độ giếng hoặc vị trí lớp tiếp nhận chưa chắc chắn.
Giảm k một nửa làm tₚ tăng từ 4,5 lên 9 giờ nhưng không đổi đỉnh 43,370246 kPa của xung có cùng V, h, φ, cₜ. Ngược lại, tăng cₜ gấp đôi cũng cho thời gian 9 giờ nhưng đỉnh giảm còn 21,685123 kPa. Hai trường hợp giống nhau về thời gian cực đại mà khác về biên độ. Đây là lý do nên khớp toàn bộ đường đáp ứng và kiểm thể tích tác động, thay vì suy một tham số chỉ từ một thời điểm.
Ngay cả khi biết D, phương trình D = k/(φμcₜ) chỉ xác định một tổ hợp tham số. Muốn suy k riêng cần thông tin về các đại lượng còn lại. Nếu V thực vào vỉa khác thể tích đo ở mặt đất do tích trữ trong giếng, biên độ cũng không đủ để suy cₜ mà không xử lý ranh giới đó. Các tham số có thể bù trừ nhau trong một đường khớp nhìn rất đẹp.
6. Phạm vi của phép diễn giải
Một đợt ép kéo dài, thay đổi lưu lượng nhiều lần hoặc dừng không đồng bộ với đồng hồ đều làm đáp ứng khác xung lý tưởng. Khi thời lượng tác động đáng kể so với tₚ, phải dùng chồng chất theo thời gian. Bộ lọc tín hiệu và quán tính cảm biến cũng có thể dịch đỉnh quan sát; đó là đặc tính hệ đo, cần tách khỏi khuếch tán của vỉa.
Biên kín, biên áp suất cố định, tính dị hướng, nhiều lớp, khe nứt và dòng nhiều pha có thể đổi dạng đường cong. Do đó kết quả 4,5 giờ trong ví dụ không phải thời gian chờ thử giếng được khuyến nghị. Nó là kết quả của một hệ giả định có thể tính lại, dùng để kiểm sự khác biệt giữa thời gian không thứ nguyên, cực đại, ngưỡng phát hiện và vận chuyển chất lưu. Báo cáo một “độ trễ” chỉ có nghĩa kỹ thuật khi định nghĩa, lịch sử tác động và mô hình đi cùng nó được công bố.
Nguồn tham khảo
[1] Lecture notes in well-testing — sections 3.3–3.4, print date 23 April 2026 — Per Arne Slotte và Carl Fredrik Berg, NTNU. Truy cập ngày 09/11/2026.
[2] Groundwater Storage in Confined Aquifers — 5 Diffusion Equation — The Groundwater Project. Truy cập ngày 09/11/2026.
[3] Radial Flow Pressure Model — ORION documentation — National Risk Assessment Partnership, ORION. Truy cập ngày 09/11/2026.
Rà soát nội dung: 09/11/2026.
Khi giếng quan sát đổi áp suất sau một tác động ở giếng khác, thời gian ghi nhận có thể được định nghĩa theo lúc vượt ngưỡng, lúc đạo hàm đổi dấu hoặc lúc áp suất đạt cực đại. Ba cách đó không nhất thiết cho cùng kết quả. Giáo trình NTNU dùng biến không thứ nguyên để đưa thời gian và khoảng cách về cùng thang khuếch tán, đồng thời chỉ rõ bán kính điều tra là một quy ước mô hình [1]. Bài này xây dựng một xung thể tích giả định để làm rõ ý nghĩa của từng thời điểm.
Xét chất lỏng một pha hơi nén được trong lớp đồng nhất, đẳng hướng, chiều dày không đổi, mở rộng vô hạn theo phương ngang. Bỏ qua thay đổi tính chất theo áp suất, trọng lực trong mặt phẳng và ảnh hưởng biên. Phương trình tuyến tính là ∂p/∂t = D[(1/r)∂/∂r(r∂p/∂r)], với D = k/(φμcₜ). k tính bằng m², φ là độ rỗng, μ dùng Pa·s và cₜ dùng Pa⁻¹; do đó D có đơn vị m²/s.
Cấu trúc tỷ số giữa khả năng dẫn và khả năng tích trữ cũng xuất hiện trong phương trình nước ngầm [2]. Nó cho thấy tăng k làm áp lực phân bố nhanh hơn, còn tăng tích trữ làm cần nhiều chất lưu hơn để tạo cùng biến đổi áp suất. D không phải vận tốc của một phần tử nước; đơn vị m²/s đã khác m/s.
2. Chọn chiều dài đặc trưng trước khi tính thời gian
Nếu lấy bán kính giếng rᵥ làm chiều dài đặc trưng, đặt rD = r/rᵥ và tD = Dt/rᵥ². Đây là một cách quy đổi dùng trong phân tích thử giếng [1]. Nếu quan tâm một giếng quan sát cách R, có thể dùng θ = Dt/R². Hai biến liên hệ bởi θ = tD/rD² tại r = R. Một giá trị tD lớn tại bán kính giếng không đồng nghĩa đáp ứng ở một khoảng cách rất xa đã đủ rõ.
Ví dụ giả định k = 5×10⁻¹⁴ m², φ = 0,18, μ = 1,2×10⁻³ Pa·s và cₜ = 1,5×10⁻⁹ Pa⁻¹ cho D = 0,154320988 m²/s. Với rᵥ = 0,10 m, sau một giờ tD = 55.555,555556. Nhưng tại R = 100 m, θ chỉ bằng 0,055556. Hai con số mô tả cùng một trạng thái, được chuẩn hóa bằng hai chiều dài khác nhau.
Thang chiều dài √(4Dt) ở một giờ bằng 47,140452 m. Đây là đường mức của nhóm r²/(4Dt), chưa phải một mặt phân cách giữa vùng “đã tới” và “chưa tới”. Nghiệm khuếch tán toán học có đuôi rất nhỏ ngoài đường mức đó. Mô hình tuyến tính này không dùng để giải quá trình truyền sóng áp suất âm học ở thời gian cực ngắn.
3. Xung thể tích cho một ví dụ có cực đại rõ ràng
Giả sử tại t = 0 có thể tích V được bổ sung tức thời vào nguồn đường xuyên suốt lớp dày h. Phân bố áp suất tăng tại r > 0, t > 0 là Δp(r,t) = V exp[−r²/(4Dt)]/(4πDthφcₜ). Ta có thể kiểm hệ số bằng tích phân tích trữ: ∫₀∞hφcₜΔp·2πr dr = V. Như vậy tổng thể tích bổ sung được bảo toàn trong mô hình, dù áp suất tại từng điểm tăng rồi giảm.
Dạng hàm mũ chia thời gian này cũng thu được khi lấy đạo hàm thời gian của đáp ứng lưu lượng bậc thang. Tài liệu ORION trình bày nghiệm Theis, nhóm r²/(4Dt) và phép chồng chất các bước lưu lượng [3]. Ở đây xung tức thời chỉ là mô hình giới hạn; một đợt ép có độ dài hữu hạn cần tích phân các đáp ứng theo lịch sử lưu lượng thực.
Lấy V = 20 m³ và h = 20 m, giữ D cùng các tính chất ở trên. Tại R = 100 m, đạo hàm của lnΔp theo t là −1/t + R²/(4Dt²). Cho bằng không được tₚ = R²/(4D) = 16.200 s = 4,5 giờ. Đạo hàm dương trước tₚ và âm sau tₚ, nên đây là cực đại duy nhất của đáp ứng tại điểm đó. Không cần giả định một vận tốc lan truyền cố định để tìm thời gian này.
Thay tₚ vào nghiệm được Δpₚ = V/(πR²hφcₜe) = 43.370,245574 Pa, tức 43,370246 kPa. Trong cùng mô hình, giếng ở 200 m có tₚ = 18 giờ và đỉnh chỉ còn 10,842561 kPa. Gấp đôi khoảng cách làm thời gian cực đại tăng bốn lần và biên độ đỉnh giảm bốn lần. Đây là dấu hiệu khác với chuyển động đối lưu đều, vốn tạo thời gian tỷ lệ bậc nhất với khoảng cách.
4. Thời điểm phát hiện phụ thuộc ngưỡng
Đặt x = t/tₚ thì Δp/Δpₚ = exp(1−1/x)/x. Tại x bằng 0,10; 0,25; 0,50; 1; 2; 4, áp suất ở 100 m lần lượt bằng 0,053523; 8,637110; 31,910043; 43,370246; 35,752723; 22,953703 kPa. Chuỗi này cho thấy tại cùng một vị trí, đáp ứng không bắt đầu bằng một bước nhảy từ không lên cực đại.
Giả sử chỉ để phân tích toán học, ngưỡng phát hiện là 10 kPa trên nền ổn định. Giải exp(1−1/x)/x = 10/43,370246 trên nhánh tăng được x = 0,262952295 và thời gian 1,183285 giờ. Nếu ngưỡng là 1 kPa, thời điểm thành 0,674973 giờ; nếu là 20 kPa, thành 1,603883 giờ. D vẫn không đổi. Vì vậy lấy “lúc đầu nhìn thấy tín hiệu” làm thời gian truyền mà không ghi ngưỡng có thể tạo ra các ước lượng khác nhau từ cùng một mô hình.
Phương trình ngưỡng còn có nghiệm trên nhánh giảm: với 10 kPa, tín hiệu xuống dưới ngưỡng sau khoảng 48,335520 giờ. Một thuật toán tìm nghiệm phải xác định đang tìm lần vượt lên hay lần đi xuống. Nếu ngưỡng cao hơn đỉnh, hoàn toàn không có nghiệm dù hai điểm vẫn liên thông trong mô hình. Không phát hiện trong cửa sổ đo vì thế chưa đủ chứng minh không có liên thông.
5. Độ nhạy giúp nhận ra tham số bị lẫn
Từ tₚ = R²φμcₜ/(4k), biến thiên nhỏ thỏa Δtₚ/tₚ ≈ 2ΔR/R + Δφ/φ + Δμ/μ + Δcₜ/cₜ − Δk/k. Nếu khoảng cách thực là 105 m thay vì 100 m, tₚ bằng 4,961250 giờ, tăng 10,25%; xấp xỉ tuyến tính cho 10%. Độ nhạy bình phương theo khoảng cách cần được giữ khi tọa độ giếng hoặc vị trí lớp tiếp nhận chưa chắc chắn.
Giảm k một nửa làm tₚ tăng từ 4,5 lên 9 giờ nhưng không đổi đỉnh 43,370246 kPa của xung có cùng V, h, φ, cₜ. Ngược lại, tăng cₜ gấp đôi cũng cho thời gian 9 giờ nhưng đỉnh giảm còn 21,685123 kPa. Hai trường hợp giống nhau về thời gian cực đại mà khác về biên độ. Đây là lý do nên khớp toàn bộ đường đáp ứng và kiểm thể tích tác động, thay vì suy một tham số chỉ từ một thời điểm.
Ngay cả khi biết D, phương trình D = k/(φμcₜ) chỉ xác định một tổ hợp tham số. Muốn suy k riêng cần thông tin về các đại lượng còn lại. Nếu V thực vào vỉa khác thể tích đo ở mặt đất do tích trữ trong giếng, biên độ cũng không đủ để suy cₜ mà không xử lý ranh giới đó. Các tham số có thể bù trừ nhau trong một đường khớp nhìn rất đẹp.
6. Phạm vi của phép diễn giải
Một đợt ép kéo dài, thay đổi lưu lượng nhiều lần hoặc dừng không đồng bộ với đồng hồ đều làm đáp ứng khác xung lý tưởng. Khi thời lượng tác động đáng kể so với tₚ, phải dùng chồng chất theo thời gian. Bộ lọc tín hiệu và quán tính cảm biến cũng có thể dịch đỉnh quan sát; đó là đặc tính hệ đo, cần tách khỏi khuếch tán của vỉa.
Biên kín, biên áp suất cố định, tính dị hướng, nhiều lớp, khe nứt và dòng nhiều pha có thể đổi dạng đường cong. Do đó kết quả 4,5 giờ trong ví dụ không phải thời gian chờ thử giếng được khuyến nghị. Nó là kết quả của một hệ giả định có thể tính lại, dùng để kiểm sự khác biệt giữa thời gian không thứ nguyên, cực đại, ngưỡng phát hiện và vận chuyển chất lưu. Báo cáo một “độ trễ” chỉ có nghĩa kỹ thuật khi định nghĩa, lịch sử tác động và mô hình đi cùng nó được công bố.
Nguồn tham khảo
[1] Lecture notes in well-testing — sections 3.3–3.4, print date 23 April 2026 — Per Arne Slotte và Carl Fredrik Berg, NTNU. Truy cập ngày 09/11/2026.
[2] Groundwater Storage in Confined Aquifers — 5 Diffusion Equation — The Groundwater Project. Truy cập ngày 09/11/2026.
[3] Radial Flow Pressure Model — ORION documentation — National Risk Assessment Partnership, ORION. Truy cập ngày 09/11/2026.
Rà soát nội dung: 09/11/2026.
