Hình 2.2: Dòng thấm không ổn định khi mực nước rút xuống đột ngột

18



c. Dòng thấm không ổn định trong trường hợp môi trường thấm thay đổi theo

thời gian. Môi trường thấm có thể thay đổi theo thời gian do sự nén chặt của đất sét

dưới ảnh hưởng của lực hoặc do sự làm chặt của đất cát rời sau khi chúng được xáo

trộn, ví dụ như do một vài tải trọng động.

Dòng thấm ổn định hoặc không ổn định được tính toán xác định thông qua các

phương trình vi phân tổng quát về thấm:

+ Trường hợp dòng thấm ổn định:

∂  ∂H

 kx

∂x  ∂x



∂θ

 ∂  ∂H 

+Q =

 +  ky

∂t

 ∂y  ∂y 



(2.1)



+ Trường hợp dòng thấm không ổn định:

∂  ∂H

 kx

∂x  ∂x



 ∂  ∂H

 +  ky

 ∂y  ∂y





∂H

mw γ w

+Q =

∂t





(2.2)



Trong đó:

H



: Tổng cột nước;



kx



: Hệ số thấm theo phương x;



ky



: Hệ số thấm theo phương y;



Q



: Lưu lượng biên;



t



: thời gian;



mw



: Hệ số góc đường cong lượng chứa nước với áp lực lỗ rỗng;



γw



: Trọng lượng đơn vị của nước.



2.8. Phương pháp tính thấm

Hiện nay có nhiều phương pháp tính thấm thông dụng có phương pháp phân

tích lý luận, đồ giải và thí nghiệm.

Phương pháp phân tích lý luận bao gồm phương pháp cơ học chất lỏng,

phương pháp phần tử hữu hạn và phương pháp thủy lực [15], [16].



2.8.2. Cơ sở lý thuyết

Chuyển động dòng thấm chính là chuyển động của các dòng nước trong đất có

chứa các lỗ rỗng, phụ thuộc vào đất nền và chênh lệch cột nước thượng hạ lưu.

Quy luật cơ bản về sự chuyển động của dòng thấm biểu thị định luật Darcy:



19



∂h



=

−

Vx

k

.



∂x



V = kJ hay 

∂

Vy = − k . h



∂y





(2.3)



trong đó

V – lưu tốc thấm (m/s) ;

K – hệ số thấm của môi trường (m/s) ;

J – gradient thấm (độ dốc thủy lực).

Lưu tốc trung bình dòng thấm trong lỗ rỗng của đất (n là độ rỗng của môi

trường):

v' =



v

n



(2.4)



Với môi trường thấm ổn định áp lực thấm không phụ thuộc vào thời gian:



∂ 2ϕ ∂ 2ϕ ∂ 2ϕ

0

+

+

=

∂x 2 ∂y 2 ∂z 2



(2.5)



Với φ = -k.h hàm thế, trong đó φ và h là những hàm điều hòa.

Thông thường trong tính toán thấm cho bài toán 2 chiều:

 ∂Vx ∂Vy

 ∂x + ∂y = 0



Phương trình liên tục của bài toán phẳng :  2

2

∂ h + ∂ h = 0

2



∂y 2

 ∂x



∂ 2ϕ ∂ 2ϕ

+

=

0

∂x 2 ∂y 2



(2.6)



2.8.3. Tính thấm qua đập đất theo phương pháp thủy lực

Phương pháp thủy lực cho kết quả tính toán là những đặc trưng trung bình của

dòng thấm.

Ưu điểm: Có thể giải được những trường hợp phức tạp trong thực tế. Mặt

khác, độ chính xác đáp ứng đuợc yêu cầu kỹ thuật, nên được áp dụng rộng rãi.

Nhược điểm: Do phải dựa vào một số tiên đề, giả thiết nhất định. Nên phương

pháp này kém chính xác



20



2.8.4. Tính thấm qua đập đất theo phương pháp cơ học chất lỏng

Phương pháp này dùng toán học làm công cụ để xác định những đặc trưng của

dòng thấm như lưu lượng, lưu tốc, gradient, áp lực, đường bão hòa... tại bất kỳ một

vị trí nào trong môi trường thấm. Do khi tính toán không đưa nhiều vào những giả

