Bài tập lớn môn học: Thuỷ văn công trình.

TÀI LIỆU THAM KHẢO - BÀI TẬP LỚN MÔN HỌC THỦY VĂN CÔNG TRÌNH CỦA TRƯỜNG CĐ THUỶ LỢI BẮC BỘ

NHIỆM VỤ BÀI TẬP LỚN

Môn học: Thuỷ văn công trình

Sinh viên: Phạm Văn Khôi

Ngày giao đề: 20/10/2007

Ngày nộp đề: 10/11/2007

Đề tài:Vẽ đường tần suất thực nghiệm theo công thức vọng số và công thức trung bình

Vẽ đường TSTN Piếc Sơn III, đường TSTN KrisKy-MenKen(K-M) và đường TSTN theo phương pháp 3 điểm:P1=10%, P2 =50%, P =90%

1)SỐ LIỆU BAN ĐẦU

1.1)Tài liệu về mực nước:

1 Nhập số liệu về mực nước và số liệu liên quan

Vẽ đường TSTN

a)Cơ sở lý luận

-Khái niệm về đường TSTN

2 -Xây dựng theo công thức : P= n

+Phân cấp tài liệu

+Xác định tần số rơi vào từng cấp(mi)

Không an toàn? giải thích, cho ví dụ cụ thể

d)Giải thích vấn đề kéo dài và hiệu chỉnh

+Phân cấp tài liệu

+Xác định hệ số tách rời  (khoảng lệch tung

độ)

+)Tính hệ số biến suất KP=  *Cv +1

+)Xác địnhcác trị số HP =H *KP (kết qủa tính

toán lập thành bảng)

+)Dựa vào các cặp trị số HP ~ P vẽ đường tần

suất lý luận Piếc Sơn III

-Khái niệm tương quan

-Phân loại tương quan

-Tương quan đường thẳng(lập phương trình

đường hồi quy

2.2.1 Nêu ngắn gọn lý thuyết áp dụng trước khi tính

2.2.2 Đường TSTN theo công thức P = 1



n

m

*100% vẽ màu đỏ, các đườngTSTN khác vẽ màu xanh

2.23 Hai đường TSTN được vẽ trên 1 trang giấy xác suất

2.2.4)Vẽ 3 đường TSLL (PIII,K-M,3điểm) trên 1 trang giấy trong đó đương TSLL theo phương pháp 3 điểm vẽ màu đỏ

Vẽ hai đường TSLL(PIII, 3 điểm) và đường TSLL theo công thứcP=

*100% vẽ màu đỏ

2.2.5 Các hình vẽ minh hoạ, bảng biểu , đồ thị phải có tên , đánh số

thứ tự

2.2.6 Các công thức phải được đánh số thứ tự

2.2.7 Nếu áp dụng tin học trong tính toán, phải đưa vào phụ lục

2.2.8 Thuyết minh khổ A4, bìa Nilon, các đồ thị vẽ trên giấy kẻ

ly(khuyến khích làm bằng vi tính) bao gồm các trình tự như sau:

-Bìa ngoài

-Nhiệm vụ bài tập lớn

-Muc lục

-Nội dung tính toán

-Phụ lục tính toán(nếu áp dụng tin học)

- Tài liệu tham khảo

2.3 Thưởng , phạt

2.3.1 Thưởng

-Ap dụng tin học………

-Nộp sớm:………

-Trình bày đẹp , đúng quy cách:………

-Lý do khác:………

2.3.2 Phạt: -Không áp dụng tin học

-Chậm tiến độ

-Trình bày xấu, không đúng quy cách ………

-Lý do khác ………

2.3.3 Đánh giá của giáo viên hướng dẫn: ………

1.Tài liệu về mực nước

Bảng 1.1: Số liệu mực nước đỉnh triều cao thang 6 năm 2002

Tuy nhiên, ta chưa thể dùng được số liệu này để tính toán Do một số chỗ

trong bảng số liệu trên có sai sót.Nếu ta bỏ đi thì số liệu của ta không đủ,

nên ta quyết định là sửa lại số liệu theo xu thế của mực nước là đang lên hay

đang xuống dựa vào các số liệu xung quanh!

Bảng số liệu này sẽ được chỉnh sửa lại cho các khoá sau!

