2.2.3 Xỏc định khoảng giỏ trị và hàm phõn bố của cỏc thụng số
Trong phƣơng phỏp GLUE, việc xỏc định khoảng giỏ trị của cỏc tham số là cần thiết. Độ rộng của dải giỏ trị phải phự hợp. Nếu rộng quỏ sẽ dẫn đến những mụ phỏng khụng cần thiết, ngƣợc lại sẽ bỏ qua nhiều giỏ trị của tham số. Trong bƣớc này cần phải chỳ ý:
- Lựa chọn tham số: cần xem xột cỏc tham số nào cú ảnh hƣởng thực sự đối với điều kiện thực tế ỏp dụng. Cú thể dựng phõn tớch độ nhạy để lựa chọn cỏc tham số này.
- Xỏc định khoảng giỏ trị và hàm phõn bố cỏc tham số: thụng thƣờng cỏc khoảng giỏ trị của cỏc tham số đƣợc xỏc định từ cỏc nghiờn cứu trƣớc đõy và đối với hàm phõn bố thƣờng đƣợc lấy là hàm phõn bố đều khi ta khụng biết nhiều về giỏ trị của chỳng.
- Phƣơng phỏp chọn mẫu: cú hai phƣơng phỏp chớnh là chọn mẫu ngẫu nhiờn (Monte Carlo) và chọn mẫu theo phƣơng phỏp siờu lập phƣơng Latin (Latin Hypercube Sampling - LHS).
2.2.4 Thiết lập quy trỡnh sử dụng chỉ tiờu đỏnh giỏ độ phự hợp để tớnh toỏn khoảng bất định khoảng bất định
Sau khi xỏc định chỉ tiờu đỏnh giỏ và khoảng giỏ trị ban đầu của tham số,
sử dụng phƣơng phỏp phõn tớch Monte Carlo tớnh toỏn với nhiều bộ thụng số. Trong thực tế nếu thời gian mụ phỏng dài (đối với cỏc mụ hỡnh phõn phối), ngƣời ta thƣờng sử dụng phƣơng phỏp chọn LHS để tăng hiệu quả của quỏ trỡnh tớnh toỏn. Thực chất LHS là phƣơng phỏp Monte Carlo cải tiến. Trong luận văn này, phƣơng phỏp LHS đƣợc sử dụng và số mẫu mụ phỏng là 200.
Với mỗi bộ thụng số đƣợc tạo ra bởi phƣơng phỏp LHS, mụ hỡnh WetSpa sẽ tớnh toỏn đƣợc lƣu lƣợng dũng chảy ra. Từ đú giỏ trị của chỉ tiờu đỏnh giỏ đó chọn đƣợc tớnh toỏn. Cú thể sử dụng một trong ba chỉ tiờu đỏnh giỏ nờu trờn. Trong nội dung thực hiện luõn văn cả ba chỉ tiờu đỏnh giỏ đƣợc thử nghiệm. Do quy trỡnh tớnh
48
toỏn chỉ ỏp dụng với một chỉ tiờu nờn sau đõy là mụ tả cỏc bƣớc tiếp theo khi sử dụng một chỉ tiờu duy nhất.
Khi cỏc giỏ trị của chỉ tiờu đỏnh giỏ đƣợc xỏc định, cỏc mụ phỏng đƣợc chấp nhận (NSi > 0.7) đƣợc giữ lại để tớnh toỏn khoảng bất định bằng việc sử dụng giỏ trị ngƣỡng loại bỏ. Giỏ trị của cỏc chỉ tiờu của cỏc mụ phỏng đƣợc chấp nhận sau đú đƣợc biến đổi sao cho tổng của chỳng bằng một đơn vị theo cụng thức sau:
RLi = Li/ (L1 + L2 +… + Ln) (2.35)
Trong đú:
RLi là giỏ trị biến đổi chỉ số độ phự hợp của mụ phỏng thứ i,
Li giỏ trị của chỉ số phự hợp của mụ phỏng thứ i,
L1 and L2 tƣơng ứng là giỏ trị của cỏc chỉ số phự hợp của cỏc mụ phỏng
đƣợc chấp nhận thứ nhất và thứ hai
LN là giỏ trị của chỉ số của mụ phỏng cuối cựng đƣợc đỏnh giỏ là phự
hợp khi dựng giỏ trị ngƣỡng loại bỏ.
Ở mỗi bƣớc thời gian, giỏ trị lƣu lƣợng ứng với 5% và 95% của hàm phõn bố lũy tớch cỏc chỉ số phự hợp đƣợc sử dụng làm khoảng bất định của giỏ trị dự bỏo.
Lƣu lƣợng Q5% và Q95% đƣợc xỏc định bởi cụng thức (5). ) ( % % nna nnb nnb nna nnb n nnb n Q Q CL CL CL CL Q Q (2.36) Trong đú:
Qn% là lƣu lƣợng tƣơng ứng với n% của hàm phõn bố lũy tớch cỏc chỉ số phự
hợp;
CLn% là giỏ trị của chỉ số phự hợp tƣơng ứng với n% của hàm phõn bố lũy
tớch cỏc chỉ số phự hợp;
CLnna, CLnnb tƣơng ứng là giỏ trị của cỏc chỉ số phự hợp ngay trờn và dƣới
giỏ trị CLn%;
49