Luận án tiến sĩ: Phương pháp phân tích ngược để ước tính tính chất thủy lực cho cát không bão hòa

Luận án tiến sĩ nghiên cứu inverse analysic to estimate hydraulic properties for unsaturated sand, phát triển phương pháp mới, đánh giá hiệu quả ứng dụng trong lĩnh vực tại Việt

• Tác giả: To Viet Nam
• Người hướng dẫn: Professor Min Tuk-Ki
• Trường học: University of Ulsan
• Chuyên ngành: Civil and Environmental Engineering
• Đề tài: Inverse Analysis to Estimate Hydraulic Properties for Unsaturated Sand
• Loại tài liệu: dissertation
• Năm xuất bản: 2013
• Địa điểm: Ulsan

Mục lục chi tiết

• Lời cảm ơn
• Tóm tắt
• Danh mục hình vẽ
• Danh mục bảng
• Danh mục ký hiệu
• 1. Chương 1: Giới thiệu
• 1.1. Bối cảnh nghiên cứu
• 1.2. Mục tiêu và đóng góp của nghiên cứu
• 1.2.1. Mục tiêu nghiên cứu
• 1.2.2. Những đóng góp mới của đề tài
• 1.3. Bố cục của luận án
• 2. Chương 2: Đặc tính thủy lực của đất không bão hòa
• 2.1. Đặc tính thủy lực của đất không bão hòa
• 2.2. Các nghiên cứu điển hình trước đây về đặc tính thủy lực của đất
• 2.3. Xác định các đặc tính thủy lực của đất không bão hòa
• 3. Chương 3: Phân tích nghịch đảo
• 3.1. Phân tích nghịch đảo đối với đặc tính thủy lực của đất không bão hòa
• 3.2. Mô hình hóa dòng chảy của nước
• 3.3. Tối ưu hóa thông số
• 4. Chương 4: Phân tích nghịch đảo để ước tính đặc tính thủy lực của cát đồng nhất
• 5. Chương 5: Phân tích nghịch đảo để ước tính đặc tính thủy lực của cát phân tầng
• 5.1. Vật liệu và phương pháp nghiên cứu
• 6. Chương 6: Kết luận
• Tài liệu tham khảo

Tóm tắt

I. Phương pháp phân tích ngược

Phương pháp phân tích ngược là một kỹ thuật quan trọng trong việc ước tính các tính chất thủy lực của đất không bão hòa. Phương pháp này sử dụng dữ liệu đầu vào từ các thí nghiệm thoát nước một chiều (1-D) để xác định các thông số thủy lực thông qua mô hình hóa và tối ưu hóa. Phân tích ngược cho phép kết hợp giữa đo lường, mô hình hóa và tối ưu hóa, giúp cải thiện độ chính xác và tốc độ xác định các tính chất thủy lực. Phương pháp này đã được áp dụng thành công trong các thí nghiệm thoát nước một bước và nhiều bước, cho thấy tiềm năng lớn trong việc ứng dụng vào thực tế.

1.1. Cơ sở lý thuyết

Phân tích ngược dựa trên việc tối ưu hóa hàm mục tiêu, được định nghĩa bởi sự khác biệt giữa dữ liệu quan sát và dữ liệu dự đoán. Một phương pháp tối ưu hóa phi tuyến kết hợp với mã phần tử hữu hạn được sử dụng để ước tính các thông số thủy lực trong mô hình van Genuchten. Các thông số này bao gồm áp suất hút đất và tốc độ thoát nước tại đáy cột cát. Kết quả cho thấy sự phù hợp cao giữa dữ liệu dự đoán và dữ liệu thực nghiệm, đặc biệt trong thí nghiệm thoát nước nhiều bước.

1.2. Ứng dụng thực tế

Phương pháp phân tích ngược đã được áp dụng để xác định các tính chất thủy lực của cát không bão hòa trong các thí nghiệm thoát nước một chiều. Kết quả cho thấy phương pháp này không chỉ đơn giản mà còn có độ chính xác cao, đặc biệt khi so sánh với các phương pháp truyền thống. Ứng dụng của phương pháp này trong thực tế giúp giảm thiểu thời gian và chi phí, đồng thời cải thiện độ tin cậy của các kết quả thí nghiệm.

II. Ước tính tính chất thủy lực

Ước tính tính chất thủy lực là một bước quan trọng trong việc hiểu và mô hình hóa hành vi của đất không bão hòa. Các tính chất thủy lực bao gồm độ dẫn thủy lực không bão hòa và đường cong đặc tính nước-đất (SWCC). Việc ước tính chính xác các thông số này giúp cải thiện độ chính xác của các mô hình thủy lực và ứng dụng trong các bài toán kỹ thuật địa chất.

2.1. Phương pháp thí nghiệm

Các thí nghiệm thoát nước một chiều được thực hiện để đo lường áp suất hút, độ bão hòa và tốc độ thoát nước theo thời gian. Các dữ liệu này được sử dụng làm đầu vào cho phương pháp phân tích ngược để ước tính các tính chất thủy lực. Kết quả cho thấy sự phù hợp cao giữa dữ liệu dự đoán và dữ liệu thực nghiệm, đặc biệt trong thí nghiệm thoát nước nhiều bước.

