• 제목/요약/키워드: 구조사상 모형

검색결과 108건 처리시간 0.043초

HEC-RAS를 적용한 아라천 홍수위산정모형 구축 (Development of Flood Routing Model in the Ara Waterway Applying HEC-RAS)

  • 노준우;이상진;박명기;심명근
    • 한국수자원학회:학술대회논문집
    • /
    • 한국수자원학회 2012년도 학술발표회
    • /
    • pp.409-413
    • /
    • 2012
  • 아라뱃길은 다양한 수리구조물의 조작을 통하여 선박운항 및 홍수처리가 이루어지며 유입홍수량, 한강수위, 그리고 서해조위 변동과 연계하여 귤현보, 서해배수문, 체절수문, 그리고 배수펌프장 등의 대표적인 구조물 운영을 지원할 수 있는 홍수위 계산모형이 필요하다. HEC-RAS모형은 이러한 수리구조물 운영을 효과적으로 경계조건으로 반영할 수 있는 장점이 있다. 특히 서해배수문의 경우 조위변동과 내수위를 동시에 고려한 수문운영이 필요하여 Rule-script기능을 적용하여 이러한 특성을 반영하였다. 모형의 검보정을 위하여 2010년 9월 홍수사상을 적용하였으며 모형으로부터 계산된 수위값과 관측된 수위값이 서로 잘 일치함을 확인하였다.

  • PDF

인공신경망과 장단기메모리 모형의 유출량 모의 성능 분석 (Comparing the Performance of Artificial Neural Networks and Long Short-Term Memory Networks for Rainfall-runoff Analysis)

  • 김지혜;강문성;김석현
    • 한국수자원학회:학술대회논문집
    • /
    • 한국수자원학회 2019년도 학술발표회
    • /
    • pp.320-320
    • /
    • 2019
  • 유역의 수문 자료를 정확하게 분석하는 것은 수리 구조물을 효율적으로 운영하기 위한 중요한 요소이다. 인공신경망(Artificial Neural Networks, ANNs) 모형은 입 출력 자료의 비선형적인 관계를 해석할 수 있는 모형으로 강우-유출 해석 등 수문 분야에 다양하게 적용되어 왔다. 이후 기존의 인공신경망 모형을 연속적인(sequential) 자료의 분석에 더 적합하도록 개선한 회귀신경망(Recurrent Neural Networks, RNNs) 모형과 회귀신경망 모형의 '장기 의존성 문제'를 개선한 장단기메모리(Long Short-Term Memory Networks, 이하 LSTM)가 차례로 제안되었다. LSTM은 최근에 주목받는 딥 러닝(Deep learning) 기법의 하나로 수문 자료와 같은 시계열 자료의 분석에 뛰어난 성능을 보일 것으로 예상되며, 수문 분야에서 이에 대한 적용성 평가가 요구되고 있다. 본 연구에서는 인공신경망 모형과 LSTM 모형으로 유출량을 모의하여 두 모형의 성능을 비교하고 향후 LSTM 모형의 활용 가능성을 검토하고자 하였다. 나주 수위관측소의 수위 자료와 인접한 기상관측소의 강우량 자료로 모형의 입 출력 자료를 구성하여 강우 사상에 대한 시간별 유출량을 모의하였다. 연구 결과, 1시간 후의 유출량에 대해서는 두 모형 모두 뛰어난 모의 능력을 보였으나, 선행 시간이 길어질수록 LSTM의 정확성은 유지되는 반면 인공신경망 모형의 정확성은 점차 떨어지는 것으로 나타났다. 앞으로의 연구에서 유역 내 다양한 수리 구조물에 의한 유 출입량을 추가로 고려한다면 LSTM 모형의 활용성을 보다 더 확장할 수 있을 것이다.

  • PDF

