1. Công thức tuyến tính cần đúng hai mốc
Quan hệ T(z) = T₀ + Gz rất tiện để ước tính nhiệt độ theo độ sâu. SLB mô tả gradient địa nhiệt là mức tăng nhiệt độ trên một đơn vị độ sâu và nêu phép tính tuyến tính này khi dùng một gradient đại diện [1]. Tuy nhiên, T₀ và z phải cùng một mốc không gian, còn nhiệt độ dùng khớp phải đại diện trạng thái đang nghiên cứu. Một công thức đúng về đại số vẫn cho kết quả sai nếu lấy nhiệt độ ngoài trời làm mốc cho chiều sâu tính từ đáy biển.
Trong bài này, z dương xuống dưới, T dùng °C và G dùng K/m hoặc °C/m. Độ chênh một kelvin bằng độ chênh một độ Celsius, nên có thể dùng G = 0,030 K/m cùng T₀ = 24°C. Khi z = 2.800 m, mô hình cho T = 108°C. Giá trị 30 K/km là giả định của ví dụ, không phải gradient được xác nhận cho một mỏ Việt Nam.
Cần phân biệt độ sâu theo đường giếng, độ sâu đứng và mốc cao độ. Một nhiệt độ đo ở đáy khi đang tuần hoàn cũng không tự bằng nhiệt độ đá nguyên trạng tại cùng vị trí. Tệp dữ liệu hữu ích phải ghi cả tọa độ, trạng thái giếng và thời điểm đo, thay vì chỉ có hai cột “độ sâu” và “nhiệt độ”.
2. Độ sâu theo đường giếng làm tăng giả tạo nhiệt độ dự báo
Xét riêng một giếng có 1.000 m đầu thẳng đứng, tiếp theo 800 m nghiêng 30° so với phương đứng và 1.200 m nghiêng 60°. Bỏ qua các đoạn chuyển tiếp cong để minh họa, chiều dài đường giếng là 3.000 m, nhưng độ sâu đứng là z = 1.000 + 800cos30° + 1.200cos60° = 2.292,820323 m. Với T₀ = 24°C và G = 0,030 K/m, nhiệt độ nền dự báo bằng 92,784610°C.
Nếu nhầm dùng 3.000 m, kết quả thành 114°C, cao hơn 21,215390°C. Sai khác đến từ hình học, không phải từ sai số cảm biến. Trong một công trình đo gradient thực tế, USGS cũng phải hiệu chỉnh do giếng lệch khỏi phương đứng [2]. Với đường giếng cong, độ sâu đứng được tính từ khảo sát quỹ đạo; không nên áp một góc trung bình tùy ý cho cả chiều dài.
Một kịch bản ngoài khơi khác có đáy biển sâu 200 m, nhiệt độ tại đáy biển giả định 4°C và điểm xét nằm 2.500 m dưới đáy biển. Nếu G = 0,030 K/m chỉ áp dụng trong trầm tích, nhiệt độ nền tại điểm đó là 79°C. Dùng nhiệt độ mặt biển 24°C cộng gradient cho toàn bộ 2.700 m cho 105°C, sai 26°C. Cột nước cần mô tả nhiệt riêng; gradient trong đá không được tự kéo qua cột nước.
3. Nhiệt độ đo sau tuần hoàn cần một mô hình phục hồi
Colorado Geological Survey lưu ý khoan và bùn tuần hoàn làm xáo trộn nhiệt độ, còn mức làm nguội và thời gian phục hồi thay đổi theo từng giếng [4]. USGS xây dựng các hiệu chỉnh BHT thực nghiệm dựa trên dữ liệu và hiệu chuẩn theo vùng [3]. Hai nhận xét đó không cung cấp một số cộng bù cố định cho mọi trường hợp; chúng nhắc rằng trạng thái đo là một phần của bài toán.
Để thấy độ nhạy của ngoại suy, tự đặt một mô hình gộp thử nghiệm T_m(t) = T∞ − A exp(−t/τ). T∞ là nhiệt độ giới hạn, A có đơn vị K, τ là hằng số thời gian dùng giờ. Đây chỉ là hàm giả định cho ví dụ số, không phải phương pháp hiệu chỉnh USGS và không khẳng định quá trình dẫn nhiệt xuyên tâm của giếng luôn là một hàm mũ đơn.
Giả sử cùng một điểm sâu 2.800 m có ba số đo tại 6, 12 và 18 giờ sau một mốc dừng tuần hoàn: 94; 101; 104,5°C. Hai mức tăng liên tiếp là 7 và 3,5 K, nên tỷ số r = 3,5/7 = 0,5. Vì các lần đo cách nhau Δ = 6 giờ, r = exp(−Δ/τ), cho τ = −6/ln0,5 = 8,656170 giờ. Nhiệt độ giới hạn suy ra T∞ = 94 + 7/(1−r) = 108°C.
Từ đó A = 28 K tại t = 0 và mô hình dự báo T_m(36 giờ) = 107,5625°C. Phép tính không biến 36 giờ thành thời gian chờ đo tiêu chuẩn. Nếu dùng ngay số 94°C làm nhiệt độ nền, gradient từ mốc 24°C chỉ là 25 K/km; dùng số cuối 104,5°C cho 28,75 K/km; dùng giới hạn mô hình 108°C cho 30 K/km. Cần ghi rõ kết quả nào là đo và kết quả nào là ngoại suy.
4. Ba điểm khớp được chưa chứng minh mô hình đúng
Giữ hai số đầu 94 và 101°C, chỉ đổi số cuối từ 104,5 lên 105°C. Khi ấy r = 4/7 và nhiệt độ giới hạn suy ra tăng thành 110,333333°C. Nếu số cuối là 104°C, r = 3/7 và giới hạn giảm còn 106,25°C. Nhiễu ±0,5°C ở một điểm đã tạo khoảng ngoại suy rộng hơn nhiều. Công thức có mẫu 1−r, nên càng gần r = 1, phép suy giới hạn càng nhạy.
Một chuỗi có mức tăng không giảm đều, có dòng mới vào giếng hoặc thay đổi vị trí đầu đo có thể không phù hợp với mô hình mũ. Khi r không nằm trong khoảng 0 đến 1, không nên ép công thức tạo một τ dương bằng cách đổi dấu dữ liệu. Cần kiểm lịch sử nhiệt, nhiều thời điểm hơn và cơ chế truyền nhiệt trước khi chấp nhận giới hạn.
Các hiệu chỉnh thực nghiệm theo khu vực trong nghiên cứu USGS [3] dựa vào một tập dữ liệu riêng. Chuyển chúng sang một khu vực khác cần đánh giá phạm vi độ sâu, lịch sử tuần hoàn và đặc điểm địa chất. Giá trị tính được bằng một công thức đã xuất bản vẫn có thể là ngoại suy ngoài phạm vi hiệu chuẩn; tên phương pháp không tự bảo đảm tính phù hợp.
5. Cùng nhiệt độ đáy có thể có đường nhiệt khác nhau
Một gradient trung bình từ hai mốc không quyết định nhiệt độ ở mọi độ sâu ở giữa. Colorado Geological Survey mô tả gradient thay đổi theo tính dẫn nhiệt của đá và cơ chế dòng nước [4]. Xét mô hình dẫn nhiệt một chiều ổn định, không có nguồn nhiệt trong các lớp và không có đối lưu. Gọi q↑ là mật độ dòng nhiệt hướng lên, đơn vị W/m². Với z dương xuống, q↑ = λdT/dz; λ là hệ số dẫn nhiệt W/(m·K).
Giả sử lớp trên dày 1.000 m có λ₁ = 2 W/(m·K), lớp dưới dày 1.800 m có λ₂ = 3 W/(m·K). Muốn giữ T₀ = 24°C và T_đáy = 108°C, dòng nhiệt phải thỏa 84 = q↑(1.000/2 + 1.800/3). Suy ra q↑ = 0,076363636 W/m², tức 76,363636 mW/m².
Gradient trong lớp trên bằng 38,181818 K/km, còn lớp dưới bằng 25,454545 K/km. Nhiệt độ tại 1.000 m là 62,181818°C và tại 2.000 m là 87,636364°C. Đường tuyến tính 30 K/km lại cho lần lượt 54 và 84°C. Hai mô hình gặp nhau tại đáy 108°C nhưng lệch đáng kể ở giữa; dùng một BHT duy nhất không phân biệt được chúng.
Đồng thời, gradient lớn không tự đồng nghĩa dòng nhiệt lớn hơn nếu λ khác. Dòng nhiệt có đơn vị công suất trên diện tích, còn gradient là độ chênh nhiệt trên chiều dài. Không thể thay một đại lượng cho đại lượng kia trong đánh giá nhiệt chỉ vì cả hai cùng tăng ở một bộ dữ liệu.
6. Độ không đảm bảo và mục đích sử dụng
Với G = (T∞−T₀)/z, độ nhạy theo nhiệt độ đáy là 1/z, theo nhiệt độ mốc là −1/z và theo độ sâu là −G/z. Nếu các độ không đảm bảo chuẩn độc lập lần lượt là 2 K, 1 K và 10 m, tại z = 2.800 m và G = 0,030 K/m, ta có uG = √(2²+1²+0,030²×10²)/2.800 = 0,000805751 K/m, tức 0,805751 K/km. Nếu các nhiệt độ có sai lệch hiệu chuẩn chung, cần thêm hiệp phương sai trước khi dùng phép cộng bình phương này.
Đây mới là lan truyền độ không đảm bảo của các đầu vào trong mô hình đã chọn. Sai lệch do dùng sai trạng thái nhiệt, thiếu lớp dẫn nhiệt hoặc có dòng đối lưu có thể lớn hơn và không được khắc phục bằng cách ghi nhiều chữ số thập phân. USGS ghi nhận dòng nước có thể che lấp một phần gradient tự nhiên [2], nên một đoạn bất thường cần đọc cùng dữ liệu dòng và cấu tạo giếng.
Nhiệt độ nền địa chất còn khác nhiệt độ chất lưu trong giếng đang sản xuất hoặc bơm ép. Muốn đánh giá thiết bị tại một chế độ vận hành, cần mô hình truyền nhiệt và lịch sử dòng phù hợp. Mô hình gradient tuyến tính vẫn hữu ích ở giai đoạn đầu, miễn là báo cáo tách rõ mốc nhiệt, độ sâu đứng, số đo, phần hiệu chỉnh và phạm vi độ sâu mà dữ liệu thực sự hỗ trợ.
Nguồn tham khảo
[1] geothermal gradient — Energy Glossary — SLB. Truy cập ngày 09/11/2026.
[2] Geophysical logging of bedrock wells for geothermal gradient characterization in New Hampshire, 2013 — Data Series 823 (2014) — James R. Degnan, Gregory Barker, Neil Olson và Leland Wilder, U.S. Geological Survey. Truy cập ngày 09/11/2026.
[3] New method for correcting bottomhole temperatures acquired from wireline logging measurements and calibrated for the onshore Gulf of Mexico Basin, U.S.A. — Open-File Report 2019-1143 — Lauri A. Burke, Ofori N. Pearson và Scott A. Kinney, U.S. Geological Survey. Truy cập ngày 09/11/2026.
[4] Geothermal — section Geothermal Gradient — Colorado Geological Survey. Truy cập ngày 09/11/2026.
Rà soát nội dung: 09/11/2026.
Quan hệ T(z) = T₀ + Gz rất tiện để ước tính nhiệt độ theo độ sâu. SLB mô tả gradient địa nhiệt là mức tăng nhiệt độ trên một đơn vị độ sâu và nêu phép tính tuyến tính này khi dùng một gradient đại diện [1]. Tuy nhiên, T₀ và z phải cùng một mốc không gian, còn nhiệt độ dùng khớp phải đại diện trạng thái đang nghiên cứu. Một công thức đúng về đại số vẫn cho kết quả sai nếu lấy nhiệt độ ngoài trời làm mốc cho chiều sâu tính từ đáy biển.
Trong bài này, z dương xuống dưới, T dùng °C và G dùng K/m hoặc °C/m. Độ chênh một kelvin bằng độ chênh một độ Celsius, nên có thể dùng G = 0,030 K/m cùng T₀ = 24°C. Khi z = 2.800 m, mô hình cho T = 108°C. Giá trị 30 K/km là giả định của ví dụ, không phải gradient được xác nhận cho một mỏ Việt Nam.
Cần phân biệt độ sâu theo đường giếng, độ sâu đứng và mốc cao độ. Một nhiệt độ đo ở đáy khi đang tuần hoàn cũng không tự bằng nhiệt độ đá nguyên trạng tại cùng vị trí. Tệp dữ liệu hữu ích phải ghi cả tọa độ, trạng thái giếng và thời điểm đo, thay vì chỉ có hai cột “độ sâu” và “nhiệt độ”.
2. Độ sâu theo đường giếng làm tăng giả tạo nhiệt độ dự báo
Xét riêng một giếng có 1.000 m đầu thẳng đứng, tiếp theo 800 m nghiêng 30° so với phương đứng và 1.200 m nghiêng 60°. Bỏ qua các đoạn chuyển tiếp cong để minh họa, chiều dài đường giếng là 3.000 m, nhưng độ sâu đứng là z = 1.000 + 800cos30° + 1.200cos60° = 2.292,820323 m. Với T₀ = 24°C và G = 0,030 K/m, nhiệt độ nền dự báo bằng 92,784610°C.
Nếu nhầm dùng 3.000 m, kết quả thành 114°C, cao hơn 21,215390°C. Sai khác đến từ hình học, không phải từ sai số cảm biến. Trong một công trình đo gradient thực tế, USGS cũng phải hiệu chỉnh do giếng lệch khỏi phương đứng [2]. Với đường giếng cong, độ sâu đứng được tính từ khảo sát quỹ đạo; không nên áp một góc trung bình tùy ý cho cả chiều dài.
Một kịch bản ngoài khơi khác có đáy biển sâu 200 m, nhiệt độ tại đáy biển giả định 4°C và điểm xét nằm 2.500 m dưới đáy biển. Nếu G = 0,030 K/m chỉ áp dụng trong trầm tích, nhiệt độ nền tại điểm đó là 79°C. Dùng nhiệt độ mặt biển 24°C cộng gradient cho toàn bộ 2.700 m cho 105°C, sai 26°C. Cột nước cần mô tả nhiệt riêng; gradient trong đá không được tự kéo qua cột nước.
3. Nhiệt độ đo sau tuần hoàn cần một mô hình phục hồi
Colorado Geological Survey lưu ý khoan và bùn tuần hoàn làm xáo trộn nhiệt độ, còn mức làm nguội và thời gian phục hồi thay đổi theo từng giếng [4]. USGS xây dựng các hiệu chỉnh BHT thực nghiệm dựa trên dữ liệu và hiệu chuẩn theo vùng [3]. Hai nhận xét đó không cung cấp một số cộng bù cố định cho mọi trường hợp; chúng nhắc rằng trạng thái đo là một phần của bài toán.
Để thấy độ nhạy của ngoại suy, tự đặt một mô hình gộp thử nghiệm T_m(t) = T∞ − A exp(−t/τ). T∞ là nhiệt độ giới hạn, A có đơn vị K, τ là hằng số thời gian dùng giờ. Đây chỉ là hàm giả định cho ví dụ số, không phải phương pháp hiệu chỉnh USGS và không khẳng định quá trình dẫn nhiệt xuyên tâm của giếng luôn là một hàm mũ đơn.
Giả sử cùng một điểm sâu 2.800 m có ba số đo tại 6, 12 và 18 giờ sau một mốc dừng tuần hoàn: 94; 101; 104,5°C. Hai mức tăng liên tiếp là 7 và 3,5 K, nên tỷ số r = 3,5/7 = 0,5. Vì các lần đo cách nhau Δ = 6 giờ, r = exp(−Δ/τ), cho τ = −6/ln0,5 = 8,656170 giờ. Nhiệt độ giới hạn suy ra T∞ = 94 + 7/(1−r) = 108°C.
Từ đó A = 28 K tại t = 0 và mô hình dự báo T_m(36 giờ) = 107,5625°C. Phép tính không biến 36 giờ thành thời gian chờ đo tiêu chuẩn. Nếu dùng ngay số 94°C làm nhiệt độ nền, gradient từ mốc 24°C chỉ là 25 K/km; dùng số cuối 104,5°C cho 28,75 K/km; dùng giới hạn mô hình 108°C cho 30 K/km. Cần ghi rõ kết quả nào là đo và kết quả nào là ngoại suy.
4. Ba điểm khớp được chưa chứng minh mô hình đúng
Giữ hai số đầu 94 và 101°C, chỉ đổi số cuối từ 104,5 lên 105°C. Khi ấy r = 4/7 và nhiệt độ giới hạn suy ra tăng thành 110,333333°C. Nếu số cuối là 104°C, r = 3/7 và giới hạn giảm còn 106,25°C. Nhiễu ±0,5°C ở một điểm đã tạo khoảng ngoại suy rộng hơn nhiều. Công thức có mẫu 1−r, nên càng gần r = 1, phép suy giới hạn càng nhạy.
Một chuỗi có mức tăng không giảm đều, có dòng mới vào giếng hoặc thay đổi vị trí đầu đo có thể không phù hợp với mô hình mũ. Khi r không nằm trong khoảng 0 đến 1, không nên ép công thức tạo một τ dương bằng cách đổi dấu dữ liệu. Cần kiểm lịch sử nhiệt, nhiều thời điểm hơn và cơ chế truyền nhiệt trước khi chấp nhận giới hạn.
Các hiệu chỉnh thực nghiệm theo khu vực trong nghiên cứu USGS [3] dựa vào một tập dữ liệu riêng. Chuyển chúng sang một khu vực khác cần đánh giá phạm vi độ sâu, lịch sử tuần hoàn và đặc điểm địa chất. Giá trị tính được bằng một công thức đã xuất bản vẫn có thể là ngoại suy ngoài phạm vi hiệu chuẩn; tên phương pháp không tự bảo đảm tính phù hợp.
5. Cùng nhiệt độ đáy có thể có đường nhiệt khác nhau
Một gradient trung bình từ hai mốc không quyết định nhiệt độ ở mọi độ sâu ở giữa. Colorado Geological Survey mô tả gradient thay đổi theo tính dẫn nhiệt của đá và cơ chế dòng nước [4]. Xét mô hình dẫn nhiệt một chiều ổn định, không có nguồn nhiệt trong các lớp và không có đối lưu. Gọi q↑ là mật độ dòng nhiệt hướng lên, đơn vị W/m². Với z dương xuống, q↑ = λdT/dz; λ là hệ số dẫn nhiệt W/(m·K).
Giả sử lớp trên dày 1.000 m có λ₁ = 2 W/(m·K), lớp dưới dày 1.800 m có λ₂ = 3 W/(m·K). Muốn giữ T₀ = 24°C và T_đáy = 108°C, dòng nhiệt phải thỏa 84 = q↑(1.000/2 + 1.800/3). Suy ra q↑ = 0,076363636 W/m², tức 76,363636 mW/m².
Gradient trong lớp trên bằng 38,181818 K/km, còn lớp dưới bằng 25,454545 K/km. Nhiệt độ tại 1.000 m là 62,181818°C và tại 2.000 m là 87,636364°C. Đường tuyến tính 30 K/km lại cho lần lượt 54 và 84°C. Hai mô hình gặp nhau tại đáy 108°C nhưng lệch đáng kể ở giữa; dùng một BHT duy nhất không phân biệt được chúng.
Đồng thời, gradient lớn không tự đồng nghĩa dòng nhiệt lớn hơn nếu λ khác. Dòng nhiệt có đơn vị công suất trên diện tích, còn gradient là độ chênh nhiệt trên chiều dài. Không thể thay một đại lượng cho đại lượng kia trong đánh giá nhiệt chỉ vì cả hai cùng tăng ở một bộ dữ liệu.
6. Độ không đảm bảo và mục đích sử dụng
Với G = (T∞−T₀)/z, độ nhạy theo nhiệt độ đáy là 1/z, theo nhiệt độ mốc là −1/z và theo độ sâu là −G/z. Nếu các độ không đảm bảo chuẩn độc lập lần lượt là 2 K, 1 K và 10 m, tại z = 2.800 m và G = 0,030 K/m, ta có uG = √(2²+1²+0,030²×10²)/2.800 = 0,000805751 K/m, tức 0,805751 K/km. Nếu các nhiệt độ có sai lệch hiệu chuẩn chung, cần thêm hiệp phương sai trước khi dùng phép cộng bình phương này.
Đây mới là lan truyền độ không đảm bảo của các đầu vào trong mô hình đã chọn. Sai lệch do dùng sai trạng thái nhiệt, thiếu lớp dẫn nhiệt hoặc có dòng đối lưu có thể lớn hơn và không được khắc phục bằng cách ghi nhiều chữ số thập phân. USGS ghi nhận dòng nước có thể che lấp một phần gradient tự nhiên [2], nên một đoạn bất thường cần đọc cùng dữ liệu dòng và cấu tạo giếng.
Nhiệt độ nền địa chất còn khác nhiệt độ chất lưu trong giếng đang sản xuất hoặc bơm ép. Muốn đánh giá thiết bị tại một chế độ vận hành, cần mô hình truyền nhiệt và lịch sử dòng phù hợp. Mô hình gradient tuyến tính vẫn hữu ích ở giai đoạn đầu, miễn là báo cáo tách rõ mốc nhiệt, độ sâu đứng, số đo, phần hiệu chỉnh và phạm vi độ sâu mà dữ liệu thực sự hỗ trợ.
Nguồn tham khảo
[1] geothermal gradient — Energy Glossary — SLB. Truy cập ngày 09/11/2026.
[2] Geophysical logging of bedrock wells for geothermal gradient characterization in New Hampshire, 2013 — Data Series 823 (2014) — James R. Degnan, Gregory Barker, Neil Olson và Leland Wilder, U.S. Geological Survey. Truy cập ngày 09/11/2026.
[3] New method for correcting bottomhole temperatures acquired from wireline logging measurements and calibrated for the onshore Gulf of Mexico Basin, U.S.A. — Open-File Report 2019-1143 — Lauri A. Burke, Ofori N. Pearson và Scott A. Kinney, U.S. Geological Survey. Truy cập ngày 09/11/2026.
[4] Geothermal — section Geothermal Gradient — Colorado Geological Survey. Truy cập ngày 09/11/2026.
Rà soát nội dung: 09/11/2026.