2.2. So sánh kết quả

So sánh giữa các kết quả dự đoán và đo lường cho thấy phương pháp phân tích ngược có độ chính xác cao trong việc ước tính các tính chất thủy lực. Đặc biệt, thí nghiệm thoát nước nhiều bước cho kết quả chính xác hơn so với thí nghiệm thoát nước một bước. Điều này khẳng định tính ứng dụng và hiệu quả của phương pháp trong thực tế.

III. Cát không bão hòa

Cát không bão hòa là một loại đất phổ biến trong các bài toán kỹ thuật địa chất. Việc hiểu và xác định các tính chất thủy lực của cát không bão hòa là rất quan trọng trong việc thiết kế và xây dựng các công trình. Phương pháp phân tích ngược đã được áp dụng để ước tính các thông số thủy lực của cát không bão hòa, cho thấy tiềm năng lớn trong việc ứng dụng vào thực tế.

3.1. Đặc tính của cát không bão hòa

Cát không bão hòa có các đặc tính thủy lực phức tạp, bao gồm độ dẫn thủy lực không bão hòa và đường cong đặc tính nước-đất. Việc xác định chính xác các thông số này giúp cải thiện độ chính xác của các mô hình thủy lực và ứng dụng trong các bài toán kỹ thuật địa chất.

3.2. Ứng dụng trong thực tế

Phương pháp phân tích ngược đã được áp dụng để ước tính các tính chất thủy lực của cát không bão hòa trong các thí nghiệm thoát nước một chiều. Kết quả cho thấy phương pháp này không chỉ đơn giản mà còn có độ chính xác cao, đặc biệt khi so sánh với các phương pháp truyền thống. Ứng dụng của phương pháp này trong thực tế giúp giảm thiểu thời gian và chi phí, đồng thời cải thiện độ tin cậy của các kết quả thí nghiệm.

Nội dung chính

Tổng quan luận án

Luận án tiến sĩ kỹ thuật "Inverse Analysis to Estimate Hydraulic Properties for Unsaturated Sand" do nghiên cứu sinh To Viet Nam thực hiện tại Khoa Kỹ thuật Xây dựng và Môi trường, Trường Sau đại học thuộc Đại học Ulsan (Hàn Quốc), dưới sự hướng dẫn của Giáo sư Min Tuk-ki (bảo vệ tháng 11 năm 2013). Luận án thuộc chuyên ngành Kỹ thuật Địa kỹ thuật - Môi trường, tập trung vào phương pháp phân tích nghịch (inverse analysis) để xác định các đặc tính thủy lực của đất cát không bão hòa.

Về tính cấp thiết và khoảng trống nghiên cứu, tác giả chỉ ra rằng khoảng một phần ba bề mặt Trái Đất phân bố ở vùng khí hậu khô hạn và bán khô hạn, nơi điều kiện đất không bão hòa chiếm ưu thế. Các đặc tính thủy lực của đất không bão hòa—bao gồm đường đặc trưng nước - đất (Soil-Water Characteristic Curve - SWCC) và hàm độ dẫn thủy lực không bão hòa—quy định các quá trình thấm, ổn định mái dốc, sức chịu tải và ô nhiễm nguồn nước. Tuy nhiên, các phương pháp thí nghiệm hiện trường và phòng thí nghiệm truyền thống đòi hỏi điều kiện cân bằng tĩnh hoặc dòng dừng (steady-state), dẫn đến quy trình kéo dài, chi phí cao, phức tạp và dễ phát sinh sai số. Do đó, việc ứng dụng bài toán nghịch dựa trên các thí nghiệm dòng không dừng (transient flow) là một hướng tiếp cận cần thiết để khắc phục các hạn chế này.

Mục tiêu nghiên cứu tổng quát của luận án là đánh giá khả năng áp dụng của kỹ thuật phân tích nghịch trong việc ước tính các đặc tính thủy lực không bão hòa cho đất hạt rời từ dữ liệu đo đạc của thí nghiệm thoát nước một chiều (1-D outflow experiment). Tác giả xác định 5 mục tiêu nghiên cứu cụ thể được đánh số trong văn bản như sau:

1. Thực hiện thí nghiệm khử bão hòa thẳng đứng một chiều (1-D vertical de-saturation experiment) cho cột cát đồng nhất và cột cát phân tầng trong điều kiện thoát nước một bước (one-step) và nhiều bước (multi-step).
2. Chế tạo và sử dụng các đầu đo điện trở suất đơn giản kết hợp đường chuẩn quan hệ giữa điện trở suất và độ bão hòa để xác định độ bão hòa trong mẫu đất.
3. Đo đạc độ hút dính (suction), độ bão hòa và lượng nước thoát lũy tích theo thời gian trong toàn bộ quá trình thí nghiệm.
4. Xây dựng mã nguồn phần tử hữu hạn và mã nguồn phân tích nghịch để ước tính các thông số thủy lực không bão hòa bằng cách sử dụng dữ liệu đo độ hút dính và tốc độ thoát nước làm dữ liệu đầu vào.
5. So sánh các đặc tính thủy lực không bão hòa giữa giá trị tính toán và giá trị đo đạc, cũng như so sánh giữa thí nghiệm thoát nước một bước và thí nghiệm thoát nước nhiều bước.

