TT N u P và Q thì n và k 1 Th p Th p Trung bình Th p 2 Trung bình Th p Trung Bình Th p 3 Cao Th p Trung bình Th p 4 R t cao Th p Cao Th p 5 Th p Trung bình Trung Bình Th p
6 Trung bình Trung bình Cao Th p
7 Cao Trung bình Cao Th p
8 R t cao Trung bình Th p Th p
9 Th p Cao Trung bình Trung bình
10 Trung Bình Cao Trung bình R t cao
11 Cao Cao Cao Cao
12 R t cao Cao R t cao Cao
Các lu t m là c s đ ra quy t đ nh cho bi n ngõ ra n và k v y nên thi t l p lu t m t t thì k t qu đ a ra s có đ chính xác cao. T các 3 t p m y u t áp l c P, 3 t p m y u t l u l ng Q thì t ng s lu t m là 4 x 3 = 12 lu t m . Thi t l p lu t m trong công c FLD d ng nguyên nhân - k t qu “N u…và….Thì” nh B ng 3.8. M nh đ k t qu n và k đ c xác đ nh d a trên k t qu phân tích h i quy h s dịng ch y rò r k theo ba y u t P, Q và n.
65
Xây d ng mơ hình HSM cho DMA0203
Mơ hình HSM xây d ng b ng cơng c SIMULINK trong Hình 3.31 g m hai kh i ngu n
đ u vào (P,Q), ph n nh n tín hi u Mux (hai ngõ vào và m t ngõ ra) đ truy n thông tin t i lõi x lý logic m (FLC - Fuzzy Logic Controller) và m t ngõ ra g m hai bi n n và
k đ c xu t tr c ti p ra môi tr ng làm vi c c a MATLAB.
Hình 3.31. Mơ hình HSM trong SIMULINK – MATLAB
Tr c khi m SIMULINK c n nh p hai t p d li u áp l c P và l u l ng n c rò r Q c a 30 đi m b vào trong mơi tr ng MATLAB. Qtrình làm vi c c a mơ hình HSM nh sau:
- Kh i ngu n s đ c tr c ti p hai t p d li u P và Q t môi tr ng làm vi c MATLAB.
- Nh p tên c a t p logic m đư l p trình trong cơng c FLD vào lõi x lỦ FLC.
- Thi t l p s l ng d li u đ u ra s đ c l u l i sau m i l n ch y mơ hình. - t tên t p đ u ra là nk đ l u l i trên môi tr ng làm vi c c a MATLAB.
- Ch y mơ hình và chuy n qua mơi tr ng làm vi c c a MATLAB đ l y k t qu n, k Xu t k t qu c a mơ hình HSM ra file exel và so sánh v i s li u th c t n u c n thi t có th đi u ch nh l i các thông s trong t p m đ có k t qu phù h p h n v i s li u kh o sát.
66
xu t h s dòng ch y rò r và th o lu n
Mơ hình HSM đ a ra h s m n và h s dòng ch y rò r k c a 30 đi m b trên MLCN
DMA0203. Giá tr c a n n m trong kho ng t 0,91 đ n 2,25 v i trung v là 0,99, trung
v ph n ánh s đông c a t p s li u nên nghiên c u l a ch n n = 0,99 đ đ a vào mơ hình th y l c MLCN. S l ng h s k t ng ng v i 30 đi m b trên MLCN, mơ hình đ a ra k n m trong kho ng t 0,0094 đ n 0,0452.
Hình 3.32. K t qu mơ hình HSM cho MLCN DMA0203
Các đi m rò r đ c coi nh các nút trong mơ hình WaterGEMS, v trí c a 30 đi m rò r đ c thi t l p trên mơ hình nh Hình 3.33. phân bi t v icác nút tiêu th (J-..) trên
MLCN thì nhưn nút (Label) c a các đi m rò r đ c đánh s theo tuy n đ ng xu t hi n đi m b (1.1; 1.2;….). 0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5 5.0 1 2 3 4 5 6 7 8 9 101112131415161718192021222324252627282930 0.0 0.5 1.0 1.5 2.0 2.5 3.0 k x 10-2 k n n
67 Hình 3.33. Mơ hình MLCN DMA0203 mơ ph ng rị r
68
H s dòng ch y rò r k t ng ng v i 30 nút rò r đ c thi t l p trong mơ hình b ng cách ch n Tab View, vào m c Flex Tables và nh n chu t trái vào m c Junction Table thì màn hình s hi n th B ng nh Hình 3.34. Trong c t “Emitter Coefficient (L/s.mH2O^n)” có m u xanh nh trong Hình 3.34, nh p k t qu c a mơ hình HSM v h s k c a 30 nút rò r . Các nút tiêu th trên MLCN (kí hi u là J-…) thì h s này b ng
khơng.
Trong mơ hình WaterGEMS, giá tr h s m rò r n đ c cài đ t b ng cách l a ch n m c Calculation Options trên Tab Analysis, nh p chu t vào ph n Base Calculation và cho h s Emitter Exponent (là h s m rò r n) b ng 0,99 nh Hình 3.35.
69
So sánh giá tr l u l ng n c rò r c a mơ hình và k t qu l u l ng kh o sát cho k t qu đ l ch ∆Q nh B ng 3.10 n m trong kho ng t 90% - 100%. Theo khái ni m t ng quan c a A.G. Asuero và c ng s thì giá tr t ng quan t 0,9 đ n 1 có th xem s li u có m i liên h r t ch t ch (A.G.Asuero và c ng s , 2006), nh v y k t qu mơ hình thu l c là hoàn toàn phù h p. i u này có ngh a là các h s (n, k) đ a vào mơ hình thu l c đư mô ph ng đúng giá tr l u l ng n c rò r trên m ng l i.