터빈 특성을 고려한 부유식 조류발전장치의 운동성능 고찰

Dynamic Behavior of Floating Tidal Current Power Device Considering Turbine Specifications

Article information

J. Ocean Eng. Technol. 2018;32(6):427-432
Publication date (electronic) : 2018 December 27
doi : https://doi.org/10.26748/KSOE.2018.32.6.427
*Department of Naval Architecture and Ocean Engineering, Inha Univ., Incheon, Korea
조철희*, 황수진*, 박홍재*, 김명주,*orcid_icon
*인하대학교 조선해양공학과
Corresponding author Myeong-Joo Kim: +82-32-860-7342, kmj1149@hanmail.net
It is noted that this paper is revised edition based on proceedings of KSOE 2017 in Geoje.
Received 2018 January 22; Revised 2018 November 6; Accepted 2018 November 30.

Trans Abstract

Tidal current power is one of the energy sources of the ocean. Electricity can be generated by converting the flow energy of the current into the rotational energy of a turbine. Unlike tidal barrage, tidal current power does not require dams, which have a severe environmental impact. A floating-type tidal current power device can reduce the expensive support and installation cost, which usually account for approximately 41% of the total cost. It can also be deployed in relatively deep water using tensioned wires. The dynamic behavior of a floater and turbine force are coupled because the thrust and moment of the turbine affect the floater excursion, and the motion of the floater can affect the incoming speed of the flow into the turbine. To maximize the power generation and stabilize the system, the coupled motion of the floater and turbine must be extensively analyzed. However, unlike pile-fixed devices, there have been few studies involving the motion analysis of a moored-type tidal current power device.

In this study, the commercial program OrcaFlex 10.1a was used for a time domain motion analysis. In addition, in-house code was used for an iterative calculation to solve the coupled problems. As a result, it was found that the maximum mooring load of 200 kN and the floater excursion of 5.5 m were increased by the turbine effect. The load that occurred on the mooring system satisfied the safety factor of 1.67 suggested by API. The optimum mooring system for the floating tidal current power device was suggested to maximize the power generation and stability of the floater.

1. 서 론

조류발전은 조수간만 차로 발생하는 강한 유동 에너지를 터빈의 기계적 회전에너지로 변환하여 전력을 생산하는 해양 재생에너지의 한 종류이다. 정확한 발전량을 예측이 용이하며 댐과 같은 대형 구조물이 필요 없어 기존 생태 환경에 미치는 영향이 매우 적다는 장점이 있다.

기존의 조류발전장치는 대부분 파일 및 자중 고정식을 사용하였다. 하지만 전통적인 방식의 조류발전 시스템은 지지구조물과 설치비용이 전체의 약 41%를 차지해 실증 연구와 상용화에 큰 걸림돌이 되고 있다(Fig. 1). 또한 지지방식의 특성상 상대적으로 깊은 수심에서 설치 및 유지보수 비용은 매우 높게 증가함으로 현재까지 중대수심에서는 조류에너지의 활용이 제한되었다.

Fig. 1.

Costs for tidal current power system component (Vazquez and Lhlesias, 2016)

Fu and Johnstone(2017)는 수치해석을 통해 터빈 익형 특성에 따른 다자유도 운동 해석에 관한 연구를 수행하며 다양한 해상 상태에서 계류라인에 작용하는 인장력 및 터빈의 궤적을 분석하였다. Jeffcoate et al.(2016)는 직접 개발한 ProteusDS의 수치해석 프로그램을 이용하여 계류라인에 작용하는 인장력 및 플랫폼의 거동을 분석하였으며 실험을 통해 이를 검증하였다. 뿐만 아니라 계류식 조류발전장치의 실험적 운동성능 고찰에 관한 연구(Jo et al., 2017), 일점계류를 적용한 수평축 조류발전장치에 관한 연구(Jo et al., 2015a) 및 덕트 형상 변화에 따른 부유식 조류발전장치의 운동성능 분석(Jo et al., 2015b) 등 최근 설치 및 유지보수 비용 절감을 위해서 최근 다양한 계류 시스템이 적용된 조류발전장치의 운동특성에 관한 연구가 다수 진행되고 있다.