Đối tượng nghiên cứu của đề tài là đất hạt rời, cụ thể là cát tiêu chuẩn Jumunjin của Hàn Quốc (cát trung) và hệ đất cát phân tầng gồm cát mịn và cát trung. Phạm vi nghiên cứu giới hạn trong điều kiện thí nghiệm mô hình cột một chiều thẳng đứng trong phòng thí nghiệm, nghiên cứu nhánh sấy khô (drying path) theo chế độ thoát nước một bước và nhiều bước.

Tổng quan tài liệu và vị trí của luận án

Trong chương tổng quan, tác giả hệ thống hóa các khái niệm nền tảng về cơ học đất không bão hòa. Đất không bão hòa được mô tả theo mô hình ba pha cơ bản (pha rắn, lỏng, khí) theo Lambe và Whitman (1969), đồng thời xét đến pha thứ tư là màng co rút khí - nước (contractile skin) theo Paddy (1969) cùng Fredlund và Morgenstern (1977). Trạng thái ứng suất của đất không bão hòa được chi phối bởi hai biến độc lập là ứng suất pháp thực $(\sigma - u_a)$ và độ hút dính ma trận $(u_a - u_w)$ theo Jennings và Burland (1962) cùng Fredlund và Morgenstern (1977).

Về mô hình hóa toán học cho SWCC và hàm độ dẫn thủy lực, luận án điểm lại các mô hình tiêu biểu:

• Mô hình Brooks và Corey (1964): Dạng hàm lũy thừa dựa trên áp lực khí vào $\psi_{ae}$ và thông số đặc tính $\lambda$. Mô hình này phù hợp với khoảng hút dính hẹp nhưng thiếu điểm uốn nên biểu diễn kém ở dải hút dính rộng.
• Mô hình van Genuchten (1980): Hàm giải tích đóng dạng chữ S mượt mà gồm ba tham số $\alpha$, $n$, $m$ (với $m = 1 - 1/n$), mô tả hiệu quả toàn bộ dải hút dính và liên kết với độ dẫn thủy lực tương đối $k_r$.
• Mô hình Fredlund và Xing (1994): Xây dựng dựa trên phân bố kích thước lỗ rỗng, tích hợp hệ số hiệu chỉnh $C(\psi)$ đưa độ hút dính về giá trị $10^6\text{ kPa}$ tại độ ẩm bằng không, kết hợp lý thuyết Childs và Collis-George (1950) để suy biến hàm $k_r$.

Về phương pháp thực nghiệm, tác giả phân loại các kỹ thuật đo đạc phòng thí nghiệm: phương pháp cột dài (long column method), phương pháp cột nước treo dùng phễu Buchner (Sharma và Mohamed, 2003), phương pháp ly tâm (Lu và Likos, 2004), phương pháp cột nước không đổi (constant-head method) và phương pháp thoát nước (outflow method). Phương pháp thoát nước bao gồm kỹ thuật một bước của Gardner (1958), nhiều bước của van Dam và cộng sự (1990) và dòng liên tục của Durner và cộng sự (1999b).

Về bài toán nghịch trong thủy lực đất, kỹ thuật này được Gardner (1956) áp dụng ban đầu, sau đó Dane và Hruska (1983) ứng dụng cho dữ liệu hiện trường để tối ưu hóa thông số van Genuchten từ quá trình thoát nước không dừng. Kool và cộng sự (1985a, 1985b) và Kool và Parker (1988) áp dụng nghiệm nghịch cho thí nghiệm một bước chỉ dựa vào dữ liệu thoát nước lũy tích. Tuy nhiên, Levasseur và cộng sự (2008) chỉ ra rằng bài toán nghịch thường đối mặt với vấn đề không duy nhất nghiệm (non-uniqueness). Để giải quyết hiện tượng này, nhiều nghiên cứu khuyến nghị bổ sung dữ liệu độ ẩm hoặc áp lực nước lỗ rỗng theo thời gian vào hàm mục tiêu (Hudson và cộng sự, 1991; van Dam và cộng sự, 1992; Toorman và cộng sự, 1992; Bohne và cộng sự, 1993; Eching và Hopmans, 1993; Simunek và cộng sự, 1998; Durner và cộng sự, 1999b; Hopmans và cộng sự, 2002).

Khoảng trống nghiên cứu được luận án lựa chọn là: thiết lập hệ thống thực nghiệm cột cát 1-D tích hợp cảm biến điện trở suất tự chế để đo độ bão hòa theo thời gian, kết hợp đầu đo tensiometer đa tầng và lưu lượng thoát nước; giải bài toán nghịch ước tính đồng thời các thông số SWCC và hàm độ dẫn thủy lực của mô hình van Genuchten cho cát tiêu chuẩn Jumunjin; đồng thời mở rộng quy trình phân tích nghịch từ môi trường cát đồng nhất sang môi trường phân tầng thẳng đứng phức tạp.