Nesting 기법을 이용한 돌발홍수 예측모형 개발 (Development of flash flood forecasting model using method)

  • 지희숙;배덕효
    • 한국수자원학회:학술대회논문집
    • /
    • 한국수자원학회 2012년도 학술발표회
    • /
    • pp.403-403
    • /
    • 2012
  • 최근 단시간 동안에 특정지역에 집중되는 국지적 호우에 의한 돌발홍수가 빈번히 발생하고 있으며, 이에 따른 위험과 손실이 증가하고 있는 추세이다. 현재 국내에서는 이러한 피해를 최소화하고자 돌발홍수 예측모형을 개발하고 예 경보 시스템을 구축하여 다양한 비구조적 대책을 마련하고 있다. 그러나 활용되는 예측모형의 경우 개념적 유출량인 한계유출량으로부터 돌발홍수능(Flash Flood Guidance, FFG)을 결정하여 예측 강우와 상대적인 대소 비교를 통해 돌발홍수의 발생가능성 유무를 판단하게 되는데, 문제는 산정되는 한계유출량은 개념적이기 때문에 검증이 어렵고 산정방법도 다양하여 불확실성이 높다는 단점이 있다. 이에 본 연구에서는 기존의 돌발홍수 예측 방법이 아닌, 수문모형 Nesting 기법을 이용한 돌발 홍수 예측 방법을 개발하였다. 저해상도의 대유역 기반의 유출량이 큰 영역의 경계값이 되고, 대유역을 이루고 있는 소유역을 고해상도의 작은 영역이라 할 때, 경계값인 대유역의 기반의 유출량을 참고 유출량으로 하여 소유역의 유출을 물리적 혹은 개념적으로 보다 타당하게 모의하는 방법이 수문모형 Nesting 기법이다. 이러한 기법에 필요한 강우-유출 모형으로는 대유역의 경우, SURR 모형(Sejong University Rainfall-Runoff model)을 선택하였으며, 대유역을 이루는 소유역의 유출모의는 물리적 기반의 분포형 모형인 CASC2D 모형을 이용하였다. 또한 실시간 활용을 위해서는 CASC2D 모형의 매개변수를 자동으로 추정하는 기술이 요구되며, 본 연구에서는 매개변수 전역 최적화 방법인 SCE-UA(The Shuffled Complex Evolution, University of Arizona) 기법을 활용하였다. 본 연구에서 사용한 수문모형의 적용성을 평가한 결과 대상유역에 대한 적용성이 높은 것으로 나타났으며, 연계된 두 모형의 유출거동이 유사하게 나타난 것으로 확인되었다. 본 연구에서는 Nesting 기법을 이용하여 0.5m 하천 수위의 상승 여부에 따라 돌발홍수의 발생 가능성을 예측하는 기법을 제안하였으며, 돌발홍수 사례와 일반호우사상으로부터 이 방법의 적용성을 평가하였다. 실제 돌발홍수가 발생한 유역을 선정하고 연계된 두 모형을 대상 유역에 적용한 결과 Nesting 기반의 돌발홍수 예측방법은 기존의 한계유출량 산정 방법에서 반영하지 못한 사상을 적절히 반영한 것으로 나타났다. 본 연구에서 개발한 Nesting 기법을 이용한 돌발홍수 예측모형은 일반적인 강우량 비교의 돌발홍수 예측방법에서 벗어나 새로운 돌발홍수 예측방법을 제안한 측면에서 큰 의미가 있다고 사료되며, 이러한 연구 결과는 실시간 돌발홍수 예측 시스템의 기본 모형으로 활용이 가능할 것으로 판단된다.

  • PDF

거대 구조물의 유체-구조 연계 해석을 위한 효과적인 보간기법에 대한 연구 (A Study on the Effective Interpolation Methods to the Fluid-Structure Interaction Analysis for Large-Scale Structure)

  • 이기두;이영신;김동수;이대열
    • 한국항공우주학회지
    • /
    • 제37권5호
    • /
    • pp.433-441
    • /
    • 2009
  • 대부분의 자연현상은 다학제 특성을 갖고 표현된다. 유체-구조 연계(FSI) 문제의 경우 기존에 검증된 전산유체 해석 프로그램 및 구조해석 프로그램을 그대로 사용할 수 있다는 장점 때문에 약결합 방식이 일반적으로 이용된다. 그러나 약결합을 이용하여 해석을 수행하기 위해서는 서로 다른 특성을 갖는 격자시스템으로 발생되는 자료의 교환을 위해서 보간 및 사상이 필수적이다. 본 연구에서는 전역지지 및 국부지지 방사기저함수(RBF)를 이용한 보간 및 가상일의 원리를 적용한 사상의 성능을 단순 3차원 형상에 적용하여 검토하였다. 국부지지 RBF에 공간분할 트리의 일종으로 빠른 공간 탐색을 가능하게 해주는 kd-tree를 사용하는 경우 효과적으로 거대 구조물의 FSI에도 보간 및 사상이 적용 가능함을 여객기 형상의 항공기 모형을 이용하여 제시하였다.

