KR101902730B1 - 해기사의 심리적 특성을 이용한 선박의 충돌 회피 방법 - Google Patents

해기사의 심리적 특성을 이용한 선박의 충돌 회피 방법 Download PDF

Info

Publication number
KR101902730B1
KR101902730B1 KR1020180016715A KR20180016715A KR101902730B1 KR 101902730 B1 KR101902730 B1 KR 101902730B1 KR 1020180016715 A KR1020180016715 A KR 1020180016715A KR 20180016715 A KR20180016715 A KR 20180016715A KR 101902730 B1 KR101902730 B1 KR 101902730B1
Authority
KR
South Korea
Prior art keywords
collision
ship
reference value
distance
relative distance
Prior art date
Application number
KR1020180016715A
Other languages
English (en)
Inventor
정재룡
임정빈
박득진
Original Assignee
목포해양대학교 산학협력단
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by 목포해양대학교 산학협력단 filed Critical 목포해양대학교 산학협력단
Priority to KR1020180016715A priority Critical patent/KR101902730B1/ko
Priority to JP2018125553A priority patent/JP6679665B2/ja
Priority to CN201810694738.0A priority patent/CN110164184A/zh
Priority to US16/026,025 priority patent/US20190251849A1/en
Application granted granted Critical
Publication of KR101902730B1 publication Critical patent/KR101902730B1/ko

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/16Devices for psychotechnics; Testing reaction times ; Devices for evaluating the psychological state
    • A61B5/165Evaluating the state of mind, e.g. depression, anxiety
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B63SHIPS OR OTHER WATERBORNE VESSELS; RELATED EQUIPMENT
    • B63BSHIPS OR OTHER WATERBORNE VESSELS; EQUIPMENT FOR SHIPPING 
    • B63B43/00Improving safety of vessels, e.g. damage control, not otherwise provided for
    • B63B43/18Improving safety of vessels, e.g. damage control, not otherwise provided for preventing collision or grounding; reducing collision damage
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B63SHIPS OR OTHER WATERBORNE VESSELS; RELATED EQUIPMENT
    • B63BSHIPS OR OTHER WATERBORNE VESSELS; EQUIPMENT FOR SHIPPING 
    • B63B43/00Improving safety of vessels, e.g. damage control, not otherwise provided for
    • B63B43/18Improving safety of vessels, e.g. damage control, not otherwise provided for preventing collision or grounding; reducing collision damage
    • B63B43/20Feelers
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B63SHIPS OR OTHER WATERBORNE VESSELS; RELATED EQUIPMENT
    • B63BSHIPS OR OTHER WATERBORNE VESSELS; EQUIPMENT FOR SHIPPING 
    • B63B49/00Arrangements of nautical instruments or navigational aids
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B63SHIPS OR OTHER WATERBORNE VESSELS; RELATED EQUIPMENT
    • B63BSHIPS OR OTHER WATERBORNE VESSELS; EQUIPMENT FOR SHIPPING 
    • B63B79/00Monitoring properties or operating parameters of vessels in operation
    • B63B79/20Monitoring properties or operating parameters of vessels in operation using models or simulation, e.g. statistical models or stochastic models
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F16/00Information retrieval; Database structures therefor; File system structures therefor
    • G06F16/20Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
    • G06F16/29Geographical information databases
    • GPHYSICS
    • G08SIGNALLING
    • G08BSIGNALLING OR CALLING SYSTEMS; ORDER TELEGRAPHS; ALARM SYSTEMS
    • G08B21/00Alarms responsive to a single specified undesired or abnormal condition and not otherwise provided for
    • G08B21/18Status alarms
    • G08B21/182Level alarms, e.g. alarms responsive to variables exceeding a threshold
    • GPHYSICS
    • G08SIGNALLING
    • G08BSIGNALLING OR CALLING SYSTEMS; ORDER TELEGRAPHS; ALARM SYSTEMS
    • G08B21/00Alarms responsive to a single specified undesired or abnormal condition and not otherwise provided for
    • G08B21/18Status alarms
    • G08B21/24Reminder alarms, e.g. anti-loss alarms
    • GPHYSICS
    • G08SIGNALLING
    • G08GTRAFFIC CONTROL SYSTEMS
    • G08G3/00Traffic control systems for marine craft
    • G08G3/02Anti-collision systems
    • B63B2732/00