Cơ sở lý thuyết và phương pháp nghiên cứu

Luận án dựa trên lý thuyết dòng chảy không dừng trong môi trường xốp không bão hòa, được mô tả bằng phương trình Richards một chiều:

$$\frac{\partial \theta}{\partial t} = \frac{\partial}{\partial z} \left[ k(h) \left( \frac{\partial h}{\partial z} + 1 \right) \right]$$

trong đó $\theta$ là độ ẩm thể tích, $h$ là áp lực cột nước hút dính, $k(h)$ là hàm độ dẫn thủy lực không bão hòa, $z$ là tọa độ không gian thẳng đứng và $t$ là thời gian. Khung phân tích thủy lực dựa trên mô hình van Genuchten (1980):

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

$$k = k_s S_e^{1/2} \left[ 1 - (1 - S_e^{1/m})^m \right]^2$$

với $S_e$ là độ bão hòa hiệu dụng, $\theta_s$ và $\theta_r$ lần lượt là độ ẩm bão hòa và độ ẩm dư, $\psi$ là độ hút dính ma trận, $k_s$ là hệ số thấm bão hòa, $\alpha$, $n$, $m$ là các hệ số thực nghiệm ($m = 1 - 1/n$).

Phương pháp nghiên cứu kết hợp giữa đo đạc thực nghiệm định lượng và mô phỏng số giải bài toán nghịch:

• Thiết lập thí nghiệm: Cột thí nghiệm hình trụ bằng vật liệu acrylic trong suốt, đường kính trong 10 cm, chiều cao 100 cm. Đáy cột lắp màng xốp hoặc đĩa gốm bão hòa nước có áp lực khí vào phù hợp. Quá trình thoát nước được kiểm soát theo chế độ hạ thế nước một bước (one-step) hoặc nhiều bước (multi-step).
• Hệ thống quan trắc: Áp lực nước lỗ rỗng (suction) được đo bằng các đầu đo tensiometer bố trí dọc theo chiều cao thân cột cách nhau 10 cm (L2, L4, L5, L8...). Độ bão hòa nước của đất được xác định liên tục thông qua hệ thống đầu đo điện trở suất tự chế, chuyển đổi qua đường chuẩn quan hệ điện trở suất - độ bão hòa. Lượng nước thoát ra đáy cột được đo bằng cân điện tử hoặc ống đo chuẩn để xác định lưu lượng lũy tích $Q(t)$ và tốc độ thoát nước $q(t)$.
• Quy trình phân tích nghịch: Tác giả phát triển một mã nguồn phần tử hữu hạn (FEM code) để mô phỏng dòng chảy thuận (forward water flow modeling) giải phương trình Richards, kết hợp thuật toán tối ưu hóa phi tuyến (nonlinear optimization / inversion code) dựa trên phương pháp bình phương tối thiểu để cực tiểu hóa hàm mục tiêu sai số dư:

$$\Phi(b) = \sum_{i=1}^{N} w_i \left[ y_i^* - y_i(b) \right]^2$$

trong đó $y_i^*$ là giá trị đo đạc thực nghiệm (độ hút dính, tốc độ thoát nước), $y_i(b)$ là giá trị dự tính từ mô hình dòng chảy thuận với vector tham số $b = (\alpha, n, k_s, \theta_r, \theta_s)$, và $w_i$ là thừa số trọng số. Mức độ phù hợp giữa mô hình và thực nghiệm được đánh giá thông qua chỉ số sai số căn bậc hai trung bình (Root Mean Square Error - RMSE).

Nguồn dữ liệu sử dụng là số liệu quan trắc thực nghiệm trực tiếp từ các kịch bản thí nghiệm trên cát đồng nhất Jumunjin và các cột cát phân tầng hai lớp.

Nội dung chính theo từng chương

Chapter 1: Introduction

Chương mở đầu trình bày bối cảnh địa kỹ thuật về sự phân bố của đất không bão hòa trên thế giới, mối liên hệ giữa điều kiện khí hậu và mực nước ngầm tạo nên tầng vadose. Tác giả nhấn mạnh những khó khăn của phương pháp đo trực tiếp và tính ưu việt của phân tích nghịch. Chương này công bố mục tiêu nghiên cứu tổng quát, 5 mục tiêu cụ thể, 4 đóng góp nguyên bản của luận án và tóm tắt bố cục toàn văn luận án từ Chương 1 đến Chương 6.

Chapter 2: Unsaturated Soil Hydraulic Properties

Chương 2 trình bày chi tiết cơ sở lý thuyết về môi trường đất không bão hòa, hệ thống ứng suất hai biến độc lập $(\sigma - u_a)$ và $(u_a - u_w)$, cùng bản chất nhiệt động học của độ hút dính đất. Tác giả định nghĩa các vùng đặc trưng trên đường SWCC sấy khô: vùng bão hòa mao dẫn (kết thúc tại giá trị khí vào $\psi_a$), vùng khử bão hòa (desaturation zone) và vùng bão hòa dư (residual saturation zone). Chương này phân tích cấu trúc toán học của các mô hình SWCC kinh điển (Brooks và Corey 1964; van Genuchten 1980; Fredlund và Xing 1994) và khảo sát ưu, nhược điểm của các phương pháp đo thí nghiệm: cột dài, cột nước treo, ly tâm, cột nước không đổi và phương pháp thoát nước (một bước, nhiều bước, dòng liên tục).

