지진은 지각내 단층운동과 함께 수십억 년 전부터 발생하여 왔다. 1960년대부터는 미국에서는 지진과 지하수위 의 연관성 연구를 본격적으로 시작하였으나, 국내에서는 2010년 경부터 지진과 지하수위 및 수리지화학적 연관성에 대 한 연구를 시작하였다. 본 연구에서는 국내학자들이 과거부터 2021년까지 연도별로 지진과 지하수의 관련성을 연구한 논문을 Web of Science에서 검색하고, 분야별(지하수위, 수리지화학, 지하수위와 수리지화학 병행, 그 외 분야) 연구 특 성을 검토하였다. 국내학술지에 게재된 지진과 지하수 관련성 연구 논문을 보면, 연도별 논문 편수는 2011년에 동일본 대지진, 2016년에 경주지진, 2017년에 포항지진과 발생과 관련되며, 이에 따라 2011년, 2018년, 2019년, 2020년에 국내 및 국제학술지 게재 논문수가 증가하였다. 대부분의 지하수위와 지진의 관련성 연구는 지진과 동시기의 지하수위 변화 에 관한 연구이며, 지진 전조와 관련한 연구는 거의 없다. 지하수위 관측자료와 함께 여러가지 수리지화학적 정보와 미 생물은 지진에 의한 기반암내 지하수의 유동과 화학적인 반응을 보다 상세하게 이해하는데 도움을 줄 수 있다. 지진 감시 및 예측을 위해서는 지진감시를 위한 지하수관측공 네트워크를 전국적으로 구축할 필요가 있다.

상⋅하수도 관로 시설은 대표적인 국가기반 시설물로서 국민생활에 있어서 매우 중요한 라이프 라인이다. 지진으로 인한 관로시설의 손상은 급수의 차단을 유발하여 심각한 피해를 초래할 가능성이 매우 크다. 그러므로 상⋅하수도 관로 시설은 지진으로부터 반듯이 안전하게 보호되어야 할 필요가 있다. 지진 및 지반침하로 인한 설계변위를 초과하는 상대변위는 구조물과 연결되는 배관의 이음부에서 손상을 발생시킨다. 벨로우즈형 신축이음관은 온도차에 의한 배관의 팽창 및 변형을 흡수하고 기계진동에 의한 배관의 손상을 막기 위한 장치이다. 본 연구에서는 매설된 상⋅하수도 관로를 보호하고자 지진변위와 지반침하 대응을 목적으로 벨로우즈를 적용하였다. 적층형 하이드로포밍 메탈 벨로우즈는 기존의 벨로우즈와 달리 지진과 같은 저주 기피로하중에 대한 내구성이 우수하다. 따라서 3ply 벨로우즈형 신축이음관을 대상으로 반복가력 굽힘시험을 수행하고 지진안전성과 내침하성능을 평가하였다. 그 결과, 3ply 벨로우즈형 신축이음관은 8.8°이상의 횡방향 변형각에 대응할 수 있는 것으로 나타났다.

The Gyeongju and Pohang earthquakes caused damages to many cultural properties; particularly, stone pagoda structures were significantly damaged among masonry cultural properties. To preserve these structures, it is necessary to understand their dynamic behavior characteristics under earthquakes. Analyses on such areas as deformation, frequency, maximum acceleration, permanent displacement, sliding, and rocking have to be performed. Although many analytical studies have already been conducted, dynamic behavior studies based on experiments are insufficient. Therefore, this study analyzed dynamic behavior characteristics by performing a shaking table experiment on a three-story stone pagoda structure at the Cheollongsa temple site damaged by the Gyeongju earthquake. As a result of the experiment, the displacements of stylobates did not occur significantly, but the tower body parts rotated. In particular, the rotation of the 1F main body stone was relatively larger than that of the other chief body stones because the 1F main body stone is relatively more slender than the other parts. In addition, the decorative top was identified as the component most vulnerable to sliding. This study found that the 1F main body stone is vulnerable to rocking, and the parts located on the upper part are more vulnerable to sliding.