_Đường TSTN là đường TS được xác định từ mẫu số liệu thực nghiệm(quan trắc)của một đại lượng ngẫu nhiên biểu thị mối quan hệ hàm số giữbiến ngẫu nhiên nghiên cứu và xác suất tích luỹ tương ứng, đồng thời thểhiện cụ thể quy luật thống kê của tập hợp mẫu

- Xây dựng ĐTSTN theo công thức P = m n *100%

+Thiết lập tập hợp mẫu cần nghiên cứu

+Xắp sếp tập hợp mẫu theo thứ tự từ lớn đến nhỏ

+Phân cấp liẹt tài liệu nếu có

+Tìm tần số mi

+Tìm tần suất tích luỹ Pi(x xi)

+Vẽ đường TSTN dựa vào các điểm thực nghiệm của biến ngẫunhiên Đường TSTN được vẽ trên giấy xác suất để tránh tình trạng quá dốc ởhai đầu đường TS

VD : Vẽ đường tần suất thực nghiệm

Tài liệu của những dòng trung bình tại một chạm thuỷ văn trong 20 năm ;176; 212; 234; 147; 288; 215; 262; 250; 192; 167; 284; 264; 275; 213; 188;221; 242; 189; 245; 196

STT

xắp

sếp mi Pi(%) Pi(%) STT

xắpsếp mi Pi(%) Pi(%)

Nhận xét : trong 20 năm xuất hiện một lần Qmax  288m3/s P(Q  288)

=5% Trong 20 năm xuất hiện 20 lần chỉ số Q  147 m3/s, P ( x 147) =

100%

Vẽ đường TSTN

-Các công thức tính đường TSTN

Ta thấy tính theo công thức P= n

m

*100% thì ứng với trị số bé nhấtcủa mẫu bao giờ cũng xuất hiện 100% nghĩa là sau này không có trị số nào

bé hơn nó nữa nên chỉ ứng dụng khi n  Để khắc phục nhược điểm đó

người ta đưa ra một số công thức sau

11 183 min max

m

Vậy m =5

+)Tần suất rơi vào từng cấp (mi) được xác định trong Bảng 5:

BẢNG 5 TÍNH TOÁN ĐTSTN

m n



) *100%

(xxi)

Dựa vào đường tần suất thực nghiệm ta thấy ,đường tần suất thực nghiêm tính theo công thức vọng số an toàn hơn cho việc tính toán các công trình chống lũ vì :ở cùng một mực nước cao thì tần suất xuất hiện của mực nước tính theo công thức vọng số lớn hơn tính theo công thức trung bình, đối với công trình chống hạn thì sử dụng công thức trung bình an toàn hơn vì : ở cùng một mực nước thấp thì tần suất xuất hiện của mực nước tính theo côngthức vọng số nhỏ hơn tính theo công thức công thức trung bình.

4) GIẢI THÍCH VỀ VẤN ĐỀ KÉO DÀI VÀ HIỆU CHỈNH ĐƯỜNG TSTN

Sau khi vẽ song đường TSTN ta không thể dùng ngay đẻ tính toán các đặctrưng thuỷ văn vì đường TSTN còn co một số nhược điểm sau :

- Bản thân đường TSTN dựa trên cơ sở thí nghiệm của tạp hợp mẫu , màtập hợp mẫu chưa thể đại diẹn cho quy luật ngẫu nhiên của tổng thể

- Đường TSTN vẽ trên cơ sở thí nghiệm thực nghiệm mang nhiều sai sốtrong quá trình đo đạc, quan trắc và sai số chủ quan trong khi vẽ đồ thị

- Đường TSTN bị hạn chế ở hai đầu đường cong vì ở đó ít hoặc không có

số liệu thực đo ở hai đầu này lại rơi vào những điểm có TS rất nhỏ hoặc rấtlớn mà thực tế lại rất cần

VD: Các công trình chống lũ cầ xác định lưu lượng với TS nhỏ ( 1%,0.1%,0.01% ).Còn các công trình giao thông , cấp nước , phát điện phải xácđịnh với các đặc trưng dòng chảy với TS lớn(95%;98% ) Đường TSTN vẽtren giấy thì hai đầu rất dốc , nếu ta tiến hành kéo dài một cách trực quan thìsai số chủ quan lớn ảnh hưởng đến quy mô , kích thước công trình về đầu tư, an toàn cung như vấn đề khai thác Chính vì những lí do trên nên ta phảikéo dài và hiệu chỉnh đường TSTN theo nguyên tắc sau :