thiết cho nên phương pháp cơ học chất lỏng cho kết quả chính xác.

Ưu điểm: Phương pháp cơ học chất lỏng chủ yếu có tầm quan trọng về mặt lý

thuyết, trên cơ sở đó người ta có thể đưa ra các giải quyết gần đúng. Ứng dụng

những lời giải của cơ học chất lỏng ta có thể lập được những biểu đồ tính toán để

dùng trong thực tế.

Nhược điểm: Phương pháp này chỉ sử dụng được trong trường hợp bài toán có

sơ đồ đơn giản. Khi gặp những sơ đồ phức tạp (điều kiện ban đầu và điều kiện biên

phức tạp) thì cách giải này gặp nhiều khó khăn về mặt toán học và trong nhiều

trường hợp gần như bế tắc. Do vậy trong thực tế thiết kế tính toán thấm, phương

pháp này ứng dụng rất hạn chế.



2.8.5. Tính thấm qua đập đất theo phương pháp thực nghiệm

Nghiên cứu thực nghiệm là dùng mô hình để xác định những đặc trưng của

dòng thấm. Phương pháp này bao gồm những loại chính sau:

- Thí nghiệm bằng máng kính: Việc thí nghiệm được tiến hành trên mô hình

trong máng kính để xác định lưu lượng, áp lực, đường bão hòa....

- Thí nghiệm khe hẹp: Dựa vào sự tương tự giữa chuyển động của chất lỏng có

độ nhớt lớn (các loại dầu nhờn) trong khe hẹp vơi sự chuyển động của nước trong

môi trường rỗng (chảy tầng) mà xác định những đặc trưng dòng thấm trên cơ sở đo

đạc những đặc trưng của dòng chảy của dầu trong khe hẹp.

- Phương pháp tương tự điện – thủy động (EGDA): Dùng tính tương tự giữa

dòng điện và dòng thấm để giải quyết những bài toán thấm trên mô hình dòng điện.

Ưu điểm: Được ứng dụng rộng rãi vì khá tiện lợi, có thể giải quyết được

những sơ đồ thấm phức tạp.

Nhược điểm: Giá thành thí nghiệm cao, kết quả thí nghiệm phụ thuộc vào điều

kiện thí nghiệm cũng như chủ quan của người làm thí nghiệm.



21



2.8.6. Tính thấm qua đập đất bằng phương pháp phần tử hữu hạn

Phương pháp phần tử hữu hạn (PTHH) sử dụng sự gần đúng về mặt vật lý của

môi trường. Việc giải bài toán chính xác trong môi trường liên tục được xấp xỉ bằng

tập hợp hữu hạn các phần tử liên kết với nhau ở các điểm nút. Theo phương pháp

này, miền thấm được chia thành những miền con gọi là các phần tử. Các phần tử

thường được sử dụng trong bài toán phẳng là phần tử tam giác và phần tử tứ giác.

Trong miền thấm, các phần tử có kích thước và hình dạng khác nhau để phù hợp với

các điều kiện biên hình học và tính chất của dòng thấm tại các khu vực khác nhau

(khu vực có cường độ thấm mạnh thì chia thành các phần tử có kích thước nhỏ và

ngược lại). Giá trị cột nước thấm tại một vị trí bất kỳ ở trong phần tử được xác định

dựa vào giá trị cột nước tại các điểm nút của phần tử.



2.8.6.1.



Trình tự giải bài toán bằng phương pháp PTHH



+ Chia miền tính toán thành nhiều miền nhỏ gọi là các phần tử, các phần tử

được nối với nhau bằng một số hữu hạn các điểm nút. Các nút này có thể là đỉnh

của các phần tử, cũng có thể là một số điểm đuợc quy ước trên cạnh của phần tử.

Trình tự tính toán của phương pháp PTHH được nghiên cứu đầy đủ trong bài

toán thấm ổn định phẳng. Cần phải xác định sự phân bố cột nước áp lực dưới đập,

gradient thấm và tổng lưu lượng dòng thấm. Ở đây cột nuớc áp lực được xem là

hàm chưa biết H. Mặt trơn của hàm H được xấp xỉ bằng một bộ các mảnh mặt