소양호 탁수거동 해석을 위한 2차원 수치 모형의 보정 및 검증 (Validation and Calibration of 2D Numerical Model for Analysis of Turbidity Current Regimes in Soyang Reservoir)

  • 류인구;정세웅;윤성완
    • 한국수자원학회:학술대회논문집
    • /
    • 한국수자원학회 2010년도 학술발표회
    • /
    • pp.253-257
    • /
    • 2010
  • 최근 들어 이상 강우 현상으로 인한 기록적인 집중 강우와 더불어 토지 이용 변화로 인해 탁수의 발생 빈도가 증가하고 있다. 강우 유출로 유입한 탁수는 저수지 내 장기간 체류하면서 하류 하천의 수질 및 수생태계 뿐만 아니라 저수지 내부의 영양단계에도 많은 영향을 준다. 특히, 성층화된 저수지의 경우 높은 영양염류 농도를 포함한 탁수는 밀도류 거동 특성을 보이면서 수평 및 수직 혼합 과정을 거치면서 국부적인 부영양화 현상과 조류의 수화 현상의 원인이 되고 있다. 따라서 대형 저수지의 수질관리에 있어 하천 유입 탁수의 밀도류 해석, 저수지의 수온 성층 구조 변화, 부유입자의 동력학적 해석이 중요한 요소로 부각되어 왔다. 본 연구에서는 소양호를 연구 대상 지역으로 선정하여 2005년과 2007년 수문 사상을 바탕으로 2차원 횡방향 수치 모형을 구축하였다. 수치모형을 통해 수온 성층 구조의 재현성을 확인하였으며, 다양한 탁수 거동 모형을 구축하여 적용성을 평가하였다. 유입수의 SS(Suspended Solid)를 단일 입경으로 가정한 TM-1 모형, SS의 입경분포에 따라 3개의 그룹(SSi)으로 구분한 TM-2 모형, 3개 그룹을 포함하면서 저수지내 탁수 장기화로 인한 탁수 저감 효과를 1차 반응상수로 매개 변수화(유기물 함량($a_0$) ${\times}$ 분해속도(${\lambda}_a$))하여 수정된 지배방정식을 적용한 TM-3 모형을 사용하였다. 각각의 탁수 거동 모형은 2005년과 2007년 수문 조건에서 수온 성층 구조를 잘 재현하였다. TM-1 모형과 TM-2 모형을 비교해보면, 탁수 중심축의 최고 탁도에 대한 예측 성능은 TM-2 모형이 우수한 결과를 나타냈었다. 하지만, 장기 탁수 모의 시 저수지 수중 잔류 SS가 지속적으로 높게 나타나 중층 탁도를 과대평가하는 경향을 보였다. TM-3 모형이 TM-2 모형에 비해 수심별 탁도 분포에 대한 중심축 탁도가 저평가되는 경향을 보였지만 저수지 내 잔류 탁도에 대한 영향 부분에서는 개선된 결과를 나타내었다. 본 연구 결과는 저수지 탁수 밀도류 해석 및 운영 시스템에 활용 될 수 있으며, 선택 취수 설비 등의 수리 구조물의 영향 평가에 활용할 수 있다.

  • PDF