하지만 이와 같은 연구들은 터빈에 입사하는 유속의 상대속도 변화 등을 고려하지 않아, 계류삭에 작용하는 하중을 정확히 판단할 수 없다. 따라서 본 연구에서는 터빈, 부유체의 연성 효과를 고려하기 위해 부유체-터빈간의 연성 해석을 수행하였다. 조류발전은 현재 일부 선진국을 중심으로 실증연구가 이루어지고 있으며 본 연구를 통해 다루어진 부유식 조류발전장치의 경우 비교할 수 있는 관련연구가 부족하나, 풍력발전의 사례에서 부유식 고정방식을 적용함으로써 설치 및 유지보수 비용을 저감하고 발전 장치의 경제성을 확보한 사례가 있다(Kim et al., 2015). 계류식 지지구조 방식을 적용한다면 구조물의 제작 및 설치는 물론 유지보수 비용을 상당부분 절감할 수 있기에 조류발전 개발에 경제성을 높일 수 있다고 기대되고 있다.

본 연구에서는 터빈과 부유체의 연성해석을 위해 상대속도에 의해 추가적으로 발생하는 터빈의 추력과 모멘트를 도출하였으며, 유의파고, 첨두주기, 유속, 유향 등을 달리하여 총 18개의 Case 설정하여 각각의 환경조건에 대한 부유체의 운동 및 계류라인 안정성을 검토하였다. 또한 터빈의 동작 유무에 따른 부유체의 운동 및 계류라인에 인장력에 미치는 영향에 대해서 분석하였다.

2. 운동해석

2.1 운동 방정식

계류 해석에는 상용 프로그램인 OrcaFlex 10.1a가 사용되었으며, 시간영역 해석을 위한 부유체 운동 방정식은 식 (1)과 같다(Cummins, 1962).

(1) (M+a)x¨+Bx˙+Cx+th(tτ)x˙(τ)dτ=F(t)

M은 질량, a는 부가질량 계수, B는 감쇠 계수, C는 복원성 계수이며, F는 구조물에 작용하는 각종 외력들을 나타낸다. 각 유체력 계수들은 포텐셜 유동 기반의 방사/산란 문제를 통해 도출할 수 있다. 본 연구에서는 선형 포텐셜 이론 기반의 상용 프로그램인 Wadam으로부터 구조물의 유체력 계수를 도출하였다. 또한 h(t-τ)는 전달함수로서 Convolution 적분을 통해 시간 영역에서의 부유체의 거동 기억 효과를 고려할 수 있다.

2.2 해석 조건

BEMT(Blade element momentum theory)는 블레이드 국부 단면에 작용하는 힘의 관계와 운동량 이론이 서로 결합된 이론으로써, 터빈의 형상 설계 및 하중을 계산할 때 주로 사용된다. Fig. 2와 같은 블레이드 익형 주위의 유동을 보면, 날개요소의 회전 방향 속도(Vw)에 의해 익형으로 입사되는 상대유속(Vr)의 크기와 방향이 결정되며, 상대유속에 의해 수직방향으로 양력, 수평 방향으로 항력이 발생되며 두 힘은 다시 회전방향의 토크와 회전면에 수직인 추력 성분으로 나눌 수 있다.

Fig. 2.

Mechanism around blade element

터빈과 부유체 운동의 연성효과를 고려하기 위해 Matlab을 이용한 OrcaFlex 구동 알고리즘을 구성하였다(Fig. 3). V0는 조류와 파도에 의한 터빈으로 입사되는 초기 유속이며, 해양 환경력이 같은 방향으로 입사하는 보수적인 조건을 고려하였으므로 조류속(Uc)와 파입자의 수평속도(Uw)를 단순 합산하였다.

Fig. 3.

Flow chart of iterative motion analysis

BEMT코드를 통해 파랑 및 조류환경과 더불어 부유체 운동에 영향을 미치는 외력으로 산정되고, 운동 해석에 적용된다. 하지만 부유체의 운동으로 인해 터빈으로 입사하는 유속은 달라지고 발생하는 추력과 모멘트 역시 변화한다. 이를 고려하기 위해 해석이 진행된 후, 부유체 거동으로 인해 발생한 터빈에서의 상대유속(UR)을 도출하고 상응하는 터빈 하중을 계산하여 초기값을 통해 계산된 하중과의 오차를 파악한다. 오차율을 만족하지 못하면 마지막 해석에서 터빈 하중을 발생시킨 유속을 이용하여 재해석이 진행되며 본 연구에서는 속도의 시간 평균을 활용하여 2%의 오차율을 적용하였다.