Chapter 3: Inverse Analysis

Chương 3 xây dựng khung lý thuyết toán học cho bài toán nghịch trong địa kỹ thuật. Tác giả phân biệt bài toán thuận (xác định đáp ứng từ nguyên nhân và tham số đã biết) và bài toán nghịch (ước tính tham số hệ thống từ đáp ứng quan trắc). Cấu trúc của phân tích nghịch xác định thông số thủy lực đất không bão hòa gồm ba thành phần tương hỗ:

1. Thí nghiệm dòng không dừng được kiểm soát điều kiện biên và đo đạc các biến dòng chảy.
2. Mô hình dòng chảy số (phương trình Richards) mô phỏng chế độ dòng chảy không dừng từ ước lượng tham số ban đầu.
3. Thuật toán tối ưu hóa phi tuyến lặp để cực tiểu hóa hàm mục tiêu sai số giữa đại lượng mô phỏng và thực nghiệm.

Chương này phân tích nguyên nhân gây ra nghiệm không duy nhất (non-uniqueness) và chứng minh sự cần thiết của việc kết hợp đồng thời áp lực hút dính, tốc độ thoát nước và độ ẩm vào hàm mục tiêu.

Chapter 4: Inverse Analysis to Estimate Unsaturated Hydraulic Properties for Homogeneous Sand

Chương 4 thực hiện quy trình nghiên cứu trên cột cát đồng nhất Jumunjin. Tác giả mô tả thiết bị cột cát, quy trình chế tạo đầu đo điện trở suất chi phí thấp và thiết lập đường chuẩn điện trở suất - độ bão hòa. Hai kịch bản thoát nước được triển khai: thoát nước một bước (one-step) và thoát nước nhiều bước (multi-step).

Đại lượng / Thông số	Thí nghiệm One-step Outflow	Thí nghiệm Multi-step Outflow
Vật liệu nghiên cứu	Cát tiêu chuẩn Jumunjin (cát trung)	Cát tiêu chuẩn Jumunjin (cát trung)
Chế độ áp lực biên	Hạ áp lực một bước duy nhất	Hạ áp lực theo nhiều bước nhỏ gián đoạn
Dữ liệu đầu vào nghịch	Áp lực nước lỗ rỗng (tensiometer) và tốc độ thoát nước $q(t)$	Áp lực nước lỗ rỗng (tensiometer) và tốc độ thoát nước $q(t)$
Vị trí kiểm chứng áp lực	L2, L5, L8 (khoảng cách 10 cm dọc cột)	L2, L5, L8 (khoảng cách 10 cm dọc cột)
Kết quả đường SWCC	Phù hợp với số liệu đo, xuất hiện tương quan tham số	Mức độ tương quan tham số giảm, độ chính xác cao hơn
Đánh giá sai số (RMSE)	Bảng 4.2: Sai số thời gian giữa dự tính và thực tế	Bảng 4.3: Sai số RMSE thấp hơn thí nghiệm one-step

Kết quả so sánh giữa số liệu dự tính từ mô hình nghịch và số liệu quan trắc thực tế cho thấy sự tương đồng cao về biến thiên áp lực nước lỗ rỗng, độ bão hòa theo chiều cao và đường cong thoát nước. Thí nghiệm nhiều bước cho thấy tính ưu việt hơn thí nghiệm một bước do giảm thiểu sự tương quan phụ thuộc giữa các tham số tối ưu hóa.

Chapter 5: Inverse Analysis to Estimate Unsaturated Hydraulic Properties for Layered Sand

Chương 5 mở rộng bài toán nghịch cho cột cát phân tầng thẳng đứng gồm hai lớp cát có cấp phối hạt khác nhau: cát mịn và cát trung. Đặc tính cơ lý cơ bản của hai loại cát được xác định độc lập (Bảng 5.1). Thí nghiệm thoát nước một bước được tiến hành trên hai cấu hình phân tầng:

• Trường hợp A (Case A): Cấu hình lớp cát trung nằm phía trên lớp cát mịn (medium sand overlies fine sand).
• Trường hợp B (Case B): Cấu hình hoán đổi vị trí, lớp cát mịn nằm phía trên lớp cát trung (fine sand overlies medium sand).

Thông số / Đặc điểm	Cấu hình Case A	Cấu hình Case B
Thứ tự lớp đất	Lớp trên: Cát trung; Lớp dưới: Cát mịn	Lớp trên: Cát mịn; Lớp dưới: Cát trung
Phản ứng thoát nước	Thoát nước nhanh ở tầng cát trung phía trên	Tầng cát mịn giữ nước, làm chậm quá trình thoát nước
Điểm quan trắc kiểm chứng	Vị trí L2, L4, L5, L8 dọc cột	Vị trí L2, L5, L8 dọc cột
Đánh giá sai số (RMSE)	Bảng 5.2: RMSE áp lực và lưu lượng miền thời gian	Bảng 5.3: RMSE áp lực và lưu lượng miền thời gian
Đánh giá mô hình nghịch	Tái tạo chính xác cấu trúc SWCC từng lớp	Dự tính phù hợp với số liệu đo thực tế