+Vẫn tôn trọng quy luật tự nhiên khách quan dựa vào các số các điểm cần

đo cảu tập hợp mẫu được biểu hiện qua đường TSTN

+ Hạn chế sai số do quan trắc đo đạc và sai số chủ quan do xây dựngđường TSTN đồng thì phải thể hiện nhiều quy luật chung của hiện tượngngẫu nhiên được biểu thị qua đường TSLL

+ Công việc tính toán phải tiến hành thuậ n lợi , dễ sử dụng Xuất phát từnhững nguyên tắc trên , trong tính toán thuỷ văn người ta áp dụng biện phápsau :

Tìm được quy luật thống kê của đại lượng thuỷ văn , kéo dài hiệu chỉnhđường TSTN bằng một đường TS lí luận tương ứng Kết quả đáp ứng được

ta vẽ đường TS tính toán phù hợp với quy luật biến đổi chung của hiệntượng thuỷ văn , đáp ứng nhu cầu xây dựng và sử dụng các công trình thuỷ

CHƯƠNG II:ĐƯỜNG TẦN SUẤT LÍ LUẬN

I-VẼ ĐƯỜNG TSLL PIÊCSON III

1) Cơ sở lí luận

- Quy luật xuất hiện của các đại lượng ngẫu nhiên được biểu thị toán họcqua hàm mật độ xác suất f(xi) hay hàm tích phân phân bố xác suất F/(xi) Nhà sinh vật học Piecson đã đưa ra 13 đường mật độ TS để biểu diễn cácquy luật khác nhau của hiện tượng ngẫu nhiên Các quy luật này nói chung

là không đối xứng, các đại lượng thuỷ văn cũng là các đại lượng ngẫu nhiênkhông đối xứng và phù hợp với dạng đường thứ 3 trong họ13 đường chonên gọi là đường PIII Họ đường cong PIII có dạng dx dy =

xb x b bo

y d x

2

21

)(

*Đặc điểm của PIII đầu trái có giới hạn tồn tại vị trí nhỏ nhất

là xmin, đầu phải là vô hạn( không có giá trị giới hạn là xmax)

2) ứng dụng

- Khi biết 3 tham số yo, a, d thì đường cong TS hòan toàn được xácđịnh Qua phân tích toán học thống kê rút ra được quan hệ giữa 3 tham số

trên với các tham số thống kê như sau: d =

2

) 4

C C

C

C y

S

C S V

C S

4(

)1

4(2

2 1

4 2

4 2

- Để thuận tiện cho việc thực hành toán toan tính toán Poster và Rupkin đãlập thành bảng tra sẵn để vẽ đường TS luỹ tích

+ Đặt hệ số tách rời  ( khoảng lệch tung độ phụ thuộc vào TS P và Cs,

(cs,p) được tra trong bảng phụ lục 1 của đường TSLL PIII

=Kp 1 Cv =f(Cs,P)

+Xác định hệ số biến suất Kp: Kp=Cv+1 (II.6)

+ Tính giá trị tung độ của đường TS Xp: Xp = X *Kp = (Cv+1)X

nhưng phải biến đổi p(cs<0) =- (100-p)(cs >0)

Giới hạn của Cs:2Cv  Cs  1 2K Cvmin

+ Khi các tham số trong liệt là lẻ thì số trung vị là số ở giữa

 Khoảng lệch quân phương  : trong tập hợp mẫu có n giá trị từ x1 xn

số trung bình là X thì ta gọi (Xi -X )là khoảng lệch của trị số Xi với X

khoảng lệch quân phương  được xác định theo công thức sau :

 phản ánh mức độ dao động của các trị số quanh trị số trung bình

 nhỏ thì liệt phân bố tập trung

 lớn thì liệt phân bố phân tán

X Xi



i=1 ,n (II.12)

BẢNG 6 TÍNH TOÁN ĐƯỜNG TẦN SUẤT LÝ LUẬN

1

2

) 1

) 1 (Ki 





i

m i

Ki

1

3

) 1 (

Xác định tham số thống kê: 1( 1)2 78.9541

0.22 360

n

i v

Ki C

Ki C

Từ các kết quả trên, ta lập bảng thống kê và vẽ đường TSLL PEAR III

BẢNG 8 VẼ ĐƯỜNG TSLL PEARSON III

1 

Vẽ đường TSLL PEARSON III trên giấy xác suất

Ki

1

2

) 1