시간 영역 운동해석은 80m 수심을 대상으로 이루어졌으며 유의파고 최대 6m의 ISSC(International ship and offshore structure congress) 스펙트럼 및 Power law 조류 프로파일이 사용되었다(Fig. 4). 경제적으로 조류발전을 적용할 수 있는 최소 유속은 통상 약 1.5m/s로, 국내에선 서남해안에서 제한적으로 발전 가능한 최소 유속이 발생하고 있다. 이는 높은 조수간만의 차와 병목현상을 유발하는 섬지형 등이 요구되기 때문인데, 이러한 환경의 특성을 고려하고 보다 실질적인 연구결과를 도출하기 위해 천해 조건 및 Wave spectrum을 적용하여 다양한 주기의 파도가 복합적으로 작용하는 해석 환경을 부여하였다. 불규칙한 파고와 파주기 및 다양한 유속의 조류가 복합적으로 작용하는 환경에서 터빈과 부유체의 연성 운동 경향을 파악하기 위해 유의파고, 첨두주기, 해양 외력의 작용 방향 등을 구분하여 아래와 같이 18개의 Load case를 산정하였다(Table 1).

Fig. 4.

ISSC wave spectrum

Load cases for motion analysis

1MW급으로 가정한 터빈 2기와 반잠수식 구조물을 이용한 부유식 조류시스템이 구성되었다(Fig. 5). 반잠수식 구조물은 외력에 대해 바지형 형태보다 안정적인 자세 유지성능으로 터빈으로 입사하는 유동의 방향을 유지시킬 수 있으며, Ballasting이 용이하여 수평유지는 물론 계류시스템에 장력을 제어할 수 있는 장점이 있고 이동 및 설치에 용이하다. 본 연구에서는 부유식 조류발전장치의 적용 가능성을 확인하는 차원에서 다양한 해양환경을 가정하였으며 특정 환경에 구조물을 최적화 시키는 것이 아닌, 보편적으로 사용되는 반잠수식 해양구조물을 대상으로 운동성능에 대한 고찰이 이루어졌다. Table 2를 통해 해석에 사용된 반잠수식 구조물의 제원을 나타내었다. OrcaFlex 상에서 제공하는 6D Buoy 모델을 이용하여 터빈을 구현하였으며, Matlab 코드를 이용하여 산정된 추력 및 모멘트를 작용 외력으로 사용하므로 항력과 양력 계수는 고려하지 않았다.

Fig. 5.

Designed floating type tidal current power device

Specification of floating structure

2.3 해석 결과

파도와 조류를 포함하는 유속에 부유체 운동으로 발생하는 상대 속도를 고려한 터빈 추력과 모멘트 평균값을 Table 3에 나타냈다. 터빈발생하중은 작용하는 유속의 제곱에 비례하여 증가하는 경향을 확인할 수 있었으며 최대 추력 및 모멘트는 864.4kN 및 564.1kNm로 강한 조류와 파도가 터빈과 수직으로 작용하는 Case 9에서 발생하였다.

Turbine loads analysis result

또한 Fig. 6Table 4를 통해 부유식 조류발전장치의 계류삭 명칭에 대한 정의와 Load case에 따른 계류삭에서 발생한 장력의 최댓값을 나타내었다. Fig. 7은 최대 계류삭 장력이 발생한 Case 18에서의 계류 하중과 터빈의 동적 영향성을 보여준다. 계류삭에 발생한 최대 하중은 약 1,903kN으로 터빈을 고려하지 않았을 때 보다 약 200kN 증가한 것으로 나타났다.

Fig. 6.

Mooring line layout

Mooring tension with a various load cases with turbine

Fig. 7.

Mooring tension with and without turbine

거동 반경 해석 결과는 Fig. 8Fig. 9에 나타내었다. 최대 거동반경은 11.7m로 Case 9에서 발생하였다. 터빈으로 입사하는 유속의 영향으로 인해 발생 추력과 모멘트의 크기가 증가한 영향으로 확인되며, 이 때 터빈 하중으로 인한 거동 반경 증가량은 약 5.5m로, 터빈을 고려하지 않았을 때보다 약 53%의 거동이 증가했다.

Fig. 8.

Excursion of structure with a turbine effect

Fig. 9.

Excursion of structure without a turbine effect

3. 계류 안정성 평가

설계된 조류발전장치의 운동해석을 통해 얻어진 계류 체인 발생 하중을 API RP 2SK(Design and analysis of station-keeping for floating structures)를 이용하여 검증하였다(API, 2015). 설계규정에서 제시하는 최소 계류 체인 안전계수는 Table 5와 같이 Intact condition에서 Dynamic 해석 방법을 사용했을 때 1.67로 규정된다.

Design criteria (API, 2005)

계류 체인은 Nominal diameter 64mm의 Studless R4 Grade chain으로 구성되었으며 3867kN의 MBL(Minimum breaking load)을 가진다. Table 6에 계류시스템 구성에 사용된 체인의 사양이 나타내었다, 해양 구조물 수명동안의 체인 부식을 고려하여 실제 체인 단면의 직경 90%만을 사용한 MBL의 사용을 제안하고 있으며(DNV, 2012) 본 연구에서도 MBL 감소율 19%를 가정한 Corroded MBL을 설계 기준으로 사용하였다.