phẳng tam giác dạng i’, j’, k’ như hình II-15 xác định trên miền con tam giác i, j, k

thuộc miền nghiên cứu trong mặt phẳng xy.

Vị trí mặt phẳng trong không gian được xác định đơn trị bằng 3 điểm không

nằm trên một đường thẳng. Hiển nhiên là để xấp xỉ hàm trơn H bằng các mảnh mặt

phẳng thì các PTHH i, j, k phải là phần tử tam giác.



22



Hình 2.3:



Minh họa hàm xấp xỉ H của phần tử



Độ sai lệch của mặt phân mảnh xấp xỉ so với mặt trơn thực tế sẽ càng lớn khi

độ cong của mặt trơn càng lớn và kích thước phần tử càng lớn, từ đó rút ra quy tắc

cơ bản xây dựng lưới các phần tử là làm dày đặc lưới tại những nơi có gradient hàm

cần tìm cao, chẳng hạn như các cột áp.

Công cụ toán học của phương pháp PTHH là đảm bảo đưa ra bài toán tích

phân phương trình vi phân song điều hòa về phép giải hệ thống các phương trình

đại số tuyến tính, trong đó giá trị cột áp tại các nút phần tử hiện diện như là các ẩn

số.

Do giới hạn lời giải của bài toán thấm dừng trên cơ sở định luật Darcy

V = −K.



Trong đó:



∆H

∆L



(2.7)



- K: Ma trận hệ số thấm.

- V: Véc tơ tốc độ thấm.

- ∆H: Chênh lệch cột nước giữa 2 mặt cắt tính toán cách nhau một



khoảng DL. Hay có thể viết dưới dạng vi phân

V = −K.



dH

ds



(2.8)



Cơ sở của phương pháp PTHH để giải bài toán thấm dựa trên nguyên lý biến

phân mà phiến hàm được chọn là:



23



 ∂H 

1    ∂H 



L = ∫∫∫   K x .

 + K y 

2    ∂x 

∂y 



v



2



2





  dV

 



(2.9)



Trong đó Kx, Ky lần lượt là các hệ số thấm theo phương x,y.

+ Phạm vi của mỗi phần tử giả thiết một dạng phân bố xác định nào đó của

hàm cần tìm. Đối với bài toán thấm thì hàm xấp xỉ có thể là hàm cột nước.

Thường giả thiết hàm xấp xỉ là những đa thức nguyên mà hệ số của nó được

gọi là các thông số. Trong phương pháp PTHH, các thông số này được biểu diễn

qua các trị số của hàm và có thể là trị số của đạo hàm của nó tại các điểm nút của

phần tử. Dạng đa thức nguyên của hàm xấp xỉ phải chọn đảm bảo để bài toán hội tụ.

Nghĩa là khi tăng số phần tử lên khá lớn thì kết quả tính toán sẽ tiệm cận đến kết

quả chính xác.

Ngoài ra hàm xấp xỉ của bài toán thấm đuợc chọn phải được đảm bảo tuân

theo định luật Darcy. Song để thỏa mãn chặt chẽ tất cả các yêu cầu thì sẽ gặp nhiều

khó khăn trong việc lựa chọn mô hình và lập thuật toán giải. Do đó trong thực tế

người ta phải giảm bớt một số yêu cầu nào đó nhưng vẫn đảm bảo được nghiệm đạt

độ chính xác yêu cầu.

+ Dựa vào phương trình cơ bản của bài toán nghiên cứu để tìm các đại lượng

khác. Đối với bài toán sử dụng định luật Darcy để tìm trường lưu tốc thấm, trường

gradien thấm, lưu lượng thấm, vv...



2.8.6.2.



Giải bài toán thấm bằng phương pháp PTHH



Do đặc thù của đập đất thường có chiều dài lớn so với chiều cao. Do vậy việc

nghiên cứu bài toán thấm trong phạm vi luận văn chỉ dừng lại ở mô hình bài toán

phẳng, thấm ổn định. Lúc đầu có thể cắt một đoạn đập dài bằng đơn vị để nghiên

cứu.

- Biên S 1 : H = H 1 .

- Biên S 2 : H = H 2 .

- Biên S 3 :



∂H

= 0.

∂n