Kết quả so sánh giữa đường cong SWCC dự tính từ bài toán nghịch và số liệu thí nghiệm độc lập cho thấy mô hình nghịch phản ánh chính xác quy luật thấm và phân bố độ ẩm của hệ đất phân tầng dị tính, khẳng định tính khả thi của phương pháp khi áp dụng cho các cấu trúc địa tầng phức tạp.

Chapter 6: Conclusions

Chương kết luận tổng hợp toàn bộ các phát hiện nghiên cứu từ Chương 4 và Chương 5. Tác giả khẳng định tính chính xác và độ tin cậy của thuật toán phân tích nghịch kết hợp mô hình phần tử hữu hạn dòng chảy một chiều. Đồng thời, chương này vạch ra các định hướng nghiên cứu tiếp theo nhằm nâng cao phạm vi ứng dụng của kỹ thuật phân tích nghịch trong địa kỹ thuật.

Kết quả và những đóng góp mới

Theo tuyên bố của tác giả tại Mục 1.2, luận án đem lại 4 đóng góp nguyên bản (original contributions) cho lĩnh vực nghiên cứu:

1. Cải tiến thiết bị và kỹ thuật thí nghiệm cột đất: Chế tạo thành công các đầu đo điện trở suất mới có cấu tạo đơn giản, chi phí thấp, cho phép xác định độ bão hòa nước trong mẫu đất nhanh chóng và liên tục trong suốt quá trình thí nghiệm. Độ bão hòa suy ra từ giá trị điện trở suất đạt độ chính xác cao và phù hợp tốt với các phép đo độc lập.
2. Ứng dụng phân tích nghịch cho cát tiêu chuẩn Hàn Quốc: Luận án là công trình nghiên cứu áp dụng kỹ thuật phân tích nghịch để xác định các đặc tính thủy lực không bão hòa cho cát tiêu chuẩn Jumunjin, tạo lập cơ sở tham chiếu và khung phương pháp luận để áp dụng cho các loại đất khác.
3. Xây dựng công cụ tính toán: Phát triển thành công mã nguồn phần tử hữu hạn giải phương trình dòng chảy thuận và mã nguồn tối ưu hóa phi tuyến giải bài toán nghịch, tạo thành công cụ phục vụ tính toán các thông số thủy lực của đất.
4. Đóng góp về mặt phương pháp học: Xác nhận các kết luận của các nghiên cứu trước đây rằng phân tích nghịch giúp cải thiện tốc độ và độ chính xác trong việc xác định thông số thủy lực đất; đồng thời chứng minh bằng thực nghiệm rằng thí nghiệm thoát nước nhiều bước (multi-step outflow) đạt độ chính xác cao hơn thí nghiệm thoát nước một bước (one-step outflow) do giảm tương quan tham số.

Về mặt kiến nghị và giải pháp, luận án đề xuất áp dụng quy trình phân tích nghịch đa biến (kết hợp đồng thời dữ liệu áp lực hút dính theo chiều cao và lưu lượng thoát nước đáy) để giải quyết triệt để bài toán không duy nhất nghiệm. Đối với các bài toán thực hành xác định thông số đất không bão hòa, tác giả khuyến nghị ưu tiên áp dụng chế độ thoát nước nhiều bước nhằm đảm bảo tính ổn định và hội tụ của thuật toán tối ưu.

Hạn chế và hướng nghiên cứu tiếp

Về mặt hạn chế, văn bản luận án phản ánh các điểm giới hạn sau:

• Phạm vi đối tượng vật liệu: Nghiên cứu mới dừng lại ở đất hạt rời thô (cát trung Jumunjin và cát mịn phân tầng); chưa tiến hành trên các loại đất hạt mịn dính (sét, bùn) vốn có độ hút dính lớn vượt quá áp lực bọt khí của các thiết bị đĩa xốp và tensiometer thông thường.
• Quy mô hình học và điều kiện biên: Thực nghiệm và mô hình hóa chỉ giới hạn trong bài toán dòng chảy một chiều (1-D) thẳng đứng trong phòng thí nghiệm; chưa xét đến dòng chảy hai chiều hoặc ba chiều.
• Nhánh quá trình: Thí nghiệm chỉ tập trung vào nhánh thoát nước/sấy khô (drying branch), chưa mô phỏng quá trình làm ướt (wetting branch) để khảo sát hiện tượng trễ thủy lực (hysteresis).

Về hướng nghiên cứu tiếp theo (stated in Chapter 6), luận án định hướng:

• Mở rộng phân tích nghịch cho các điều kiện biên và cấu trúc địa tầng có độ phức tạp cao hơn.
• Nghiên cứu ứng dụng trên các loại đất có tính chất cơ lý khác nhau, bao gồm đất dính và đất lẫn nhiều thành phần hạt.
• Phát triển mô hình bài toán nghịch cho không gian dòng chảy nhiều chiều (2D/3D) phục vụ các bài toán công trình thực tế ngoài hiện trường.