동적수자원평가모형을 이용한 순천동천 유역의 홍수량 산정 (Estimation of flood in Suncheon Dongcheon watershed using dynamic water resources assessment Tool)

  • 김덕환;김현준;장철희
    • 한국수자원학회:학술대회논문집
    • /
    • 한국수자원학회 2022년도 학술발표회
    • /
    • pp.285-285
    • /
    • 2022
  • 기후변화가 현실화되면서 수자원평가 (Water Resources Assessment)에 대한 관심과 중요성이 높아지고 있다. 본 연구에서는 '주민공감 문제기획리빙랩' 대상지인 순천은 동천을 중심으로 홍수량을 정량적으로 분석하였다. 순천시의 가장 시급한 사안 중 하나인 범람 및 침수 문제로, 최근 3년(2018~2020)간의 집중호우로 인한 내수배제로 주택 및 도로침수, 산사태 등의 피해를 겪었다. 시기마다 고질적으로 반복되는 동천 인근 지역의 침수문제를 사전에 예방하고 피해의 빈도나 규모를 줄이기 위하여 분석을 수행하였다. 이에 본 연구에서는 환경부의 지원을 받아 한강홍수통제소와 한국건설기술연구원이 공동으로 개발한 동적수자원평가모형(DWAT, Dynamic Water Resources Assessment Tool)을 이용하여 정량적으로 홍수량 산정을 하고자 한다. 본 모형은 전 세계가 무료로 이용할 수 있는 수자원평가도구로 사용자의 편의를 위해 GIS전처리 기능을 포함하고 있어, 자동으로 유역 매개변수 및 면적 평균강우량을 Thiessen method를 사용하여 산정할 수 있다. 또한, 물의 순환과정을 투수 및 불투수지역으로 구분되며, 투수지역은 1개의 토양층과 1개의 불압대수층으로 구성되고, 유출기여역과 함양역으로 유역을 분할하여 적용할 수 있으며, 대수층을 통하여 지하수의 흐름을 산정할 수 있다. 기상청에서 제공하는 기상자료를 분석하여 과거 관측 강우사상 3개를 선정하여 검·보정을 수행하였으며, 그 결과 모형 효율계수(Nash-Sutcliffe efficiency) 및 결정계수(Coefficient of Determination)가 0.78~0.94, 0.82~0.94로 우수한 모의 결과를 산정할 수 있었다. 빈도별 확률강우량을 Huff 4분위법을 사용하여 확률홍수량을 산정하였다. 미래 홍수량 증감량 산정을 위하여 RCP(Representative Concentration Pathways) 기후변화 시나리오를 사용하였다. 관측값과 모의값의 누적확률분포 이용하여 모의값의 확률분포를 관측값의 확률분포에 사상시키는 방법인 분위사상법(Quantile Mapping)을 사용하여 시나리오자료를 보정하였다. 본 연구에서 산정한 홍수량을 바탕으로 침수피해를 막기 위한 구조적 및 비구조적 방안을 위한 기초자료로 사용될 것으로 판단된다.

  • PDF

한강하류부 지천의 유입량이 한강수위에 미치는 영향 분석 (Effect of Tributary Inflow on Downstream Water Level of the Han River)

  • 이정규;장기환;양희성
    • 한국수자원학회:학술대회논문집
    • /
    • 한국수자원학회 2004년도 학술발표회
    • /
    • pp.1366-1370
    • /
    • 2004
  • 한강은 남한 제1의 하천으로 지리적으로 한반도의 중심부에 위치하여 동에서 서로 흐르는 대하천이다. 매년 발생하는 홍수에 대한 한강하류부의 홍수관리는 대단히 중요하다. 본 연구는 한강의 지천유입량이 한강본류 수위에 미치는 영향을 수리학적 홍수추적 모형을 이용하여 분석하였다. 수리학적 홍수추적 모형은 하천의 지형적 특성과 여러인자들을 반영하여 하천의 흐름을 해석하기 때문에 하천에 설치되어 있는 교량, 수중보, 지천유입량, 조도계수 등을 적절히 반영하였을 때 보다 정확한 결과를 얻을 수 있다. 본 연구에서는 1차원 부정류 흐름으로 해석하였으며 미국 NWS(National Weather Service)에서 기존의 DWOPER 모형과 댐파괴 모형인 DAMBRK 모형을 통합하여 개선한 FLDWAV 모형을 사용하였다. 이 모형의 지배방정식은 연속방정식과 운동량방정식으로 구성되어 확장된 Saint-Venant 방정식이고, 수치기법으로는 가중4점음해법을 이용하려 비선형 연립방정식을 Newton-Raphson 방법으로 해석한다. 모형의 상류경계는 팔당댐 시간 방류량을 하류경계로는 전류수위표 지점의 수위를 이용하였다. 한강하류부에는 많은 지천들이 흘러들어 유입되고 있는데 지천유입량이 한강수위에 미치는 영향을 분석하기 위하여 비교적 유역면적이 큰 왕숙천, 탄천, 중랑천, 안양천 등의 지천유입량을 고려하였다. 2000년과 2001년에 발생한 2개의 실측사상을 FLDWAV 모형을 이용하여 계산된 수위값과 검증해 본 결과 실측수위를 잘 재현하고 있는 것을 알 수 있었다. 검증된 모형을 가지고 동일한 사상에 대해서 지천 유입량을 고려하지 않고 모의를 한 결과 지천유입량을 고려하였을 때와 고려하지 않았을 때의 지천유입량에 대한 수위상승량은 $1\~5cm$ 정도로 나타났다. 지천유입량에 의한 수위상승이 어느정도 일어나고 있지만 홍수기간에 이 정도의 수위상승은 한강의 흐름이나 한강하류부에 설치된 구조물의 침수피해에 크게 영향을 주지않을 것으로 판단되었다. 한강하류부의 수위는 지천유입량의 영향보다는 팔당댐의 방류량에 의해서 좌우되고 있고 홍수시에는 한강의 수위상승으로 인하여 지천에 미치는 배수영향이 훨씬 더 클 것으로 판단된다.

  • PDF

