Journal Search Engine

Download PDF Export Citation Korean Bibliography
ISSN : 1226-525X(Print)
ISSN : 2234-1099(Online)
Journal of the Earthquake Engineering Society of Korea Vol.20 No.6 pp.413-424
DOI : https://doi.org/10.5000/EESK.2016.20.6.413

Development of Korean Standard Vertical Design Spectrum Based on the Domestic and Overseas Intra-plate Earthquake Records

Jae Kwan Kim1) , Jung Han Kim2) , Jin Ho Lee3) , Tae Min Heo1)*
1)Department of Civil and Environmental Engineering, Seoul National University,
2)Integrated Safety Assessment Division, Korea Atomic Energy Research Institute,
3)Department of Ocean Engineering, Pukyong National University
Corresponding author: Heo, Tae Min intm2314@snu.ac.kr
September 7, 2016 October 14, 2016 October 14, 2016

Abstract

The vertical design spectrum for Korea, which is known to belong to an intra-plate region, is developed from the ground motion records of the earthquakes occurred in Korea and overseas intra-plate regions. From the statistical analysis of the vertical response spectra, a mean plus one standard deviation spectrum in lognormal distribution is obtained. Regression analysis is performed on this curve to determine the shape of spectrum including transition periods. The developed design spectrum is valid for the estimation both spectral acceleration and displacement. The ratio of vertical to horizontal response spectrum for each record is calculated. Statistical analysis of the ratios rendered the vertical to horizontal ratio (V/H ratio). Subsequently the ratio between the peak vertical ground acceleration to the horizontal one is obtained.


국내외 판내부 지진기록을 이용한 한국 표준수직설계스 펙트럼의 개발

김 재관1) , 김 정한2) , 이 진호3) , 허 태민1)*
1)서울대학교 건설환경공학부,
2)한국원자력연구원 종합안전평가부,
3)부경대학교 해양공학과