Giá trị tham khảo

Luận án là tài liệu tham khảo trong các lĩnh vực sau:

• Đối với nghiên cứu sinh và học viên cao học ngành Địa kỹ thuật, Thủy văn đất và Kỹ thuật Môi trường: Luận án cung cấp hệ thống lý thuyết về cơ học đất không bão hòa, phương trình Richards và quy trình giải bài toán nghịch tối ưu hóa phi tuyến.
• Đối với các phòng thí nghiệm cơ học đất: Phần thiết kế và hiệu chuẩn đầu đo điện trở suất tự chế (Chương 4) là giải pháp thực nghiệm có giá trị ứng dụng cao để quan trắc độ bão hòa liên tục với chi phí thấp mà không phá hủy mẫu.
• Đối với kỹ sư thiết kế và chuyên gia mô phỏng địa kỹ thuật: Khung phân tích nghịch kết hợp mã nguồn phần tử hữu hạn và thuật toán tối ưu (Chương 3 và Chương 4) cung cấp phương pháp luận để xác định nhanh các thông số van Genuchten ($\alpha$, $n$, $k_s$, $\theta_r$, $\theta_s$), phục vụ phân tích thấm qua đập, ổn định mái dốc ven sông và lan truyền chất ô nhiễm trong tầng vadose.

Câu hỏi thường gặp

1. Mô hình toán học nào được lựa chọn làm cơ sở để ước tính các đặc tính thủy lực trong luận án?

Luận án sử dụng mô hình van Genuchten (1980) dạng giải tích đóng ba tham số: $$S_e = \left[ 1 + (\alpha \psi)^n \right]^{-m}$$ với điều kiện $m = 1 - 1/n$. Hàm độ dẫn thủy lực tương đối được liên kết thông qua công thức: $$k = k_s S_e^{1/2} \left[ 1 - (1 - S_e^{1/m})^m \right]^2$$ Mô hình này được lựa chọn nhờ khả năng mô tả dạng đường cong chữ S mượt mà, bao quát toàn bộ dải hút dính và liên kết giải tích trực tiếp với hàm độ dẫn thủy lực.

2. Vì sao thí nghiệm thoát nước nhiều bước (multi-step) được đánh giá phù hợp hơn thí nghiệm một bước (one-step)?

Theo phân tích tại Chương 2, Chương 4 và Mục 1.2, trong thí nghiệm một bước, áp lực thay đổi đột ngột tạo ra sự sụt giảm độ ẩm tức thời, dẫn đến hiện tượng tương quan mạnh giữa các tham số tối ưu hóa (parameter correlation). Thí nghiệm nhiều bước hạ áp lực theo từng khoảng nhỏ giúp quá trình thoát nước diễn ra tuần tự, các biến dòng chảy phân bố đều hơn trong miền thời gian, làm giảm sự phụ thuộc lẫn nhau giữa các tham số, từ đó tăng độ chính xác và tính ổn định của nghiệm nghịch.

3. Phương pháp đo độ bão hòa của đất trong suốt quá trình thí nghiệm được thực hiện như thế nào?

Tác giả đã tự chế tạo các đầu đo điện trở suất (electrical resistivity probes) kích thước nhỏ, giá thành thấp đặt tại các cao độ khác nhau trong cột cát. Độ bão hòa được xác định thông qua đường chuẩn quan hệ thực nghiệm giữa điện trở suất và độ bão hòa (calibration relationship between electrical resistivity and saturation) được thiết lập từ các phép thử độc lập trước khi tiến hành thí nghiệm cột.

4. Vấn đề không duy nhất nghiệm (non-uniqueness) trong bài toán nghịch được tác giả khắc phục bằng giải pháp nào?

Tác giả giải quyết hiện tượng không duy nhất nghiệm bằng cách tích hợp đồng thời nhiều loại dữ liệu đo thực nghiệm vào hàm mục tiêu tối ưu hóa, bao gồm: áp lực nước lỗ rỗng (độ hút dính) theo thời gian đo tại nhiều vị trí bằng tensiometer, tốc độ thoát nước $q(t)$ hoặc lưu lượng lũy tích $Q(t)$ ở đáy cột, và dữ liệu độ bão hòa thu được từ đầu đo điện trở suất. Việc bổ sung đa nguồn dữ liệu giúp thu hẹp không gian tìm kiếm tham số và đảm bảo nghiệm tối ưu hội tụ về giá trị thực.

5. Cấu trúc địa tầng của hai cột cát phân tầng trong Chương 5 được bố trí như thế nào?

Thí nghiệm cột phân tầng được thực hiện với hai trường hợp hoán đổi vị trí của hai loại cát (cát mịn và cát trung):

• Case A: Lớp cát trung (medium sand) bố trí ở nửa trên, lớp cát mịn (fine sand) bố trí ở nửa dưới.
• Case B: Lớp cát mịn (fine sand) bố trí ở nửa trên, lớp cát trung (medium sand) bố trí ở nửa dưới. Cả hai trường hợp đều được khảo sát theo quy trình thoát nước một bước để đánh giá khả năng nhận diện thông số từng lớp của thuật toán phân tích nghịch.

