Đường cong giữ nước của đất phản ánh điều gì?

Đường cong giữ nước phản ánh mối quan hệ giữa thế ma dẫn và lượng nước liên kết trong cấu trúc lỗ rỗng của đất không bão hòa. Đường quan hệ này giúp đánh giá khả năng tích trữ và giải phóng nước của đất trong các điều kiện khí hậu và tưới tiêu khác nhau.

Dung tích ẩm đồng ruộng và điểm héo vĩnh viễn khác nhau thế nào?

Dung tích ẩm đồng ruộng là lượng nước đất giữ lại sau khi nước trọng lực tự do đã thoát hết, thường ở thế ma dẫn từ -10 kPa đến -33 kPa. Điểm héo vĩnh viễn là ngưỡng thế ma dẫn -1500 kPa mà tại đó rễ cây không thể thắng được lực hút của đất để lấy nước, khiến thực vật bị khô héo không thể phục hồi.

Hiện tượng trễ trên đường cong giữ nước bắt nguồn từ đâu?

Hiện tượng trễ giữa quá trình làm khô và làm ẩm sinh ra do hiệu ứng cổ chai lọ của lỗ rỗng, sự chênh lệch góc thấm ướt tiếp xúc của màng nước và bọt khí bị bẫy trong đất. Do đó, ở cùng một mức áp suất ma trận, độ ẩm của đất trong quá trình giải hấp luôn lớn hơn quá trình hấp phụ cấp ẩm.

Mô hình van Genuchten (1980) có vai trò gì trong vật lý đất?

Mô hình van Genuchten cung cấp phương trình giải tích dạng đóng để tham số hóa đường cong giữ nước và dự báo độ dẫn thủy lực không bão hòa của đất. Phương trình này là thành phần cốt lõi trong các phần mềm mô phỏng dòng thấm nước và chuyển hóa chất ô nhiễm trong tầng đất ngầm.

Đường cong giữ nước là gì? Bản chất và mô hình toán

Đường cong giữ nước (soil water retention curve, SWRC hoặc soil water characteristic curve, SWCC) là hàm số thể hiện mối quan hệ cân bằng giữa hàm lượng nước chứa trong đất và thế ma dẫn của nước trong môi trường đất không bão hòa. Đường quan hệ này phản ánh năng lượng liên kết giữa các phân tử nước với bề mặt khoáng sét và hệ thống vi quản mao dẫn, đóng vai trò nền tảng trong vật lý đất, thủy văn học ngầm và địa kỹ thuật công trình.

Bản chất vật lý và năng lượng liên kết nước trong đất

Trong môi trường đất không bão hòa, nước chịu tác động đồng thời của trọng lực, lực hấp phụ phân tử bề mặt và lực căng bề mặt tại mặt phân cách pha khí và pha lỏng. Tổng năng lượng tự do trên một đơn vị thể tích nước đất so với trạng thái nước tự do ở cùng nhiệt độ và áp suất khí quyển được xác định thông qua khái niệm thế năng nước của đất (soil water potential).

Thành phần chiếm ưu thế quyết định khả năng giữ nước trong đất không bão hòa là thế ma dẫn hay áp suất ma trận (matric potential), ký hiệu là $\psi_m$ hoặc cột áp hút $h$ . Do nước bị giữ lại bởi lực căng mao dẫn và lực dính bám của hạt đất, thế ma dẫn luôn mang giá trị âm. Khi biểu diễn dưới dạng cột áp hút nước (suction head) hoặc độ hút ẩm ma dẫn, giá trị này được lấy dương: $h = -\psi_m / (\rho_w g)$ , trong đó $\rho_w$ là khối lượng riêng của nước và $g$ là gia tốc trọng trường.

Mối quan hệ giữa kích thước lỗ rỗng và độ hút ẩm mao dẫn được quy định bởi phương trình Young và Laplace. Đối với một lỗ rỗng hình trụ có bán kính hiệu dụng $r$ , áp suất hút ma dẫn giới hạn để đẩy nước ra khỏi lỗ rỗng được xác định theo biểu thức:

h = \frac{2 \sigma \cos \gamma}{\rho_w g r}