초록


    Ministry of Public Safety and Security
    MPSS-자연-2013-71

    1.서 론

    지진원에서 만들어진 지진파가 전파하면 지반이 진동하게 된다. 어떤 특 정 관측점에서의 지반운동은 직교하는 수평 2축 방향 성분과 수직 방향 성 분을 갖고 있다. 그러므로 내진설계에서는 원칙적으로 이 3가지 병진 운동 성분에 의한 영향을 고려해야 한다. 우리나라는 그 동안 1997년에 발간된 “내진설계기준연구” 보고서에 제시된 수평 지반운동에 대한 스펙트럼을 설계스펙트럼으로 사용했다[1]. 그리고 수직 지반운동을 고려해야 할 경우 에는 수평설계스펙트럼을 2/3배로 줄여서 사용하도록 권고하고 있다. 이것 은 수직 지반운동의 스펙트럼이 수평스펙트럼과 같고 단지 수직 지반가속 도만 차이가 있다는 것을 전제한다. 도로교 설계기준 케이블교량편에서는 (1) 수직 지반운동 PGA를 수평 지반운동 PGA의 1/2~2/3으로 취할 수 있 고 (2) 수직설계스펙트럼의 형상을 수평설계스펙트럼과 같게 취할 수 있도 록 허용하고 있다[2]. 그러나 실제 지진기록에서는 수직 지반운동은 수평 지반운동과 세기에서 차이가 있을 뿐 아니라 스펙트럼 형상에도 차이가 있 다. USNRC는 원자력 시설의 내진설계에서 수직 지반운동을 고려하도록 규정하고 있고 수직설계스펙트럼을 별도로 제시하고 있다[3]. Eurocode [4]와 최근에는 NEHRP[5]에서도 수직 설계스펙트럼을 별도로 제시하고 있다. 이러한 세계적 추세에 비추어 보면 우리나라도 수직설계스펙트럼을 별도로 규정할 필요가 있다.

    우리나라의 설계지진하중은 우리나라에서 발생한 지진의 기록을 사용 하여 구하는 것이 가장 이상적일 것이다. 그러나 우리나라는 큰 지진이 자주 발생하지 않고 지진계로 관측을 시작한 시간도 오래되지 않았다. 최근 기상 청이 지진관측망을 확대 구축하는 사업을 추진하여 최근에 발생한 지진의 기록은 많아지고 있지만 강진 기록의 수는 여전히 제한적이다. 우리나라 지 진학자들은 우리나라에서 예상되는 최대지진의 크기를 규모 6.8로 추정하 고 있다. 강진국가에서 지진위험도 해석과 설계스펙트럼 도출에는 일반적 으로 규모 5 이상의 지진에 의한 기록이 사용되고 있다. 말하자면 고려하는 지진 크기의 하한이 규모 5인 것이다[6]. 그 이유는 규모 5 이하의 지진에 의 해서는 유의미한 구조적 피해가 거의 발생하지 않기 때문이다[7]. 그러나 규모의 하한을 5.0으로 정하면 스펙트럼을 구하는 데 사용될 수 있는 우리 나라에서 발생한 지진의 수는 극히 적다. 그리고 규모 5이하의 지진에 의해 서 경미한 피해가 발생할 가능성이 있다. 그런데 USGS의 자료에 따르면 규 모 4.5 이하의 지진에 의해서는 지진 피해가 거의 발생하지 않는다고 한다. 그래서 이 연구에서는 고려하는 지진의 크기의 하한을 규모 4.5로 결정하였 다. 그러므로 우리나라 설계스펙트럼을 구하기 위해서 규모가 4.5에서 7 정 도인 우리나라에서 발생한 지진의 기록이 필요하였다.

    우리나라에서 기록된 최대지진은 규모 5.2로 평가되는 울진지진이다.1) 만약 우리나라 지진기록만 사용한다면 규모 5.2에서 규모 7까지 범위의 지진 에 대한 데이터가 없다. 심각한 피해는 주로 이 범위에 속하는 지진에 의해서 발생하는 것을 감안할 때 다른 대안을 사용해서 필요한 정보를 확보할 필요가 있다. 가능한 방법으로는 참고문헌[8]에서 서술한 바와 같은 경험적 그린함수 법(Empirical Green’s Function method, EGF)[9-15], 추계학적 방법 등 이 있다[16-20], 그리고 우리나라와 지진환경이 유사한 해외 판내부지역에 서 발생한 지진 기록을 활용하는 방법이 사용될 수 있을 것이다.

    수직설계스펙트럼은 실제 수직지반운동 기록으로부터 도출할 필요가 있 다. 추계학적 방법에 따르면 수평지반운동과 수직지반운동의 스펙트럼이 같게 얻어진다. EGF 방법은 실제 적용하였을 때 장주기 영역에서 오차가 확 대되는 문제가 있다는 것을 발견하였다. 이 문제에 대해서는 별도의 논문에 서 보고할 계획이다. 이러한 문제 때문에 또 수평설계스펙트럼과 도출방법 에 일관성이 있어야 하기 때문에 이 연구에서는 최종적으로 국내외 판내부 지진 기록을 사용하여 수직설계스펙트럼을 구하는 방안을 채택하였다.

    수직하중을 규정할 때 수직스펙트럼만으로는 충분하지 않다. 수평스펙 트럼에 대한 수직스펙트럼의 비에 대한 정보가 필요하다. 최소한 최대수평 지반가속도에 대한 최대수직지반가속도의 비가 제공되어야 한다. 이 연구 에서는 각 지진기록마다 기하평균스펙트럼으로 정의한 수평스펙트럼과 수 직스펙트럼을 계산하고 그 비를 구하였다. 이 비의 중앙값을 계산하고 주기 0.01초에서 값을 최대수평지반가속도에 대한 최대수직지반가속도의 비로 정의하였다.

    판내부 지진 기록은 구하기가 쉽지 않았다. 관측 기록 중에서 암반이나 단단한 지반에서 관측된 지진 기록은 더욱 수 제한적이었다. 그러나 한정된 개수의 기록으로 분석하였지만 통계적으로 유의한 결과가 얻어졌다. 이 결 과를 우리나라 설계스펙트럼으로 채택하는 것은 관계 전문가들이 의견을 수렴해서 결정할 사안이라고 생각된다.

    연구를 수행하면서 고무적인 사항과 개선되어야 할 사항들을 인지하게 되었다. 기상청이 지진관측망을 잘 구축하였기 때문에 비록 규모는 작지만 좋은 지진기록을 확보할 수 있었다. 그런데 국외 판내부 지진 기록을 구하는 문제는 앞으로 지속적으로 노력해야 할 사항이라고 생각된다.

    2.수직설계스펙트럼에 관한 외국 기준

    2.1.미국 원자력 시설의 설계 기준에서의 관련 규정

    USNRC Reg. Guide 1.6[3]에는 수직설계스펙트럼이 수평설계스펙트 럼과는 별도로 제시되어 있다. 3.5 Hz이상의 단주기 영역에서는 수직설계 스펙트럼은 수평설계스펙트럼의 값과 같아지고, 2.5 Hz 이하의 장주기 영 역에서는 수직설계스펙트럼은 수평설계스펙트럼의 2/3로 그 값이 주어진 다. 말하자면 V/H 비를 단주기 영역에서는 1, 장주기 영역에서는 2/3로 규 정하고 있다. 가속도에 민감한 구간인 T0TTS에서 증폭계수는 2.61부 터 2.98까지 변동한다.

    미국 ASCE의 4-98[21]의 규정도 Reg. Guide 1.6과 유사하다. 설계 지 반운동이 원거리 지진에 의해서 발생한다면 수직설계스펙트럼은 전체 주 기 영역에서 수평설계스펙트럼의 2/3로 취하도록 규정하고 있다. 그런데 만약 가까운 거리에서 발생한 지진에 의해서 설계 지반운동이 발생한다면 5 Hz 이상의 (0.2초 이하의) 단주기 영역에서는 수평설계스펙트럼과 같은 값 을 취하고 3 Hz 이하의 (0.33초 이상의) 장주기 영역에는 수평설계스펙트 럼의 2/3로 그 값을 취하도록 규정하고 있으며 5 Hz와 3 Hz 사이에서는 그 값이 1에서 2/3로 변동하도록 규정하고 있다.

    2.2.미국 NEHRP(2009)[5]에서의 규정

    미국 NEHRP에서 제시한 수직설계스펙트럼은 Fig. 2와 같다.

    Fig. 2에서 SDS는 가속도에 민감한 구간의 수평설계스펙트럼으로부터 정의되는데 이 값은 지반의 영향이 반영되어 있다. 그리고 Cv는 vertical coefficient로 정의되는 값이며 다음 Table 1에 주어져 있다.

    SS = 0.5이고 표준 지반인 SB의 경우에는 Cv = 0.9로 주어진다. 수직가 속도스펙트럼의 최대값은 SA = 0.8CvSDS이고 ZPA(Zero Period Acceleration) 즉, 최대수직지반가속도는 PGAv = 0.3CvSDS로 주어진다. 따라서 최대수 직지반가속도에 대한 최대가속도스펙트럼의 비는 0.8/0.3 = 2.67로 계산 되면 수평설계스펙트럼에서의 비인 2.5에 비하여 조금 큰 값을 갖는다. 최 대수평지반가속도 PGAh는 다음과 같이 계산될 수 있다.

    PGA h = 1 2.5 S DS = 0.4 S DS
    (1)

    그러면 여기서 우리는 최대수직지반가속도 PGAv와 최대수평지반가속 도 PGAh의 비를 추정할 수 있다.

    PGA v PGA h = 0.3 C v S DS 0.4 S DS = 0.3 × 0.9 0.4 = 0.675
    (2)

    수직설계스펙트럼은 주기 T ≥ 0.15초 이상에서는 T0.75에 반비례하여 감소하며 T ≥ 2초 이상인 장주기 영역에서는 수평설계스펙트럼의 1/2로 취할 수 있다. 그런데 주기가 0.05초≤ T ≤ 0.15초인 구간에서는 수평설계 스펙트럼에 대한 수직설계스펙트럼의 비는 0.8 × 0.9 = 0.7이다.

    이 결과는 비록 원자력 시설의 결과와 차이는 있지만 전체적으로 수직설 계스펙트럼의 주기에 따른 변동은 같은 경향을 갖고 있다고 서술하고 있다.

    2.3.Eurocode 8[4]에서의 규정

    앞서 언급하였듯이 Eurocode 8은 수직설계스펙트럼을 명확하게 별도 로 제시하고 있다. 수평설계스펙트럼과는 달리 설계지진규모에 무관하게 Fig. 3에 묘사된 하나의 설계스펙트럼이 사용된다.

    Eurocode 8에서는 최대수직지반가속도 PGAv의 최대수평지반가속도 PGAh에 대한 비를 다음과 같이 정의하고 있다.

    PGA v PGA h = 0.9 , 지진 규모M 5.5
    (3)

    또는

    PGA v PGA h = 0.45 , 지진 규모M < 5.5.
    (4)

    설계스펙트럼의 전이 주기 TB(T0), TC(TS), TD(TL)는 다음과 같이 규정 되고 있다.

    TB = 0.05초, TC = 0.15초, TD = 1.0초 가속도설계스펙트럼은 TBTTC에 일정하고 수직 지반 가속도에 대한 증폭 비율은 3η 로 주어지는데 이 때 η 는 감쇠비에 따라서 결정되는 parameter이다. 감쇠비 ξ = 5%이면 η = 1.0로 주어진다. 그러므로 증폭 비율은 3이 된다. TCTTD 구간에서 가속도설계스펙트럼은 1/T에 비례하여 감소하고 변위에 민감한 T > TD 구 간에서는 1/T2에 비례하여 감소한다.

    3.지진 목록과 분류

    분석에는 총 16개 지진에 의해 진앙거리 200 km 이하의 관측소에 관측 된 49개 지진기록이 사용되었다. 그 중에서 국내 지진 기록은 5개 지진에서 관측된 24개 기록으로 구성되어있고 국외 지진 기록은 11개 지진에서 관측 된 25개 기록으로 구성되어있다. 참고문헌[8]의 수평지진기록의 개수보다 적은데 그 이유는 외국의 2개 지진의 경우 수평 기록만 있고 수직성분의 기 록은 없기 때문이다. 규모 4.5 이상인 국내 지진의 개수는 더 많다. 그러나 진앙거리가 200 km 이하이고 암반에서 관측된 기록으로 제한하면 5개 지 진으로 줄어든다[22].

    Table 2에 국내 지진 기록을 수록하였다. 이 지진기록의 규모와 거리에 따른 히스토그램을 Fig. 4(a)와 Fig. 4(B)에 각각 제시하였다. 지진 기록은 모두 규모 5.0 근처에 모여 있는 것을 알 수 있다.

    국외 지진 기록은 Table 3에 수록하였다. 이 지진기록의 규모와 거리에 따른 히스토그램을 Fig. 5(a)와 Fig. 5(b)에 각각 제시하였다. 국내 기록에 비해서 비교적 규모가 큰 지진에 의한 기록이 많고 규모가 넓은 범위에 펼쳐 져 분포되어 있는 것을 알 수 있다.

    국외 지진기록의 지역별 분포를 Fig. 6에 정리하였다. 지진은 대부분 북 미 대륙과 유럽, 이란 그리고 인도에서 발생한 것을 알 수 있다. 아직 이 목록 으로는 세계 전체의 판내부지역을 망라하지는 못하고 있다는 것을 또한 인 지할 수 있다.

    국내, 국외 지진 기록을 합친 기록의 규모와 거리에 따른 분포를 Fig. 7(a)와 Fig. 7(b)에 히스토그램으로 나타냈다.

    비록 규모 5 근처 기록이 상대적으로 많지만 규모 6 이상의 지진 기록도 상당수 포함되어 있어서 전체적으로 균형을 맞추어주고 있다. 거리에 따라 서는 비교적 고른 분포를 보이고 있다.

    4.분석 방법과 절차

    4.1.V/H 비 계산

    한 관측점에서의 지진기록은 직교하는 수평 두 축 방향으로의 2개의 가 속도 성분과 하나의 수직 방향 가속도 성분으로 구성된다. V/H 비를 계산하 기 위해서 주기 T에서 수평스펙트럼과 수직스펙트럼이 필요하다. 수평스 펙트럼은 기하평균스펙트럼인 GMRotI50으로 정의하였다[8]. 수직스펙 트럼을 계산한 다음 그 값을 수평스펙트럼의 값으로 나누면 V/H 비가 얻어 진다. 즉, 한 관측점에서의 기록으로부터 하나의 V/H 비의 곡선이 얻어진 다. 같은 계산을 모든 지진기록에 대해서 실시한다. 그 결과로 우리는 49개 곡선을 얻게 된다. 각 주기마다 이 49개 곡선의 중앙값을 계산하면 우리는 V/H 비 중앙값 곡선을 얻게 된다. V/H 비는 주기에 따라서 변동한다.

    설계지진의 세기는 최대수평지반가속도로 정의되고 있다. 그러므로 수 직 지진하중을 완전하게 정의하기 위해서는 서론에서 언급하였듯이 최대 수평지반가속도에 대한 최대수직지반가속도의 비가 주어져야 한다. 이 연 구에서는 그 비를 주기 T = 0.01초에서 스펙트럼의 비 즉 V/H 비로 정의하 였다.

    4.2.수직설계스펙트럼을 구하기 위한 통계분석

    V/H 비를 구하기 위해서 계산된 수직스펙트럼을 T = 0.01 sec, 즉 100 Hz에서 가속도스펙트럼 값이 1.0 g가 되도록 정규화하였다. 여기서 g는 중 력가속도를 의미한다. 이렇게 하여 주어진 감쇠비에 대해서 각 관측점에서 하나의 정규화된 스펙트럼 곡선이 얻어진다. 모든 관측점에서 수직스펙트 럼을 구하면 우리는 일련의 스펙트럼 곡선을 얻게 된다. 개별 스펙트럼 곡선 은 주기에 대해서 들쭉날쭉한 모양을 갖고 특정 주기에서는 일련의 스펙트 럼 값은 어떤 확률 분포를 갖게 된다. 응답스펙트럼은 대수정규분포를 갖는 것으로 알려져 있다[23]. 대수정규분포의 특성은 평균 mean과 분포의 폭 파라미터인 표준편차 σ 에 의해서 서술될 수 있다. 원하는 신뢰도에서 스펙 트럼은 mean + 형태로 정의될 수 있다. 여기서 t는 표준편차의 배수를 의미하고 그 값은 신뢰도에 따라서 결정될 수 있다. 각 주기별로 mean + 에 해당되는 값을 연결하면 mean + 이 얻어진다. 만약 t = 0이면 초과 확 률이 50%인 평균스펙트럼이 얻어지고 t = 1이면 초과 확률이 약 16%인 mean + 1.0σ 스펙트럼이 얻어진다.

    원하는 신뢰도 수준에 따라서 t의 값을 결정하면 해당하는 스펙트럼 곡 선을 얻을 수 있다. 이 곡선은 개별 지진 기록의 스펙트럼보다는 매끈하지만 설계스펙트럼으로 가공하기 위해서는 회귀분석이 필요하다. Fig. 1, 2, 3에 예시한 세계 여러 나라 설계기준에서 채택하는 가속도설계스펙트럼은 대 체로 Fig. 8에 예시한 조각적 연속함수의 형태로 주어진다.

    Fig. 8의 설계스펙트럼은 4개의 구간으로 구분된다. 각 구간의 절점주기 혹은 전이주기인 T0, TS, TL을 경계로 구분된다. 원점은 T = 0.01초를 나타 낸다. 이 주기에서 가속도스펙트럼은 최대지반가속도(Peak Ground Acceleration, PGA)와 같아진다. 그런데 어떤 지진 기록에서는 가속도스 펙트럼이 PGA와 같아지는 주기 Tz가 0.01초보다 더 길 수 있다. T0는 가속 도스펙트럼이 평탄한 구간이 시작되는 주기, TS는 가속도에 민감한 구간과 속도에 민감한 구간의 전이주기, TL은 속도에 민감한 구간과 변위에 민감한 구간의 전이주기이다. 설계스펙트럼의 형상을 결정하기 위해서는 각 구간 에서 곡선의 함수 형태가 미리 정해져야 한다. 설계스펙트럼 곡선은 원점에 서 T0까지 구간에서는 직선으로 증가하고, T0에서 TS까지 구간에서는 수평 직선, TS에서 TL까지 속도에 민감한 구간에서는 1/T에 비례하여 감소하는 곡선, TL 이상에서는 1/T2 비례하여 감소하는 곡선의 형상을 갖는 것으로 알 려져 있다. 최종적으로 가속도가 평탄한 구간에서 가속도스펙트럼의 최대 지반가속도(PGA)에 대한 비인 증폭계수 αA를 알아야 한다.

    증폭계수 αA와 전이주기 T0, TS, TL은 통계분석 결과로 얻은 mean + 가속도스펙트럼 곡선에 회귀분석을 실시해서 구할 수 있다. 그런데 장주기 영역에서 가속도스펙트럼의 값은 작아진다. 따라서 TL이 정확하게 식별되 지 않을 우려가 있다. 그러나 변위스펙트럼에서는 장주기 영역에서 곡선의 변화가 더 분명하게 드러날 수 있다. 변위스펙트럼 SD와 가속도스펙트럼 SA의 관계는 다음과 같다.

    SD ω = SA ω ω 2 = T 2 4 π 2 SA ω
    (5)

    이 식에 따르면 SD는 속도에 민감한 구간에서는 주기 T에 비례하여 직 선적으로 증가하고 변위에 민감한 구간에서는 일정한 상수이다. 변위스펙 트럼의 전형적인 형상을 Fig. 9에 예시하였다. 변위스펙트럼은 TTL 구간 에서 일정한 값을 유지하다가 다시 감소하고 T가 계속 길어지면 다시 일정 한 값으로 수렴하는 경향을 보인다. 이 두 가지 전이주기를 TL1, TL2로 명명 하기로 하자. 만약 변위설계스펙트럼을 별도로 제시해야 한다면 전이주기 TL1, TL2도 식별되어야 할 필요가 있다. TL은 변위스펙트럼에서 더 명확하게 식별될 수 있다는 것을 Fig. 9는 보여준다.

    이 연구에서 TL, TL1, TL2은 변위스펙트럼에서 식별하고, TL이 고정된 조 건에서 가속도스펙트럼으로부터 증폭계수 αA와 전이주기 T0, TS를 식별하 는 방안을 채택하였다. 그리고 회귀분석에서 제곱근제곱합(SRSS) 오차를 최소화하는 방법을 사용하여 목표 함수의 파라미터를 구했다. 변위스펙트 럼에서 그 값이 최대가 되는 TLTTL1에서의 값을 βD라고 하자. βD는 회귀분석에서 TL 등과 같이 식별될 수 있고 동시에 변위설계스펙트럼도 구 할 수 있다.

    5.분석 결과

    5.1.V/H 비

    각 지진기록의 수직스펙트럼의 값을 기하평균(GMRotI50)으로 정의한 수평스펙트럼으로 나누어서 주기마다 비를 구하였다. 이때 감쇠비 ξ = 5%이다.

    V / H T = SA v T SA h T
    (6)

    여기서 SAv(T) 는 주기 T에서 수직스펙트럼이고 SAh(T) 는 같은 주기에서 수평스펙트럼이다. 이리하며 각 지진기록마다 하나의 곡선이 얻어지고 모 두 49개의 곡선이 구해진다. 이 곡선들을 Fig. 10에 도시하였다. Fig. 11에 는 각 주기마다 중앙값을 계산하여 연결한 곡선이 주어져 있다.

    Fig. 11에서 mean을 취하지 않고 median을 취한 이유는 SAv(T)와 SAh(T)가 대수정규분포를 갖고 따라서 V/H(T)도 대수정규분포 함수이기 때문이다. Fig. 11에서 우리는 중앙값 곡선이 주기 T에 따라서 변동의 폭이 아주 작고 거의 일정하다는 것을 볼 수 있다. 이 현상은 수직스펙트럼과 수 평스펙트럼의 형상이 유사할 것이라는 것을 시사한다고 하겠다.

    일반적으로 수직성분 지진하중은 PGAv = 1.0 g에 정규화된 스펙트럼 과 PGAv / PGAh 비로 규정된다. 이 비와 V/H 비는 다음 식에 의해서 주어 진다.

    PGA v PGA h = V / H 0.01 = SA v 0.01 SA h 0.01
    (7)

    여기서 T = 0.01초로 가정한다. Fig. 11에서 T = 0.01초에서 V/H 비를 읽 으면 PGAv / PGAh = 0.77 으로 얻어지며 설계스펙트럼에서 주로 관심을 갖는 TL = 0.3초 까지 평균값 역시 0.77이다. 따라서 이는 넓은 범위의 주기 에 걸쳐 상당히 합리적인 값임을 알 수 있다. 이 값을 국내 및 국외 설계기준 에서의 값과 Table 4에서 비교하였다.

    지금까지 우리나라에서 사용되어 왔던 2/3에 비해서 큰 값이지만 Eurocode 8에서의 0.9에 비하면 작은 편이다.

    5.2.수직설계스펙트럼의 형상

    국내외 49개 지진기록의 수직응답스펙트럼을 Fig. 12(a)에 도시하였다. 그림에 mean, mean + 1.0σ 스펙트럼도 함께 도시하였다. Fig. 12(b)에는 이 두 가지 스펙트럼만 별도로 도시하였다.

    회귀분석을 실시하기 전에 표본의 적정성을 검정하였다. 즉 Type 3 오 류가 발생하는지 Cross Validation 방법으로 확인하였다. 49개 지진 기록 에서 각 데이터 1개씩을 제외한 48개 데이터로 구성된 49개의 표본을 작성 하였다. 이 49개의 표본의 각각에서 mean, mean + 1.0σ 를 구하고 변동성 을 조사하였고 그 결과를 Fig. 13(a), (b)에 각각 도시하였다. 이 결과들에서 표본에 따른 변동이 아주 작다는 것을 볼 수 있다. 결과적으로 49개 지진 기 록으로 구성된 표본이 적정하다는 것이 검정되었다고 하겠다.

    두 가지 신뢰도에 대한 스펙트럼에 회귀분석을 실시하여 설계스펙트럼 을 구하였다. 앞서 설명한 바와 같이 변위스펙트럼에서 TL, TL1, TL2를 식별 하여 Fig. 14에 그 결과를 제시하였다.

    mean 스펙트럼에서 TL = 3초, βD = 0.0095s2 , mean + 1.0σ 스펙트럼 에서 TL = 3초, βD = 0.056s2로 식별되었다. TL을 고정시키고 TL에서의 SA(TL) = (4π2/TL2)SD(TL) = (4π2/TL2)βD × PGA로 고정시킨 다음 가속 도스펙트럼에 회귀분석을 실시하였다. 그 결과로 얻어진 설계스펙트럼과 전이주기 및 증폭계수를 Fig. 15에 제시하였다. Fig. 15에서 제시한 수직설 계스펙트럼을 참고문헌[8]에서 제시한 수평설계스펙트럼의 형상과 Fig. 16에서 비교하였고 주요 파라미터를 Table 5에서 비교하였다.

    가속도에 민감한 구간에서 수직스펙트럼의 최대증폭계수 αA의 값은 수 평스펙트럼의 증폭계수에 비해서 약간 큰 값을 갖는다. 그러나 그 차는 1% 보다도 작다. 그리고 전이주기는 실질적으로 같은 값을 갖는 것을 확인할 수 있다. 이 사실은 수평설계스펙트럼의 형상을 수직설계스펙트럼의 형상으 로 사용할 수 있다는 것을 시사한다.

    앞서 살펴본 외국의 수직설계스펙트럼은 대응하는 수평설계스펙트럼 과 비교하여 단주기 영역으로 이동하고 가속도에 민감한 구간이 좁아지는 것을 보여주고 있다. 그런데 본 연구에서는 그렇지 않다는 결론이 도출되었 다. 그 원인을 알아보기 위해서 국내지진 기록만 사용하여 얻은 수평/수직 스펙트럼과 국외지진기록만을 사용한 수평/수직 스펙트럼을 계산하여 비 교하였다. Fig. 17(a)에는 국외지진의 스펙트럼이, Fig. 17(b)에는 국내지 진의 스펙트럼이 비교되어 있다.

    국외지진의 경우에는 수직스펙트럼은 수평스펙트럼에 비해서 단주기 로 이동하는 것을 Fig. 17(a)는 분명하게 보여주고 있다. 그런데 국내지진 의 경우에는 수직스펙트럼의 전이주기가 수평스펙트럼의 전이주기와 거의 같다는 것을 Fig. 17(b)에서 알 수 있다. 즉, Fig. 17(b)는 국내지진의 특성 이 외국지진과 다르다는 것을 시사한다. 다른 하나의 가능성은 국내지진의 규모는 대체로 5.0 정도인데 비해서 국외지진은 큰 지진도 포함되어 있다. 규모가 6.0보다 작은 국외 3개 지진의 기록에 의한 수직 및 수평 mean + 1.0σ 스펙트럼을 Fig. 18에서 비교하였다. 여기서 우리는 국외 지진의 경 우에도 규모가 작으면 수직스펙트럼의 전이주기가 수평스펙트럼의 전이주 기와 유사하다는 것을 볼 수 있다.

    따라서 Fig. 18에서 수직스펙트럼의 형상이 수평스펙트럼의 형상과 거 의 같게 도출된 것은 국내지진과 규모가 작은 국외지진의 영향이 크게 반영 된 결과라는 것을 이해할 수 있다. 이러한 사실들은 본 연구에서 얻은 수직 스펙트럼의 형상이 타당하다는 것을 보여주고 있다고 하겠다.

    5.3.수직설계스펙트럼

    5.1과 5.2절의 결과를 종합하면 수직성분 지진하중 즉 수직설계스펙트 럼이 얻어진다. 지진계수는 수평지진가속도로 표현되므로 수직설계스펙 트럼도 최대수평지반가속도 PGAh = 1.0 g에 대해서 정규화하였다. Fig. 19에는 수직설계스펙트럼을 도시하였고 Table 6에는 전이주기를 비롯한 파라미터들이 수록되어 있다.

    6.국내외 기준과의 비교

    6.1.국외기준과의 비교

    앞서 언급하였듯이 설계기준에서는 보수적으로 mean + 1.0σ 스펙트 럼을 채택하고 있다. 이 연구에서 얻은 스펙트럼을 USNRC 스펙트럼, NEHRP 스펙트럼 및 판내부 지역에 속하는 유럽의 설계스펙트럼과 비교 하였다. Fig. 20에서 비교한 스펙트럼들은 최대수평지반가속도 PGAh = 1.0 g에 대해서 정규화하였다. Table 7에는 증폭계수와 전이주기들이 비교 되어 있다.

    Fig. 20과 Table 7에서 Eurocode 8과 NEHRP 스펙트럼은 단주기 영역 이 잘 발달되어 있고 USNRC 스펙트럼은 상대적으로 장주기 영역이 잘 발 달되어 있는 것을 볼 수 있다. 본 연구에서 제시한 스펙트럼은 단주기 영역 과 장주기 영역의 폭넓은 범위에 걸쳐서 발달되어 있는 것을 알 수 있다. 어 떻게 보면 보수적인 결과라고 생각할 수 있다. 가속도스펙트럼의 최대값은 Eurocode 8과 USNRC 스펙트럼보다는 작고 NEHRP 스펙트럼보다는 크 다는 것을 또한 알 수 있다.

    6.2.국내 기준과의 비교

    본 연구 결과로 제시되는 수직설계스펙트럼과 기존 우리나라 내진설계 기준II[1]의 설계스펙트럼의 형상을 Fig. 21에 그리고 증폭계수와 전이주 기는 Table 8에 비교하였다. 내진설계기준II[1]에서는 최대수직지반가속 도는 최대수평지반가속도의 2/3로 취하도록 규정하고 있다. 두 가지 설계 스펙트럼을 최대수평지반가속도 PGAh = 1.0 g에 대해서 정규화하였다.

    Fig. 21과 Table 8는 가속도에 민감한 구간에서는 기존의 스펙트럼에 비해서 새로운 스펙트럼의 값이 대폭 증가하는 것을 보여주고 있다. 그러나 속도에 민감한 구간에서는 거의 값이 같고 변위에 민감한 구간에서는 오히 려 새로운 스펙트럼이 작아지게 된다. 지진하중의 수직성분은 단주기 구조 물보다는 장주기 구조물 특히 사장교나 현수교와 같은 장대 교량에서 중요 할 수 있다. 이러한 측면에서 본다면 본 연구에서 제시하는 수직설계스펙트 럼은 상당히 합리적이라고 말할 수 있다.

    7.결 론

    국내외 판내부 지진기록에 근거하여 V/H 비와 수직설계스펙트럼을 구 하였다. 규모 4.5 이상의 국내외 판내부 지진 기록 49개를 수집하여 통계 분 석하였다. 비록 기록의 개수는 많지는 않지만 통계적으로 유의한 결과를 얻 을 수 있었다.

    V/H 비로부터 PGAv / PGAh = 0.77으로 얻어졌다. 스펙트럼 분석에서 가속도설계스펙트럼뿐 아니라 변위설계스펙트럼도 같이 얻어졌고 형상을 결정하는 파라미터도 회귀분석에 의해서 αA = 2.8, T0 = 0.06초, TS = 0.3 초, TL = 3초로 구해졌다. 이 결과들로부터 최대수평지반가속도 PGAh = 1.0에 대해서 정규화된 수직설계스펙트럼을 제시하였다.

    외국의 설계스펙트럼에 비교하여 새로 제시한 수직설계스펙트럼은 단 주기에서 장주기 영역에 걸쳐서 잘 발달되어 있고 그렇기 때문에 보수적이 라고 판단된다. 증폭계수도 외국의 기준에 비해서 비교적 합리적임을 확인 하였다.

    내진설계기준연구II에 제시된 설계스펙트럼과 비교하였을 때 본 연구 에서 제시한 설계스펙트럼을 적용한다면 주기가 0.39초보다 짧은 단주기 구조물의 경우에는 설계하중이 증가하고 주기가 0.39초보다 긴 장주기 구 조물의 경우에는 설계하중이 감소하게 된다. 특히 전이주기 TL이 명확하게 제시되었기 때문에 변위에 민감한 영역에서 설계변위는 크게 감소할 것으 로 예상된다.

    서론에서 언급하였듯이 본 연구의 결과는 규모 4.5 이상의 국내 지진과 국외 판내부 지진 기록을 사용하여 얻었다. 여러 가지 제약 조건으로 인해서 더 많은 수의 기록을 확보하지 못하였고 지진학적 관점에서 더 자세하게 분 석하지는 못했다. 실제 지진기록을 이용한 우리나라 지진하중을 구하는 본 격적인 연구는 이제 시작되었다고 해도 지나치지 않을 것이다. 앞으로 후속 연구를 통해서 본 연구의 한계를 극복하고 좀 더 완전한 결과가 제시되기를 기대한다.

    감사의 글

    본 연구는 정부(국민안전처)의 재원으로 재난안전기술개발사업단의 지 원을 받아 수행된 ‘지진가속도 분석 등을 통한 설계지반운동 표시방법 고도 화기술 개발’ [MPSS-자연-2013-71] 과제의 지원을 받아 수행한 연구입니 다. 이에 감사드립니다. 특히 저자들은 연구의 주요 고비마다 해결의 실마 리를 제공해준 그리스 Patras 대학교의 Apostolos S. Papageorgiou 교수 께 큰 감사를 드리고, 구하기 어려운 판내부 지진 데이터를 흔쾌히 공유한 미국 Notre Dame 대학교의 George P. Mavroedis 교수께도 깊은 감사의 뜻을 전합니다.

    Figure

    EESK-20-6-413_F1.gif

    Vertical design spectrum in USNRC Regulation Guide 1.6[3]

    EESK-20-6-413_F2.gif

    Vertical design spectrum in NEHRP[5]

    EESK-20-6-413_F3.gif

    Vertical design spectrum in Eurocode 8[4]

    EESK-20-6-413_F4.gif

    Distribution of domestic earthquake records: (a) Frequency histogram w.r.t. magnitude, (b) Frequency histogram w.r.t. distance

    EESK-20-6-413_F5.gif

    Distribution of overseas intra-plate earthquake records: (a) Frequency histogram w.r.t. magnitude, (b) Frequency histogram w.r.t. distance

    EESK-20-6-413_F6.gif

    Regional distribution of overseas intra-plate earthquake records

    EESK-20-6-413_F7.gif

    Distribution of domestic and overseas intra-plate earthquake records: (a) Frequency histogram w.r.t. magnitude, (b) Frequency histogram w.r.t. distance

    EESK-20-6-413_F8.gif

    The shape of a typical acceleration design spectrum

    EESK-20-6-413_F9.gif

    The shape of a typical displacement design spectrum

    EESK-20-6-413_F10.gif

    V/H ratios of earthquake records both domestic and overseas

    EESK-20-6-413_F11.gif

    The median value curve of V/H ratios

    EESK-20-6-413_F12.gif

    Geometric mean spectra: (a) spectra for 49 records, (b) mean and mean + 1.0σ spectra

    EESK-20-6-413_F13.gif

    Cross Validation results: (a) mean spectrum, (b) mean + 1.0σ spectrum

    EESK-20-6-413_F14.gif

    Displacement spectra and their parameters: (a) mean spectrum, (b) mean + 1.0σ spectrum

    EESK-20-6-413_F15.gif

    Acceleration spectra and their parameters: (a) mean spectrum, (b) mean + 1.0σ spectrum

    EESK-20-6-413_F16.gif

    Comparison between horizontal and vertical spectra

    EESK-20-6-413_F17.gif

    Comparison of vertical spectrum with horizontal one: (a) Overseas earthquakes; (b) Domestic earthquakes

    EESK-20-6-413_F18.gif

    Comparison of vertical spectrum with horizontal one for overseas earthquakes of magnitude less than 6.

    EESK-20-6-413_F19.gif

    Vertical design spectrum

    EESK-20-6-413_F20.gif

    Comparison of mean + 1.0σ spectrum with the design spectra for various overseas intra-plate regions

    EESK-20-6-413_F21.gif

    Comparison of mean + 1.0σ spectrum with the design spectra for 1997 design code

    Table

    Values of vertical coefficient Cv

    List of domestic earthquake records[22]

    List of overseas intra-plate earthquake records

    Comparison of the PGAv / PGAh ratios

    Comparison of amplification factor and other parameters

    Parameters for vertical design spectrum

    Comparison of parameters of mean + 1.0σ spectrum with those of design spectra for various overseas intra-plate regions

    Comparison of parameters of mean + 1.0σ spectrum with those of the design spectra for 1997 design code

    Reference

    1. EESK, KICT (1997) The Research of Seismic Design Code , Ministry of Construction and Transportation, ; pp.37
    2. Korea Road Association (2014) Korean Highway Bridge Design Code (Limit state design) Cable Bridge part , Ministry of Construction and Transportation, ; pp.9-19
    3. US Nuclear Regulatory Commission (NRC) (2014) Design Response Spectra for Seismic Design of Nuclear Power Plants (Regulatory Guide 1.60) , Washington DC. USA,
    4. Comité EuropEen de Normalisation (2004) Eurocode 8:Design of Structures for Earthquake Resistance Part 1 General Rules Seismic Actions and Rules for Buildings , Comit Europ en de Normalisation c2004, ; pp.40
    5. BSSC, NEHRP (2009) Recommended Seismic Provisions for New Buildings and Other Structures (FEMA P-750) , Washington DC. USA Federal Emergency Management Agency,
    6. McGuire RK (2004) Seismic hazard and Risk Analysis EERI Monograph MNO-10 , Oakland California USA Earthquake Engineering Research Institute, ; pp.221
    7. Nakamura Y (1996) Real-Time Information Systems for Hazard Mitigation, Paper No. 2134
    8. Kim JK , Kim JH , Lee JH , Heo TM (2016) Development of Korean standard horizontal design spectrum based on the domestic and overseas intra-plate earthquake records , Journal of Earthquake Engineering, Vol.20 (6) ; pp.369-378
    9. Hartzell SH (1978) Earthquake Aftershocks as Green's Functions , Geophysical Research Letters, Vol.5 (1) ; pp.1-4
    10. Kamae K , Irikura K , Pitarka A (1998) A Technique for Simulating Strong Ground Motion Using Hybrid Green’s Functions , Bulletin of the Seismological Society of America, Vol.88 (2) ; pp.357-367
    11. Mueller CS (1985) Source Pulse Enhancement by Deconvolution of an Empirical Green’s Function , Geophysical Research Letters, Vol.12 (1) ; pp.33-36
    12. Kanamori H (1979) A Semi-Empirical Approach to Prediction of Long-Period Ground Motions from Great Earthquakes , Bulletin of the Seismological Society of America, Vol.69 (6) ; pp.1645-1668
    13. Irikura K (1983) Semi-empirical Estimation of Strong Ground Motions During Large Earthquakes , Bull. Disast. Prev. Res. Inst. Kyoto Univ, Vol.33-2 (298) ; pp.63-104
    14. Imagawa K , Mikami N , Mikumo T (1984) Analytical and Semi-Empirical Synthesis of Near-Field Seismic Waveforms for Investigating the Rupture Mechanism of Major Earthquakes , J. Phys. Earth, Vol.32 ; pp.317-38
    15. Mikumo T , Miyatake T BOLT B.A (1987) Numerical Modeling of Realistic Fault Rupture Processes , Chapter 3 in Seismic Strong Motion Synthetics, ACADEMIC PRESS INC,
    16. Boore DM (1983) Stochastic Simulation of High-frequency Ground Motions Based on Seismological Models of the Radiated Spectra , Bulletin of the Seismological Society of America, Vol.73 (6A) ; pp.1865-1894
    17. Boore DM (2003) Simulation of Ground Motion Using the Stochastic Method , Pure and Applied Geophysics, Vol.60 ; pp.635-676
    18. Jo ND , Baag CE (2001) Stochastic Prediction of Strong Ground Motions in Southern Korea , Earthquake Engineering Society of Korea, Vol.5 (4) ; pp.17-26
    19. Jo ND , Baag CE (2003) Estimation of Spectrum Decay Parameter and Stochastic , Prediction of Strong Ground Motions in Southern Korea, Vol.7 (6) ; pp.59-70
    20. Junn JG , Jo ND , Baag CE (2002) Stochastic prediction of ground motions in southern Korea , Geosciences Journal, Vol.6 (3) ; pp.203-214
    21. American Society of Civil Engineers (ASCE) (1999) Seismic Analysis of Safety-Related Nuclear Structures and Commentary , ASCE 4-98, ; pp.2-6
    22. KMADomestic Earthquake List. http://www.kma.go.kr/weather/earthquake_volcano/domesticlist.jsp,
    23. Davenport WB Jr (1970) Probability and random processes an introduction for applied scientists and engineers, McGraw-Hill, ; pp.542
    Journal Abbreviation J. Earthq. Eng. Soc. Korea
    Frequency Bimonthly
    Doi Prefix 10.5000/EESK
    Year of Launching 1997
    Publisher Earthquake Engineering Society of Korea
    Indexed/Tracked/Covered By