Current seismic fragility functions for buildings were developed by defining damage state threshold based on story drift concerning foreign references and using the capacity spectrum method based on spectral displacement. In this study, insufficient details and dependence on the core location of piloti-type buildings were not reflected in the fragility function because it was developed before the Pohang earthquake. In order to develop an improved one for piloti-type buildings, several types of core were selected, damage state threshold was determined based on the capacity of structural members, and three-dimensional analyses were utilized. As a result, seismic fragility functions based on spectral acceleration were developed for various core locations and different shear strengths of the column stirrup. The fragility of piloti-type buildings significantly varied according to core location, an additional single wall, and whether the contribution of column stirrup was included or not. To estimate fragility more reasonably, it is necessary to prepare the parameters to reflect actual state well.

The stochastic method is applied to simulate strong ground motions at seismic stations of seven metropolises in South Korea, creating an earthquake scenario based on the causative fault of the 2016 Gyeongju earthquake. Input parameters are established according to what has been revealed so far for the causative fault of the Gyeongju earthquake, while the ratio of differences in response spectra between observed and simulated strong ground motions is assumed to be an adjustment factor. The calculations confirm the applicability and reproducibility of strong ground motion simulations based on the relatively small bias in response spectra between observed and simulated strong ground motions. Based on this result, strong ground motions by a scenario earthquake on the causative fault of the Gyeongju earthquake with moment magnitude 6.5 are simulated, assuming that the ratios of its fault length to width are 2:1, 3:1, and 4:1. The results are similar to those of the empirical Green’s function method. Although actual site response factors of seismic stations should be supplemented later, the simulated strong ground motions can be used as input data for developing ground motion prediction equations and input data for calculating the design response spectra of major facilities in South Korea.

The spatial variation characteristics of seismic motions at the nuclear power plant's site and structures were analyzed using earthquake records obtained at the Fukushima nuclear power plant during the Great East Japan Earthquake. The ground responses amplified as they approached the soil surface from the lower rock surface, and the amplification occurred intensively at about 50 m near the ground. Due to the soil layer's nonlinear characteristics caused by the strong seismic motion, the ground's natural frequency derived from the response spectrum ratio appeared to be smaller than that calculated from the shear wave velocity profile. The spatial variation of the peak ground acceleration at the ground surface of the power plant site showed a significant difference of about 0.6 g at the maximum. As a result of comparing the response spectrums at the basement of the structure with the design response spectrum, there was a large variability by each power plant unit. The difference was more significant in the Fukushima Daiichi site record, which showed larger peak ground acceleration at the surface. The earthquake motions input to the basement of the structure amplified according to the structure's height. The natural frequency obtained from the recorded results was lower than that indicated in the previous research. Also, the floor response spectrum change according to the location at the same height was investigated. The vertical response on the foundation surface showed a significant difference in spectral acceleration depending on the location. The amplified response in the structure showed a different variability depending on the type of structure and the target frequency.

Existing old reinforced concrete buildings could be vulnerable to earthquakes because they were constructed without satisfying seismic design and detail requirements. In current seismic design standards, the target collapse probability for a given Maximum Considered Earthquake (MCE) ground-shaking hazard is defined as 10% for ordinary buildings. This study aims to estimate the collapse probabilities of a three-story, old, reinforced concrete building designed by only considering gravity loads. Four different seismic design categories (SDC), A, B, C, and D, are considered. This study reveals that the RC building located in the SDC A region satisfies the target collapse probability. However, buildings located in SDC B, C, and D regions do not meet the target collapse probability. Since the degree of exceedance of the target probability increases with an increase in the SDC level, it is imminent to retrofit non-ductile RC buildings similar to the model building. It can be confirmed that repair and reinforcement of old reinforced concrete buildings are required.

본 연구에서는 비틀림비정형성과 수직비정형성을 가진 RC 필로티 건축물의 지진동에 대한 거동을 층강성을 적용하여 간단하게 모델링하는 선형 동적해석 프로그램을 개발하고자 한다. 개발된 동적 해석 프로그램을 적용하여 필로티 건축물의 동적 거동 및 필로티층 각 기둥의 전단력을 분석하고, 필로티층에 전단벽 또는 가새를 보강하였을 때 보강효과를 평가하고자 한다. 모서리코어가 있는 필로티 건축물에서 필로티층의 코어 반대편 모서리를 전단벽이나 K형 가새로 보강하였을 때 변위와 기둥 전단력이 크게 감소하는 것으로 나타났으며, 모서리 양면을 K형 가새로 보강하는 것보다 한 면을 전단벽으로 보강하는 것이 보강효과가 큰 것으로 나타났다.