Trong công thức này, $\sigma$ là hệ số căng bề mặt của nước tiếp xúc với không khí, và $\gamma$ là góc thấm ướt tiếp xúc giữa màng nước và thành hạt khoáng. Biểu thức chỉ ra rằng lỗ rỗng có kích thước càng nhỏ thì lực hút giữ nước càng lớn, đòi hỏi áp suất ma trận âm càng cao để làm rỗng thể tích nước bên trong.

Các điểm ẩm đặc trưng trên đường cong giữ nước

Trên đường cong giữ nước, có một số mốc trạng thái độ ẩm quan trọng được tiêu chuẩn hóa trong khoa học đất và nông học để đánh giá khả năng cấp nước cho cây trồng cũng như tính chất cơ lý của đất.

Điểm ẩm đặc trưng	Ký hiệu	Thế ma dẫn quy ước	Ý nghĩa vật lý và sinh học
Độ ẩm bão hòa	$\theta_s$	0 kPa	Toàn bộ hệ thống lỗ rỗng trong đất chứa đầy nước, không còn pha khí.
Dung tích ẩm đồng ruộng	$\theta_{FC}$	-10 kPa đến -33 kPa	Lượng nước đất giữ lại được sau khi nước trọng lực tự do đã thoát hết.
Điểm héo vĩnh viễn	$\theta_{PWP}$	-1500 kPa	Ngưỡng hút nước mà rễ cây trồng không thể thắng được lực giữ của đất.
Độ ẩm dư	$\theta_r$	Hút ẩm cực đại	Lượng nước liên kết màng mỏng khó di chuyển ở điều kiện áp suất hút rất cao.

Theo tiêu chuẩn quốc gia TCVN 6651:2021 , dung tích ẩm đồng ruộng đối với đất cát thường tương ứng với thế ma dẫn -10 kPa, trong khi đối với đất thịt mịn hoặc đất sét, giá trị này nằm ở mức -33 kPa. Theo quy ước nông học quốc tế (Vereecken và cs., 2010), điểm héo vĩnh viễn của thực vật tương ứng với thế ma dẫn -1500 kPa (tương đương 15 bar ). Khoảng chênh lệch giữa dung tích ẩm đồng ruộng và điểm héo vĩnh viễn tạo nên dải nước hữu dụng cho cây trồng (plant available water), quyết định chu kỳ tưới tiêu và khả năng chống chịu hạn hán của thảm thực vật.

Các mô hình toán học giải tích mô tả đường cong

Để đưa dữ liệu thực nghiệm vào các mô hình mô phỏng chuyển động của dòng thấm và chất tan (như mô hình HYDRUS hoặc Richards), các nhà khoa học đã xây dựng nhiều phương trình dạng đóng (closed-form equations).

Mô hình van Genuchten (1980)

Mô hình được áp dụng rộng rãi nhất trong thủy văn học và vật lý đất hiện đại là phương trình của nhà khoa học van Genuchten (1980). Phương trình biểu diễn mối quan hệ giữa hàm lượng nước không thứ nguyên $\Theta$ và cột áp hút $h$ theo dạng:

\Theta = \frac{\theta - \theta_r}{\theta_s - \theta_r} = \left[ \frac{1}{1 + (\alpha h)^n} \right]^m

Trong hệ thức này, $\theta$ là hàm lượng nước thể tích tại cột áp hút $h$ ; $\theta_s$ là hàm lượng nước bão hòa; $\theta_r$ là hàm lượng nước dư; $\alpha$ là thông số thực nghiệm liên quan nghịch đảo tới điểm áp suất khí xâm nhập (air-entry value); $n$ là hệ số phân bố kích thước lỗ rỗng; và $m$ là tham số thực nghiệm thường được liên kết với $n$ theo điều kiện $m = 1 - 1/n$ khi kết hợp với mô hình độ dẫn thủy lực của Mualem.

Mô hình Fredlund và Xing (1994)

Trong lĩnh vực cơ học đất không bão hòa và địa kỹ thuật công trình, mô hình Fredlund và Xing (1994) được đánh giá cao nhờ khả năng biểu diễn chính xác trên toàn bộ dải thế ma dẫn từ 0 đến 1 000 000 kPa (tương ứng trạng thái sấy khô hoàn toàn). Phương trình có dạng:

\theta(h) = C(h) \frac{\theta_s}{\left[ \ln\left( e + (h/a)^n \right) \right]^m}