기계학습 기법을 이용한 다지점 일강수량 모의 모형 개발 (Deelopment of a Multisite Daily Rainfall Simulation Model Using a Machine Learning)

  • 소병진;권현한
    • 한국수자원학회:학술대회논문집
    • /
    • 한국수자원학회 2017년도 학술발표회
    • /
    • pp.83-83
    • /
    • 2017
  • 수자원공학에서 일강수량 모의기법은 다양한 목적으로 활용되고 있지만, 일반적으로 홍수와 가뭄의 영향을 고려할 수 있는 수공구조물의 위험도 및 신뢰성 평가 및 수자원 계획을 수립하기 위한 입력 자료생성을 목적으로 활용된다. 유역 단위의 분석시 단일 지점에 대한 강수 모의 기법을 적용할 경우 각각의 지점에서 관측된 강수 자료의 시계열 및 통계치 특성이 효과적으로 재현되지만 공간적으로 발생하는 즉, 지점 간의 종속관계를 재현하지 못하는 문제가 발생한다. 이러한 이유로 공간적인 전이 특성이 있는 가뭄 분석 및 유역내 유출량의 공간적 변동 특성 분석에 단일지점별 모의 결과를 이용할 경우 관측 자료와 상반된 공간적 변동성으로 인하여 잘못된 가뭄 및 유출 분석 결과가 도출되는 문제점이 있다. 따라서, 실제적으로 발생하는 강수 특성을 반영한 유역 단위의 홍수 및 가뭄 등의 수문 분석을 위해서는 지점간의 종속성을 반영할 수 있는 다지점 강수 모의 모형의 적용이 필수적이다. 본 연구에서는 다지점 모의에 있어서, Wilks 모형의 지점별 시변동 특성과 공간상관성 재현 능력, HMM 모형이 갖는 강수 사상별로 분포된 양적 분포 패턴 재현 능력을 복합적으로 나타낼 수 있는 새로운 다지점 일강수량 모의 모형인 기계학습 기반 범주화 기법을 이용한 다지점 일강수량 모의 모형(ML-MRS)을 개발하였다. 또한, 지점별 강수량에 적용되는 확률분포모형은 Gamma 분포로 구성된 혼합모형을 적용하여 단일 확률 분포 모형의 자료 적합 문제를 개선하였다. 모의를 통한 일강수량 시계열 자료는 일 강수자료의 통계량을 효과적으로 모의하였으며, 다지점 모형의 모의 결과를 적용한 가뭄 모의 결과 관측 자료에서 나타나는 공간적 패턴이 재현되었다. 본 모형은 시 공간적 사상을 효과적으로 재현함으로서 지역의 변동특성을 반영한 가뭄, 홍수, 기상 현상 분석 등 활용도가 매우 높을 것으로 판단된다.

  • PDF