A shake table test is conducted for the three-story reinforced concrete building structure using 0.28 g, 0.5 g, 0.75 g, and 1.0 g of seismic input motions based on the Gyeongju earthquake. Computational efforts are made in parallel to explore the mechanical details in the structure. For engineering practice, the elastic modulus of concrete and rebar in the dynamic analysis is reduced to 38% and 50%, respectively, to calibrate the structure's natural frequencies. The engineering approach to the reduced modulus of elasticity is believed to be due to the inability to specify the flexibility of the actual boundary conditions. This aspect may lead to disadvantages of nonlinear dynamic analysis that can distort local stress and strain relationships. The initial elastic modulus can be applied directly without the so-called engineering adjustment with infinite element models with spring and spring-dashpot boundary conditions. This has the advantage of imposing the system flexibility of the structure on the sub-boundary conditions of springs and damping devices to control its sensitivity in a serial arrangement. This can reflect the flexibility of realistic boundary conditions and the effects of system damping (such as the gap between a concrete footing and shake table, loosening of steel anchors, etc.) in scalar quantities. However, these spring and dashpot coefficients can only be coordinated based on experimental results, making it challenging to select the coefficients in-prior to perform an experimental test.

This paper presents a real-time, false-pick filter based on deep learning to reduce false alarms of an onsite Earthquake Early Warning (EEW) system. Most onsite EEW systems use P-wave to predict S-wave. Therefore, it is essential to properly distinguish P-waves from noises or other seismic phases to avoid false alarms. To reduce false-picks causing false alarms, this study made the EEWNet Part 1 'False-Pick Filter' model based on Convolutional Neural Network (CNN). Specifically, it modified the Pick_FP (Lomax et al.) to generate input data such as the amplitude, velocity, and displacement of three components from 2 seconds ahead and 2 seconds after the P-wave arrival following one-second time steps. This model extracts log-mel power spectrum features from this input data, then classifies P-waves and others using these features. The dataset consisted of 3,189,583 samples: 81,394 samples from event data (727 events in the Korean Peninsula, 103 teleseismic events, and 1,734 events in Taiwan) and 3,108,189 samples from continuous data (recorded by seismic stations in South Korea for 27 months from 2018 to 2020). This model was trained with 1,826,357 samples through balancing, then tested on continuous data samples of the year 2019, filtering more than 99% of strong false-picks that could trigger false alarms. This model was developed as a module for USGS Earthworm and is written in C language to operate with minimal computing resources.

Historical records of earthquakes are generally used as a basis to extrapolate the instrumental earthquake catalog in time and space during the probabilistic seismic hazard analysis (PSHA). However, the historical catalogs’ input parameters determined through historical descriptions rather than any quantitative measurements are accompanied by considerable uncertainty in PSHA. Therefore, quantitative assessment to verify the historical earthquake parameters is essential for refining the reliability of PSHA. This study presents an approach and its application to constrain reliable ranges of the magnitude and corresponding epicenter of historical earthquakes. First, ranges rather than specific values of ground motion intensities are estimated at multiple locations with distances between each other for selected historical earthquakes by reviewing observed co-seismic natural phenomena, structural damage levels, or felt areas described in their historical records. Based on specific objective criteria, this study selects only one earthquake (July 24, 1643), which is potentially one of the largest historical earthquakes. Then, ground motion simulations are performed for sufficiently broadly distributed epicenters, with a regular grid to prevent one from relying on strong assumptions. Calculated peak ground accelerations and velocities in areas with the historical descriptions on corresponding earthquakes are converted to intensities with an empirical ground motion-intensity conversion equation to compare them with historical descriptions. For the ground motion simulation, ground motion prediction equations and a frequencywavenumber method are used to consider the effects of possible source mechanisms and stress drop. From these quantitative calculations, reliable ranges of epicenters and magnitudes and the trade-off between them are inferred for the earthquake that can conservatively match the upper and lower boundaries of intensity values from historical descriptions.