Kết luận

Luận án tiến sĩ của To Viet Nam đã thiết lập thành công khung phương pháp luận phân tích nghịch kết hợp mô hình số phần tử hữu hạn để xác định đồng thời đường đặc trưng nước - đất (SWCC) và hàm độ dẫn thủy lực không bão hòa cho đất cát. Thông qua các thí nghiệm cột 1-D trên cát tiêu chuẩn Jumunjin và hệ cát phân tầng, nghiên cứu đã chứng minh tính khả thi, độ tin cậy và sự phù hợp của mô hình nghịch khi sử dụng dữ liệu kết hợp từ tensiometer và lưu lượng thoát nước. Công trình đã khẳng định ưu thế giảm tương quan tham số của thí nghiệm nhiều bước so với thí nghiệm một bước, đồng thời cung cấp giải pháp đo đạc độ bão hòa bằng đầu đo điện trở suất chi phí thấp phục vụ công tác nghiên cứu thực nghiệm địa kỹ thuật.

Trích đoạn nội dung tài liệu

Doctor of Philosophy Inverse Analysis to Estimate Hydraulic Properties for Unsaturated Sand The Graduate School of the University of Ulsan Department of Civil and Environmental Engineering To Viet Nam 1 Doctor of Philosophy Inverse Analysis to Estimate Hydraulic Properties for Unsaturated Sand The Graduate School of the University of Ulsan Department of Civil and Environmental Engineering To Viet Nam 2 Inverse Analysis to Estimate Hydraulic Properties for Unsaturated Sand Supervisor: Professor Min Tuk-ki A Dissertation Submited to The Graduate School of the University of Ulsan in Partial Fulfillment of the Requirements for the Degree of Doctor of Philosophy By To Viet Nam Department of Civil and Environmental Engineering Ulsan, Korea November 2013 3 4 ACKNOWLEDGMENTS This is my great pleasure to have this opportunity to thank people who have helped me during my dissertation and my study at Laboratory of Geotechnics, Department of Civil and Environmental Engineering, University of Ulsan, Korea. First of all, I would like to sincerely thank my supervisor Professor Min Tuk-Ki for his brilliant supervising, guidance, patience, consistent supports and encouragements throughout my Ph. I would also like to thank and gratitude to Professor Shin Ho-sung for accepting me into his group, giving me freedom in my work, supporting my attendance at various conferences, engaging me in new ideas, and demanding a high quality of work in all my attemps. Additionally, I would like to acknowledge the committee members Professor Min Tuk-Ki, Professor Lee Jong Seok, Professor Yoo Ji Hyeong, Professor Kang Suk-Bong and Professor Shin Ho-Sung for their comments, advises and suggestions in my study.

I also express my appreciation to all members of the Geotechnical Lab. Every result described in this thesis was accomplished with the help and support of fellow labmates and collaborators. Thank are also due to professors and staff members of Department of civil and Environmental Engineering for their kindly helps. I would like to thank my parents, my elder sister and my elder brother, who provided endless help, patience and good advice to have me through my program.

Special thanks and many loves to my wife, Nguyen Thi Thuy Linh and my daughter, To Nguyen Quynh Huong for endless love, great cares, encouragements and sufferings of a separation during the time I am away from home. Finally, I would like to thank for the support from all my friends from Union of Vietnamese Students in University of Ulsan. With warmest regards Ulsan 2013 To Viet Nam 5 AN INVERSE ANALYSIS TO ESTIMATE HYDRAULIC PROPERTIES FOR UNSATURATED SAND By To Viet Nam Department of civil and Environmental Engineering University of Ulsan ABSTRACT Since the term “unsaturated soils” was defined, and relation of unsaturated soils to many engineering problems was proven, the determination of unsaturated hydraulic properties, which is one of most important parameters in unsaturated soils, is a hot topic in geotechnical engineering. Realizing importance of unsaturated soil hydraulic properties, for a long time, many methods for determining the unsaturated hydraulic properties of soils have been proposed by many researchers.

However, these methods usually encounter the conditions of tediousness, time and money consumption, complexity, inaccuracy, narrow applicability, etc. Therefore, up to now, understanding and discovering a simple, appropriate and accurate method for determination of unsaturated hydraulic properties of soils is still received many interests from the scientists. In recent decades, the inverse analysis method becomes a useful technique in determining the unsaturated hydraulic properties of soils. This method, in many cases, has been shown to be an exciting new tool that allows for soil hydraulic characterization using a wide range of transient laboratory and field experiments.

The objective of this dissertation is to evaluate the applicability of inverse analysis technique in estimating the unsaturated hydraulic properties for grain soils using the input data measured from the one-dimension (1-D) outflow experiments. The applicability of inverse analysis method in this dissertation is evaluated through two following parts. 6 In the first part, the 1-D vertical draining test for homogeneous medium sand was carried out in one-step and multi-step outflow conditions. An inverse analysis method was conducted to determine the hydraulic properties of unsaturated soil.