Trong đó, $e$ là cơ số logarit tự nhiên; $a$ , $n$ , $m$ là các tham số điều chỉnh thực nghiệm đặc trưng cho đất; và $C(h)$ là hàm số hiệu chỉnh đưa hàm lượng ẩm về đúng giá trị không tại mốc 1 000 000 kPa.

Mô hình phân bố hai bậc rỗng của Durner (1994)

Đối với đất có cấu trúc phân tầng tự nhiên, đất canh tác có nhiều khe nứt sinh học do giun đất hoặc rễ cây mục nát, mô hình đơn đỉnh thường đánh giá thấp lượng nước thoát ở áp suất thấp. Durner (1994) đã đề xuất mô hình hai bậc rỗng (bimodal) biểu diễn đường cong giữ nước dưới dạng tổng hợp tuyến tính của hai hàm thành phần:

\Theta = \sum_{i=1}^{k} w_i \left[ \frac{1}{1 + (\alpha_i h)^{n_i}} \right]^{m_i}

Với $w_i$ là trọng số tỷ lệ của từng hệ thống lỗ rỗng thỏa mãn tổng các trọng số bằng một. Mô hình này phản ánh rõ nét sự khác biệt giữa hệ thống lỗ rỗng đại thể liên kết nhanh và hệ thống mao dẫn vi mô trong ma trận hạt đất.

Hiện tượng trễ hút ẩm và giải hấp

Một đặc tính phức tạp của đường cong giữ nước là hiện tượng trễ (hysteresis). Khi tiến hành đo đạc quá trình làm khô đất (desorption hoặc drying curve) và quá trình làm ẩm đất (wetting hoặc sorption curve), hai đường cong thực nghiệm không trùng nhau. Tại cùng một mức thế ma dẫn, hàm lượng nước trên đường làm khô luôn cao hơn hàm lượng nước trên đường làm ẩm.

Hiện tượng trễ bắt nguồn từ ba cơ chế vật lý cơ bản:

Hiệu ứng cổ chai lọ (ink-bottle effect): Lỗ rỗng trong đất có hình dạng bất đối xứng với phần thân rộng nhưng cổ hẹp. Khi làm khô, nước chỉ thoát ra khi lực hút vượt qua kích thước cổ hẹp; ngược lại khi làm ẩm, nước chỉ làm đầy khoang thân khi mực nước đạt đến bán kính lớn hơn.
Sự khác biệt về góc thấm ướt tiếp xúc: Góc tiếp xúc của màng nước đang tiến (advancing contact angle) trong quá trình làm ẩm luôn lớn hơn góc tiếp xúc đang lùi (receding contact angle) trong quá trình giải hấp.
Khí bị bẫy cục bộ: Trong quá trình cấp ẩm nhanh, một phần bọt khí bị giữ lại trong các khoang rỗng cụt, làm giảm độ bão hòa thực tế so với mẫu đất bão hòa ban đầu.

Phương pháp thực nghiệm xác định đặc tính giữ nước

Các phương pháp phòng thử nghiệm được chuẩn hóa trong TCVN 6651:2021 bao gồm bàn hút cát (từ 0 đến -50 kPa) và buồng nén màng áp lực (từ -5 kPa đến -1500 kPa). Để mở rộng dải đo ngoài hiện trường hoặc ở vùng áp suất hút cực cao, các phương pháp bổ trợ khác cũng được áp dụng:

Áp kế tensiometer: Đo trực tiếp thế ma dẫn tại hiện trường hoặc trong phòng thí nghiệm. Thiết bị cấu tạo gồm ống sứ xốp nối với cảm biến áp suất. Do hiện tượng xâm thực sôi bọt khí chân không của nước ở áp suất âm lớn, tensiometer chỉ hoạt động chính xác trong dải từ 0 đến 85 kPa.
Bàn hút cát và buồng nén màng áp lực (pressure plate extractor): Sử dụng màng gốm có ngưỡng khí xâm nhập cao để tác động áp suất khí nén lên mẫu đất bão hòa, ép nước thoát ra ngoài cho đến khi đạt cân bằng thế năng. Phương pháp này áp dụng cho dải áp suất từ 10 kPa đến 1500 kPa.
Thiết bị gương lạnh đo điểm sương: Đo độ ẩm tương đối của pha khí cân bằng với mẫu đất để suy ra tổng thế năng nước thông qua phương trình Kelvin, đặc biệt hiệu quả ở dải hút ẩm cực cao.