WRF 예측강우를 활용한 홍수 및 침수예측에 관한 연구 (A Study on the Flood and Inundation Prediction using Forecasted Rainfall of the WRF Model)

  • 윤성심;트란앙푸옹;배덕효
    • 한국수자원학회:학술대회논문집
    • /
    • 한국수자원학회 2010년도 학술발표회
    • /
    • pp.263-267
    • /
    • 2010
  • 최근 지구온난화, 엘니뇨 및 라니냐 등 지구환경 변화에 따른 기후변화의 영향으로 지구상의 많은 지역에서 집중호우가 발생하고 있으며 우리나라도 예외 없이 매년 되풀이되고 있다. 이로 인해 발생하는 홍수피해를 경감하기 위해서 홍수조절용 다목적 댐 건설과 같은 구조적 대책과 홍수를 사전에 예측할 수 있는 홍수예경보 시스템 구축과 같은 비구조적 대책의 마련이 필요하다. 일반적인 홍수예경보 시스템은 강우 관측치를 강우-유출 모형 및 수리해석 모형의 입력 자료로 하여 홍수량 및 홍수위를 계산하고 그 결과를 이용하여 운영된다. 그러나 집중호우와 같은 악기상 조건에서는 관측강우자료를 이용한 유출해석 결과로 홍수예경보 시스템을 운영 할 경우 예방 대응시간의 부족으로 인해 방재 효율성이 떨어지게 된다. 따라서 미래에 발생할 강우를 사전에 예측하고, 이를 효율적으로 유출 모형과 연계하여 홍수발생 이전에 홍수발생 가능성을 예측할 수 있는 홍수 모의시스템을 구축하는 것이 필요하다. 이를 위해 본 연구에서는 중규모 수치예보모형인 WRF 모형(Weather Research and Forecasting model)으로 모의된 2007년 태풍 '나리' 사상의 예측강우를 이용하여 유역평균강우를 산정하였으며, 산정된 예측강우를 도시유역유출모형인 SWMM과 2차원 침수모의가 가능하도록 개선한 CASC2D 모형에 활용하여 침수현상을 모의하였다. 실제 침수흔적과 모의된 결과의 비교를 통해 예측강우를 이용한 침수예측 및 홍수예보의 가능성을 평가한 결과, 과소추정된 예측강우의 영향으로 인해 모의된 침수심이 실제보다 작게 발생하였으나 침수발생 위치는 대체적으로 정확하게 모의하는 것으로 나타났다.

  • PDF

레이더 반사도 자료의 오차구조 및 거리오차 보정 (The Error Structure of the Radar Reflectivity and the Correction of the Range Dependent Error)

  • 윤정수;유철상;김정호;전창현
    • 한국수자원학회:학술대회논문집
    • /
    • 한국수자원학회 2012년도 학술발표회
    • /
    • pp.468-468
    • /
    • 2012
  • 레이더 강우를 적극적으로 활용하기 위해서는 레이더 강우에 포함된 각 오차들에 대한 특성을 파악하고 정량화하는 것이 무엇보다 중요하다. 본 연구에서는 레이더 반사도의 오차구조를 파악하고 그러한 오차구조를 갖는 반사도로부터 표출한 CAPPI의 거리오차를 보정하였다. 이러한 거리오차 파악을 위해서는 참 값으로 가정할 수 있는 기준 반사도가 필요하며 본 연구에서는 VPR 모형으로부터 기준 반사도인 지상 반사도를 추정하였다. 그 결과 일정한 VPR 모형을 적용하게 되면 거리와 상관없이 오차는 일정하고 오직 고도에 의해서만 영향을 받는다. 그러나 일정 고도에서의 반사도 표출 방법인 CAPPI는 지구곡률효과로 인해 실제로 거리가 멀어 질수록 관측 고도가 높아진다. 이에 따라서 오차는 거리가 멀어질수록 커지게 된다. 이는 실제 호우사상에 적용한 결과에서도 유사하게 나타났다. 강릉 기상 레이더의 경우 1.5km CAPPI는 약 100km까지 1.5km 고도를 유지하다 그 이상부터 고도가 점점 높아진다. CAPPI의 오차를 거리에 따라 분포시킨 결과에서도 100km까지는 어느 정도 일정한 오차를 보이다 그 이상부터 오차가 점점 증가하는 것으로 나타났다. CAPPI의 오차를 2차원 평면으로 나타낸 결과에서도 호우가 전반적으로 퍼져있는 시점부터 원거리에서 큰 오차를 보이고 있다. 이는 오차의 평균에서 더욱 명확히 나타났다. 이와 같이 CAPPI는 원거리 자료에서 오차가 크게 나타나고 있다. 이에 CAPPI에 포함된 거리오차를 VPR 모형을 이용하여 보정하였다. 그 결과 원거리에서의 오차가 감소하였음을 확인하였다.

  • PDF