Ki

1

3

) 1 (

n

Ki C

n

i v

3 1

3

) 1 (

v

n

i s

II- VẼ ĐƯỜNG TSLL KRISKY - MEKEN (K - M)

1-Cở sở lí thuyết

Trong thực tế tính toán thuỷ văn vẫn tồn tại một số trường hợp Cs < 2Cv lúc đó

đường PIII có giá trị âm nên không phù hợp.Xuất phát từ đườngPIII với Cs = 2Cv và

một số điều kiện toán học phù hợp với thực tế ,Krisky-Melken đề nghị dùng đường

mật độtần suất K-M để dùng cho từng trường hợp Cs<2Cv Hàm mật độ K-M có

III-VẼ ĐƯỜNG TẦN SUẤT LÍ LUẬN THEO PHƯƠNG PHÁP 3 ĐIỂM

1-Cở sở lí thuyết

Để giảm bớt khối lượng tính toán và thử được nhanh hơn, Alếchxâyép đề

nghị phương pháp đồ giải 3 điểm để xác định ĐTTT ĐTSLL được dùng là

, 3 ( ) , 1 (

) , 2 ( ) , 3 ( ) , 1 ( 2

1

2 3

Cs P Cs

P Cs P

Cs P Cs

P Cs P P

P

P P

X X

X X

Khi P đã chọn thì S chỉ phụ thuộc Cs và  Từ S tra bảng tìm đợc Cs

- Từ Cs ta tìm được X và  như sau:

3 1 )

, 3 ( ) , 1 (

3 1

P P

P P Cs

P Cs P

+ Nếu nó phù hợp với đường TSTN thì kéo dài ĐTSTN theo xu thế của

đường tần suất lý luận ta được một đường TSTT Dựa vào TSTK ta được giá

trị của biến ngẫu nhiên tương ứng với Xp thiết kế

2-Trình tự tính toán

vẽ đường TSLL theo phương pháp 3 điểm(phương pháp đồ giải):P1 =3%,

P2=50%, P3=97%,

Dựa vào đường tần suất TN ta xác định được 3 điểm:

70 155

119

* 2 70 155

% 90

70 155

% 90

% 10

% 90

% 10

81 33

Từ Cs,Cv, , H ta lập được bảng theo PIII

 0.118 -0.142 -0.431 -0.592 -0.791 -1.33 -1.908 -2.147 -2.813 -4.086

1 

* Khái niệm: Khi nghiên cứu các đại lượng thuỷ văn nhiều trường hợpliệt tài liệu mẫu quá ngắn dẫn đến sai số lấy mẫu lớn, không cho phép sửdụng các kết quả trong thực tế vì vậy cần tìm cách phát triển dung lượng củamẫu Công việc đó còn gọi là kéo dài số liệu.

Trong thực nghiệm ta nhận thấy sự biến đổi của đại lượng X nào đó cótương quan đến sự biến đổi của một hoặc một số đại lợng khác Y, Z Việcphân tích xác lập về mặt số lượng của đại lợng X với các đại lợng Y, Z trong

lý thuyết thống kê được gọi là phân tích tương quan

*, Phân loại:

- Tương quan hàm số: X, Y có tương quan chặt chẽ, với mỗi giá trị của X

sẽ có một hay nhiều giá trị của Y hoặc ngược lại

- Không tương quan: quan hệ X, Y rời rạc ta không nhận thấy sự liên hệgiữa chúng hoặc sự thay đổi của đại lượng này không làm thay đổi đại lượngkia

- Tương quan thống kê: mỗi giá trị của X có thể có hay không giá trị của

Y mà ta không thể khảng định trước được Vì Y không những phụ thuộc và

X mà còn phụ thuộc nhiều các giá trị khác nữa

Tuy vậy qua nhiều số liệu thống kê ta có thể nhận thấy một xu thế rõ rệtnào đó gọi là tương quan thống kê

b Tương quan đường thẳng

* Phương pháp giải tích:

- Phương trình biểu diễn quan hệ thống kê giữa hai đại lượng ngẫu nhiên

X và Y được gọi là phương trình đường hồi quy

- Giả sử X và Y là cặp trị số đối ứng nhau quan trắc được của hai biến sốngẫu nhiên X và Y Trong tập hợp mẫu nghiên cứu có n số hạng như vậy ta

có n cặp trị số (X,Y) đối ứng nhau từng đôi một đánh dấu các cặp điểm(X,Y) tạo thành một băng điểm hẹp trên đồ thị Đờng thẳng đi qua trọng tâmbăng điểm đó và phối hợp tốt nhất các điểm (X,Y) trên đồ thị gọi là đườnghồi quy

Đường hồi quy có dạng tổng quát của phương trình đưường thẳng y=ax+b(V.1)

y i  i 



=> y i  yy i  (axb)

) (ax b y

y i  i  



Khoảng lệch này có thể mang dấu âm hoặc dương tuỳ theo các điểm thực

đo nằm trên hoặc dưới đường hồi quy

Theo nguyên lý bình phương tối thiểu, đường thẳng nào đó có tổng bình ương là min sẽ là đường phối hợp tốt nhất các điểm thực đo Nó chính là đ-ường hồi quy

ph-min )]