Hàm truyền thổ nhưỡng và ước tính tham số

Việc đo đạc trực tiếp đường cong giữ nước trong phòng thí nghiệm đòi hỏi thời gian kéo dài từ vài tuần đến vài tháng cho mỗi mẫu đất. Để khắc phục rào cản chi phí và thời gian trong các dự án quy mô lớn, các nhà khoa học sử dụng hàm truyền thổ nhưỡng (pedotransfer functions, PTF).

Tổng quan nghiên cứu của Vereecken và cs. (2010) chỉ ra rằng các hàm truyền thổ nhưỡng khai thác các chỉ tiêu cơ lý cơ bản dễ đo lường—như tỷ lệ thành phần hạt cát, limon, sét, khối lượng thể tích khô và hàm lượng carbon hữu cơ—để dự báo các tham số của mô hình van Genuchten. Các phần mềm nổi tiếng như Rosetta của Bộ Nông nghiệp Hoa Kỳ (USDA) sử dụng mạng nơ-ron nhân tạo huấn luyện trên hàng nghìn bộ dữ liệu đất toàn cầu để đưa ra ước tính tham số thủy lực với độ tin cậy chấp nhận được cho phân tích mô phỏng diện rộng.

Ứng dụng thực tiễn trong nông nghiệp và địa kỹ thuật tại Việt Nam

Tại Việt Nam, nghiên cứu đường cong giữ nước có vai trò đặc biệt sống còn trong hai lĩnh vực chính:

Quy hoạch tưới tiết kiệm và quản lý đất nông nghiệp: Tại các vùng khô hạn như Ninh Thuận, Bình Thuận và vùng đất đỏ bazan Tây Nguyên, xác định chính xác đường cong giữ nước giúp kỹ sư xác định ngưỡng độ ẩm tối ưu để kích hoạt hệ thống tưới nhỏ giọt. Duy trì thế ma dẫn của đất ở ngưỡng ẩm tối ưu giúp rễ cây hút nước dễ dàng mà không gây thất thoát phân bón hay rửa trôi dinh dưỡng.
Đánh giá độ ổn định mái dốc và sạt lở đất: Ở khu vực miền núi phía Bắc và dọc tuyến đường Trường Sơn, các khối trượt mái dốc thường xảy ra trong tầng đất không bão hòa khi gặp mưa kéo dài. Sức kháng cắt của đất phụ thuộc trực tiếp vào lực hút ma dẫn thông qua ứng suất hữu hiệu Bishop. Khi nước mưa ngấm vào làm giảm thế ma dẫn về gần 0, lực dính biểu kiến mất đi khiến mái dốc sạt lở. Việc kết hợp đường cong giữ nước vào các mô hình ổn định mái dốc được các chuyên gia địa kỹ thuật khuyến nghị áp dụng trong các mô hình phân tích trượt lở nâng cao.

Hạn chế và thách thức nghiên cứu còn mở

Mặc dù lý thuyết về đường cong giữ nước đã hoàn thiện tương đối vững chắc, các nhà nghiên cứu vẫn đối mặt với một số giới hạn phương pháp luận:

Hiện tượng co ngót và trương nở: Các mô hình kinh điển như van Genuchten ngầm định thể tích khung đất là cố định. Đối với các loại đất sét trương nở giàu khoáng montmorillonite (phổ biến ở một số vùng trũng tại Đồng bằng sông Cửu Long), thể tích lỗ rỗng biến đổi mạnh theo hàm lượng nước, đòi hỏi các mô hình tích hợp biến dạng thể tích phức tạp hơn.
Ảnh hưởng của nhiệt độ và nồng độ muối: Hệ số căng bề mặt và độ nhớt của nước biến thiên theo nhiệt độ môi trường, trong khi thế thẩm thấu do nồng độ ion muối hòa tan tác động trực tiếp đến khả năng hút nước của rễ cây, đòi hỏi phân tách rạch ròi giữa thế ma dẫn và tổng thế năng nước.
Biến động không gian: Cấu trúc lỗ rỗng của đất có tính dị hướng và biến động không gian rất lớn trong tự nhiên, khiến việc chuyển đổi từ dữ liệu mẫu đất kích thước nhỏ trong phòng thí nghiệm sang quy mô lưu vực thủy văn luôn đi kèm sai số không nhỏ.