본 연구에서는 2001년부터 2020년까지 지진과 관련된 지하수위, 수질, 라돈, 원격탐사, 전기비저항, 중력, 지자기 분야의 세계적으로 학술지에 게재된 논문 편수를 Web of Science에서 검색하여 그 경향성을 분석하였다. 그리고 논문 편수와 Mw 5.0 이상, Mw 6.0 이상, Mw 7.0 이상, Mw 8.0 이상, Mw 9.0 이상 지진 발생 건수를 비교 분석하였다. 지진과 관련한 중력, 라돈, 지하수(지하수위, 수질), 전기비저항, 지자기분야 논문 편수는 장기적으로 증가하는 추세를 보인다. 이는 원격탐사 기술의 발달, 측정 장비의 고도화, 빅데이터 분석 등을 통한 종합적인 자료 해석이 가능해지면 서 여러 분야에서 지진 전조 및 지진 현상 연구가 활발해지고 있기 때문이다. Mann-Kendall과 Sen 추세 검정에 의하면, 중력 관련 논문의 경우 1.30편/년의 증가추세를 보이고, 라돈 0.60편/년, 지하수 0.70편/년, 전기비저항 0.25편/년, 원격탐사 0.67편/년의 증가추세를 보인다. Mw 5.0 이상, Mw 6.0 이상, Mw 7.0 이상, Mw 8.0 이상, Mw 9.0 이상의 지진 발생 건수와 경향성을 제거한 분야별 논문 편수 간의 교차상관분석에 의하면, 라돈과 원격탐사 분야의 교차상관성이 높으며, 지연시간은 1년이다. 또한 2004년과 2005년 수마트라 지진, 2008년 쓰촨성 지진, 2010년 아이티 지진, 2010 칠레 지진 등의 큰 규모의 지진 발생이 논문 편수 증가와 관련되는 것으로 추정된다.

2004년 5월 29일 발생한 울진해역지진(Mw 5.1)과 관련된 대기 인프라사운드 신호가 철원(진앙 거리 321 km) 및 대전(256 km) 관측소에 기록되었다. 신호의 지속시간은 수 분 이상이며, 음원 방향을 지시하는 후방-방위각은 28 o 이상의 큰 변화를 보였다. 역-투사 방법과 신호 감쇄 보정 결과, 인프라사운드 신호는 삼척-울진-포항까지 연결되는 약 4,600 km2 면적의 지반운동으로 발생하였으며, 음원 최대 크기(BSP)는 11.1 Pa로 계산되었다. 이 결과는 최대지반가속 도(PGA) 자료로 계산한 음원 최대 크기(PSP)와도 부합하고 있으며, 지진 발생 당시 인프라사운드 신호 탐지를 가능케 했던 최소 지반운동은 ~3.0 cm s−2 이상으로 확인되었다. 울진해역지진이 비록 동해 해역에서 발생하였지만, 진앙과 가 까운 강원도 남부-경상북도의 고지대를 따라 전파한 표면파의 지반운동으로 회절 인프라사운드가 효과적으로 발생한 것으로 해석된다. 인프라사운드 관측을 통한 원거리 지진 지반운동 특성 추정 방법은 지진관측망이 설치되어 있지 않거나 관측소 수가 적은 지역을 대상으로 활용이 가능할 것이다.

유연한 지반에 놓인 액체저장탱크의 지진 거동은 유체-구조물-지반 상호작용에 의해 복잡하게 나타나므로, 이 시스템의 지진응답 과 피해를 정확하게 예측하기 위해서는 이를 엄밀히 고려하여야 한다. 이 연구에서는 유체-구조물-지반 상호작용을 엄밀히 고려하여 유연한 지반에 놓인 직사각형 액체저장탱크의 지진응답 해석을 수행하고 그 응답 특성을 분석하고자 한다. 이를 위해 지진하중 작용 시 발생하는 유체의 동수압력 및 지반과 구조물 간의 상호작용력을 유한요소 기법을 사용하여 산정한다. 이때, 반무한 지반에서의 에너지 방사를 고려하기 위해 mid-point integrated finite element와 점성 감쇠기를 사용하여 지반 원역의 거동을 모사한다. 이와 같이 산정된 동수압력과 지반-구조물 상호작용력을 구조물의 유한요소에 작용시킨다. 자유장 해석을 통하여 입사 지진파에 의한 유효 지진력을 산정한다. 유연한 지반에 놓인 직사각형 액체저장탱크의 지진응답 해석을 통하여 지반-구조물 상호작용의 효과로 인해 시스템 응답의 변화가 다양하게 나타남을 확인할 수 있다. 그러므로, 유연한 지반에 놓인 직사각형 액체저장탱크의 내진설계를 수행하거나 내진성능을 검토할 때는 유체-구조물-지반 상호작용을 엄밀히 고려하여야 할 것이다.