Landscapes

  • Engineering & Computer Science (AREA)
  • Ocean & Marine Engineering (AREA)
  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Mechanical Engineering (AREA)
  • Chemical & Material Sciences (AREA)
  • Combustion & Propulsion (AREA)
  • Remote Sensing (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Public Health (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Emergency Management (AREA)
  • Business, Economics & Management (AREA)
  • Databases & Information Systems (AREA)
  • Theoretical Computer Science (AREA)
  • Psychiatry (AREA)
  • Developmental Disabilities (AREA)
  • Pathology (AREA)
  • Medical Informatics (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Biomedical Technology (AREA)
  • Veterinary Medicine (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Biophysics (AREA)
  • Social Psychology (AREA)
  • Psychology (AREA)
  • Data Mining & Analysis (AREA)
  • General Engineering & Computer Science (AREA)
  • Hospice & Palliative Care (AREA)
  • Educational Technology (AREA)
  • Child & Adolescent Psychology (AREA)
  • Probability & Statistics with Applications (AREA)
  • Traffic Control Systems (AREA)

Abstract

본 발명은 해양사고 예방을 위해 해기사의 충돌 위기에 대한 심리적 특성을 이용하여 선박의 충돌을 회피하도록 지원해 주는 해기사의 심리적 특성을 이용한 선박의 충돌 회피 방법에 관한 것으로, 상대선박의 정보와 본선 정보를 이용하여 두 선박 간의 상대거리(RD)와 상대방위(RB)를 계산하는 제1단계와, 특정 상황을 샘플링(sampling)한 이산적(discrete) 데이터 형태인 해기사의 충돌위기(CR) 지각데이터를 이용하여 상기 상대거리(RD)와 상기 상대방위(RB)에 대한 충돌위기수준(CL)을 추정하고, 그 추정된 결과를 3차원 좌표 데이터로 변환하는 충돌위기예측모델(CREM)을 모델링하는 제2단계와, 선박영역의 거리(DSD)와 상기 모델링된 충돌위기예측모델(CREM)을 이용한 상기 충돌위기수준(CL)을 계산하는 제3단계와, 두 선박 간에 충돌위험이 있다고 판단하기 위한 기준 거리에 해당하는 공간적인 측면의 기준 값과 두 선박 간에 충돌위험이 있다고 판단하기 위한 충돌위기수준(CL)에 해당하는 심리적인 측면의 기준 값을 정하는 제4단계와, 상기 상대거리(RD)와 상기 공간적인 측면의 기준 값을 비교하여 충돌 가능 거리임을 알리는 선박영역(SD) 경보를 발생하거나, 상기 상대거리(RD)에 대한 충돌위기수준(CL(RD))과 상기 심리적인 측면의 기준 값을 비교하여 충돌위기의 수준임을 알리는 충돌위기(CR) 경보를 발생하는 제5단계로 이루어지는 것이 특징인 발명이다.

Description

해기사의 심리적 특성을 이용한 선박의 충돌 회피 방법{Ship collision avoidance method using psychological character of ship officer}
본 발명은 해양사고 예방을 위해 해기사의 충돌 위기에 대한 심리적 특성을 이용하여 선박의 충돌을 회피하도록 지원해 주는 해기사의 심리적 특성을 이용한 선박의 충돌 회피 방법에 관한 것이다.
현재 선박에는 ARPA/Radar(Automatic Radar Plotting Aids/Radar), ECDIS(Electric Chart Display and Information System), AIS(Automatic Identification System) 등 다양한 전자항해장비를 이용하여 선박 충돌을 탐지하고 있다.
선박 사이의 충돌 위험성 평가에는 선박영역(Ship Domain) 이론을 바탕으로, 선박 사이의 최근접거리(Distance at Closest Point of Approach, DCPA)와 최근접시간(Time to the Closest Point of Approach, TCPA)을 이용하고 있다.
선박영역 이론과 DCPA 및 TCPA을 이용한 선박 충돌회피 방법에는 다수의 문제가 남아 있어서 아직도 해결을 못하고 있는 실정이다.
아울러 충돌예방을 위한 다양한 항법장비가 사용되고 있음에도 불구하고 아직도 충돌사고가 빈번히 발생하고 있는데, 이러한 충돌사고가 발생하는 이유는 선박을 조종하는 해기사(OOW: Officer on the Watch)의 심리적인 특성(특히, 선박 충돌상황 시에 느끼는 충돌의 위기감)을 고려하지 않기 때문이다.
일반적으로 위험한 상황이 발생하거나 발생할 것으로 예상되는 경우 인간은 위기를 지각하게 되는데, 위기의 지각 정도는 사람마다 다르다.
선박이 서로 충돌하는 상황이 발생하면, OOW는 충돌위기(Collision Risk: CR)를 지각하는 것으로 알려져 있다.
선박이 서로 충돌 가능한 조우 상황이 발생하는 경우, 선박을 조종하는 OOW는 국제해상충돌예방규칙(Convention on the International Regulations for Preventing Collisions at Sea: COLREG)에 의거하여 충돌을 피하기 위한 일련의 동작을 취해야 한다.
한편, 세계적으로 해양사고의 70% 이상이 OOW의 인적오류에 기인한 것으로 보고되고 있는데, 국제해사기구(International Maritime Organization: IMO)에서도 인적오류의 심각성을 인지하여 인적오류 예방을 위한 다양한 활동을 전개하고 있다.
선박 충돌 상황이 발생하는 경우, OOW가 지각한 CR은 인적오류 예방에 중요하다. 그 이유는 CR을 분석하면 다양한 충돌상황에 대해서 OOW 개인 또는 특정 OOW 그룹이 지각하는 특징을 도출하여 인적오류로 인한 해양사고 예방 방안을 모색할 수 있기 때문이다
그러나 종래 기술에서 OOW의 인적오류에 관한 연구는 주로 인적오류의 원인과 분류에 관한 연구가 다수를 차지하고 있을 뿐 실제 해상에서 항해 중인 선박에 적용하기 위한 연구는 없는 실정이다. 그 이유는, 실제 해상에서 항해 중인 선박과 OOW를 대상으로 하는 실험은 대단히 위험하고 어려우며 고가의 실험 비용이 소요되기 때문이다. 이러한 이유로 OOW가 지각한 CR을 실제 선박의 충돌회피 또는 충돌예방에 적용한 연구는 전무한 실정이다.
본 발명의 목적은 상기한 점을 감안하여 안출한 것으로, 특히 선박 충돌상황 시에 OOW가 느끼는 충돌 위기감과 선박영역 이론을 이용하여 선박의 충돌을 회피할 수 있도록 지원해주는 해기사의 심리적 특성을 이용한 선박의 충돌 회피 방법을 제공하는 데 있다.
상기한 목적을 달성하기 위한 본 발명에 따른 해기사의 심리적 특성을 이용한 선박의 충돌 회피 방법의 특징은, 해기사의 충돌 위기에 대한 심리적 특성을 이용하여 선박의 충돌을 회피하도록 지원해주는 방법으로써, 상대선박의 정보와 본선 정보를 이용하여 두 선박 간의 상대거리(RD)와 상대방위(RB)를 계산하는 제1단계와, 특정 상황을 샘플링(sampling)한 이산적(discrete) 데이터 형태인 해기사의 충돌위기(CR) 지각데이터를 이용하여 상기 상대거리(RD)와 상기 상대방위(RB)에 대한 충돌위기수준(CL)을 추정하고, 그 추정된 결과를 3차원 좌표 데이터로 변환하는 충돌위기예측모델(CREM)을 모델링하는 제2단계와, 선박영역의 거리(DSD)와 상기 모델링된 충돌위기예측모델(CREM)을 이용한 상기 충돌위기수준(CL)을 계산하는 제3단계와, 두 선박 간에 충돌위험이 있다고 판단하기 위한 기준 거리에 해당하는 공간적인 측면의 기준 값과 두 선박 간에 충돌위험이 있다고 판단하기 위한 충돌위기수준(CL)에 해당하는 심리적인 측면의 기준 값을 정하는 제4단계와, 상기 상대거리(RD)와 상기 공간적인 측면의 기준 값을 비교하여 충돌 가능 거리임을 알리는 선박영역(SD) 경보를 발생하거나, 상기 상대거리(RD)에 대한 충돌위기수준(CL(RD))과 상기 심리적인 측면의 기준 값을 비교하여 충돌위기의 수준임을 알리는 충돌위기(CR) 경보를 발생하는 제5단계로 이루어지는 것이다.
바람직하게, 상기 제4단계는 상기 상대방위(RB)에 대한 선박영역의 거리(DSD)를 상기 공간적인 측면의 기준 값으로 정하고, 상기 충돌위기예측모델(CREM)에서 추정한 충돌위기수준(CL) 중에서 상기 상대방위(RB)에 대한 선박영역의 거리(DSD)에 해당하는 충돌위기수준(CL(DSD))을 상기 심리적인 측면의 기준 값으로 정할 수 있다.
바람직하게, 상기 제5단계는 상기 상대거리(RD)와 상기 공간적인 측면의 기준 값을 비교하여, 상기 상대거리(RD)가 상기 공간적인 측면의 기준 값 보다 크면 상기 선박영역(SD) 경보를 발생할 수 있다.
바람직하게, 상기 상대거리(RD)에 대한 충돌위기수준(CL(RD))과 상기 심리적인 측면의 기준 값을 비교하여, 상기 심리적인 측면의 기준 값이 상기 상대거리(RD)에 대한 충돌위기수준(CL(RD)) 보다 크면 상기 충돌위기(CR) 경보를 발생할 수 있다.
바람직하게, 상기 제2단계에서 상기 충돌위기예측모델(CREM)의 입력 변수는 상기 상대방위(RB)와 상기 상대거리(RD)를 포함하고, 상기 충돌위기예측모델(CREM)의 출력 변수는 연속된 상대방위(RB)와 상대거리(RD)에 대해서 추정한 충돌위기(CR)의 상기 충돌위기수준(CL)을 포함할 수 있다.
보다 바람직하게, 상기 출력 변수는 상기 충돌위기수준(CL)을 3차원 복합지도에 표시하기 위한 좌표 값을 더 포함할 수 있다.
보다 바람직하게, 상기 제2단계에서 상기 충돌위기예측모델(CREM)은 극대치(GEV) 분포의 확률밀도함수(pdf)를 이용하여 상기 입력 변수에 대한 상기 충돌위기수준(CL)을 추정할 수 있다.
본 발명에 따르면, 선박영역 이론과 OOW의 심리적인 충돌위기수준(Collision Level: CL)을 복합적으로 이용하기 때문에, 기존에 선박영역 이론만을 적용하는 방법과 비교하여 충돌회피 동작 실행에 대한 불이행의 확률을 감소시킬 수 있어서 충돌회피 또는 충돌예방 가능성을 증대시킬 수 있다.
공간적 및 심리적인 선박 충돌 영역을 이용한 복수의 경보 기능을 통해서 OOW의 충돌회피 동작의 실행에 대한 신뢰성을 향상 시킬 수 있다.
선박 사이의 상대방위와 상대거리에 따라서 충돌 경보가 필요 충분하게 그리고 서로 상호적으로 발생하기 때문에 OOW의 주의력 집중과 충돌회피 동작이 서로 복합된 고도의 충돌예방이 가능하다.
OOW의 심리적인 충돌 위기감을 충돌회피 방식에 적용하기 때문에 인적오류를 적극 예방할 수 있다.
충돌회피 동작에 필요한 최소 거리를 알 수 있고, 평균적인 OOW들이 인지하는 상대방위와 상대거리에 대한 충돌의 위험성을 알 수 있다.
물리적인 공간 확보와 함께 심리적인 공간확보가 가능하여 차세대 항법 시스템 개발에 적용 가능하다.
육상의 해상교통관제센터(Vessel Traffic System, VTS)의 관제사에게 적용하여 이들 관제사들의 주의력 결핍에 의한 사고예방에 기여할 수 있다. 특히, 다수의 선박을 동시에 관리해야 하는 VTS에 적용하는 경우에는 자동으로 선박 사이의 충돌위기를 관제사에게 통보할 수 있어서 적극적인 관제 지원이 가능하다.
도 1은 본 발명의 일 실시 예에 따른 해기사의 심리적 특성을 이용한 선박의 충돌 회피 방법에서 CREM의 입출력 구조를 도시한 다이어그램이고,
도 2는 본 발명의 일 실시 예에 따른 해기사의 심리적 특성을 이용한 선박의 충돌 회피 방법에서 CREM에 사용된 GEV 분포의 pdf 예를 도시한 그래프이고,
도 3은
Figure 112018014873253-pat00001
로 나타낸 SD 영역을 도시한 그래프이고,
도 4는 샘플링된 다섯 가지 선박충돌 조우 상황(S1 ~ S5)에 대한 ES(
Figure 112018014873253-pat00002
데이터를 도시한 그래프이고,
도 5는 다섯 가지 선박충돌 조우 상황 (S1 ~ S5)에 대해서 측정한 PCR 데이터를 도시한 그래프이고,
도 6은 연속된 상대방위
Figure 112018014873253-pat00003
에 대해서 추정한 결과를 도시한 그래프이고,
도 7은 ES 데이터와 PCR 데이터를 CREM에 적용하여 가시화한 CR 영역을 도시한 그래프이고,
도 8은 SD 영역과 추정된 CR 영역을 비교하기 위하여 SD 영역과 CR 영역을 겹쳐서 도시한 그래프이고,
도 9는 DSD의 상대거리와 CL(DSD)의 계산결과를 도시한 그래프이고,
도 10은 본 발명에 따른 해기사의 심리적 특성을 이용한 선박의 충돌 회피 방법에서 DSD와 CL(DSD)을 이용하여 충돌위험을 알리기 위한 경보(Warning)를 설정하는 방법을 설명하기 위한 그래프이고,
도 11은 본 발명의 일 실시 예에 따른 해기사의 심리적 특성을 이용한 선박의 충돌 회피 절차를 도시한 다이어그램이다.
본 발명의 다른 목적, 특징 및 이점들은 첨부한 도면을 참조한 실시 예들의 상세한 설명을 통해 명백해질 것이다.
이하, 첨부된 도면을 참조하여 본 발명의 실시 예의 구성과 그 작용을 설명하며, 도면에 도시되고 또 이것에 의해서 설명되는 본 발명의 구성과 작용은 적어도 하나의 실시 예로서 설명되는 것이며, 이것에 의해서 상기한 본 발명의 기술적 사상과 그 핵심 구성 및 작용이 제한되지는 않는다.
이하, 첨부한 도면을 참조하여 본 발명에 따른 해기사의 심리적 특성을 이용한 선박의 충돌 회피 방법의 바람직한 실시 예를 자세히 설명한다.
본 발명은 선박이 서로 충돌할 수 있는 조우 상황에서 당직 OOW가 느끼는 충돌위기(Collision Risk: CR)를 이용한다.
CR을 추정하기 위해 충돌위기예측모델(Collision Risk Estimation Model: CREM)을 사용하며, CREM은 일반화된 극대치(Generalized Extremely Value: GEV) 분포의 확률밀도함수(probability density function: pdf)의 형상 자체를 이용하여 구성한다.
이어, CREM을 이용하여 선박 사이의 상대방위와 상대거리에 대한 충돌위기수준(Collision Level: CL)을 추정한다.
이어, 2차원 직교좌표와 극-좌표 및 3차원 콘토맵(contour map)을 혼합한 3차원 복합지도에 추정된 CL에 대한 충돌위기 영역(CR domain)을 표시한다.
충돌위기 영역(CR domain)은 선박영역(Ship Domain: SD)과 비교하고, 그 두 영역 사이의 차이를 산출한다. 그 산출된 차이 값을 충돌위기 판단을 위한 기준 값으로 설정한다. 기준 값의 설정에 대해서는 다음에 상세한다.
본 발명에서 충돌의 가능성은 두 척의 선박이 서로 접근하는 경우 OOW가 느끼게 될 CL과 SD 거리를 이용하여 판단한다.
여기서, 선박이 서로 충돌할 수 있는 조우 상황에 근접할 때 발생하는 OOW의 CR을 CL로 추정하고, 3차원 복합지도에 표시된 SD와 CL을 통해 충돌의 가능성을 공간적으로 나타낸다.
이와 같이 본 발명은 물리적인 요소(거리, 속력, 방위 등) 뿐만 아니라 사람의 인지적인 지각요소를 공간화한 인자까지 충돌을 회피하는데 사용한다.
선박의 충돌회피를 위한 개념 중 하나인 선박영역(Ship Domain: SD) 이론은 선박이 충돌을 회피하기 위하여 필요 충분한 공간적인 영역을 확보하기 위한 물리적인 개념이다.
SD 이론을 위해, Radar, ARPA(Automatic Radar Plotting Aids), ECDIS(Electric Chart Display and Information System) 등의 충돌회피를 지원하는 전자항해장비를 적용한다.
본 발명은 선박 사이의 최근접거리(Distance at Closest Point of Approach: DCPA)와 최근접시간(Time to the Closest Point of Approach: TCPA) 등의 공간적인 영역을 산출하기 위한 SD 이론에 OOW가 지각한 CR을 통해 산출되는 충돌에 대한 상황인식(Situation Awareness: SA)을 더 적용한다. 즉, 본 발명에서는 SD 이론과 CR과 SA의 관계를 통해 선박의 충돌을 회피할 수 있도록 지원해 준다.
OOW가 지각하는 CR을 예측하는 것이 중요하다. CR을 예측하기 위해 CREM을 사용한다.
OOW가 지각한 CR은 설문지를 이용하여 조사하거나 OOW에게 특정 장치를 부착하여 위험한 상황에 대한 심장 박동수와 혈압 등을 측정하여 수집한다.
수집된 CR 데이터는 특정 상황을 샘플링(sampling)한 이산적(discrete)인 데이터 형태이므로, 연속된 입력에 대해서 원하는 데이터 형태로 출력할 수 있는 CREM을 이용한다.
본 발명에서는 연속된 입력에 대해 이산적인 CR 데이터를 이용하여 연속된 충돌위기수준(Collision Level: CL)을 추정할 수 있고, 그 추정한 결과를 3차원 좌표 데이터로 변환하기 위한 CREM을 모델링한다.
도 1은 본 발명의 일 실시 예에 따른 해기사의 심리적 특성을 이용한 선박의 충돌 회피 방법에서 CREM의 입출력 구조를 도시한 다이어그램이다.
도 1을 참조하면, CREM의 입력 변수는 선박 사이에 발생하는 상대방위(Relative Bearing: RB)와 상대거리(Relative Distance: RD)이고, CREM의 출력 변수는 연속된 상대방위와 상대거리에 대해서 추정한 CR의 충돌위기수준(Collision risk Level: CL)과 CL의 가시화를 위한 X 좌표 값과 Y 좌표 값이다.
특히, CREM은 미리 추정된 GEV(Generalized Extreme Value) 분포의 파라미터를 이용하여 입력 변수에 대한 CL을 추정한다. 이때, 추정된 결과는 상대방위와 상대거리에 대해서 일정한 양의 오차(또는 잡음)을 가질 수 있다.
다음은 CREM을 모델링하는 절차를 보다 구체적으로 설명한다.
본 발명에서는 샘플링한 CR 데이터의 특징 때문에 GEV 분포의 확률밀도함수(probability density function: pdf)를 이용한다.
GEV 분포의 pdf는 식(1)과 같이 주어진 표준정규(standard normal) 데이터
Figure 112018014873253-pat00004
에 대해서 형상 파라미터(shape parameter)
Figure 112018014873253-pat00005
, 위치 파라미터(location parameter)
Figure 112018014873253-pat00006
, 스케일 파라미터(scale parameter)
Figure 112018014873253-pat00007
등으로 정의할 수 있고, 식(1)의 GEV 분포의 pdf는 식(2)과 같이
Figure 112018014873253-pat00008
로 간단하게 나타낼 수 있다.
Figure 112018014873253-pat00009
Figure 112018014873253-pat00010
(1)
Figure 112018014873253-pat00011
(2)
도 2는 식(2)을 이용하여 계산한 GEV 함수의 pdf의 예를 도시한 것으로,
Figure 112018014873253-pat00012
,
Figure 112018014873253-pat00013
,
Figure 112018014873253-pat00014
,
Figure 112018014873253-pat00015
등을 식(2)에 적용한 경우의
Figure 112018014873253-pat00016
이다.
도 2에서 사각형 박스는 GEV 분포의 pdf(
Figure 112018014873253-pat00017
)의 좌측한계부터 최대 확률을 나타내는 우측한계까지만을 도시한 것으로, GEV 분포의 pdf를 CREM에 적용하는 방법을 설명하기 위한 것이다. 도 2에서 x축은 표준정규 데이터
Figure 112018014873253-pat00018
이고, y축은 GEV 분포의 pdf의 확률밀도 값이다.
도 2의 사각형 박스에 표시된 A, B, C의 의미는 다음과 같다. A는
Figure 112018014873253-pat00019
부터
Figure 112018014873253-pat00020
가 최대가 되는
Figure 112018014873253-pat00021
까지의 표준정규 데이터이고, B는
Figure 112018014873253-pat00022
의 기울기이며, C는
Figure 112018014873253-pat00023
에 대한
Figure 112018014873253-pat00024
값이다. 본 발명에서는 GEV 분포의 pdf에서 전술된 세 가지 특징(A, B, C)을 이용하여 CREM을 모델링한다.
CREM의 모델링은 다음과 같은 4단계 절차를 거친다.
우선,
Figure 112018014873253-pat00025
의 상대방위(RB) 각각에 대한 j
Figure 112018014873253-pat00026
의 상대거리(RD)에서 측정한 충돌위기를
Figure 112018014873253-pat00027
로 나타내면,
Figure 112018014873253-pat00028
는 I-by-J의 차원을 갖는 매트릭스 형태가 된다. CREM의 수식 전개를 간단히 하기 위하여
Figure 112018014873253-pat00029
에서 임의의 i에 대한 CR을
Figure 112018014873253-pat00030
로 정의한다.
- 1단계(곡선근사)
수열의 길이 j
Figure 112018014873253-pat00031
을 갖는 상대거리(RD)
Figure 112018014873253-pat00032
에서 측정한 샘플 데이터
Figure 112018014873253-pat00033
에 최적 근사한 b차의 다항계수
Figure 112018014873253-pat00034
는 식(3)과 같고, 식(4)와 같이
Figure 112018014873253-pat00035
이용하여 수열의 길이 n
Figure 112018014873253-pat00036
를 갖는
Figure 112018014873253-pat00037
을 추정할 수 있다.
Figure 112018014873253-pat00038
(3)
Figure 112018014873253-pat00039
(4)
여기서,
Figure 112018014873253-pat00040
Figure 112018014873253-pat00041
에 대한 상대거리(RD)를 나타낸 것으로,
Figure 112018014873253-pat00042
이고,
Figure 112018014873253-pat00043
으로 계산한 것이다. rA는 측정시작 지점에서의 상대거리이고, rZ는 측정종료 지점(또는 최대 CR 값이 나타난 지점)에서의 상대거리이다.
- 2단계(GEV 분포에 대한 pdf의 파라미터 추정)
Figure 112018014873253-pat00044
에 최적인 GEV 분포에 대한 pdf의 파라미터 세트
Figure 112018014873253-pat00045
을 표 1의 탐색 범위에서 추정한다. 표 1은 수열의 길이 w
Figure 112018014873253-pat00046
를 갖는 표준정규 데이터
Figure 112018014873253-pat00047
Figure 112018014873253-pat00048
의 탐색 범위(좌측 한계 값은 아래 첨자 LT로 표시하고, 우측 한계 값은 아래 첨자 RT로 표시)를 나타낸다.
표 1에 나타낸
Figure 112018014873253-pat00049
Figure 112018014873253-pat00050
을 식(2)에 적용하면 수열의 길이 w를 갖는 GEV 분포에 대한 pdf의
Figure 112018014873253-pat00051
를 획득할 수 있다. 여기서, 모델에 적용할 GEV 분포에 대한 pdf는 도 2에서 설명한 바와 같이 N개의 수열의 길이를 갖는
Figure 112018014873253-pat00052
에 최적 근사할 수 있는
Figure 112018014873253-pat00053
의 일부분을 추정하는 것이기 때문에 W
Figure 112018014873253-pat00054
의 조건이 필요하다. 표 1의 탐색범위는 이러한 조건들을 사전에 고려하여 결정한 것으로,
Figure 112018014873253-pat00055
으로 고정시켜서
Figure 112018014873253-pat00056
을 탐색한다.
Standard normal data
Figure 112018014873253-pat00057
Shape parameter
Figure 112018014873253-pat00058
Scale parameter
Figure 112018014873253-pat00059
Position parameter
Figure 112018014873253-pat00060
Left limits
Right limits
Figure 112018014873253-pat00061

Figure 112018014873253-pat00062
Figure 112018014873253-pat00063

Figure 112018014873253-pat00064
Figure 112018014873253-pat00065

Figure 112018014873253-pat00066
Figure 112018014873253-pat00067
Spaces
Figure 112018014873253-pat00068
Figure 112018014873253-pat00069
Figure 112018014873253-pat00070
-
Index w = 1, 2, 3,
Figure 112018014873253-pat00071
, W
q = 1, 2, 3,
Figure 112018014873253-pat00072
, Q
v = 1, 2, 3,
Figure 112018014873253-pat00073
, V
-
먼저, 수열의 길이를 w = W로 고정한 정규표준 데이터
Figure 112018014873253-pat00074
에 대해서
Figure 112018014873253-pat00075
을 표 1의 범위에서 변화시켰을 때의 GEV 분포에 대한 pdf,
Figure 112018014873253-pat00076
을 식(5)으로 계산한다.
Figure 112018014873253-pat00077
(5)
그리고
Figure 112018014873253-pat00078
의 최대 값
Figure 112018014873253-pat00079
을 구하고,
Figure 112018014873253-pat00080
의 좌측한계
Figure 112018014873253-pat00081
부터
Figure 112018014873253-pat00082
까지의 수열의 길이
Figure 112018014873253-pat00083
(여기서,
Figure 112018014873253-pat00084
는 표 1에 나타낸
Figure 112018014873253-pat00085
의 수열의 간격)을 구한 후,
Figure 112018014873253-pat00086
에 해당하는
Figure 112018014873253-pat00087
을 식(6)으로 발췌한다.
Figure 112018014873253-pat00088
(6)
이어서 식(4)의
Figure 112018014873253-pat00089
에서
Figure 112018014873253-pat00090
에 해당하는
Figure 112018014873253-pat00091
(식(7))을 구하고,
Figure 112018014873253-pat00092
의 최대 값
Figure 112018014873253-pat00093
을 구한 후,
Figure 112018014873253-pat00094
과 동일한 최대 값을 갖는
Figure 112018014873253-pat00095
을 식(8)로 구한다.
Figure 112018014873253-pat00096
(7)
Figure 112018014873253-pat00097
(8)
그리고 식(4)의
Figure 112018014873253-pat00098
에서
Figure 112018014873253-pat00099
일 때의
Figure 112018014873253-pat00100
의 길이에 해당하는 최소상대거리
Figure 112018014873253-pat00101
을 식(9)로 구한다.
Figure 112018014873253-pat00102
(9)
식(6)의
Figure 112018014873253-pat00103
와 식(8)의
Figure 112018014873253-pat00104
사이의 평균오차
Figure 112018014873253-pat00105
을 식(10)으로 구하고,
Figure 112018014873253-pat00106
가 최소가 될 때의
Figure 112018014873253-pat00107
을 식(11)로 구한다.
Figure 112018014873253-pat00108
(10)
Figure 112018014873253-pat00109
(11)
즉, 식(3)부터 식(11)를 사용하여 상대거리
Figure 112018014873253-pat00110
에서 측정한 샘플 데이터
Figure 112018014873253-pat00111
에 대한 GEV 분포 pdf의 파라미터를 추정한다.
다음에는 충돌위기
Figure 112018014873253-pat00112
전체에 대해서 전술된 파라미터 추정 과정을 반복하여 i에 대한
Figure 112018014873253-pat00113
을 식(11)로 추정하고,
Figure 112018014873253-pat00114
을 이용하여 GEV 형상 파라미터
Figure 112018014873253-pat00115
와 스케일 파라미터
Figure 112018014873253-pat00116
그리고 수열의 길이
Figure 112018014873253-pat00117
을 구한다. 그리고 이 결과를 이용하여 식(12)부터 식(16)까지를 계산한다.
Figure 112018014873253-pat00118
(12)
Figure 112018014873253-pat00119
(13)
Figure 112018014873253-pat00120
(14)
Figure 112018014873253-pat00121
(15)
Figure 112018014873253-pat00122
(16)
- 3단계(추정한 파라미터 보간)
전술된 2단계에서 추정한 결과는 이산적인 상대방위
Figure 112018014873253-pat00123
에 대한 것이기 때문에 연속된 상대방위
Figure 112018014873253-pat00124
에 대한 모델 파라미터가 필요하다.
본 발명에서는 다음과 같은 보간을 이용하여 연속된 파라미터를 추정한다.
보간은 Matlab 코드 'Outputs
Figure 112018014873253-pat00125
를 이용할 수 있다.
'interp1'은 1차원 보간을 의미하고, Var1과 Var2 그리고 Var3는 입력 변수들이다. 그리고 'pchip'은 Piecewise Cubic Hermit (PCH) 보간을 의미하는 것으로, 본 발명에서는 주어진 데이터의 특성을 최대한 유지하기 위하여 PHC 보간을 적용한다.
Matlab codes Meanings
Figure 112018014873253-pat00126
Interpolated GEV shape parameter
Figure 112018014873253-pat00127
Interpolated GEV scale parameter
Figure 112018014873253-pat00128
Interpolated maximum DP
Figure 112018014873253-pat00129
Interpolated minimum RD
표 2는 연속된 모델 파라미터를 획득하기 위한 PHC 보간에 적용된 Matlab 코드를 나타낸 것이다.
- 4단계(연속된 방위에 대한 계산)
표 2의 보간 결과를 이용하여 상대방위
Figure 112018014873253-pat00130
에 대해서 추정하면 다음과 같다.
먼저, w개의 수열의 길이를 갖는 정규표준 데이터
Figure 112018014873253-pat00131
에 대한 GEV 분포의 pdf
Figure 112018014873253-pat00132
(
Figure 112018014873253-pat00133
을 식(17)로 구하고,
Figure 112018014873253-pat00134
의 최대 값
Figure 112018014873253-pat00135
을 식(18)로 구한다. 그리고
Figure 112018014873253-pat00136
의 좌측한계
Figure 112018014873253-pat00137
부터
Figure 112018014873253-pat00138
인 지점의
Figure 112018014873253-pat00139
까지의 수열의 길이
Figure 112018014873253-pat00140
를 식(19)로 구한다.
Figure 112018014873253-pat00141
(17)
Figure 112018014873253-pat00142
(18)
Figure 112018014873253-pat00143
(19)
식(17)의
Figure 112018014873253-pat00144
에서 w
Figure 112018014873253-pat00145
에 해당하는
Figure 112018014873253-pat00146
을 구하고, 다음 식(20)을 이용하여 최대 값이 1.0인
Figure 112018014873253-pat00147
을 구한다.
Figure 112018014873253-pat00148
(20)
여기서,
Figure 112018014873253-pat00149
Figure 112018014873253-pat00150
의 최대 값이고,
Figure 112018014873253-pat00151
은 원래 데이터에서 충돌위기수준의 최대 값이다.
식(20)의
Figure 112018014873253-pat00152
이 CR을 가시화하는데 적용할 충돌위기수준(CL)이 된다. 마지막으로 식(21)로 계산한 좌표 값과
Figure 112018014873253-pat00153
을 이용하여 CR 영역을 형성하기 위한 3차원 좌표의 변수 세트
Figure 112018014873253-pat00154
을 산출한다.
Figure 112018014873253-pat00155
(21)
여기서
Figure 112018014873253-pat00156
(radian)으로 x-y 직교좌표계에서 방위 90도(즉,
Figure 112018014873253-pat00157
)을 기준 방위 0도로 정한 것이고,
Figure 112018014873253-pat00158
이다.
전술된 CREM 모델링에 적용된 SD(Ship Domain) 데이터와, ES(Emotional Sensitivity) 데이터와 PCR(Perceived CR) 데이터에 대해 설명한다.
SD 데이터의 선박의 조우 상황에 따라서 충돌 회피에 필요한 여유공간을 정의하기 위한 것으로 이론적인 개념을 컴퓨터로 계산한 거리 데이터이다. ES 데이터는 제한된 수역에서 선박 OOW들의 선박 조종의 어려움을 측정하기 위한 것으로 선박조종시뮬레이터를 이용하여 선박의 충돌 조우 상황에 대하여 측정한 데이터이다. PCR 데이터는 OOW들이 충돌 상황에 대하여 지각하는 위기감을 측정하기 위한 것이다.
SD 데이터는 SD 스케일을 측정한 후 연속된 방위에 대한 거리로 변환시킨 것으로, 반경 Cr을 갖는 원의 중심에 가상 선박(Phantom ship)을 두고, 실제 선박(Real ship)을 원의 중심에서 x-축과 y-축으로 일정한 거리만큼 편이(offset)시킨 형태이다.
SD 데이터를 계산하는 절차는 다음과 같다.
먼저, 극-좌표를 직교좌표로 변환시키는 매트랩 코드 '
Figure 112018014873253-pat00159
= pol2cart(
Figure 112018014873253-pat00160
,
Figure 112018014873253-pat00161
)'을 이용하여 n개의 직교좌표 데이터 세트
Figure 112018014873253-pat00162
을 계산한다. 여기서,
Figure 112018014873253-pat00163
은 방위
Figure 112018014873253-pat00164
에 대해서 식(22)을 이용하여 360도 방위표시법으로 계산한 라디안(radian) 단위의 방위이다.
Figure 112018014873253-pat00165
(22)
여기서
Figure 112018014873253-pat00166
은 0도부터 360도까지를 N개로 구분한 방위이고,
Figure 112018014873253-pat00167
이다.
SD 데이터 세트
Figure 112018014873253-pat00168
는 다음과 같이 계산할 수 있다.
Figure 112018014873253-pat00169
,
Figure 112018014873253-pat00170
(23)
여기서,
Figure 112018014873253-pat00171
은 x-축의 편이 값(287.06 m)이고,
Figure 112018014873253-pat00172
는 y-축의 편이 값(864.27 m)이다.
도 3은
Figure 112018014873253-pat00173
로 나타낸 SD 영역을 도시한 것이다. x-축과 y-축은 1852m의 단위로 표시한 거리이고, 좌표 중심(O)으로부터 최대 거리(maxR)가 되는 18.4도 방향으로 편이가 된 원까지의 공간이 SD 영역이 된다. 도 3에서 본선의 선수방위를 0도 정하여 시계방향으로 360도까지 나타내는 360도 방위표시법을 적용하고, 국제해상마일(international nautical mile, NM)의 1.0 NM에 해당하는 1,852 meter을 거리의 단위를 적용한 것이다. 본 발명에서는 0도부터 180도까지 굵은 선으로 나타낸 SD의 우측 반원을 CR 영역과의 비교 분석에 적용한다.
ES(Emotional Sensitivity) 데이터는 선박이 충돌할 수 있는 조우 상황에 대해서 OOW가 지각하는 위험의 정도를 선박조종 시뮬레이터를 이용하여 측정한 데이터이다.
도 4는 샘플링된 다섯 가지 선박충돌 조우 상황(S1 ~ S5)에 대한 ES(
Figure 112018014873253-pat00174
데이터를 도시한 것이다.
도 4에서 S1부터 S5까지는 두 척의 선박이 서로 상대방위 0도, 45도, 90도, 135도, 180도 등으로 조우하여 충돌하는 상황을 각각 나타낸다. 도 4의 데이터는 일정한 간격으로 구분된 상대거리에 대응하는 ES 값을 산출하되 ES 값의 변화가 큰 상대거리는 간격을 줄여서 샘플링한 것이다.
도 4에서 상대거리 0m부터
Figure 112018014873253-pat00175
m까지는 다른 조우 상황과 비교하여 S4에 대한 ES 값이 크다가 그 이후부터는 S5의 ES 값이 크게 된다. 그리고 최대 ES 값에 도달하는 상대거리가 다섯 가지 조우 상황마다 각각 다르게 나타난다. 도 4에 도시된 바와 같이 ES 데이터는 세 가지 특징(A, B, C)으로 대표할 수 있다. 즉, A는 최대 ES 값이 나타나는 최소 상대거리이고, B는 상대거리에 따른 ES 값의 증가와 감소의 변화 양상이고, C는 상대거리에 대한 ES 값이다.
PCR 데이터는 실제 함정에서 측정한 CRPI를 참고하여 획득한 데이터이다. 도 5는 다섯 가지 선박충돌 조우 상황 (S1 ~ S5)에 대해서 측정한 PCR 데이터를 도시한 것이다. 여기서, S1부터 S5까지는 두 척의 선박이 서로 상대방위 0도, 45도, 90도, 135도, 180도 등으로 조우하여 충돌하는 상황을 각각 나타낸다. 도 4의 ES 데이터와 마찬가지로 PCR 데이터 역시 세 가지 특징(A, B, C)으로 대표할 수 있다.
CREM의 모델 변수에 대한 추정 결과에 대해 설명하면, 표 3은 ES 데이터와 PCR 데이터를 CREM에 입력하여 추정한
Figure 112018014873253-pat00176
Figure 112018014873253-pat00177
을 나타낸 것으로 식(12)부터 식(16)까지를 이용하여 추정한 이산적인 상대방위
Figure 112018014873253-pat00178
에 대한 결과이다.
Data Parameters S1 (i=1) S2 (i=2) S3 (i=3) S4 (i=4) S5 (i=5)
ES
Figure 112018014873253-pat00179
-1.0 -1.0 -1.0 -0.9 -0.8
Figure 112018014873253-pat00180
2.6 3.1 5.0 12.7 4.3
PCR
Figure 112018014873253-pat00181
-0.2 -0.6 -0.6 -0.1 -0.3
Figure 112018014873253-pat00182
9.7 9.3 9.2 12.6 0.3
도 6은 연속된 상대방위
Figure 112018014873253-pat00183
에 대해서 추정한 결과를 도시한 것으로,
Figure 112018014873253-pat00184
,
Figure 112018014873253-pat00185
,
Figure 112018014873253-pat00186
의 계산결과이다. 도 6(A)는 ES 데이터에 대한 보간 결과이고, 도 6(B)는 PCR 데이터에 대한 보간 결과이다.
도 6(A)와 6(B)의 상단 박스를 보면, 모두 상대방위 135도 부근에서
Figure 112018014873253-pat00188
의 변화가 크게 나타난다. 도 6(A)의 하단 박스에서 상대방위 0도부터 상대방위 180도까지
Figure 112018014873253-pat00189
가 증가함을 나타낸다. 즉, 상대방위 증가에 따라서 최대 ES 값 1,000이 나타나는 상대거리가 점차 감소됨을 의미한다. 도 6(B)의 하단 박스에서, 최대 PCR 값이 나타나는
Figure 112018014873253-pat00190
가 상대방위 45도 부근에서 약간 증가하고 전반적으로 일정한 값을 나타낸다.
한편, 모델링한 CREM에서 추정한 CR 영역에 대한 분석 결과에 대해 설명하면, 도 7은 ES 데이터와 PCR 데이터를 CREM에 적용하여 가시화한 CR 영역을 도시한 것이다. 등고선은 0도부터 360도까지의 상대방위에 대해서 추정한
Figure 112018014873253-pat00191
(
Figure 112018014873253-pat00192
)을 나타내고,
Figure 112018014873253-pat00193
은 가장 낮은 충돌위기수준을 의미하며
Figure 112018014873253-pat00194
은 가장 높은 충돌위기수준을 나타낸다. 도 7은 좌표의 중심을 본선과 상대선박이 가장 먼 RD =
Figure 112018014873253-pat00195
)m로 나타내서 좌표의 중심으로부터 멀어질수록 상대거리가 감소되도록 나타낸 것이다
도 7에서 좌측 그래프는 ES 데이터를 CREM에 적용하여 나타낸 영역이고, 우측 그래프는 PCR 데이터를 CREM에 적용하여 나타낸 영역이다.
도 7의 좌측 그래프에 도시된 ES 데이터를 보면, 상대방위 135도 부근에서 동일한 크기의 충돌위기수준이 다른 상대방위와 비교하여 더 긴 상대거리에서 나타난다. 예를 들어, 상대방위 0도에서
Figure 112018014873253-pat00196
이 약
Figure 112018014873253-pat00197
) m에서 나타나는데, 135도 부근에서는
Figure 112018014873253-pat00198
이 약 3(
Figure 112018014873253-pat00199
) m에서 나타난다. 즉, 본선과 상대선박 사이의 상대거리가 먼데도 불구하고 135도 부근의 선박 조우상황을 다른 선박 조우상황과 비교하여 더 일찍 OOW들이 느끼고 있음을 나타낸다. 그리고 135도 부근에서는
Figure 112018014873253-pat00200
부터
Figure 112018014873253-pat00201
까지 충돌위기수준의 증가율이 대략 일정한 비율로 발생하는 상태를 나타낸다. 이러한 상태는 도 7의 우측 그래프에 도시된 PCR 데이터에서도 동일하게 나타난다. 예를 들어, 상대방위 90도 부근의 선박조우상황은
Figure 112018014873253-pat00202
가 약
Figure 112018014873253-pat00203
m에서 나타나고 있는데, 135도 부근의 경우는
Figure 112018014873253-pat00204
가 약
Figure 112018014873253-pat00205
) m에서 나타난다.
따라서, ES 데이터와 PCR 데이터 모두 상대방위 135도 부근에서 발생하는 선박 조우상황이 다른 상대방위와 비교하여 동일한 충돌위기수준이 더 빠르게 나타난다.
도 7의 좌측 그래프에 도시된 ES 데이터를 보면, 상대방위 135도 부근을 제외하고 어느 정도 상대거리가 감소된 이후부터 급격하게 충돌위기수준이 증가하고 있다. 예를 들어, 상대방위 0도부터 90도 부근 사이에서는 상대거리 약 1.0(
Figure 112018014873253-pat00206
)에서
Figure 112018014873253-pat00207
이 나타나고,
Figure 112018014873253-pat00208
부터
Figure 112018014873253-pat00209
까지 급격하게 충돌위기수준이 증가하고 있다. 그리고 상대방위 135도 부근에서는 충돌위기수준이 어느 정도 균일하게 나타나고, 상대방위 180도 부근에서는 상대거리
Figure 112018014873253-pat00210
) m에서
Figure 112018014873253-pat00211
이 나타나고 1
Figure 112018014873253-pat00212
m 부근에서
Figure 112018014873253-pat00213
이 나타난다.
반면, 도 7의 우측 그래프에 도시된 PCR 데이터의 경우는 상대거리 감소에 따라서 충돌위기수준이 균일하게 증가하고 있다. 따라서, ES 데이터의 경우, 상대방위 135도 부근을 제외하고 상대거리가 어느 정도 감소한 지점부터 충돌위기수준이 급격하게 증가함을 나타내고, PCR 데이터의 경우는 상대거리에 대해서 대략 일정한 크기로 충돌위기수준이 나타난다.
도 8은 SD 영역과 추정된 CR 영역을 비교하기 위하여 SD 영역과 CR 영역을 겹쳐서 나타낸 그래프이다
도 7은 좌표의 중심을 3
Figure 112018014873253-pat00214
)m로 나타내었으나 도 8은 좌표의 중심을 0m로 나타낸 것이다. 이와 같이, 도 8의 좌표중심을 0m로 정하여 SD 영역을 CR의 영역과 겹쳐서 나타낸다. 도 8에 나타낸 SD 영역은 CR 영역과 시각적으로 쉽게 비교하기 위하여 SD의 반경을 2배로 증가시켜서 나타낸 것이다.
도 8에서 x-축은 0.0m부터 3
Figure 112018014873253-pat00215
)m까지의 상대거리이고, y-축은 0.0m을 중심으로 상하로 나타낸 상대거리이다. 등고선은 상대방위 0도부터 180도까지에 대해서 동일한
Figure 112018014873253-pat00216
(
Figure 112018014873253-pat00217
)을 나타낸다. 등고선의 제일 바깥쪽 반원은 상대거리 0.0m를 임의로 연결한 것이다.
도 8의 좌측 그래프는 SD 영역과 ES 데이터에 대한 CR 영역을 도시한 것으로, 충돌회피를 위한 공간적인 측면에서 SD의 경우는 상대방위 0도부터 90도 사이의 공간이 90도부터 180도 사이의 공간과 비교하여 더 넓다. 그리고 SD 데이터에서 조사한 결과 상대방위 30도 부근이 가장 상대거리가 길게 나타난다. 이에 반하여 ES 데이터에 대한 CR 영역은 상대방위 135도 부근이 다른 상대방위와 비교하여 더 먼 상대거리에서부터
Figure 112018014873253-pat00218
이 나타나고, 충돌위기수준이 나타나는 공간이 더 넓게 형성된다. 결국, SD는 상대방위 30도 부근의 선박 조우상황에 더 많이 공간이 할당되어 있으나, ES 데이터에 대한 CR 영역은 상대방위 135도 부근의 선박 조우상황에 대해서 충돌위기수준이 나타나는 공간이 더 넓게 나타난다.
도 8의 우측 그래프는 SD 영역과 PCR 데이터에 대한 CR 영역을 도시한 것으로, CR 영역의 경우 상대방위 135도 부근에서 다른 상대방위와 비교하여 더 먼 상대거리에서
Figure 112018014873253-pat00219
가 나타난다. 그러나 상대방위 135도 부근의 충돌위기수준이 차지하는 공간은 다른 상대방위와 비교하여 현저한 특징을 나타내지 않는다. 반면, 상대방위 30도 부근에서 넓은 공간을 차지하는 SD와 비교하면 CR 영역은 상대방위 135도 부근에서 넓게 분포된다.
이상에서 설명된 SD 영역과 ES 및 PCT 데이터를 이용한 CR 영역과의 비교 결과를 요약하면, SD와 ES 데이터에 대한 CR 영역의 경우에 SD는 상대방위 30도 부근에 더 넓은 공간이 할당되나 CR 영역의 경우는 상대방위 135도 부근으로 더 넓은 충돌위기수준 공간이 나타난다. SD와 PCR 데이터에 대한 CR 영역의 경우에도 CR 영역이 상대방위 135도 부근으로 충돌위기수준의 분포 형태가 편중되어 나타난다. ES 데이터에 대한 CR 영역과 PCR 데이터에 대한 CR 영역의 경우는 모두 상대방위 135도 부근에서 동일한 크기의 충돌위기수준이 더 먼 거리에서 나타난다. 따라서 SD에서는 상대방위 30도 부근에 더 넓은 공간이 할당되고, ES와 PCR에 대한 CR 영역의 경우는 상대방위 135도 부근 쪽으로 더 넓게 충돌위기수준이 분포된다.
이상에서 모델링된 CREM을 이용하여 ES 데이터와 PCR 데이터에 대한 CR 영역을 3차원 좌표에 표시하면, SD와 비교하여 서로 반대되는 현상이 나타난다. 즉, SD는 상대방위 30도 부근으로 넓은 공간이 형성되고, CR 영역의 경우는 상대방위 135도 부근에 더 넓은 지각공간이 형성된다. 이는 실제 선박충돌 회피를 위한 SD의 지리학적인 공간과 OOW가 지각하는 CR의 심리적인 공간과는 서로 다르다는 것을 의미하며, 선박을 조종하는 OOW의 심리적인 공간과 실제 선박 충돌회피를 위한 지리학적인 공간을 함께 고려해야 한다는 것을 의미한다.
다음은 충돌위기 판단을 위한 기준 값의 설정과 그 기준 값을 이용한 충돌 위험의 판단 예에 대해 설명한다. 본 발명에서는 공간적인 측면의 기준 값과 심리적인 측면의 기준 값을 포함하는 두 종류의 기준 값을 사용한다.
공간적인 측면의 기준 값을 설정하기 위해서는 선박 사이에 충돌의 위험이 있다고 판단하기 위한 기준 거리를 정해야 하는데 본 발명에서는 상대방위에 대한 SD의 거리를 기준 거리(DSD)로 정한다. 그리고 선박 간의 상대거리(RD)와 그 기준 값을 비교하여 RD > DSD인 경우 충돌의 위험이 있다고 판단한다.
심리적인 측면의 기준 값을 설정하기 위해서는 충돌의 위험이 있다고 판단하기 위한 충돌위기수준(Collision Level: CL)의 기준을 정해야 하는데 본 발명에서는 전술된 CREM에서 추정한 CL 중에서 상대방위에 대한 SD의 거리를 기준 거리(DSD)에 해당하는 CL을 충돌위기수준의 기준(CL(DSD))으로 정한다. 그리고 심리적인 측면의 기준 값을 상대거리(RD)에 대한 CL 즉, CL(RD)과 비교하여 CL(DSD) > CL(RD)인 경우 충돌의 수준이 기준 값을 초월했다고 판단한다.
도 9는 DSD의 상대거리와 CL(DSD)의 계산결과를 도시한 그래프이다. 도 9에서 DSD와 CL(DSD)은 상대방위 약 85도 부근을 기준으로 서로 반대로 변하고 있다. SD는 상대방위 18.4도에서 상대거리
Figure 112018014873253-pat00220
로 최대가 되고, 상대방위 180도에서 상대거리가
Figure 112018014873253-pat00221
로 최소가 된다. CR 영역의 경우는 상대방위 29도에서
Figure 112018014873253-pat00222
이 0.38로 최소가 되고, 상대방위 177도에서
Figure 112018014873253-pat00223
이 0.96로 최대가 된다. 따라서 DSD는 상대방위 18.4도에서 최대 상대거리이고, CL(DSD)은 상대방위 177도에서 최대의 충돌위기수준이 된다.
다음은 기준 값을 이용하여 충돌위험에 대해 판단한 결과를 토대로 경보를 설정하는 방법을 설명한다.
도 10은 본 발명에 따른 해기사의 심리적 특성을 이용한 선박의 충돌 회피 방법에서 DSD와 CL(DSD)을 이용하여 충돌위험을 알리기 위한 경보(Warning)를 설정하는 방법을 설명하기 위한 그래프이다.
도 10을 참조하면, 본 발명에서는 세 가지의 상황(Case 1, Case 2, Case 3)에 대해 경보를 설정하는 예를 나타낸다.
Case 1: 본선과 상대선박이 상대방위 40도로 충돌 조우상황이 형성된 경우이다. 양 선박 사이의 상대거리가 감소하여 DSD을 지나가는 순간 SDW1(Case 1의 상황에 대한 SD의 Warning을 의미함)을 발생하고, 계속해서 상대거리가 감소하여 CL(DSD)의 충돌위기수준을 지나가는 순간 CRW1(Case 1의 상황에 대한 CR의 Warning 의미함)를 발생한다. 따라서 OOW는 두 가지 경보를 청취할 수 있다. 먼저, OOW는 SDW1을 통하여 양 선박 사이의 상대거리가 충돌회피 동작을 취해야 할 정도로 가까워졌음을 인지하여 충돌회피를 위한 선박조종을 실시할 수 있다. 만약 OOW가 충돌회피 동작을 실시하지 않은 경우에는 두 번째 경보(CRW1)가 발생하여 충돌에 직면했음을 재차 인지할 수 있다.
DSD 뿐만 아니라 CL(DSD)까지 적용하면, DSD만을 이용함에 따라 비록 충돌회피 동작에 실패했을지라도 CRW1이 발생하기 때문에 충돌회피 동작을 취할 수 있는 기회가 OOW에게 주어일 뿐만 아니라 OOW의 주의환기 수준을 향상시킬 수 있다. 이를 통해서 충돌사고를 예방할 수 있다.
Case 2: 본선과 상대선박이 상대방위 약 85도로 충돌 조우상황이 형성된 경우이다. 양 선박 사이의 상대거리가 감소하여 DSD와 CL(DSD)를 지나가는 순간 SDW2와 CRW2가 동시에 발생한다. 충돌회피에 필요한 공간적인 거리뿐만 아니라 위험수준임을 동시에 OOW에게 알릴 수 있기 때문에 충돌회피 동작이 필요함을 알 수 있고 OOW의 주의력도 동시에 향상시킬 수 있다.
Case 3: 본선과 상대선박이 상대방위 약 140도로 충돌 조우상황이 형성된 경우이다. 특징적으로 Case 1과 반대의 경보가 발생한다. 양 선박 사이의 상대거리가 감소하여 CL(DSD)를 지나는 순간 CRW3가 발생하여 충돌의 위험성이 있음을 알린다. 계속하여 상대거리가 감소하여 DSD을 지나는 순간 SDW3가 발생하여 충돌회피 동작이 필요한 시기임을 알릴 수 있다. 이를 통해서 충돌회피 동작을 취하기 전에 CRW3에 의해서 OOW의 주의력 향상을 제고시켜서 SDW3에서의 충돌회피 동작 불이행을 방지할 수 있다.
도 11은 본 발명의 일 실시 예에 따른 해기사의 심리적 특성을 이용한 선박의 충돌 회피 절차를 도시한 다이어그램으로, 본 발명의 방법에서는 선박 자동식별장치(Automatic Identification System: AIS)를 사용할 수 있다.
도 11을 참조하면, 선박에 설치된 AIS를 하여 상대선박의 정보(AIS 데이터)를 수신한다(Step1).
다음에는 상대선박의 정보(AIS 데이터)와 본선 정보를 이용하여 두 선박 간의 상대거리(RD)와 상대방위(RB)를 계산한다(Step2).
SD의 거리(DSD)와 모델링된 CREM을 이용한 충돌위기수준(CL)을 계산한다(Step3).
RD와 DSD를 비교하여 만약 RD가 DSD보다 크면 충돌할 수 있는 거리임을 알리는 SD 경보를 발생하고(Step4 & Step5), 이에 대한 적절한 충돌회피 동작을 행한다(Step 6).
만약 RD가 DSD보다 크지 않으면 전술된 과정들(Step1~Step3)을 반복한다.
한편, CL(DSD)와 CL(RD)를 더 비교하여 만약 CL(DSD)가 CL(RD) 보다 크면 충돌위기의 수준이 기준 값을 초과했음을 알리는 CR 경보를 발생하고(Step4 & Step5), 최대의 주의력을 집중하여 상황인식의 정도를 증가시키도록 한다(Step 6).
만약 CL(DSD)이 CL(RD)보다 크지 않으면 전술된 과정들(Step1~Step3)을 반복한다.
지금까지 본 발명의 바람직한 실시 예에 대해 설명하였으나, 본 발명이 속하는 기술분야에서 통상의 지식을 가진 자는 본 발명의 본질적인 특성을 벗어나지 않는 범위 내에서 변형된 형태로 구현할 수 있을 것이다.
그러므로 여기서 설명한 본 발명의 실시 예는 한정적인 관점이 아니라 설명적인 관점에서 고려되어야 하고, 본 발명의 범위는 상술한 설명이 아니라 특허청구범위에 나타나 있으며, 그와 동등한 범위 내에 있는 모든 차이점은 본 발명에 포함되는 것으로 해석되어야 한다.

Claims (8)

  1. 해기사의 충돌 위기에 대한 심리적 특성을 이용하여 선박의 충돌을 회피하도록 지원해주는 방법에 있어서,
    상대선박의 정보와 본선 정보를 이용하여 두 선박 간의 상대거리(RD)와 상대방위(RB)를 계산하는 제1단계;
    특정 상황을 샘플링(sampling)한 이산적(discrete) 데이터 형태인 해기사의 충돌위기(CR) 지각데이터를 이용하여 상기 상대거리(RD)와 상기 상대방위(RB)에 대한 충돌위기수준(CL)을 추정하고, 그 추정된 결과를 3차원 좌표 데이터로 변환하는 충돌위기예측모델(CREM)을 모델링하는 제2단계;
    선박영역의 거리(DSD)와 상기 모델링된 충돌위기예측모델(CREM)을 이용한 상기 충돌위기수준(CL)을 계산하는 제3단계;
    두 선박 간에 충돌위험이 있다고 판단하기 위한 기준 거리에 해당하는 공간적인 측면의 기준 값과 두 선박 간에 충돌위험이 있다고 판단하기 위한 충돌위기수준(CL)에 해당하는 심리적인 측면의 기준 값을 정하는 제4단계;
    상기 상대거리(RD)와 상기 공간적인 측면의 기준 값을 비교하여 충돌 가능 거리임을 알리는 선박영역(SD) 경보를 발생하거나, 상기 상대거리(RD)에 대한 충돌위기수준(CL(RD))과 상기 심리적인 측면의 기준 값을 비교하여 충돌위기의 수준임을 알리는 충돌위기(CR) 경보를 발생하는 제5단계로 이루어지는 것을 특징으로 하는 해기사의 심리적 특성을 이용한 선박의 충돌 회피 방법.
  2. 제 1 항에 있어서,
    상기 제4단계는,
    상기 상대방위(RB)에 대한 선박영역의 거리(DSD)를 상기 공간적인 측면의 기준 값으로 정하는 것을 특징으로 하는 해기사의 심리적 특성을 이용한 선박의 충돌 회피 방법.
  3. 제 1 항에 있어서,
    상기 제4단계는,
    상기 충돌위기예측모델(CREM)에서 추정한 충돌위기수준(CL) 중에서 상기 상대방위(RB)에 대한 선박영역의 거리(DSD)에 해당하는 충돌위기수준(CL(DSD))을 상기 심리적인 측면의 기준 값으로 정하는 것을 특징으로 하는 해기사의 심리적 특성을 이용한 선박의 충돌 회피 방법.
  4. 제 1 항에 있어서,
    상기 제5단계는,
    상기 상대거리(RD)와 상기 공간적인 측면의 기준 값을 비교하여, 상기 상대거리(RD)가 상기 공간적인 측면의 기준 값 보다 크면 상기 선박영역(SD) 경보를 발생하는 것을 특징으로 하는 해기사의 심리적 특성을 이용한 선박의 충돌 회피 방법.
  5. 제 1 항에 있어서,
    상기 제5단계는,
    상기 상대거리(RD)에 대한 충돌위기수준(CL(RD))과 상기 심리적인 측면의 기준 값을 비교하여, 상기 심리적인 측면의 기준 값이 상기 상대거리(RD)에 대한 충돌위기수준(CL(RD)) 보다 크면 상기 충돌위기(CR) 경보를 발생하는 것을 특징으로 하는 해기사의 심리적 특성을 이용한 선박의 충돌 회피 방법.
  6. 제 1 항에 있어서,
    상기 제2단계에서,
    상기 충돌위기예측모델(CREM)의 입력 변수는 상기 상대방위(RB)와 상기 상대거리(RD)를 포함하고, 상기 충돌위기예측모델(CREM)의 출력 변수는 연속된 상대방위(RB)와 상대거리(RD)에 대해서 추정한 충돌위기(CR)의 상기 충돌위기수준(CL)을 포함하는 것을 특징으로 하는 해기사의 심리적 특성을 이용한 선박의 충돌 회피 방법.
  7. 제 6 항에 있어서,
    상기 제2단계에서,
    상기 충돌위기예측모델(CREM)의 상기 출력 변수는 상기 충돌위기수준(CL)을 3차원 복합지도에 표시하기 위한 좌표 값을 더 포함하는 것을 특징으로 하는 해기사의 심리적 특성을 이용한 선박의 충돌 회피 방법.
  8. 제 6 항에 있어서,
    상기 제2단계에서,
    상기 충돌위기예측모델(CREM)은 극대치(GEV) 분포의 확률밀도함수(pdf)를 이용하여 상기 입력 변수에 대한 상기 충돌위기수준(CL)을 추정하는 것을 특징으로 하는 해기사의 심리적 특성을 이용한 선박의 충돌 회피 방법.
KR1020180016715A 2018-02-12 2018-02-12 해기사의 심리적 특성을 이용한 선박의 충돌 회피 방법 KR101902730B1 (ko)

Priority Applications (4)

Application Number Priority Date Filing Date Title
KR1020180016715A KR101902730B1 (ko) 2018-02-12 2018-02-12 해기사의 심리적 특성을 이용한 선박의 충돌 회피 방법
JP2018125553A JP6679665B2 (ja) 2018-02-12 2018-06-29 航海士の心理的特性を用いた船舶の衝突回避方法
CN201810694738.0A CN110164184A (zh) 2018-02-12 2018-06-29 利用海员心理特征避免船舶碰撞的方法
US16/026,025 US20190251849A1 (en) 2018-02-12 2018-07-02 Ship collision avoidance method using psychological character of ship officer

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
KR1020180016715A KR101902730B1 (ko) 2018-02-12 2018-02-12 해기사의 심리적 특성을 이용한 선박의 충돌 회피 방법

Publications (1)

Publication Number Publication Date
KR101902730B1 true KR101902730B1 (ko) 2018-09-28

Family

ID=63721447

Family Applications (1)

Application Number Title Priority Date Filing Date
KR1020180016715A KR101902730B1 (ko) 2018-02-12 2018-02-12 해기사의 심리적 특성을 이용한 선박의 충돌 회피 방법

Country Status (4)

Country Link
US (1) US20190251849A1 (ko)
JP (1) JP6679665B2 (ko)
KR (1) KR101902730B1 (ko)
CN (1) CN110164184A (ko)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR102135081B1 (ko) 2020-04-20 2020-07-17 한화시스템 주식회사 실시간 시공간 함정 표적 db를 이용한 차세대 함정용 충돌 회피 시스템

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10854090B2 (en) * 2018-09-27 2020-12-01 United States Of America As Represented By The Secretary Of The Navy Collision avoidance system and method for a watercraft
CN110758669B (zh) * 2019-10-30 2021-01-26 青岛科技大学 一种开阔繁忙水域的双层相对领域获取方法
CN110794843B (zh) * 2019-11-15 2022-11-01 山东交通学院 基于观测器的非线性船舶时滞动力定位船鲁棒镇定系统
WO2021149447A1 (ja) * 2020-01-20 2021-07-29 古野電気株式会社 操船支援装置
CN111899489B (zh) * 2020-06-01 2022-03-01 武汉理工大学 一种船舶会遇意图辨识方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR100967886B1 (ko) 2007-10-10 2010-07-05 삼성중공업 주식회사 선박의 충돌 방지 시스템
KR101647743B1 (ko) 2015-07-07 2016-08-11 한국해양과학기술원 시계열 그래픽 표시를 이용한 선박 충돌회피 안내시스템
KR101696615B1 (ko) 2016-08-04 2017-01-16 한국해양과학기술원 시계열 그래픽 표시를 이용한 다중 선박간 충돌회피 관제시스템

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0922500A (ja) * 1995-07-07 1997-01-21 Tokimec Inc 避航支援装置
KR101800453B1 (ko) * 2015-07-07 2017-12-20 (주)뉴월드마리타임 Ais정보를 이용한 선박 충돌방지 시스템 및 방법
CN105137971B (zh) * 2015-08-03 2018-07-06 大连海事大学 一种辅助船舶进行避碰决策的方法
CN105185162B (zh) * 2015-10-26 2017-10-17 中国电子科技集团公司第二十八研究所 一种基于ais信息的多目标防撞预警方法
JP6646477B2 (ja) * 2016-03-07 2020-02-14 本田技研工業株式会社 小型船舶の接触回避装置
JP6759673B2 (ja) * 2016-03-31 2020-09-23 富士通株式会社 衝突リスク算出プログラム、衝突リスク算出方法および衝突リスク算出装置
CN107284618B (zh) * 2017-07-03 2019-10-18 吴克 一种船舶防撞的方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR100967886B1 (ko) 2007-10-10 2010-07-05 삼성중공업 주식회사 선박의 충돌 방지 시스템
KR101647743B1 (ko) 2015-07-07 2016-08-11 한국해양과학기술원 시계열 그래픽 표시를 이용한 선박 충돌회피 안내시스템
KR101696615B1 (ko) 2016-08-04 2017-01-16 한국해양과학기술원 시계열 그래픽 표시를 이용한 다중 선박간 충돌회피 관제시스템

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR102135081B1 (ko) 2020-04-20 2020-07-17 한화시스템 주식회사 실시간 시공간 함정 표적 db를 이용한 차세대 함정용 충돌 회피 시스템

Also Published As

Publication number Publication date
US20190251849A1 (en) 2019-08-15
JP2019139731A (ja) 2019-08-22
CN110164184A (zh) 2019-08-23
JP6679665B2 (ja) 2020-04-15

Similar Documents

Publication Publication Date Title
KR101902730B1 (ko) 해기사의 심리적 특성을 이용한 선박의 충돌 회피 방법
Zhen et al. A novel analytic framework of real-time multi-vessel collision risk assessment for maritime traffic surveillance
JP6759673B2 (ja) 衝突リスク算出プログラム、衝突リスク算出方法および衝突リスク算出装置
US6898528B2 (en) Collision and injury mitigation system using fuzzy cluster tracking
CN108873799A (zh) 船载智能驾驶辅助终端
CN111536962B (zh) 智能船舶的航线规划方法及装置、存储介质、计算机设备
JP2022103443A (ja) 衝突警報装置及び衝突警報方法
JP7250462B2 (ja) 船舶用航行支援装置
JP2017166971A (ja) 物体検出装置および物体検出プログラム
KR20200133732A (ko) 선박용 항행 지원 장치
CN115983627B (zh) 船舶航行环境风险预测方法、装置、电子设备和介质
US20170261996A1 (en) Location and mapping device and method
KR101275277B1 (ko) 사각단위항로 생성기법을 이용한 가항항로내 선박 피항 항로 탐색 지원 시스템
CN111913183A (zh) 车辆侧向避障方法、装置、设备及车辆
CN113433946A (zh) 一种防撞预警方法、装置、系统、电子设备以及存储介质
CN113538734A (zh) 用于处理驾驶数据的方法、装置、电子设备及存储介质
CN112101558B (zh) 基于d-s证据理论的无人船态势感知及避碰决策方法、装置
Liu et al. Spatial logical relationship analysis model of ship encounter space
CN110937082B (zh) 一种基于随机风场和海浪的船舶倾覆风险测试方法
CN117368877A (zh) 基于生成对抗学习的雷达图像杂波抑制与目标检测方法
WO2022113610A1 (ja) 船舶監視システム、船舶監視方法、情報処理装置、及びプログラム
JP6333437B1 (ja) 物体認識処理装置、物体認識処理方法および車両制御システム
KR20160081139A (ko) 전장 상황 가시화 시스템 및 방법
US11624803B2 (en) System for detecting direction of a vehicle honk and determining location of honking vehicle
JP7215503B2 (ja) 危険箇所と時刻の算出方法、危険箇所と時刻の算出装置および危険箇所と時刻の算出プログラム

Legal Events

Date Code Title Description
E701 Decision to grant or registration of patent right
GRNT Written decision to grant