Specification of mooring chain

계류 체인에서 발생한 최대하중은 1,944kN으로, 가장 큰 유의 파고와 조류속이 적용되어 하중이 In-line condition으로 작용하는Lload case19에서 발생하였다. 이는 Corroded MBL 3,534.8kN 대비 안전계수 1.82를 확보하는 값으로 최소 요구 안전계수 1.67을 상회하는 값이다. Table 7은 사용된 계류 체인의 MBL 및 해석에서 발생한 계류 체인 최대 인장력과 안전 계수를 나타낸다.

Analysis result summary

4. 결 론

본 연구에서는 조류발전 장치개발 비용 약 41%를 차지하는 유지보수 및 지지구조 비용을 절감할 수 있는 부유식 조류발전장치 운동해석에 관한 연구가 수행되었다. 약 2만 5천톤의 반잠수식 구조물에 1MW급 조류발전 2기를 장착하여 진행된 운동 해석과 계류 설계가 이루어졌으며 부유체 거동에 따른 상대 유속 증가를 고려하기 위해 작성된 Matlab 코드가 활용되었다.

해석 결과, 터빈 하중 최대값은 조류 및 파도가 터빈과 수직으로 입사하는 Case 9에서 발생했지만, 가장 강한 계류 체인 하중은 조류와 파도가 In-line condition으로 작용하는 Case 18에서 발생하였다. 터빈 최대 추력 및 모멘트는 약 864.4kN 및 약 564.1kN.m로 계류 체인 하중 1,903kN에 대해 터빈을 고려하지 않았을 때보다 약 200kN의 인장 하중을 증가시키는 것으로 확인되었다. 또한 터빈은 부유체의 거동 반경을 약 53% 증가시켰으며 최대 거동 반경은 약 11.7m로 확인되었다. 계류시스템은 64mm Studless R4 grade 계류체인 12개를 이용하여 구현하였고, Case 18에서 1.82의 최소 안전 계수를 확보했으며 이는 API (American Petroleum Institute)에서 제시하는 설계 안전 계수 1.67을 상회하는 값이다.

본 연구를 통해 터빈의 동적 효과를 고려한 부유체의 거동 특성을 파악할 수 있었으며 API 기준을 만족하는 계류 시스템이 설계되었다. 연구 결과를 통해 부유식 조류발전장치의 설치 가능성을 확인하였고 추후 이루어질 실증 연구의 기초 자료로 활용될 것으로 기대한다. 특히 부유식 고정 방식을 적용한 조류발전장치에 관한 연구는 개발비용의 상당부분을 차지하는 설치 및 유지보수 비용을 대폭 절감시켜 발전 단가를 경제적인 수준으로 저감시키고 조류발전의 상용화에 큰 도움이 될 것이라 기대한다.

Acknowledgements

이 논문은 2018년 해양수산부 재원으로 해양수산과학기술진흥원의 지원을 받아 수행된 연구임(해양에너지 융복합 인력양성).

본 연구는 산업통상자원부(MOTIE)와 한국에너지기술평가원(KETEP)의 지원을 받아 수행한 연구 과제입니다(No. 201630300 71850).

References

American Petroleum Institute(API). 2005. Design and Analysis of Stationkeeping Systems for Floating Structure RP-2SK. 3dr Editionth ed.
Det Norske Veritas(DNV). 2012. DNV OS-C101 Design of Offshore Steel Sturctures, General. Det Norske Vertias.
Cummins W. 1962. The Impulse Response Function and Ship Motions David Taylor Model Basin. Washington DC:
Kim HJ, Han SO, Choung JM. 2015;A Study on Effect of Aerodynamic Loads on Mooring Line Responses of a Floating Offshore Wind Turbine. Journal of the Society of Naval Architects of Korea 52(1):43–51.
Fu S, Johnstone C. 2017. Evaluating the Dynamic of Tension Mooring Supported Tidal Turbines. In : Proceedings of The 12th European Wave and Tidal Energy Conference..
Jeffcoate P, Fiore F, O’Farrell E, Starzmann R, Bischof S. 2016. Comparison of Simulations of Taut-Moored Platform PLAT-O Using ProteusDS with Experiments. In : Proceedings of AWTEC..
Jo CH, Lee KH, Cho SJ, Goo CH. 2015a. Performance of the Duct Implemented in Single Point Moored TCP System with Various Diffuser Angles. In : Proceedings of Korean Association of Ocean Science and Technology Societies Conference. Jeju Korea.
Jo CH, Rho YH, Kim DY, Kim MJ. 2015b. Dynamic Response of Single Point Mooring Tidal Energy Converter with Various Duct Configurations. In : Proceedings of The Korean Society for New and Renewable Energy..
Jo CH, Park HJ, Cho BK, Kim MJ. 2017;Investigation of Motion of Single Point Moored Duct-type TCP System by Both Numerical and Experimental Method. Transactions of the Korean Hydrogen and New Energy Society 28(2):212–219.
Vazquez A, Lhlesias G. 2016;Capital Cost in Tidal Stream Energy Projects – A Spatial Approach. Energy 107:215–226.

