KR20110124525A - 하천제방 안전관리를 위한 제방붕괴 예측시스템 - Google Patents

하천제방 안전관리를 위한 제방붕괴 예측시스템 Download PDF

Info

Publication number
KR20110124525A
KR20110124525A KR1020100043949A KR20100043949A KR20110124525A KR 20110124525 A KR20110124525 A KR 20110124525A KR 1020100043949 A KR1020100043949 A KR 1020100043949A KR 20100043949 A KR20100043949 A KR 20100043949A KR 20110124525 A KR20110124525 A KR 20110124525A
Authority
KR
South Korea
Prior art keywords
displacement
information
collapse
slope
river
Prior art date
Application number
KR1020100043949A
Other languages
English (en)
Other versions
KR101178557B1 (ko
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 KR20100043949A priority Critical patent/KR101178557B1/ko
Publication of KR20110124525A publication Critical patent/KR20110124525A/ko
Application granted granted Critical
Publication of KR101178557B1 publication Critical patent/KR101178557B1/ko

Links

Images

Classifications

    • EFIXED CONSTRUCTIONS
    • E02HYDRAULIC ENGINEERING; FOUNDATIONS; SOIL SHIFTING
    • E02DFOUNDATIONS; EXCAVATIONS; EMBANKMENTS; UNDERGROUND OR UNDERWATER STRUCTURES
    • E02D1/00Investigation of foundation soil in situ
    • E02D1/02Investigation of foundation soil in situ before construction work
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B21/00Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant
    • G01B21/32Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant for measuring the deformation in a solid
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01MTESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M5/00Investigating the elasticity of structures, e.g. deflection of bridges or air-craft wings
    • G01M5/0033Investigating the elasticity of structures, e.g. deflection of bridges or air-craft wings by determining damage, crack or wear
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01MTESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M5/00Investigating the elasticity of structures, e.g. deflection of bridges or air-craft wings
    • G01M5/0041Investigating the elasticity of structures, e.g. deflection of bridges or air-craft wings by determining deflection or stress
    • G01M5/005Investigating the elasticity of structures, e.g. deflection of bridges or air-craft wings by determining deflection or stress by means of external apparatus, e.g. test benches or portable test systems
    • G01M5/0058Investigating the elasticity of structures, e.g. deflection of bridges or air-craft wings by determining deflection or stress by means of external apparatus, e.g. test benches or portable test systems of elongated objects, e.g. pipes, masts, towers or railways
    • GPHYSICS
    • G08SIGNALLING
    • G08CTRANSMISSION SYSTEMS FOR MEASURED VALUES, CONTROL OR SIMILAR SIGNALS
    • G08C17/00Arrangements for transmitting signals characterised by the use of a wireless electrical link
    • G08C17/02Arrangements for transmitting signals characterised by the use of a wireless electrical link using a radio link

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Testing Or Calibration Of Command Recording Devices (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Soil Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Mining & Mineral Resources (AREA)
  • Paleontology (AREA)
  • Civil Engineering (AREA)
  • General Engineering & Computer Science (AREA)
  • Structural Engineering (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Chemical & Material Sciences (AREA)
  • Pit Excavations, Shoring, Fill Or Stabilisation Of Slopes (AREA)

Abstract

원격지의 사면붕괴 측정센서에서 계측된 거동 데이터를 실시간으로 분석하여 제방의 안전성 여부를 판단할 수 있도록 한 하천제방 안전관리를 위한 제방붕괴 예측시스템이 개시된다.
개시된 하천제방 안전관리를 위한 제방붕괴 예측시스템은, 지중에 매설되며, 응력에 따라 탄성력이 가변하여 변위를 발생하고, 상기 응력이 정해진 기준 이상 변경될 경우 내장된 측정센서를 정해진 위치로부터 이탈시켜 거동 데이터를 발생하는 사면붕괴 측정장치를 포함하는 현장 시스템과; 현장에 설치된 상기 현장 시스템에서 전송된 거동데이터를 분석하여 하천제방의 상태를 계측하는 상태 예측수단을 구비한다.

Description

하천제방 안전관리를 위한 제방붕괴 예측시스템{Levee failure forecasting system for safety management of levee facilities}
본 발명은 하천 제방의 붕괴를 예측하는 기술에 관한 것으로서, 더욱 상세하게는 원격지의 사면붕괴 측정센서에서 계측된 거동 데이터를 실시간으로 분석하여 제방의 안전성 여부를 판단할 수 있도록 한 하천제방 안전관리를 위한 제방붕괴 예측시스템에 관한 기술이다.
우리나라는 지난 2002년도부터 2004년도까지 태풍 '루사', '매미' 등 크고 작은 태풍의 내습과 잦은 강우로 인하여 모든 유역의 대·소 하천 모두가 호안(제방) 붕괴, 유실 및 범람으로 많은 수해를 입었다.
이는 하천 호안(제방)의 시설물 설계 당시 기존의 호우사상과는 달리 특정 지역에 국지적이며 집중적(게릴라 성)인 강우가 발생하여 예측이 불가능하고, 천문학적 확률의 초과 강우로 인하여 하천의 홍수지속시간이 3일('02년도 설계개념)인 72시간을 넘어섬에 따라 제방의 제체가 포화되어 슬라이딩 등이 발생하고 파이핑 및 제체, 구조물 부위의 누수 등이 발생한 것으로 추정된다.
따라서, 효과적인 제방안전관리를 위하여 수자원(수리, 수문학적), 토질, 구조적 측면의 다각적인 내용을 고려할 수 있는 미래지향적 첨단 안전관리기법의 제시가 절실히 필요한 상황이다.
한편, 제방붕괴의 위험성이 있는 다양한 실 계측 데이터들이 토목기술자들에게 주어졌지만, 이를 해석, 가공 후 실시간 사면 거동 양상을 분석하고, 사면안정에 대한 결론을 유도하는 일은 결코 쉬운 일이 아니다.
현재 제방의 누수부위를 탐지하거나 제체 구조의 건전성을 평가하기 위한 다양한 방법이 국내외적으로 제안되고 있다.
국내에서는 자연전위탐사나 전기 비저항탐사, 저온탐사와 같은 물리 탐사법이 다양하게 적용되고 있다. 이러한 방법은 분해능을 높이기 위하여 전극 간격을 매우 작게 설정할 경우 가탐심도가 낮아지게 되므로, 일정 규모 이상의 하천에서는 전극 간격을 적정간격 이상으로 유지해야 한다.
국외에서는 계측자료에 의한 위험 예측에 관한 연구 중 수학적 해석기법과 유사한 연구로서 비선형동적모델(NDS)을 이용한 연구가 있으며, 이는 지진을 포함한 구조물에 대한 붕괴시간을 해석하는 방법을 제안한 것이다.
그러나 이러한 모든 기존의 하천제방 계측관리방법은, 효과적인 과업수행을 위해 제방별 관리기준치를 산정하여 계측결과와 설계시 산정한 값과의 비교분석을 수행하는 방식이다.
현재 토목계측분야에서는 계측자료와 비교할 관리기준치가 명확하지 않은 것으로 확인되었으며, 현재 사용되고 있는 관리기준치는 일본과 미국의 자료를 그대로 인용하고 있는데, 이 또한 구조물의 설계시 산정한 값이 아니라 경험적인 값을 일괄적으로 유사 구조물에 사용해야 하는 단점이 있다.
특히, 수리시설물은 타 구조물에서 사용하는 단순 비교 형태의 관리기준방법을 적용하기에는 너무나 많은 변수들이 존재하므로 수리시설물의 특성에 적합한 관리기준치 산정방법이 절실히 요구된다.
이에 본 발명은 상기와 같은 종래의 요구에 의해 제안된 것으로서,
본 발명이 해결하려는 과제는 원격지의 사면붕괴 측정센서에서 계측된 거동 데이터를 실시간으로 분석하여 제방의 안전성 여부를 판단할 수 있도록 한 하천제방 안전관리를 위한 제방붕괴 예측시스템을 제공하는 데 있다.
본 발명이 해결하려는 다른 과제는 실시간 감시를 통해 제방 변위에 따른 거동양상을 예측할 수 있는 하천제방 안전관리를 위한 제방붕괴 예측시스템을 제공하는 데 있다.
상기와 같은 과제를 해결하기 위한 본 발명에 따른 "하천제방 안전관리를 위한 제방붕괴 예측시스템"은,
지중에 매설되며, 응력에 따라 탄성력이 가변하여 변위를 발생하고, 상기 응력이 정해진 기준 이상 변경될 경우 내장된 측정센서를 정해진 위치로부터 이탈시켜 거동 데이터를 발생하는 사면붕괴 측정장치를 포함하는 현장 시스템과;
현장에 설치된 상기 현장 시스템에서 전송된 거동데이터를 분석하여 하천제방의 상태를 계측하는 상태 예측수단을 포함한다.
상기 사면붕괴 측정장치는,
응력에 따라 탄성력이 가변하여 변위를 발생하고, 상기 응력이 정해진 기준 이상 변경될 경우 내장된 측정센서를 정해진 위치로부터 이탈시켜 변위 정보를 발생하는 센서 변위 수단을 포함하고,
상기 현장 시스템은,
상기 센서 변위 수단에서 획득한 변위 정보를 연결된 통신망에 대응하는 데이터 포맷으로 변환하여 거동 데이터로 상기 원격에 위치한 상태 예측수단에 전송하는 변위정보 전송수단을 포함한다.
상기 사면붕괴 측정장치는,
내장된 측정센서에 습기가 유입되는 것을 방지하는 방수관을 더 포함하고,
상기 센서 변위 수단은,
응력에 따라 탄성력이 가변하여 상기 측정센서를 정해진 위치로부터 이탈시키는 복수의 압축 스프링을 구비한 것을 특징으로 한다.
상기 상태 예측수단은,
상기 변위정보 전송수단에서 전송된 복수 타깃의 변위 정보를 수신하여 처리하는 호스트 컴퓨터와;
상기 호스트 컴퓨터에서 처리한 변위 정보를 데이터베이스에 저장하고 이를 관리하기 위한 데이터베이스 서버와;
상기 데이터베이스에 저장된 복수 타깃의 변위 정보를 그룹화하고, 그 그룹화한 타깃의 변위 정보를 SAM 및 MDM 알고리즘으로 분석하여 사면 거동을 판단하고, 그 판단 결과 거동 가능한 그룹이 존재할 경우 승법 분해를 수행하여 미래 상황을 예측하며, 그 결과를 하천 호안 시설물의 미래 상황 예측 정보로 제공하는 예측수단과;
상기 예측수단에서 예측한 하천 호안 시설물의 상태 정보를 통신망을 통해 일반 사용자 및 정보 참조자의 단말기로 전송해주는 예측정보 전송수단을 포함한다.
본 발명에 따르면, 원격지의 사면붕괴 측정센서에서 계측된 거동 데이터를 실시간으로 분석하여 제방의 안전성 여부를 정확하게 판단할 수 있는 장점이 있다.
또한, 본 발명에 따르면 정확하게 사면붕괴를 예측할 수 있으므로, 하천 호안(제방) 시설물의 대규모 손상이나 붕괴로 인한 피해를 사전에 방지할 수 있도록 도모해주는 효과도 있다.
도 1은 본 발명에 따른 하천제방 안전관리를 위한 제방붕괴 예측시스템의 구성도.
도 2는 본 발명에 적용된 사면붕괴 측정장치의 개략구성도.
도 3a 및 도 3b는 사면 붕괴시 측정센서가 이탈되는 동작상태도.
이하 본 발명의 바람직한 실시 예를 첨부한 도면에 의거 상세히 설명하면 다음과 같다. 본 발명을 설명하기에 앞서 관련된 공지 기능 혹은 구성에 대한 구체적인 설명이 본 발명의 요지를 불필요하게 흐릴 수 있다고 판단되는 경우에는 그에 대한 상세한 설명은 생략한다.
도 1은 본 발명에 따른 "하천제방 안전관리를 위한 제방붕괴 예측시스템"의 전체 구성도로서, 현장 시스템(100), 하천 호안 시설물의 상태 예측수단(200), 통신망(300), 감시 시스템(400)으로 구성된다.
현장 시스템(100)은 응력에 따라 탄성력이 가변하여 변위를 발생하고, 상기 응력이 정해진 기준 이상 변경될 경우 내장된 측정센서(130)를 정해진 위치로부터 이탈시켜 변위 정보를 발생하는 센서 변위 수단(110)을 포함하는 사면붕괴 측정장치(100)와, 상기 센서 변위 수단(110)에서 획득한 변위 정보를 연결된 통신망에 대응하는 데이터 포맷으로 변환하여 거동 데이터로 상기 원격에 위치한 상태 예측수단(20)에 전송하는 변위정보 전송수단(150)을 포함한다.
사면붕괴 측정장치(100)는, 도 2, 도 3a 및 도 3b에 도시한 바와 같이, 내장된 측정센서(130)에 습기가 유입되는 것을 방지하는 방수관(120)을 더 포함하고, 상기 센서 변위 수단(110)은 응력에 따라 탄성력이 가변하여 상기 측정센서(130)를 정해진 위치로부터 이탈시키는 복수의 압축 스프링(111 ~ 113)을 구비한다. 여기서 압축 스프링은 일 예에 불과하며, 응력에 따라 길이가 가변하고, 임의의 응력에서 탄성력이 최소화되어 절단이 이루어져 튕겨나가며, 땅속에 묻혀있어도 장기간 녹슬지 않는 재질로 이루어진 모든 탄성 기기를 이용할 수 있다.
방수관(120)은 방수가 이루어지는 재질로 구현하는 것이 바람직하나, 유리섬유·탄소 섬유·케블라 등의 방향족 나일론섬유와 불포화 폴리에스테르·에폭시수지 등의 열경화성 수지가 결합하며, 철보다 강하고 알루미늄보다 가벼우며 녹슬지 않고 가공하기 쉽다는 것이 장점이 있는 유리섬유강화플라스틱(FRP: glass fiber reinforced plastic)를 사용하는 것이 더욱 바람직하다.
상기 사면붕괴 측정장치(100)는, 상기 복수의 압축 스프링(111 ~ 113)이 외부의 이물질에 의해 변성되는 것을 방지하기 위해 외부 이물질의 유입을 방지하는 커버(141)(142)를 포함한다.
변위정보 전송수단(150)은 각 변위 측정 센서(130)로부터 획득한 변위 정보를 신호 처리하여 접속되는 단말기로 유선 전송하기 위한 신호처리기를 포함할 수 있으며, 신호처리기에서 처리된 변위 정보를 무선 통신망으로 전송하기 위한 무선 데이터 포맷으로 변환하여 무선 전송하는 무선 전송기를 포함할 수도 있다. 유선 단말기를 이용할 경우 관리자가 신호처리기에 접근하여 직접 단말기를 신호처리기에 접속해야하는 불편함이 있으며, 무선 전송 방식일 경우에는 관리자가 직접 신호처리기의 위치까지 답사하는 불편함을 제거할 수 있는 장점이 있다. 그러나 유선 방식과 무선 방식은 주지한 바와 같이 다양한 장·단점이 존재하므로, 하천제방 시설물의 설치 위치나 기타 주변 환경을 고려하여 적절한 방식을 취사선택하거나 두 가지 방식을 병행하여 사용하는 것도 가능하다.
상태 예측수단(20)은 현장에 설치된 상기 현장 시스템(10)에서 전송된 거동데이터를 분석하여 하천제방의 상태를 계측하는 기능을 수행한다.
이러한 상태 예측수단(20)은, 상기 변위정보 전송수단(150)에서 전송된 복수 타깃의 변위 정보를 수신하여 처리하는 호스트 컴퓨터(210)와; 상기 호스트 컴퓨터(210)에서 처리한 변위 정보를 데이터베이스(221)에 저장하고 이를 관리하기 위한 데이터베이스 서버(220)와; 상기 데이터베이스(221)에 저장된 복수 타깃의 변위 정보를 그룹화하고, 그 그룹화한 타깃의 변위 정보를 SAM 및 MDM 알고리즘으로 분석하여 사면 거동을 판단하고, 그 판단 결과 거동 가능한 그룹이 존재할 경우 승법 분해를 수행하여 미래 상황을 예측하며, 그 결과를 하천제방의 미래 상황 예측 정보로 제공하는 예측수단(230)과; 상기 예측수단(230)에서 예측한 하천제방의 상태 정보를 통신망을 통해 일반 사용자 및 정보 참조자의 단말기로 전송해주는 예측정보 전송수단(240)을 포함한다. 여기서 예측수단(230)은 통상의 정보처리용 컴퓨터를 의미하므로, 그에 대한 구체적인 설명은 생략한다.
감시 시스템(400)은 상기 통신망(300)을 통해 변위정보 분석정보를 수신하여 사용하기 위한 일반 사용자(비 전문가) 시스템(410)과, 상기 통신망(300)을 통해 변위정보 분석정보를 수신하여 하천 호안 시설물의 대규모 손상이나 붕괴를 통보해주거나 대규모 손상이나 붕괴로부터 피해를 방지하기 위한 대책을 수립하는 분석 전문가인 정보 참조자 시스템(420)을 포함한다.
여기서 상태 예측수단(20)은 상기 현장에 설치된 현장 시스템(100)과 마찬가지로 현장에 설치할 수도 있으며, 원격에 설치하고 별도의 통신망(예를 들어, CDMA)을 통해 데이터 통신을 수행하는 방식을 이용할 수도 있다. 관리의 편리성이나 기타 주변 환경 등을 고려할 때 원격에 설치하는 것이 바람직하다.
이와 같이 구성된 본 발명에 따른 하천제방 안전관리를 위한 제방붕괴예측시스템은, 도 2와 같이 사면붕괴 측정장치(100)를 계측하고자 하는 하천 호안(제방) 등의 지중에 매설하게 된다. 여기서 매설 깊이는 기존 지중에 매설되는 센서와 같은 깊이로 매설하는 것이 바람직하다.
사면붕괴 측정장치(100)가 지중에 매설된 상태에서 동작을 시작하면, 방수관(120)에 내장된 측정센서(130)가 동작을 하여 자신의 위치 정보(변위 정보)를 실시간으로 출력을 하게 된다. 여기서 측정센서(130)는 도 3에 도시한 바와 같이, 방수관(120)에 의해 외부의 이물질 접촉이 차단되어 항상 정상적인 동작을 할 수 있다. 측정센서(130)는 자신의 위치 정보를 전송할 수 있는 GPS모듈, 무선신호를 송출할 수 있는 RFID, 변경된 위치 정보를 전송할 수 있는 변위 센서 등으로 구현할 수 있고, 측정센서가 지중에 매설된 경우 측정센서의 측정신호를 외부로 전송할 수 있는 별도의 정보전송장치를 결합하여 사용할 수도 있다.
사면붕괴 측정장치(100)의 센서 변위 수단(110)은 탄성력을 가지는 복수의 압축 스프링(111 ~ 113) 등으로 이루어지며, 그 중 압축 스프링(113)은 후속하는 측정센서(230)를 내장하기 위한 방수관(220)과 연결되며, 압축 스프링(111)은 상부 커버(141)에 연결되고, 압축 스프링(112)은 하부 커버(142)에 연결된다.
각각의 압축 스프링은 응력에 따라 탄성력이 가변하여 변위가 이루어지도록 구현하게 되는 데, 여기서 응력과 탄성력과의 관계는 실험에 의해 최적 치로 설정하는 것이 바람직하다. 이것은 지중의 상태 즉, 매설물이 흙인지 모래인지 등의 성분에 따라 응력이 달라지므로, 이러한 매설물의 종류에 따라 실험을 통해 최적 치를 산출하여 적용하는 것이 바람직하다.
바람직하게 응력이 감소하게 되면 각각의 압축 스프링의 탄성력이 감소하여 앞쪽의 방수관과 뒤쪽의 방수관이 점점 벌어지도록, 아울러 커버와 방수관이 점점 벌어지게 되도록 설계하며, 어느 정도의 응력, 예를 들어 사면붕괴가 60% 정도 발생하였을 경우에 해당하는 응력 감소시에는 압축 스프링의 탄성력 최소화로 압축 스프링이 끊어져 상·하부 커버와 방수관이 분리되고, 앞쪽의 방수관과 뒤쪽의 방수관이 분리되도록 한다.
여기서 상·하부 커버와 연결된 압축 스프링이 먼저 끊어지고 가운데 위치한 압축 스프링이 나중에 끊어지도록 구현하는 것이 바람직하며, 이러한 방법 이외에 동시에 3개의 압축 스프링이 전부 절단되도록 할 수도 있다.
이렇게 응력에 따라 압축 스프링의 탄성력이 저하되도록 설계를 하고, 지중에 매설하면, 사면붕괴가 이루어질 경우 응력이 점점 감소하게 되며, 이에 따라 압축 스프링의 탄성력이 줄어들어 사면붕괴가 시작되는 쪽에 위치한 방수관(120)이 후속 방수관(220)으로부터 점점 벌어지게 되고, 따라서 내장된 측정센서(130)의 변위가 사면붕괴에 따라 다르게 된다.
이렇게 사면붕괴 측정장치(100)로부터 송출된 계측정보(변위정보)는 실시간으로 변위정보 전송수단(150)을 통해 상태 예측수단(20)에 전송된다.
상태 예측수단(20)에서는 이러한 변위 정보를 이용하여 사면붕괴가 시작하고 있는지, 아니면 사면붕괴가 어느 정도 이루어진 상태(예를 들어, 60% 붕괴)인지 등을 판단하게 된다.
이를 좀 더 구체적으로 설명하면, 측정된 변위 정보는 각 단면의 센서위치 정보와 단면 번호 정보를 포함하며, 이렇게 얻어지는 변위 정보는 대표 단면을 설정하여 하천제방의 파괴라인 추정 즉, 깊은 파괴면 인지 아니면 얕은 파괴면으로 붕괴가 될 것인지 예측하는 데 사용된다. 분석하고자 하는 단면에서 발생 가능한 잠재적인 파괴 라인은 센서 수에 따라 비례하며, 모든 경우의 수를 포함하여 각 단계별로 분석을 수행한다.
상기에서 현장 시스템(100)은 필요에 따라 유선 전송 방식과 무선 전송 방식을 병행하여 구현하는 것도 가능하다.
하천 제방의 상태 예측수단(20)은 호스트 컴퓨터(210)에서 상기 변위정보 전송수단(150)에서 전송된 변위 정보를 수신하여 처리하게 되는 데, 호스트 컴퓨터(210)는 변위정보 전송수단(150)에 대응하는 기술적 구성을 갖고 변위 정보를 수신하게 된다. 예를 들어, 전송수단이 유선 방식(RS-485)일 경우 수신수단도 유선 방식의 구성을 갖게 되며, 전송수단이 무선 방식일 경우에는 수신수단도 무선 방식의 구성을 갖게 된다.
호스트 컴퓨터(210)에서 하천제방에 설치된 복수 타깃의 변위 정보를 수신하게 되면, 데이터베이스 서버(220)에서 이를 데이터베이스(221)에 저장하게 되고, 예측수단(230)은 상기 데이터베이스(221)에 저장된 복수 타깃의 변위 정보를 기 설정된 SAM 알고리즘 및 승법 분해 모델(MDM: Multiplicative Decomposition Model)을 적용하여 하천제방의 상태를 예측하게 된다.
이를 위해 근접성 분석과 상대변위 차 분석, 그리고 상기 근접성 분석과 상대변위 분석 결과를 이용한 유사성을 분석한다. 먼저, 상기 타깃의 변위 정보를 기반으로 근접성(Closeness)을 분석하게 된다. 여기서 근접성은 관측점 사이의 시간 열 데이터의 차이, 즉, 각 관측점의 시간 열 데이터의 상호 관련 정도를 의미한다.
사면에 설치된 각 관측점의 시간 열 데이터의 상호 관련 정도를 측정하기 위하여, 각 관측점의 변위 데이터를 비교하여, 관측점들의 연동성, 즉, 거동양상을 분석하여 거동하는 사면의 크기와 속도 등을 확인하게 된다.
계측대상사면 표면에 관측점 혹은 폴이 N개 있다고 가정을 하고, 각각의 관측점을 P1, P2,...,Pn이라고 한다. 데이터를 측정하는 시간분포를 t1, t2,...,tn이라 하고, 시간에 따른 각 관측점의 위치변화, 즉, 사면 상부의 변위를 D라고 한다. 그러므로 시간에 따른 Pi점의 변위는 Di(t1),Di(t2),...,Di(tn)으로 표현하거나, 간단히 Di(tn), k=1, 2,...,tn으로 표현할 수 있다.
같은 방법으로 Pi점의 경우, 시간 열에 따른 변위는 Dj(t1),Dj(t2),...,Dj(tn)으로 표현된다.
시간이 tk일 때, 두 관측점 Pi와 Pj 사이의 변위 차는 아래의 수식1과 같이 정의한다.
Figure pat00001
두 관측점 Pi와 Pj사이의 시간영역이 [t1, t2]일 때, 변위에 관한 시간 열 데이터, Di,j(tk)이 평균치는 아래의 수식 2와 같이 표현할 수 있다.
Figure pat00002
그러나 변위차 Di(tk)는 영역이 [-∞ ~ ∞]사이인 거리 차원이 있다.
사면의 붕괴가 사면 하단으로 붕괴가 일어나는 얕은 붕괴 같은 경우, 변위가 +∞방향으로 발생할 것이며, 깊은 붕괴와 같이 붕괴 면이 회전할 경우, 변위는 -∞방향으로 움직일 것이다. 여기서 ∞의 공학적 의미는 사면붕괴에 해당하며, +는 두 관측점 간의 변위의 증가, -는 변위의 감소를 의미한다. 사면에서 발생한 변위 차는 추후 전개될 수식 및 그래프 상에서 비교분석되어야 하므로, 변위 차를 무 차원의 간단한 상수로 바꾸어 계산하는 것이 좋을 것이다. 그러므로 변위 차가 [-∞ ~ ∞]사이에서 움직이고 있으므로, [-90° ~ 90°] 사이에서 [-∞ ~ ∞]로 움직이는 tan함수의 특성을 이용할 경우, 함수를 무 차원 상수로 변환할 수 있다. 변환된 무 차원 상수는 [-1 ~ 1] 구간에서 움직일 것이다. 이 변환방법은 아래의 수식3과 같다.
Figure pat00003
혹은
Figure pat00004
이다.
이 식은 다시 아래의 수식 4와 같이 다시 정리할 수 있다.
Figure pat00005
여기서
Figure pat00006
의 영역은 -900≤
Figure pat00007
≤900이므로,
Figure pat00008
는 [-1 ~ 1] 구간 내의 상수 값으로 계산된다. 그러므로
Figure pat00009
는 두 사면 관측점의 계측변위 데이터의 차이를 거리 차원이 아닌 무 차원 상수로 나타낼 수 있는 값이다.
그러나 사면이 얼마나 같이 연동하는가 즉, 사면의 두 관측점이 같이 움직이는가를 알면 붕괴 면의 크기를 알 수 있을 것이다. 다시 말하면, 사면 내 두 관측점의 같은 시간대 변위 차가 크면, 두 지점은 서로 다른 붕괴면 위에서 따로 움직이고 있는 것이고, 변위 차가 적거나 거의 없다면 두 관측점은 같은 붕괴면 상에서 움직이고 있으므로, 붕괴 면의 크기와 속도 등을 유추할 수 있다.
그러므로 두 관측점의 변위 데이터의 차이, 즉 관측점의 근접성에 관한 지수를 새로 도입할 필요가 있다. 근접성 지수는 전적으로 두 관측점의 변위 차와 다른 개념이며, 근접성 지수는 상수로 정리된 변위 차에서 최대값을 빼면 된다. 사인함수로 표시된 변위의 최대값의 경우, +변위의 경우 1이며, -변위의 경우 -1이다. 따라서 근접성 지수는 무 차원 상수로 정의된 변위차, 즉,
Figure pat00010
를 최대값 1에서 빼주면된다.
시간차 i, j에서 발생한 사면 내 두 관측점 사이의 근접성 지수는,
근접성 지수(
Figure pat00011
)=발생가능한 최대변위 차(-∞ ~ ∞) - 실제 변위차(
Figure pat00012
) = 1(
Figure pat00013
의 양의 최대변위 값) -
Figure pat00014
(실제변위 값)
혹은, = -1 - (-
Figure pat00015
)(음의 실제변위 값)
그러므로 붕괴 사면의 회전으로 인한 두 관측점 간의 상대적 변위가 감소하는 깊은 붕괴의 경우, 근접성 지수는 음의 값을 가지게 되며, 붕괴 사면의 병진으로 상대적 변위가 증가하는 얕은 붕괴의 경우, 양의 지수 값을 가지게 된다.
결론적으로 근접성 지수의 절대값이 크면 두 관측점의 시간 열 데이터의 차이가 작다는 것이고, 이는 두 지점의 변위량 차이가 작다는 것을 의미한다. 즉,
Figure pat00016
근접성 지수의 절대값이 작을수록 두 관측점의 변위 차가 작으므로 거동 양상은 비슷하다. 즉, 지수가 1 혹은 -1에 근접할수록 사면의 거동 양상은 비슷하며, 이는 같은 파괴면 상에 있다. 변위량이 지속적으로 증가하고 있는 일반적인 사면 거동에 대한 근접성 지수는 아래의 수식 5에 보인 바와 같으며, 만약 현저한 깊은 파괴로 인한 상대변위 감소 구간에서는 아래의 수식 5에 대해 부호만 다르게 정의된다.
Figure pat00017
다음으로, 타깃의 변위 정보를 기반으로 상대변위 차를 분석하게 된다. 여기서 상대변위 차는 데이터 시리즈의 그래프 형상의 비슷한 정도를 의미한다.
사면의 붕괴에 따른 거동양상은 긴 시간에 걸쳐 일어나므로, 두 관측점 사이의 변위 차가 시간에 따라 다른 거동을 보이는 형상을 찾는 것이 유리하다. 그러나 근접성 지수의 경우 각 관측점 사이의 전체 시간 열 데이터를 단순히 하나의 지수로 나타내므로, 각 관측점 사이의 거동차이가 일어나는 과정을 분석하기가 힘들다. 즉, 전체 거동차이를 하나의 값으로 나타내는 근접성 지수보다 과정을 알 수 있는 그래프를 이용하는 것이 유리하다. 그러므로 이러한 요구에 합당하게 유사성이라는 개념을 도입하고자 한다.
다시 말하면, 근접성 지수는 평균치라는 개념을 이용하고 있으므로, 안정적인 데이터와 불안정한 데이터 모두 평균치는 같을 수 있으나, 데이터를 그래프로 그려보면 유사한 형상을 보이는 그래프의 데이터가 같은 사면 거동양상을 보임을 쉽게 알 수 있다.
유사성의 개념을 이용하기 위해 먼저 상대 변위 차(relative displacement)를 정의한다. 상대변위 차는 두 관측점 사이의 최대변위차(
Figure pat00018
)와 최소변위차(
Figure pat00019
)와 해당 시점의 변위차를 사용하여 구할 수 있다. 두 관측점의 최대변위 차와 최소변위 차의 차이가 변위차 차이의 최대값이 될 것이며, 이는 두 관측 데이터를 이용한 변위차 그래프의 최대간격(
Figure pat00020
-
Figure pat00021
)이 될 것이다. 그리고 최대 변위차값을 해당 시점의 변위차값으로 빼면 이는 변위차 그래프의 해당 시점의 간격이 되며(
Figure pat00022
-
Figure pat00023
), 해당 간격을 최대 간격으로 나누면 이는 상대 변위차가 될 것이다. 먼저, 최대 및 최소 상대변위 차는 아래의 수식 6 및 수식 7로 정의되며, 이는 도 5와 같이 나타낼 수 있다.
Figure pat00024
Figure pat00025
Figure pat00026
두 관측점의 관측시간 i, j구간에서 어떤 시간 tk에 해당하는 시간 열 데이터에 관한 상대변위 차는 아래의 수식 8로 정의할 수 있다.
Figure pat00027
Figure pat00028
두 시간 열 데이터의 평균값을 구하면, 아래의 수식 9로 정의할 수 있다.
Figure pat00029
Figure pat00030
의 범위는 [0 ~ 1]인 무차원수이므로, 당연히 이의 평균값인
Figure pat00031
역시 범위가 [0 ~ 1] 사이인 무차원수이다.
Figure pat00032
값이 크거나 혹은 1에 근접할 때 1, i, j시간 구간의 데이터 시리즈는 유사한 모양을 갖게 될 것이다. 다시 말하면,
Figure pat00033
값이 작거나 0에 근접하면 i, j시간 구간의 데이터 시리즈는 서로 다른 모양을 갖게 될 것이다.
다음으로, 상기 근접성 분석과 상태변위 차 분석 결과를 이용하여 유사성 분석을 하고, 그 유사성 분석 결과를 기반으로 사면 거동을 판단하게 된다.
두 관측점의 시간 열 데이터의 차이, 즉, 두 지점의 변위량 차이를 수치화하기 위하여 근접성 지수(
Figure pat00034
)라는 개념을 정의하였고, 데이터 시리즈의 그래프 형상의 비슷한 정도를 나타내기 위하여 상대 변위차(
Figure pat00035
)라는 개념을 정의하였다. 상기 두 가지 개념을 각각 사용하는 것이 아니라 함께 사용한다면, 두 관측지점의 시간 열 데이터의 변위 차의 정도 및 그래프로 나타난 모양이 비슷한 정도(각 지점의 시간에 따른 데이터의 증감 정도)를 동시에 나타낼 수 있을 것이다.
두 관측점의 시간 열 i와 j에서 만들어진 근접성 지수와 상대 변위차의 곱을 이용하여 유사성을 정의한다. 이 유사성(
Figure pat00036
)은 아래의 수식 10과 같이 정의한다. 구간이 [0 ~ 1] 사이인 두 무차원수로 구성되었으므로, 유사성은 당연히 [0 ~ 1]사이의 범위 구간에 존재한다.
Figure pat00037
주지한 바와 같이, 근접성, 상대변위 차 및 유사성이라는 개념을 이용함으로써, 사면 붕괴의 위험성이 있는 다양한 실 계측 데이터들을 해석, 가공한 후 실시간 사면 거동 양상을 분석하여 사면안정에 대한 결론을 유도할 수 있게 된다.
상기 유사성을 분석한 결과 거동 가능한 사면이 존재할 경우, 승법 분해 모델(MDM)을 이용하여 향후 데이터를 예측하게 된다.
시계열분석(Time-Series Analysis)은 과거 데이터를 축적하여 이후 상황을 예측하기 위한 가장 훌륭한 방법이다. 이는 과거에 발생했던 사실을 바탕으로 미래에 일어날 가능성을 예측하는 것이다. 만약 구조물의 추후 파괴가능성을 미리 예측하고 싶다면 과거의 주기적인 데이터를 이용하여 능동적으로 예측할 수 있다.
시계열 모델에는 이동평균(Moving Averages), 이동가중평균법(Weighted Moving Averages), 지수평활법(Exponential Smoothing), 선형추세분석(Linear Trend Analysis), 계절변동분석(Seasonality Analysis), 그리고 승법분해(Multiplicative Decomposition) 등의 방법이 있다.
통계적 시계열 분석은 일반적으로 크게 2가지 방법이 있다. 그 중 가장 폭넓게 사용하는 방법으로 승법분해가 있으며, 본 발명에서도 이를 이용하여 미래 데이터를 예측하게 된다. 승법 분해는 4가지 즉, 추세변동(Trend), 계절변동(Seasonality), 순환변동(Cycles), 불규칙 변동(Random)으로 구성되어 있으며, 이는 아래의 수식 11과 같다.
Figure pat00038
승법분해모델(MDM)은 크게 2가지 성분으로 구분할 수 있는 데, 계절적인 영향 요소와 추세 및 주기의 조합으로 구성된다. 불규칙변동 요소는 우연적이며 급격한 변화에 의한 비 반복적인 요인으로서 분별할 수 있는 가능성이 없다. 다시 말하면 승법분해모델은 추세변동과 주기적인 계절변동으로 승법분해되는 분석방법이다.
추세변동분석(Trend Analysis)은 과거 측정된 데이터들로부터 추세변동을 알 수 있는 방정식이라고 할 수 있으며, 추세변동분석 그래프를 통해 중장기적인 미래의 데이터를 예측할 수 있다. 이를 위해 본 발명에서는 회귀분석법을 사용하였으며, 이는 회귀방정식으로 표현되며 아래의 수식 12에 나타내었다.
Figure pat00039
Figure pat00040
Figure pat00041
예측된 데이터의 정확도는 실제 혹은 계측된 데이터와 예측된 데이터의 비교에 의해 결정할 수 있다. 만약 시간에 따라 예측된 값을 Ft라고 하고 실제 값을 At라 하면 예측에러(Forecast error)는 아래의 수식 13과 같이 정의된다.
Figure pat00042
다음으로, 데이터 분석을 수행한다.
여기서 승법분해모델(MDM) 분석은 유사성 분석결과 사면 중 위험 가능성 사면으로 분류한 센서번호를 우선적으로 분석한다. MDM분석의 목적은 실제 값과 예측 값의 차이가 거의 없는 상황이라면 추후 계측될 데이터의 값을 추정할 수 있다는 것이다. MDM분석시 구할 수 있는 MAD, MSE, MAPE는 예측에러 정도를 파악할 수 있다.
상기 예측에러 정도를 파악한 결과, 예측에러가 적은 경우에는 추후 데이터를 예측할 가능성이 크다는 것을 의미하고, 그 분석 결과를 감시 시스템으로 전송하게 된다.
그러나 예측에러 정도를 파악한 결과 예측에러가 크다는 것은 두 그래프가 거의 일치하지 않는다는 의미이지만, 예측 에러의 원인규명 노력을 통해 사면 거동에 관한 또 다른 정보를 연구할 수 있는 가능성을 줄 수 있다. 다시 말하면 주기적인 요인이 아니라 비 예측적인(붕괴 인자 등) 요인이 지배적이라는 의미가 되므로 현장조사 등을 통해 데이터의 적합성 유무를 판단한 후 적절한 조치를 취할 수가 있다.
이렇게 판단한 하천제방의 안전관리를 위한 분석정보는 인터넷과 같은 통신망을 이용하거나 문자메시지(SMS) 같은 문자 전송 방식을 이용하여 관리자(하천제방을 관리하는 전문 관리자) 및 일반 사용자(해당 정보를 이용하는 사용자로서 비 전문가)에게 전송된다.
관리자는 전송된 하천제방의 안전관리를 위한 분석 정보를 확인하고, 하천제방이 대규모 손상이 발생할 조짐이 있거나 붕괴 우려가 있다고 판단이 되면, 경보를 발령하고, 후속 조치를 신속하게 취함으로써, 추후 하천제방의 대규모 손상이나 붕괴로 인해 발생할 수 있는 피해를 최소화한다.
본 발명은 상술한 특정의 바람직한 실시 예에 한정되지 아니하며, 청구범위에서 청구하는 본 발명의 요지를 벗어남이 없이 당해 발명이 속하는 기술분야에서 통상의 지식을 가진 자라면 누구든지 다양한 변형실시가 가능한 것은 물론이고, 그와 같은 변경은 청구범위 기재의 범위 내에 있게 된다.
10… 현장 시스템
100… 사면붕괴 측정장치
110… 센서 변위 수단
150… 변위정보 전송수단
20… 상태 예측수단
210… 호스트 컴퓨터
230… 예측수단
240… 예측정보 전송수단

Claims (5)

  1. 하천제방 안전관리를 위한 제방붕괴 예측시스템에 있어서,
    지중에 매설되며, 응력에 따라 탄성력이 가변하여 변위를 발생하고, 상기 응력이 정해진 기준 이상 변경될 경우 내장된 측정센서를 정해진 위치로부터 이탈시켜 거동 데이터를 발생하는 사면붕괴 측정장치를 포함하는 현장 시스템과;
    현장에 설치된 상기 현장 시스템에서 전송된 거동데이터를 분석하여 하천제방의 상태를 계측하는 상태 예측수단을 포함하는 것을 특징으로 하는 하천제방 안전관리를 위한 제방붕괴 예측시스템.
  2. 제1항에 있어서, 상기 사면붕괴 측정장치는,
    응력에 따라 탄성력이 가변하여 변위를 발생하고, 상기 응력이 정해진 기준 이상 변경될 경우 내장된 측정센서를 정해진 위치로부터 이탈시켜 변위 정보를 발생하는 센서 변위 수단을 포함하는 것을 특징으로 하는 하천제방 안전관리를 위한 제방붕괴 예측시스템.
  3. 제1항에 있어서, 상기 현장 시스템은,
    상기 센서 변위 수단에서 획득한 변위 정보를 연결된 통신망에 대응하는 데이터 포맷으로 변환하여 거동 데이터로 상기 원격에 위치한 상태 예측수단에 전송하는 변위정보 전송수단을 포함하는 것을 특징으로 하는 하천제방 안전관리를 위한 제방붕괴 예측시스템.
  4. 제2항에 있어서, 상기 사면붕괴 측정장치는,
    내장된 측정센서에 습기가 유입되는 것을 방지하는 방수관을 더 포함하고,
    상기 센서 변위 수단은,
    응력에 따라 탄성력이 가변하여 상기 측정센서를 정해진 위치로부터 이탈시키는 복수의 압축 스프링을 구비한 것을 특징으로 하는 하천제방 안전관리를 위한 제방붕괴 예측시스템.
  5. 제1항에 있어서, 상기 상태 예측수단은,
    상기 변위정보 전송수단에서 전송된 복수 타깃의 변위 정보를 수신하여 처리하는 호스트 컴퓨터와;
    상기 호스트 컴퓨터에서 처리한 변위 정보를 데이터베이스에 저장하고 이를 관리하기 위한 데이터베이스 서버와;
    상기 데이터베이스에 저장된 복수 타깃의 변위 정보를 그룹화하고, 그 그룹화한 타깃의 변위 정보를 SAM 및 MDM 알고리즘으로 분석하여 사면 거동을 판단하고, 그 판단 결과 거동 가능한 그룹이 존재할 경우 승법 분해를 수행하여 미래 상황을 예측하며, 그 결과를 하천 호안 시설물의 미래 상황 예측 정보로 제공하는 예측수단과;
    상기 예측수단에서 예측한 하천 호안 시설물의 상태 정보를 통신망을 통해 일반 사용자 및 정보 참조자의 단말기로 전송해주는 예측정보 전송수단을 포함하는 것을 특징으로 하는 하천제방 안전관리를 위한 제방붕괴 예측시스템.
KR20100043949A 2010-05-11 2010-05-11 하천제방 안전관리를 위한 제방붕괴 예측시스템 KR101178557B1 (ko)

Priority Applications (1)

Application Number Priority Date Filing Date Title
KR20100043949A KR101178557B1 (ko) 2010-05-11 2010-05-11 하천제방 안전관리를 위한 제방붕괴 예측시스템

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
KR20100043949A KR101178557B1 (ko) 2010-05-11 2010-05-11 하천제방 안전관리를 위한 제방붕괴 예측시스템

Publications (2)

Publication Number Publication Date
KR20110124525A true KR20110124525A (ko) 2011-11-17
KR101178557B1 KR101178557B1 (ko) 2012-08-30

Family

ID=45394257

Family Applications (1)

Application Number Title Priority Date Filing Date
KR20100043949A KR101178557B1 (ko) 2010-05-11 2010-05-11 하천제방 안전관리를 위한 제방붕괴 예측시스템

Country Status (1)

Country Link
KR (1) KR101178557B1 (ko)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20150129487A (ko) * 2014-05-12 2015-11-20 한국건설기술연구원 제방의 위험도 평가 방법 및 그에 적합한 평가 시스템
KR102615082B1 (ko) * 2022-11-29 2023-12-19 라온구조안전기술(주) 저수지 및 하천 제방 관리를 위한 매립형 복합 IoT 센서 장치를 이용한 모니터링 시스템 및 이의 실행 방법
CN117349602A (zh) * 2023-12-06 2024-01-05 江西省水投江河信息技术有限公司 一种水利设施运行状态预测方法、系统及计算机

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101796078B1 (ko) 2017-06-22 2017-11-09 주식회사 유신 항만시설용 침하 및 재해 감시 설비

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR100755474B1 (ko) 2005-12-22 2007-09-04 (주)지엠지 센서와 연동되는 붕괴 위험 위치 영상 획득 시스템

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20150129487A (ko) * 2014-05-12 2015-11-20 한국건설기술연구원 제방의 위험도 평가 방법 및 그에 적합한 평가 시스템
KR102615082B1 (ko) * 2022-11-29 2023-12-19 라온구조안전기술(주) 저수지 및 하천 제방 관리를 위한 매립형 복합 IoT 센서 장치를 이용한 모니터링 시스템 및 이의 실행 방법
CN117349602A (zh) * 2023-12-06 2024-01-05 江西省水投江河信息技术有限公司 一种水利设施运行状态预测方法、系统及计算机

Also Published As

Publication number Publication date
KR101178557B1 (ko) 2012-08-30

Similar Documents

Publication Publication Date Title
Manfreda et al. A digital elevation model based method for a rapid estimation of flood inundation depth
Garousi‐Nejad et al. Terrain analysis enhancements to the height above nearest drainage flood inundation mapping method
Wang et al. Modeling tunneling-induced ground surface settlement development using a wavelet smooth relevance vector machine
Jalbert et al. Temporal uncertainty estimation of discharges from rating curves using a variographic analysis
Wu et al. Modeling risk analysis for forecasting peak discharge during flooding prevention and warning operation
Horritt et al. Modelling the hydraulics of the Carlisle 2005 flood event
Titova et al. Lifetime of earth dams
JP7347918B2 (ja) 水位予測方法、水位予測プログラム及び水位予測装置
De Risi et al. Bayesian tsunami fragility modeling considering input data uncertainty
JP2015004245A (ja) 浸水予測システム、浸水予測方法およびプログラム
KR101178557B1 (ko) 하천제방 안전관리를 위한 제방붕괴 예측시스템
Sharafati et al. Assessment of dam overtopping reliability using SUFI based overtopping threshold curve
Liu et al. Current profile analysis and extreme value prediction in the LH11-1 oil field of the South China Sea based on prototype monitoring
Yang et al. An indirect approach for discharge estimation: a combination among micro-genetic algorithm, hydraulic model, and in situ measurement
Prieto et al. Preface: Advances in flood risk assessment and management
Pulido et al. Study of the spatiotemporal correlation between sediment-related blockage events in the sewer system in Bogotá (Colombia)
Shastry et al. Using Steady‐State Backwater Analysis to Predict Inundated Area from National Water Model Streamflow Simulations
Tang et al. An optimized weighted Naïve Bayes method for flood risk assessment
Lee et al. Construction of rating curve at high water level considering rainfall effect in a tidal river
Romstad et al. A GIS method for assessment of rock slide tsunami hazard in all Norwegian lakes and reservoirs
KR101056056B1 (ko) 하천 호안 시설물의 안전관리를 위한 분석방법
KR101056057B1 (ko) 해안 호안 시설물의 상태 계측방법
Kim et al. Probabilistic flood hazard assessment method considering local intense precipitation at NPP sites
Ling et al. Improved Kalman filter method considering multiple factors and its application in landslide prediction
KR20100113404A (ko) 하천 호안 시설물의 실시간 사면붕괴 예측방법

Legal Events

Date Code Title Description
A201 Request for examination
N231 Notification of change of applicant
E902 Notification of reason for refusal
E701 Decision to grant or registration of patent right
GRNT Written decision to grant
FPAY Annual fee payment

Payment date: 20160728

Year of fee payment: 5

FPAY Annual fee payment

Payment date: 20170727

Year of fee payment: 6

FPAY Annual fee payment

Payment date: 20180730

Year of fee payment: 7

FPAY Annual fee payment

Payment date: 20190723

Year of fee payment: 8