( [ )

(

2 1

y

(1)Muốn đạt được điều kiện trên ta phải đạo hàm nó với a và b, tức là

0

) (

i

i ;

0

) (

(.2) Giải 2 phương trình trên ta sẽ được:

) )(

(

X X

Y Y X X

a

i

i i

(3)

X X

X

Y Y X X Y

y

b

i

i i

) )(

(

.4)Thay a, b vào phương trình (V-1) ta nhận đợc phương trình đồng quy của ytheo x

) (

) (

) )(

(

X X

Y Y X X

y

y

i

i i

) )(

(

2 1

y y

y y x x n i

i

i n

x,ylà giá trị trung bình của 2 biến ngẫu nhiên X,Y

b, Hệ số tương quan: biểu thị mức độ tương quan giữa hai biến ngẫu nhiên X và Y Tương quan thống kê có thể biểu diễn bằng hai phương trình (V-5) và (V-6) Dạng tổng quát của chúng

Nếu aa1   1 thì X, Y có tương quan hàm số.

Nếu aa1   1 thì X, Y có tương quan thống kê

Nếu aa1   1 càng nhỏ thì X, Y có tương quan kém chặt chẽ

Biểu thị mức độ tương quan giữa 2 biến ngẫu nhiên X và Y Thay a, a1 vào

1

) (

( ) (

) )(

(

y yi x

xi

y y x

x i

n

i i (7)

Qua thực tế tính toán thuỷ văn điều kiện để ứng dụng: Hệ số ≥0.8 và phải >

10 năm quan trắc đồng thời giữa hai biến ngẫu nhiên X và Y

Ưu điểm: Thông qua đờng hồi quy để tính toán không mắc phải sai số chủ

quan, có tiêu chuẩn để đánh giá mức độ tương quan, tránh được việc áp

dụng tuỳ tiện

Nhược điểm: quan hệ X, Y thể hiện bằng hai đườngtương quan Nếu kết quả

của 2 đường không giống nhau thì đó là điều không hợp lý

* Phương pháp tương quan đồ giải

- Trong trường hợp tính toán sơ bộ yêu cầu tính nhanh và mức độ chính xác

không cao ta có thể dùng phương pháp đồ giải

+ Đánh dấu các cặp điểm (X,Y) trên trục toạ độ x0y Nếu các điểm phân

bố thành một băng hẹp ≥0.8 ta có thể dùng một đường thẳng đi qua tâm

băng điểm đó

+ Ta có thể dùng đường quan hệ để kéo dài số hiệu khi biết các giá trị của

X, ta có thể suy ra giá trị của Y và ngược lại

Bài toán không thoả mãn điều kiện để áp dụng phương pháp tương quan đường

thẳng nhưng ở góc độ bài tập vận dụng nên ta vẫn coi như tương quan để ta tính toánmẫu

-Phương trình đường hồi quy của x theo y:

) (

) )(

(

2 1

y y

y y x x x

x

n i i

i n

i i

Từ công thức (5) ta có: ( )

) (

) )(

(

X X

Y Y X X y

y

i

i i