In order to improve the ground-motion prediction equation, which is an important factor in seismic hazard assessment, it is essential to obtain good quality seismic data for a region. The Korean Peninsula has an environment in which it is difficult to obtain strong ground motion data. However, because digital seismic observation networks have become denser since the mid-2000s and moderate earthquake events such as the Odaesan earthquake (Jan. 20, 2007, ML 4.8), the 9.12 Gyeongju earthquake (Sep. 12, 2016, ML 5.8), and the Pohang earthquake (Nov. 15, 2017, ML 5.4) have occurred, some good empirical data on ground motion could have been accumulated. In this study, we tried to build a ground motion database that can be used for the development of the ground motion attenuation equation by collecting seismic data accumulated since the 2000s. The database was constructed in the form of a flat file with RotD50 peak ground acceleration, 5% damped pseudo-spectral acceleration, and meta information related to hypocenter, path, site, and data processing. The seismic data used were the velocity and accelerogram data for events over ML 3.0 observed between 2003 and 2019 by the Korean National Seismic Network administered by the Korea Meteorological Administration. The final flat file contains 10,795 ground motion data items for 141 events. Although this study focuses mainly on organizing earthquake ground-motion waveforms and their data processing, it is thought that the study will contribute to reducing uncertainty in evaluating seismic hazard in the Korean Peninsula if detailed information about epicenters and stations is supplemented in the future.

This study simulated strong ground motion waveforms in the southern Korean Peninsula, based on the physical earthquake modeling of the Southern California Earthquake Center (SCEC) BroadBand Platform (BBP). Characteristics of intensity attenuation were investigated for M 6.0-7.0 events, incorporating the site effects. The SCEC BBP is software generates broadband (0-10 Hz) ground-motion waveforms for earthquake scenarios. Among five available modeling methods in the v16.5 platform, we used the Song Model. Approximately 50 earthquake scenarios each were simulated for M 6.0, 6.5, and 7.0 events. Representative metrics such as peak ground acceleration (PGA) and peak ground velocity (PGV) were obtained from the synthetic waveforms that were simulated before and after the consideration of site effects (VS30). They were then empirically converted to distribution of instrumental intensity. The intensity that considers the site effects is amplified at low rather than high VS30 zones.

지진 하중에 의해 피해를 받는 교각의 보강을 위한 방법론과 다양성에 대한 연구가 활발히 진행되었지만, 단자유도에 대한 비선형 정적 분석에 대한 연구가 주로 수행되었으며, 교각을 구성하는 재료적 특성만이 변수로 간주되었다. 그러나, 교 각의 기하학적 요소와 보강 수단의 관계성에 관한 연구는 수행되지 않았다. 이에 본 연구에서는 CFRP 재킷을 사용하여 캘리포 니아에 존재하는 비 내진 상세 콘크리트 교량 교각을 보강하고 형상 변화에 따라 (교각 직경, 전단 경간 비, 보강 재킷의 길이) 교량 교각의 지진 취약성 곡선 도출하고 내진 성능을 평가하였다. 상기의 목표를 달성하기 위해 Opensees 프로그램를 사용하여 문헌에서 인용할 수 있는 실험체를 모델링하고 반복 하중을 가하여 결과 비교를 통해 해석방법의 적절성을 결정하였다.

The empirical Green’s function method is applied to the foreshock and the mainshock of the 2016 Gyeongju earthquake to simulate strong ground motions of the mainshock and scenario earthquake at seismic stations of seven metropolises in South Korea, respectively. To identify the applicability of the method in advance, the mainshock is simulated, assuming the foreshock as the empirical Green’s function. As a result of the simulation, the overall shape, the amplitude of PGA, and the duration and response spectra of the simulated seismic waveforms are similar with those of the observed seismic waveforms. Based on this result, a scenario earthquake on the causative fault of Gyeongju earthquake with a moment magnitude 6.5 is simulated, assuming that the mainshock serves as the empirical Green’s function. As a result, the amplitude of PGA and the duration of simulated seismic waveforms are significantly increased and extended, and the spectral amplitude of the low frequency band is relatively increased compared with that of the high frequency band. If the empirical Green’s function method is applied to several recent well-recorded moderate earthquakes, the simulated seismic waveforms can be used as not only input data for developing ground motion prediction equations, but also input data for creating the design response spectra of major facilities in South Korea.