Article information Continued

Fig. 2.

Mechanism around blade element

Fig. 3.

Flow chart of iterative motion analysis

Fig. 4.

ISSC wave spectrum

Fig. 5.

Designed floating type tidal current power device

Fig. 6.

Mooring line layout

Fig. 7.

Mooring tension with and without turbine

Fig. 8.

Excursion of structure with a turbine effect

Fig. 9.

Excursion of structure without a turbine effect

Table 1.

Load cases for motion analysis

Load case No. Significant wave height [m] Peak period [s] Current speed at surface [m/s] Condition
1 2 5.7 1.5 in-between
2 2 5.7 2.0 in-between
3 2 5.7 2.5 in-between
4 4 8.1 1.5 in-between
5 4 8.1 2.0 in-between
6 4 8.1 2.5 in-between
7 6 9.9 1.5 in-between
8 6 9.9 2.0 in-between
9 6 9.9 2.5 in-between
10 2 5.7 1.5 in-line
11 2 5.7 2.0 in-line
12 2 5.7 2.5 in-line
13 4 8.1 1.5 in-line
14 4 8.1 2.0 in-line
15 4 8.1 2.5 in-line
16 6 9.9 1.5 in-line
17 6 9.9 2.0 in-line
18 6 9.9 2.5 in-line

Table 2.

Specification of floating structure

Description Value
Width [m] 55
Length [m] 80
Height [m] 40
Weight [ton] 25,270
Disp. [m3] 24.7

Table 3.

Turbine loads analysis result

Load case No. Thrust [kN] Thrust dev. Moment [kNm] Moment dev.
1 309.6 6.8 338.9 4.2
2 547.4 10.0 450.7 4.9
3 852.5 13.8 562.4 5.6
4 311.3 38.9 339.2 21.3
5 549.3 52.0 451.0 21.5
6 854.6 65.3 562.7 21.7
7 319.3 92.0 340.6 49.5
8 558.2 122.1 452.4 49.7
9 864.4 152.5 564.1 49.9
10 154.8 3.4 239.6 3.0
11 273.7 5.0 318.6 3.4
12 426.2 6.9 397.6 3.9
13 155.7 19.5 239.9 15.1
14 274.7 26.1 318.9 15.3
15 427.3 32.7 397.9 15.4
16 159.6 46.1 240.8 35.1
17 279.0 61.2 319.8 35.2
18 432.1 76.5 398.9 35.4

Table 4.

Mooring tension with a various load cases with turbine

Load case No. Max. tension [kN] Line No. Std. Dev.
1 402.3338 1-1 12.5
2 558.8096 1-1 15.9
3 817.9315 1-1 24.3
4 513.5729 1-1 26.6
5 701.7993 1-1 39.8
6 1030.969 1-1 59.7
7 666.8946 1-1 46.8
8 930.9297 1-1 67.6
9 1311.669 1-1 107.4
10 659.3693 4-1 57.0
11 792.5558 4-1 57.6
12 989.2022 4-1 59.9
13 1081.043 4-1 124.8
14 1310.183 4-1 137.3
15 1641.159 4-1 143.9
16 1357.692 4-1 136.6
17 1702.308 4-1 163.2
18 1903.318 4-1 187.7

Table 5.

Design criteria (API, 2005)

Analysis method Tension limits [% of MBL] Equivalent FOS
Intact Quasi-static 50 2.0
Dynamic 60 1.67
Damaged Quasi-static 70 1.43
Dynamic 80 1.25

Table 6.

Specification of mooring chain

Specification Value
Type 64mm R4 studless
Length [m] 425.7
Normal drag diameter [m] 0.064
Mass per unit length [kg /m] 82
Axial stiffness [kN] 349,798
Minimum breaking load [kN] 3534.84
Factor of safety 1.67

Table 7.

Analysis result summary

Condition Mooring Legs

Max tension [kN] Factor of safety
Intact / Dynamic 1944.32 1.82
Minimum breaking load [kN] 3534.84 (Corroded)