7-1
<Table 7.1> 기조력관련천문상수 항목크기 태양의질량지구의질량달의질량지구의평균반경지구-달중심간평균거리지구-태양중심간평균거리달의질량 / 지구의질량태양의질량 / 지구의질량 1.99 10 33 g 5.98 10 27 g 7.34 10 25 g 6.37 10 8 cm 3.84 10 10 cm 1.49 10 13 cm 1/81.5 3.28 10 5 7-2
예 ) 9 시 10 분의조위를표척관측으로얻고자할경우 9 시 9 분 9 시 10 분 9 시 11 분 고 455.5 456.5 456.5 저 450.0 452.5 453.0 평균 452.75 454.5 454.75 9 시 10 분의조위 = (452.75+454.5+454.75)/3=454.0 7-3
7-4
<Table 7.2> 4 대주요분조 기호명칭 각속도 / 시간 ( /hr) 조화상수기호 반조차 지각 M 2 태음반일주조 28.984 H m k m S 2 태양반일주조 30.000 H s k s K 1 태양 태음합성일주조 15.041 H' k' O 1 태음일주조 13.943 H o k o ζ 기타분조 7-5
7-6
Earth Moon P 3 P 2 a P O 1 m e R l m m m
Earth
소조 지구 하현 상현 대조 만월 신월 태양
(7.1) (7.2)
(7.3) (7.4)
(7.5) α αα (7.6)
ζ ζ φ φ φ (7.8) φ (7.9)
<Table 7.3> 천문학적기본주기 주기 각진동수 ( o /hr) 기호 각속도 변화율 평균태양일 1.00 ( 태양일 ) 15.0 ω 0 C s 평균태음일 1.0351 ( 태양일 ) 14.4921 ω 1 C l 항성월 27.3217 ( 태양일 ) 0.5490 ω 2 s 회귀년 365.2422 ( 태양일 ) 0.0411 ω 3 h 달의근지점 8.85 ( 년 ) 0.0046 ω 4 p 달의공전평면주기 18.61 ( 년 ) 0.0022 ω 5 N 근일점 20942 ( 년 ) - ω 6 p' ω ω ωωωω γα ω ω γ ω α ω ω α ω ω (7.10)
ω ω ω (7.11) ωωωω ωωωωωωω ωωωωωωω
반일주조성분 분조명기호 <Table 7.4> 주요한조화성분 속력 ( 평균태양시당각도 ) 주기 ( 태양시 ) 계수비 주태음반일주조 28.98410 12.42 100.0 주태양반일주조 30.00000 12.00 46.6 주태음타원반일주조 28.43973 12.66 19.2 월일합성반일주조 30.08214 11.97 12.7 주태양타원반일주조 29.95893 12.01 2.7 부태음타원반일주조 29.52848 12.19 2.8 2 차태음타원반일주조 27.89535 12.91 2.5 주태음출차반일주조 28.51258 12.63 3.6 부태음출자반일주조 29.45563 12.22 0.7 변차반일주조 27.96821 12.87 3.1 일주조성분 일월합성일주조 15.04107 23.93 58.4 주태음일주조 13.94304 25.82 41.5 주태양일주조 14.95893 24.07 19.4 주태음타원일주조 13.39866 26.87 7.9 부태음타원일주조 14.49205 24.84 3.3 소태음타원일주조 15.58544 23.10 3.3 장주기조성분 태음반월주조 1.09803 327.67 17.2 태음월주조 0.54437 661.30 9.1 태양반년주조 0.08214 2191.43 8.0 ω ωω ω ω ω ω ω ω ω ωω ωωωω ωωωω
ω σ σ ωπ ω ωω ωωωω ω ω ωω ω ωω σ
위치노량리미조리삼천포 분조 위 경 기 도 34 56' 16''N 34 42' 30''N 34 55' 27''N 도 127 52' 38''E 128 02' 57''E 128 04' 11''E 간 1980. 6. 1. 1995. 11. 1. 2001. 11. 1. ~ 1981. 5. 31. ~1997. 10. 31. ~ 2002.10. 31. 조화상수 진폭 (H) 지각 (G) 진폭 (H) 지각 (G) 진폭 (H) 지각 (G) (cm) (deg.) (cm) (deg.) (cm) (deg.) M2 106.5 259.0 84.5 249.0 94.3 250.9 S2 45.0 291.1 36.7 273.0 41.7 272.6 K1 17.2 173.6 18.6 178.1 16.6 169.0 O1 11.10 156.3 12.8 150.2 12.2 152.7 평균해면의높이 So=Hm+Hs+H'+Ho 180.4 152.6 164.8 평균고조간격 08h 56m 08h 35m 08h 39m 평균저조간격 15h 07m 14h 47m 14h 51m 고극조위 - - - 약최고고조위 (2 So) 360.8 306.4 329.6 대조평균고조위 (So+Hm+Hs) 332.5 278.8 300.8 평균고조위 (So+Hm) 286.9 240.0 259.1 소조평균고조위 (So+Hm-Hs) 241.3 201.2 217.4 평균해면 (So) 180.4 153.2 164.8 소조평균저조위 (So-(Hm-Hs)) 119.5 105.2 112.2 평균저조위 (So-Hm) 73.9 66.4 70.5 대조평균저조위 (So-(Hm+Hs)) 28.3 27.6 28.8 약최저저조위 ( 0 ) 0.0 0.0 0.0 저극조위 - - - 대조차 2(Hm+Hs) 304.2 251.2 272.0 평균조차 2 Hm 213.0 173.6 188.6 소조차 2(Hm-Hs) 121.8 96.0 105.2
조위도 삼천포 400 CM 350 329.6 APPROX. H.H.W 300.8 300 H. W. O. S. T 259.1 H. W. O. M. T 250 217.4 H. W. O. N. T 200 SP. RANGE 272.0 MN. RANGE 188.6 NP. RANGE 105.2 164.8 112.2 150 M. S. L L. W. O. N. T 100 70.5 L. W. O. M. T 50 28.8 L. W. O. S. T 0.0 0 APPROX. L.L.W -25
조위명최고고조위 (highest high water) 약최고고조위 (approx highest high water) 대조평균고조위 (high water ordinary spring tide) 평균고조위 (high water ordinary mean tide) 조화상수와의관계실측값 소조평균고조위 (high water ordinary neap tide) 평균해면 (mean sea level) 소조평균저조위 (low water ordinary neap tide) 평균저조위 (low water ordinary mean tide) 대조평균저조위 (low water ordinary spring tide) 약최저저조위 (approx lowest low water) 기본수준면 (datum level) (0.0) 최저조위 (lowest low water) 실측값 대조차 (spring range) 소조차 (neap range) 평균조차 (mean range)
<Figure 7.7> 국립해양조사원관할의조위검조소 ω ω ω ω ω ω ω ω ω ω
1 1 0.5 M2 const. M4 const. M2+M4 const. 0.5 M2 const. M4 const. M2+M4 const. H (m) 0-0.5 0 4 8 12 Time (hr) H (m) 0-0.5 0 4 8 12 Time (hr) -1-1 Ψ Ψ
1 1 M2 const. M2 const. M4 const. M4 const. 0.5 M2+M4 const. 0.5 M2+M4 const. H (m) 0-0.5 0 4 8 12 Time (hr) H (m) 0-0.5 0 4 8 12 Time (hr) -1-1 Ψ Ψ η η Ψψφ φ η ω φ η ω φ φ φ Ψ Ψ Ψ Ψ
1. 조위및수심추정 2. 연안구조물및지형변화의영향에대한연구 3. 폐수및열오염확산연구 4. 잔차류추정 5. 하구매몰에관한연구 과거의기상자료를이용한후측 (hindcast) 을통한모델의수립을거쳐실제예보체계 (real time prediction) 수립농지및용수확보를위한하구및반개방된연안지형을조력발전을위해체절하여방조제를건설할때변화되는조석및이상고조에대한연구수행해양및연안오염물질의확산에필수적입력자료인조류및해양순환자료를모델로부터추정모델에의해잔차류 ( 조석, 바람, 밀도분포에의한 ) 의분포를추정하여해수면의부유물질, 유류오염등의이동상황을추정연안에서의매몰토의정량적인조사는실제적으로중요한문제이며, 근년에 3 차원모델의결과를이용하는연구가활발히추진됨
ζ ζ ζ τ ρ τ τ ρ τ ζ ζ
ρ τ τ τ τ ζ ζ τ ζ τ ρ τ τ ζ ζ τ ρ τ τ ζ
ζ ζ ω ω (7.20) (7.21) ω ζ ω ω ζ ζ ω ζ
Wave Rider Pressure Wave Gage 20 T3 10 광양제철 공사중 광양군터미널너이테콘 CTMES 5 P1 ADP1 묘도 공사중 태인도 갈도 10 ADP Tide Gage CTMES 10 10 20 20 20 10 5 10 20 20 10 10 10 5 5 2 2 2 5 5 5 5 10 10 ADP2 20 20 20 20 남해도 170000 Gwangyang 165000 160000 2 여수반도 T1 20 WR 10 30 30 P2 20 10 T2 ============================================================================== TM LONGI TUDE LATI TUDE ------------------------------------------------------------------- NO. X Y DEG MI N SEC DEG MI N SEC ============================================================================== T1 149045. 717 270922. 659 127 46 41. 974 34 50 4. 189 T2 149088. 858 274155. 712 127 48 49. 232 34 50 4. 759 T3 160789. 148 279461. 098 127 52 22. 035 34 56 22. 993 ADP1 154673. 311 260975. 875 127 40 11. 961 34 53 9. 136 ADP2 153108. 098 270131. 962 127 46 12. 078 34 52 16. 215 P1 155131. 042 260819. 773 127 40 5. 933 34 53 24. 024 P2 149110. 670 271839. 098 127 47 18. 063 34 50 6. 066 CTMES 154701. 703 259774. 781 127 39 24. 665 34 53 10. 314 WR 149120. 396 271339. 824 127 46 58. 416 34 50 6. 508 ============================================================================== 0m 1000m 2000m 155000 150000 Container Terminal Myodo Yeosu Namhae 260000 270000 280000 Distance (m)
Velocity(m/sec) 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 Velocity(m/sec) 0.8 0.9 1 0 160000 160000 155000 155000 150000 150000 255000 260000 <Figure 7.19> 모델링에 의한 광양항 대조시의 최강창조 조류속 0.2 255000 265000 Distance (m) 0.1 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.5 0.6 0.7 0.8 0.9 1 265000 Distance (m) <Figure 7.20> 모델링에 의한 광양항 대조시의 최강낙조 조류속 Velocity(m/sec) 0.8 0.9 1 0 160000 160000 155000 155000 150000 150000 255000 0.4 260000 Velocity(m/sec) 0 0.3 260000 Distance (m) 255000 265000 <Figure 7.21> 광양항 3-2공사 이후 대조시의 최강창조 조류속 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 260000 0.9 1 265000 Distance (m) <Figure 7.22> 광양항 3-2공사 이후 대조시의 최강낙조 조류속 7.9.7 Internet을 통한 모델링의 발전방향 다양한 분야에서 접근하는 해안 해양의 발전은 각각의 분야에 정진하고 유기적인 연 관을 가져야 그 효과가 배가될 수 있다. 그렇지만 점차 다변화 다양화되는 현시점에서 공 통주제 없이는 이와 같은 유기성이 결여되고 그 효과 또한 회의적일 수 있다. 따라서 해 안 해양공학분야에서도 근자에 사회 전반에 매우 유용한 도구로서 인정 받는 인터넷상 의 정보를 공유하고 향유할 수 있다면 부족한 면을 상당부분 완화 시키면서 공통된 주제 를 향해 용이하게 접근할 수 있을 것이다. 현재까지 알려진 과학기술이라는 것이 바닷가의 수많은 모래알 중 불과 몇 알에 지나 지 않는 것에 비유되는데, 우리의 관심사인 해양과학 및 해안 해양공학분야는 이중에서 7-39
연구회나연구모임 ( 학회의분과운영 ) 을통한조석 조류자료, 정보의공유와제공 학회차원의 web server 구축운영으로학회지의목차와초록및조석 조류모델링성과제 공 인터넷을통한가상연구회 (virtual workshop) 혹은가상세미나 (virtual seminar) 의운영 benchmark test 등을통한조석수치모델기법의신뢰성향상제고 실측조류자료공유와모델접목을통한예측능력향상도모 해양기초자료를수집하고분석하는정부기관과의연계운영방안모색
(7.23) (7.24) (7.25) <Figure 7.23> 교호격자로공간이산화를표현한 모식도와시간단계및격자에서의계산위치
(7.31) (7.32)
C*** C DIMENSION U(20,20),UN(20,20),V(20,20),VN(20,20),H(20,20), > Z(20,20),IS(20),IE(20),IB(20),JB(20),NB(20) PI=3.141592 OPEN(5,FILE='TAPE5',STATUS='OLD') OPEN(6,FILE='TAPE6',STATUS='UNKNOWN') READ(5,*) DT,DX,CF,F,IM,JM,KB,NM,AMPL,PER DO J=1,JM-1 READ(5,*) IS(J),IE(J) ENDDO DO J=1,JM-1 DO I=IS(J)-1,IE(J)+1 READ(5,*) H(I,J) H(I,J)=15.0 ENDDO ENDDO DO I=IS(J)-1,IE(J)+1 DO I=1,20 DO J=1,20 U(I,J)=0.0 V(I,J)=0.0 Z(I,J)=0.0 UN(I,J)=0.0 VN(I,J)=0.0 ENDDO ENDDO DO K=1,KB READ(5,*) IB(K),JB(K),NB(K) ENDDO N=0 T=0.0 EK=0.0 120 N=N+1 T=T+DT DO J=2,JM-2 DO I=IS(J),IE(J) Z(I,J)=Z(I,J)-DT/2/DX*(U(I+1,J)*(H(I,J)+H(I+1,J))- > U(I,J)*(H(I,J)+H(I-1,J))+V(I,J+1)*(H(I,J+1)+H(I,J))- > V(I,J)*(H(I,J)+H(I,J-1))) ENDDO ENDDO DO J=2,JM-2 DO I=IS(J)+1,IE(J) VV=(V(I,J)+V(I-1,J)+V(I,J+1)+V(I-1,J+1))/4. HM=(H(I,J)+H(I-1,J))/2. UN(I,J)=U(I,J)-DT*(((U(I,J)+U(I+1,J))**2-(U(I,J)+ > U(I-1,J))**2)/8./DX+VV*(U(I,J+1)-U(I,J-1))/2./DX+9.81* > (Z(I,J)-Z(I-1,J))/DX-F*VV+CF*U(I,J)*SQRT(VV**2+U(I,J)**2) > /HM) ENDDO ENDDO DO J=3,JM-2 DO I=IS(J),IE(J) UU=(U(I,J)+U(I+1,J)+U(I,J-1)+U(I+1,J-1))/4. HM=(H(I,J)+H(I,J-1))/2. VN(I,J)=V(I,J)-DT*(((V(I,J+1)+V(I,J))**2-(V(I,J)+ > V(I,J-1))**2)/8./DX+UU*(V(I+1,J)-V(I-1,J))/2./DX+9.81* > (Z(I,J)-Z(I,J-1))/DX+F*UU+CF*V(I,J)*SQRT(UU**2+V(I,J)**2) > /HM) ENDDO ENDDO DO K=1,KB I=IB(K) J=JB(K) IF((NB(K)-1).EQ.1) GOTO 350 IF((NB(K)-1).EQ.2) GOTO 360 IF((NB(K)-1).EQ.3) GOTO 370 IF((NB(K)-1).EQ.4) GOTO 380 IF((NB(K)-1).EQ.5) GOTO 390 IF((NB(K)-1).EQ.6) GOTO 400 IF((NB(K)-1).EQ.7) GOTO 410 IF((NB(K)-1).EQ.8) GOTO 420 IF((NB(K)-1).EQ.9) GOTO 430 IF((NB(K)-1).EQ.10) GOTO 440 IF((NB(K)-1).EQ.11) GOTO 450 <Figure 7.24> 2 차원장파모델프로그램
350 UN(I,J)=0.0 GOTO 465 360 VN(I,J)=0.0 GOTO 465 370 UN(I,J)=0.0 VN(I,J)=0.0 GOTO 465 380 UN(I,J)=UN(I+1,J) VN(I,J)=VN(I+1,J) GOTO 460 390 UN(I+1,J)=UN(I,J) VN(I,J)=VN(I-1,J) GOTO 460 400 VN(I,J+1)=VN(I,J) UN(I,J)=UN(I,J-1) GOTO 460 410 UN(I,J)=UN(I,J+1) VN(I,J)=VN(I,J+1) GOTO 460 420 UN(I,J)=UN(I+1,J) VN(I,J)=0.0 GOTO 460 430 UN(I+1,J)=UN(I,J) VN(I,J)=0.0 GOTO 460 440 VN(I,J+1)=VN(I,J) UN(I,J)=0.0 GOTO 460 450 VN(I,J)=VN(I,J+1) UN(I,J)=0.0 GOTO 460 460 Z(I,J)=AMPL*SIN(2*PI*T/PER) 465 ENDDO DO J=1,JM DO I=1,IM U(I,J)=UN(I,J) V(I,J)=VN(I,J) ENDDO ENDDO IF(N/2.NE.(INT(N/2.))) GOTO 120 EK=0.0 DO J=2,JM-2 DO I=IS(J),IE(J) EK=EK+((U(I,J)+U(I+1,J))**2+(V(I,J)+V(I,J+1))**2)*H(I,J)/8. ENDDO ENDDO WRITE(6,*) EK WRITE(6,*) 'UU' DO J=JM-2,2,-1 DO I=2,IM-2 WRITE(6,*) (U(I,J)+U(I+1,J))/2 ENDDO ENDDO WRITE(6,*) 'VV' DO J=JM-2,2,-1 DO I=2,IM-2 WRITE(6,*) (V(I,J)+V(I,J+J))/2 ENDDO ENDDO WRITE(6,*) 'ZZ' DO J=JM-2,2,-1 DO I=2,IM-2 WRITE(6,*) (V(I,J)+V(I,J+J))/2 ENDDO ENDDO IF(N.LT.NM) GOTO 120 STOP END <Figure 7.24> 2 차원장파모델프로그램
240. 4000. 0.01 0.0001 9 9 11 200 1. 4800. 2 7 2 7 2 7 2 7 2 7 2 7 2 7 2 7 2 2 12 3 2 8 4 2 8 5 2 8 6 2 8 7 2 8 2 3 2 2 4 2 2 5 2 2 6 2 2 7 2 <Figure 7.25> 모델링을실시하기위해장치번호 5 를통해입력 되는주요변수및입력자료 <Figure 7.26> 모델적용을위한이상적인만의형상과특별 격자망의속성구분
<Figure 7.27> 모델적용을위한단순해역과특별격자의정보