Động lực học cát biển - Hướng dẫn các ứng dụng thực hành - Chương 8 pdf

113 Chơng 8. Trầm tích lơ lửng 8.1 Tổng quan Đối với vận tốc dòng chảy hoặc điều kiện sóng đáng kể trên ngỡng chuyển động, cát bị kéo lên khỏi đáy và đi vào trạng thái lơ lửng, tại đó nó đợc mang đi với cùng vận tốc dòng chảy. Khi điều này xảy ra, phần trầm tích đợc mang đi trong trạng thái lơ lửng nói chung lớn hơn nhiều phần trầm tích đáy đợc mang đi đồng thời, và do đó dòng lơ lửng là một thành phần quan trọng cho suất vận chuyển trầm tích tổng cộng. Một yếu tố quan trọng trong việc thiết kế các công trình lấy nớc làm lạnh cho các nhà máy điện là ngăn ngừa sự xâm nhập của trầm tích lơ lửng, đối với chúng đòi hỏi tính toán nồng độ và kích thớc hạt tại cao độ công trình lấy nớc. 8.2. Chỉ tiêu lơ lửng và kích thớc hạt Kiến thức Đối với các hạt đang ở trạng thái lơ lửng, vận tốc chìm lắng của chúng phải nhỏ hơn thành phần rối thẳng đứng của vận tốc, liên quan đến u * . Điều này dẫn đến chỉ tiêu đối với ngỡng lơ lửng của trầm tích, đợc cho xấp xỉ bằng quan hệ sau đây: ss* wu (96) trong đó u *s = vận tốc ma sát lớp đệm w s = vận tốc chìm lắng hạt (xem mục 8.3) Đối với trầm tích hỗn hợp, phơng trình (96) có thể áp dụng cho mỗi nhóm kích thớc hạt. Nếu vật liệu đáy có cấp phối rộng, chỉ có các nhóm mịn hơn là lơ lửng, còn các nhóm thô hơn chuyển động nh trầm tích di đáy. Quy trình tốt nhất trong trờng hợp này là chia trầm tích ra một số nhóm kích thớc hạt, mỗi nhóm gồm một dải hẹp các đờng kính hạt và xử lý riêng biệt cho từng nhóm. Một cách tiếp cận đơn giản hơn, nhng ít chính xác hơn, là chọn một kích thớc hạt đơn lẻ đại diện cho toàn bộ mẫu. Acker và White (1973) thấy rằng d 35 của trầm tích đáy cho dự báo tốt nhất về suất vận chuyển trầm tích tổng cộng trong sông. Van Rijn (1984) liên hệ đờng kính trung vị hạt lơ lửng d 50,s với đờng kính trung vị hạt đáy d 50,b thông qua tham số chọn lọc s = 0,5(d 84 /d 50 +d 50 /d 16 ) và tham số vận chuyển T s = crcros bằng quan hệ: )25)(1(011,01/ ,50,50 ssbs Tdd với 0 < s T < 25 (97a) 114 = 1 với 25 s T . (97b) Phơng trình (97) chỉ hiệu lực khi ( s -1) < [0,011(25-T s ) ] -1 ; nếu không nó sẽ dự báo d 50 , s < 0. Fredsoe và Deigaard (1992) loại bỏ tất cả các hạt có w s > 0,8u *s khỏi quá trình lơ lửng, và lấy đờng kính trung vị hạt của phần còn lại làm kích thớc hạt đại biểu cho lơ lửng. Từ đo đạc hiện trờng Whitehouse (1995) thấy rằng đờng kính trung vị hạt lơ lửng d 50,s thích hợp với các giá trị giữa d 2 và d 15 của trầm tích đáy (thô hơn đối với dòng chảy mạnh hơn), với d 10 là giá trị tiêu biểu. 8.3. Vận tốc chìm lắng Kiến thức Vận tốc chìm lắng (hoặc vận tốc rơi, vận tốc kết thúc) của hạt cát trong nớc xác định bằng đờng kính hạt, mật độ của chúng và độ nhớt của nớc. Tại phía mịn nhất của dải đờng kính hạt cát (d = 0,062mm), các hạt chìm lắng theo định luật Stokes về sức cản do nhớt; tại phía thô nhất (d = 2mm) chúng tuân theo định luật sức cản khối đứng bậc hai; và kích thớc trung bình chịu tác động hỗn hợp giữa sức cản do nhớt và sức cản khối đứng. Sức cản lên hạt cát có hình dạng tự nhiên không đều đơn giản hơn so với hạt hình cầu, bởi vì bề mặt góc cạnh và các thay đổi hình dạng giữa các hạt có xu hớng tạo ra quá trình phân tách dòng chảy dần dần hơn. Nh vậy, tốt hơn hết là không coi các hạt cát nh hình cầu đối với các tính toán nh thế. Có một vài công thức để tính toán vận tốc chìm lắng w s cho các hạt cát tách biệt trong nớc tĩnh. Một trong số chúng đòi hỏi phải tính toán kích thớc hạt phi thứ nguyên D * : 50 31 2 1 d )s(g D / * (98) trong đó g = 9,81ms -2 = gia tốc trọng trờng = độ nhớt động học của nớc d = đờng kính sàng trung vị của các hạt s = tỷ lệ giữa mật độ hạt và nớc . Công thức Gibbs và nnk (1971) đối với hạt hình cầu: d dsgd w s 074405,0011607,0 3024801,010869,319 2/1 322 (99) trong đó mọi đơn vị tính trong hệ CGS, trớc đây đợc sử dụng rộng rãi cho cát (nhng không thích hợp). Công thức của Hallermeier (1993) đối với cát tự nhiên là: 115 d D w * s 18 3 39 3 * D (100a) d D w s 6 1,2 * 43 * 1039 D (100b) d D w s 5,1 * 05,1 63 * 4 10310 D (100c) Công thức của Van Rijn (1984) đối với cát tự nhiên là: d D w * s 18 3 187,16 3 * D (101a) 101,01 10 2/1 3 * D d w s 16187187,16 3 * D (101b) d D w s 5,1 * 1,1 16187 3 * D (101c) Zanke (1977) cũng đa ra phơng trình (101b). Soulsby dẫn xuất công thức sau đây đối với cát tự nhiên, dựa trên việc tối u hoá 2 hệ số trong định luật kết hợp độ nhớt với sức cản khối đứng theo số liệu của các hạt không đều: 36,10049,136,10 2/1 3 * 2 D d w s . SC (102) Một so sánh các dự báo vận tốc chìm lắng bằng 4 công thức này (phơng trình (99-102)) đợc thực hiện trên tập hợp lớn số liệu của 115 số đo vận tốc chìm lắng của cát tự nhiên và các hạt trọng lợng nhẹ có hình dạng không đều. Số liệu đợc thu thập và lập bảng bởi Hallermeier (1981). Bảng 11 cho thấy phần trăm dự báo nằm trong khoảng 10% hoặc 20% của các quan trắc. Phơng trình (102) cho kết quả tốt nhất, và cũng là đơn giản nhất. Sự phù hợp tốt này một phần do các hệ số đợc tối u hoá với tập hợp số liệu kiểm tra. Các công thức của Hallermeier và Van Rijn hầu nh tốt, nhng phức tạp hơn, do đó phơng trình (102) đợc kiến nghị. Sự kém cỏi của công thức Gibbs là do nó đợc dự kiến cho hạt hình cầu, không phải hạt tự nhiên. Hình 26 thể hiện hình vẽ phi thứ nguyên của 116 Hình 26. Vận tốc chìm lắng của hạt cát (hình vẽ tổng hợp) w s /[(s-1) ] 1/3 theo D * đối với tập hợp số liệu này, và đờng cong ứng với phơng trình (102). Hình 27 thể hiện các đờng cong w s theo d cho trờng hợp g = 9,81ms -1 , s = 2650kgm -3 , nhiệt độ = 10 0 C, độ muối bằng 35 o/oo, tiêu biểu cho cát trong nớc biển, và đối với trờng hợp cát trong nớc ngọt ở 20 0 C. Với nồng độ cao, dòng chảy xung quanh các hạt đang chìm lắng kề nhau tác động tơng hỗ nên chịu sức cản lớn hơn so với cùng loại hạt tách biệt (chìm lắng bị cản trở). Điều này làm cho vận tốc chìm lắng bị cản trở, w sC ở nồng độ cao sẽ nhỏ hơn so với khi ở nồng độ thấp, w s . áp dụng lý giải tơng tự nh lý giải đã dẫn tới phơng 117 trình (102), nhng giờ đây bao gồm hệ số (1- C) -4.7 trong công thức sức cản hạt (do Wen và Yu nhận đợc bằng kinh nghiệm, 1966), ta có công thức sau đây đối với vận tốc chìm lắng w sC của hạt trong trạng thái lơ lửng dày đặc có nồng độ C: 36,101049,136,10 2/1 3 * 7,4 2 DC d w sC đối với mọi D * và C. SC (103) Phơng trình (103) đơn giản thành phơng trình (102) khi C 0, và nó tơng thích với phơng trình (18) đối với sự lỏng hoá. Bảng 11. So sánh các dự báo vận tốc chìm lắng Công thức Phơng trình 10% 20% Gibbs và nnk 99 35 50 Hallermeier 100 60 89 Van Rijn 101 59 90 Soulby 102 66 90 Đối với các giá trị D * nhỏ, phơng trình (103) cho thấy tỷ số w sC / w s =(1- C) -4,7 , trong khi đối với các giá trị D * lớn, phơng trình (103) cho ta w sC / w s =(1- C) 2,35 . Trong thực tế, chỉ cần xét đến hiệu ứng của vận tốc chìm lắng bị cản trở đối với nồng độ lớn hơn 0,05, thờng xảy ra chỉ trong vài mm tại đáy, vì sự khác biệt giữa w s và w sC nhỏ hơn 20% đối với các nồng độ thấp hơn. Quy trình 1. Để tính toán vận tốc chìm lắng của hạt tại nồng độ thấp, cho kích thớc hạt trung vị d 50 của mẫu cát (nên lấy từ lơ lửng) bằng phân tích sàng (xem mục 2.2). Nếu nhiệt độ và độ muối gần với 10 0 C và 35 o/oo tơng ứng, đọc giá trị w s và d = d 50 từ hình 27, hoặc nếu cần chính xác hơn: Ví dụ 8.1. Vận tốc chìm lắng - Cho nhiệt độ bằng 0 C 10 - Cho độ muối bằng o/oo 35 - Cho mật độ kích thớc hạt s 2650kgm -3 - Cho đờng kính hạt d 0,2mm - Tính toán độ nhớt động học (xem ví dụ 2.1 ) 36 x 10 -6 m 2 s -1 - Tính toán mật độ nớc (xem ví dụ 2.1 ) 1027kgm -3 - Tính toán s s 2,58 118 - Tính toán D * theo phơng trình (98) D * 4,06 - Tính toán vận tốc chìm lắng theo phơng trình (102) w s 0,0202ms -1 Để so sánh, giá trị w s dự báo theo phơng pháp Van Rijn phơng trình (101) là 0,0198 ms -1 và phơng pháp Hallermeier (phơng trình (100)) là 0,0216ms -1 . 2. Nếu nồng độ thể tích lớn hơn 0,05, sử dụng phơng trình (103) thay vào đó. Hình 27. Vận tốc chìm lắng của hạt thạch anh có đờng kính sàng d với nồng độ thấp trong nớc tĩnh 8.4. Nồng độ trầm tích lơ lửng dới tác động dòng chảy Kiến thức Trong cát lơ lửng, sự chìm lắng hạt về phía đáy đợc cân bằng tơng ứng bởi sự khuyếch tán cát ngợc lên do các chuyển động rối của nớc (kể cả thành phần thẳng đứng của vận tốc) gần đáy. Phơng trình điều khiển sự cân bằng này là: 119 dz dC KCw ss (104) trong đó w s = vận tốc chìm lắng của hạt trầm tích C = nồng độ thể tích của trầm tích tại độ cao z K s = độ khuếch tán rối của trầm tích. Khuếch tán rối phụ thuộc vào rối trong dòng chảy và vào độ cao trên đáy. Có thể giải phơng trình (104) để nhận đợc phân bố thẳng đứng của nồng độ trầm tích lơ lửng, phù hợp với các giả thiết cụ thể. Các giả thiết khác nhau về khuếch tán rối của trầm tích dẫn đến các biểu thức khác nhau đối với phân bố nồng độ. Hình dạng của phân bố phụ thuộc vào tỷ số: * ku w b s (105) trong đó b = số Rouse, hoặc tham số lơ lửng k = 0,40 = hằng số von Karman u * = vận tốc ma sát tổng cộng. Chú ý rằng sự cuốn theo của trầm tích từ đáy đợc điều khiển bởi ma sát lớp đệm 2 *sos u , trong khi khuếch tán trầm tích lên cao hơn trong cột nớc đợc điều khiển bởi ứng suất trợt tổng cộng 2 *sos u . Đó là vì sức cản hình dạng của gợn cát không tác động trực tiếp lên các hạt nằm trên bề mặt của đáy, nhng nó tạo ra rối điều khiển quá trình khuếch tán. Sự khác biệt không còn đối với điều kiện dòng trầm tích sát đáy (xem mục 72) trong đó u * = u *s . Trầm tích lơ lửng do dòng chảy trong biển cũng tơng tự nh trong sông. Nếu khuếch tán rối đợc giả thiết tăng tuyến tính theo độ cao trên đáy (K = k z u * z), phân bố nồng độ tơng ứng là phân bố theo định luật số mũ: b a a z z CzC . SC (106) Nếu khuếch tán rối đợc giả thiết biến đổi theo độ cao bằng dạng parabôn K * = ku * z[1-(z/h)], nhận đợc phân bố Rouse: b a a a zh zh z z CzC . SC (107) Nếu khuếch tán rối đợc giả thiết biến đổi với độ cao bằng dạng parabôn ở nửa dới, và không đổi ở nửa trên của cột nớc, nhận đợc phân bố Van Rijn (1984). Thấy rằng Van Rijn cũng xét đến sự khác biệt giữa khuếch tán trầm tích và động lợng chất lỏng, và sự phân tầng mật độ do trầm tích, bằng cách đa ra một dạng sửa đổi có số mũ b: 120 ' b a a a zh zh z z CzC với 2 h zz a SC (108a) )(exp ' ' 2 1 4 h z b zh z CzC b a a a với hz h 2 SC (108b) 2 1 ' B B b b SC (108c) 2 * 1 21 u w B s với 110 u w s , SC (108d) =2 với 1 * u w s 4,0 8,0 * 2 65,0 5,2 as C u w B với 1010 * , u w s . SC (108e) Trong các phơng trình (106) đến (108): z = độ cao trên đáy biển z a = độ cao tham chiếu gần đáy biển C(z) = nồng độ trầm tích tại độ cao z C a = nồng độ trầm tích tham chiếu tại độ cao z a h = độ sâu nớc b = số Rouse (phơng trình (105)). Nồng độ có thể biểu thị là thể tích/thể tích hoặc khối lợng/thể tích, thể hiện rằng C(z) và C a có cùng đơn vị (xem mục 2.3). Một so sánh hình dạng của 3 phân bố đối với trờng hợp b = 1 với z a /h = 0,01 đợc đa ra trên hình 28a. Phân bố Rouse đợc sử dụng rộng rãi nhất, đặc biệt trong sông. Các ví dụ của phân bố Rouse với các giá trị b khác nhau đợc cho trên hình 28b. Đối với các hạt mịn và dòng chảy mạnh (b nhỏ) trầm tích đợc xáo trộn mạnh trong toàn bộ cột nớc, trong khi đối với các hạt thô và dòng chảy yếu (b lớn), trầm tích tập trung chủ yếu gần đáy. Phân bố Rouse ít phù hợp hơn khi sử dụng trong biển bởi vì độ nhớt rối dạng parabôn giảm tới không tại bề mặt, làm cho nồng độ trầm tích dự báo bằng không tại mặt nớc. Điều này ngợc lại với quan trắc, đặc biệt nếu có sóng. Cả phân bố Van Rijn và theo định luật hàm mũ đều có khuếch tán khác không tại mặt nớc, và đơng nhiên là phù hợp hơn. Định luật số mũ hấp dẫn do sự đơn giản của thao tác toán học, đặc biệt bởi vì độ chính xác lớn xuất hiện ở phần thấp hơn của cột nớc, nơi nồng độ lớn nhất. Phân bố Van Rijn có lẽ phù hợp với số liệu tốt nhất, và đợc kiến nghị sử dụng chỉ đối với dòng chảy trong biển. 121 Hình 28. Phân bố nồng độ trầm tích lơ lửng (a) so sánh 3 công thức (trục lôgarít ), b) biến đổi kích thớc hạt và vận tốc dòng chảy (trục tuyến tính) Nồng độ tham chiếu C a , và độ cao tham chiếu z a phải đợc xác định để các phơng trình (106) -(108) cho dự báo nồng độ tiện lợi nhất. Một số biểu thức cho 122 chúng đang có. Garcia và Parker (1991) kiểm tra 7 trong số chúng theo một tập hợp lớn số liệu, và kết luận rằng 2 công thức tốt nhất là: 1. Smith và McLean (1977) s s a T T C 0024,01 00156,0 SC (109) tại độ cao 121 3,26 50 d sg T z scr a 2. Van Rijn (1984) 3,0 * 2/3 015,0 Dz dT C a s a SC (110) tại độ cao z a = s /2, với s cho bởi phơng trình (83) và giá trị nhỏ nhất của z a = 0,01h. 3. Biểu thức gần đây hơn cả, cũng cho kết quả tốt là của Zyserman và Fredsoe (1994): 75,1 75,1 045,0720,01 045,0331,0 s s a C SC (111) tại độ cao z a = 2d 50 . Trong các phơng trình (109) - (111): C a = nồng độ (thể tích/thể tích) tại độ cao z a z a = độ cao tham chiếu os = ứng suất trợt ma sát lớp đệm tại đáy cr = ngỡng ứng suất trợt tại đáy đối với chuyển động trầm tích T s = crcros / d 50 = đờng kính hạt trung vị h = độ sâu nớc g = gia tốc trọng trờng = mật độ nớc s = mật độ vật liệu trầm tích / s s = độ nhớt động học của nớc D * = 50 3/1 2 1 d sg 50 1 dsg os s = tham số Shields ma sát lớp đệm [...]... 1027 9 ,81 2, 58 1 12 0,00156 4,41 1 0,0024 4,41 1,30 x 1 0-3 m 6 ,81 x 1 0-3 Ca - Tính toán số Rouse 0,0202 ws u* 0,40 0,0623 b 0 ,81 1 - Chọn độ cao tại đó yêu cầu tính toán nồng độ z 1,0m - Tính toán nồng độ tại độ cao z bằng cách sử dụng phương trình (106) 6 ,81 10 3 1,0 1,30 10 3 0 ,81 1 3,1 x 1 0-5 C(z) - Nếu muốn chuyển đổi thành kgm-3, nhân với 2650kgm-3 (xem bảng 4) CM(z) 0, 082 kgm-3 2... Nếu s < 0 ,8 , đáy gợn cát (cũng xem có sóng cát không, trong ví dụ này là không gợn cát) - Tính toán hệ số ma sát tổng cộng CD trên đáy gợn cát, sử dụng phương trình (37) với z0= 0,006m (bảng 7 đối với gợn cát) CD 0,00 388 - Tính toán ma sát lớp đệm tổng cộng 1 2 u* C D U 0,0623ms-1 - Tính toán ngưỡng ứng suất trượt cr (xem ví dụ 6.3) Ts - Tính toán - Tính toán độ cao tham chiếu 123 0,176ms-1 os ... + góc giữa hướng sóng và dòng chảy T - Tính toán vận tốc chìm lắng ws 10m 1027kgm-3 1,0ms-1 0,310ms-1 6s 0o 0,0202ms-1 - Tính toán ứng suất ma sát lớp đệm tại đáy, sử dụng phương pháp DATA13 (xem ví dụ 5.1) với độ nhám liên quan đến hạt z0s=d50/12 0,0167mm - ứng suất trượt ma sát lớp đệm trung bình tại đáy ms 2,53Nm-2 - ứng suất trượt ma sát lớp đệm lớn nhất tại đáy max, s 2,95Nm-2 - Vận tốc ma... 1,74 x 1 0-3 0,10m - Tính toán nồng độ tại độ cao z phương trình (115b) C(z) 8, 82 x 1 0-4 - Chuyển đổi thành nồng độ khối lượng s C z CM(z) 2,34kgm-3 Giá trị CM(z) = 2,34kgm-3 tại z = 0,10m so với giá trị 1,19kgm-3 tại cùng độ cao trong ví dụ 8. 3 đối với cùng đặc trưng sóng không có dòng chảy Tính toán tương tự sử dụng z = 1,0m trong ví dụ 8. 4 cho ta CM(z) = 0,225kgm-3, so với giá trị 0, 082 kgm-3 tại... xét đến hiệu ứng gợn cát bằng cách sử dụng ứng suất trượt ma sát lớp đệm tại đáy khi tính toán Ca, nhưng sử dụng ứng suất trượt tổng cộng tại đáy khi tính toán bmax và bm 3 Một so sánh có thể thực hiện với các dự báo theo chương trình TRANSPOR của Van Rijn, với đầu vào tương tự như trong ví dụ 8. 4 (h =10m, U = 1,0ms-1, Hs= 1,0m, Tp= 6s, =0o, d50= 0,2 x1 0-3 m, d90= 0,23 x1 0-3 m, d50,s= 0,2 x1 0-3 m, 10oC,... gợn cát, ví dụ 7.3 cho ta r 0,0452m r 0,265m r - Tính toán theo phương trình (113d) 0 ,84 9 0,0202ms-1 - Tính toán ws (xem ví dụ 8. 1) ws - Tính toán Uw/ ws 49,5 - Tính toán l từ phương trình (113b) l 0,0633m C(z) 4, 48 x 1 0-4 - Tính toán nồng độ tại độ cao z bằng phương trình (112) - Hoặc khối lượng trên thể tích s C (z ) CM(z) 1,19kgm-3 2 Để so sánh, nếu giả thiết dòng trầm tích sát đáy trên đáy gợn cát. .. s 2, 58 (Các giá trị ví dụ như trong ví dụ 7.3) - Tính toán giá trị U w T /(5d 50 ) r 592 - Tính toán hệ số ma sát sóng phương trình (60b) fwr 0,01 18 ws 0, 183 - Tính toán biên độ tham số Shields liên quan đến hạt Nếu ws < 0 ,8, đáy gợn cát và phải sử dụng phương trình (112) và (113) Nếu ws > 0 ,8, đáy phẳng với dòng trầm tích sát đáy và phải sử dụng phương trình (114) - Trong ví dụ này ws < 0 ,8 do... kết hợp sóng và dòng chảy vẫn còn hiểu biết một cách nghèo nàn, do vậy kết quả cần được xử lý thận trọng Các kết quả phụ thuộc nhiều vào phương pháp sử dụng, các giả thiết đã thực 1 28 hiện, và các tham số đầu vào như kích thước hạt lơ lửng và nhiệt độ nước Kiến nghị khi phân tích độ nhạy phải kể đến tất cả các điều này khi thực hiện dự báo cho các ứng dụng thực tế Quy trình 1 Quy trình tính toán nồng... (114) Các giá trị thực nghiệm của b không nhất thiết tuân thủ các dự kiến lý thuyết, và tiêu biểu là b =1,7 đối với d50= 0,13mm và b = 2,1 đối với d50= 0,21mm, đối với tất cả các điều kiện sóng trong phạm vi 0,4 < Uw < 1,3ms-1 và 5 < T < 12s Nếu za= 1cm, giá trị của Ca 4 x 1 0-4 đối với Uw= 0,5ms-1 và Ca= 8 x 1 0-3 đối với Uw= 1,3ms-1 được thấy bằng thực nghiệm đối với d50= 0,21mm (Ribberink và Al-Salem,... minh hoạ bằng cách kết hợp các giá trị đầu vào của dòng chảy ổn định sử dụng trong ví dụ 8. 2 đồng tuyến ( 0 0 ) với giá trị đầu vào của sóng sử dụng trong ví dụ 8. 3 Xem các ví dụ này đối với một số bước: Ví dụ 8. 4 Nồng độ do sóng và dòng chảy kết hợp - Cho các giá trị của: + độ sâu nước + mật độ nước (10oC, 35o/oo) h + đường kính hạt trung vị d50 0,2mm + mật độ của hạt trầm tích s 2650kgm-3 + vận tốc . 2650kgm -3 - Cho đờng kính hạt d 0,2mm - Tính toán độ nhớt động học (xem ví dụ 2.1 ) 36 x 10 -6 m 2 s -1 - Tính toán mật độ nớc (xem ví dụ 2.1 ) 1027kgm -3 - Tính toán s s 2, 58. - Tính toán nồng độ tại độ cao z bằng cách sử dụng phơng trình (106) 81 1,0 3 3 1030,1 0,1 1 081 ,6 C(z) 3,1 x 10 -5 - Nếu muốn chuyển đổi thành kgm -3 , nhân với 2650kgm -3 . bảng 4) C M (z) 0, 082 kgm -3 2. Để so sánh: - sử dụng phân bố Van Rijn C a , z a 0,0616kgm -3 - sử dụng phân bố Rouse và Zyserman và Fredsoe C a , z a 0,121kgm -3 3. Các phơng pháp khác