Liquefaction Evaluation by One-Dimensional Effective Stress Analysis Using UBC3D-PLM Model

UBC3D-PLM 모델을 이용한 1차원 유효응력해석에 의한 액상화 평가

Abstract

This study compares the revised method in loose saturated sandy ground where the LNG storage tank will be installed with an evaluation method by one-dimensional effective stress analysis using the UBC3D-PLM model. Various laboratory and field tests were conducted to establish the parameters necessary for evaluation. The revised liquefaction evaluation method using the seismic response analysis result and N value from standard penetration testing evaluated the possibility of liquefaction as high, but assessment using effective stress analysis, which can consider various liquefaction resistance factors, found the site to be somewhat stable against liquefaction. One-dimensional finite element analysis using UBC3D-PLM modeling facilitated easier assessment of stability against liquefaction than the other methods and minimized the area required for reinforcement against liquefaction. In addition, it is expected that two-and three-dimensional numerical analysis considering the foundation of the LNG storage tank can identify the seismic design and behavior when liquefaction occurs.

본 연구는 LNG저장탱크가 설치될 느슨한 포화사질지반을 대상으로 개정된 액상화 평가법과 UBC3D-PLM 모델을 이용한 1차원 유효응력해석에 의한 액상화 평가법을 비교한 것이다. 이를 위해 여러가지 실내 및 현장시험을 실시하여 필요한 Parameter를 산정하였다. 검토결과, 지진응답해석결과와 SPT N 값을 이용하는 개정 액상화 평가법은 액상화 발생 가능성을 높게 평가하였지만, 다양한 액상화 저항인자를 고려할 수 있는 유효응력해석법은 액상화에 다소 안정한 것으로 분석되었다. UBC3D-PLM 모델을 이용한 1차원 유한요소해석을 할 경우 보다 간편하게 액상화 안정성 검토가 가능하였고, 액상화 보강 영역을 최소화 할 수 있었다. 또한, LNG저장탱크의 기초를 고려한 2·3차원 수치해석 시에는 액상화 발생 시 내진설계 및 거동특성을 규명하는 것에 활용될 수 있을 것으로 기대된다.

서론

LNG(Liquefied Natural Gas)는 천연가스를 -162°C 상태에서 냉각하여 액화시킨 뒤 부피를 1/600로 압축시킨 것이다. 이러한 LNG를 저장하는 시설이 LNG저장탱크(LNG Storage Tank)이다. LNG는 외국의 천연가스를 액체화시킨 뒤 LNG선으로 국내로 운반하기에 LNG저장탱크는 대부분 해안가에 위치하고 있다. 국내 해안가의 지반은 연약한 점토 및 느슨한 모래층의 출현빈도가 높아 LNG저장탱크 기초공법으로는 주로 강관말뚝이 많이 이용되고 있다. 또한 기초지반조건이 느슨한 모래층일 경우에는 반드시 액상화에 대한 안정성 검토가 필요하다. Fig. 1은 강관말뚝기초로 지지된 LNG저장탱크의 단면 및 평면도를 보여주고 있다.

Fig. 1. LNG storage tank.

액상화는 느슨하고 포화된 사질토 지반에 지진 하중이 작용하였을 때, 지반 내 간극수압 증가에 따른 유효응력 감소로 지반의 강도를 잃는 현상이다. 이러한 액상화 현상은 지진 시 유발되는 대표적인 재난적 지반공학적 현상으로 사회적, 경제적으로 많은 피해를 유발할 수 있다. 액상화로 인한 대표적 피해사례는 1964년 일본 니이가타(Niigata)에서 발생한 아파트 기울어짐과 교량 상판의 붕괴사고, 1976년 중국 탕산에서의 건물 붕괴사고(24만 명의 주민 사망), 2011년 뉴질랜드 크라이스처치(Christchurch)에서의 지반침하사고 등이 있다. 한편 국내에서는 2017년 포항시에서 공원 및 논밭 등지에 액상화가 발생하였으나, 큰 피해는 없었다. 이와 같이 세계 곳곳에서 액상화로 인한 인명 및 재산상의 손실을 발생시키고 있으며, 이에 수 많은 지반공학자들을 중심으로 하여 액상화 피해저감을 위한 연구가 지속적으로 이루어지고 있다.

이러한 액상화 현상 발생 가능성을 평가하기 위해, 2018년 한국지반공학회에서는 국가건설기준센터와 함께 국가건설 기준 “내진설계일반(KDS 17 10 00)”內 액상화 공통설계기준을 제정하였으며, 이를 반영하여 2020년 한국시설안전공단에서는 “기존 시설물(기초및지반) 내진성능 평가 요령(안)”內 액상화 평가를 개정한 바 있다. 기존 시설물(기초및지반) 내진성능 평가 요령(안)(KISTEC, 2020)에 따르면 액상화 평가는 예비평가와 본 평가 2단계로 구분되어져 있다. 액상화 예비평가는 지반조건(지하수위, 입도, 심도 등)을 고려하여 액상화 평가 생략여부를 결정하고, 액상화 본평가는 지반응답해석과 N값을 이용하여 액상화 발생 가능성을 평가하는 것이다. 하지만, 액상화 발생 영향인자는 N값, 과잉간극수압, 간극 수압, 연직유효응력, 지반의 투수계수 및 지반의 전단강도(점착력 및 내부마찰각) 등이 있으나, 국내기준에는 이러한 것들을 고려할 수 없다는 것이다. Park(2008)과 Park and Kwak(2009)에 따르면 지진응답해석과 N값을 이용한 액상화 평가방법은 액상화 발생 가능성을 과대평가하는 경향이 있으며, 보다 “합리적인 판단을 위해서는 유효응력해석이 필요하다”라고 제안하였다.

따라서, 본 연구는 2020년 개정된 액상화 예비평가 및 본 평가 결과(지진응답해석과 N값을 이용)와 유효응력해석에 의한 액상화 평가결과를 비교하고자 한다. 유효응력해석법에 의한 액상화 평가 시에는 지층모델링은 1차원으로 하였고, 예비 평가 및 본평가를 통해 확인된 액상화 발생 가능성이 있는 지층에는 UBC3D-PLM 모델을 적용하였다. 특히 UBC3D-PLM모델이 적용된 지층은 1 m간격으로, SPT조사 간격과 같게 모델링 함으로써 가급적으로 동일한 조건에서 이 둘을 비교하였다.

관련연구동향

유효응력해석에 의한 액상화 평가를 해야할 경우, 액상화 발생 시 지반의 응력-변형률 거동 및 과잉간극수압을 정확히 예측할 수 있는 신뢰성 높은 지반해석 모델의 선정이 매우 중요하다. 이렇게 액상화 평가에 이용 가능한 지반해석모델은 Iai 모델(Iai et al., 1990), Finn 모델(Byrne, 1991), PM4Sand 모델(Boulanger and Ziotopoulou, 2015), UBCSand 및 UBC 3D-PLM 모델(Beaty and Byrne, 2011) 등이 있다. 여기서 UBCSand 모델은 2차원 모델이며, UBC3D-PLM 모델은 2차원 및 3차원 적용이 가능한 모델이다.

유효응력해석에 의한 액상화 평가 관련연구는 국내외에서 다방면으로 이루어져왔다. Park(2008)은 유효응력해석과 등 가선형해석을 이용하여 액상화 평가 연구를 수행하여 등가선형해석은 지반의 전단응력을 과대평가하여 유효응력해석보다 액상화가 일어나기 쉬운 것으로 예측되었다. 또한 Park and Kwak(2009)은 간편예측법은 액상화 발생 가능성을 과대평가하는 경향이 있으며 반드시 정확한 판단을 위해서는 유효응력해석이 필요하다라고 제안하였다. 그리고 Daftari(2015), Tran et al.(2017) 및 Tung et al.(2020)에 의하면 Finn 모델은 액상화 후 발생되는 과잉간극수압을 예측할 수 없다고 평가하였다. Park and Kim(2022)은 California Imperial Valley의 Wildlife Site에서 1987년 11월 지진이 발생하였고, 이때 실측된 값(가속도와 간극수압)과 Finn 모델과 UBC3D-PLM 모델을 이용한 유효응력해석결과와 비교하였다. 그 결과 Finn 모델보다 UBC3D-PLM 모델에 의한 해석결과가 실측값과 유사한 과잉간극수압 거동을 예측하였고, 실측과 동일하게 액상화 발생여부를 평가할 수 있는 것을 확인하였다. 이외에도 액상화 발생 시 케이슨 안벽 거동 평가(Lee and Noh, 2016; Lee et al., 2018; Kim et al., 2019)에 유효응력해석법이 많이 사용되고 있다.

액상화 평가방법은 2010년 이전에는 주로 등가선형해석 및 간편예측법(액상화 본 평가)이 많이 사용되었으나, 2010년 이후에는 유효응력해석법의 사용이 점점 많아지고 있었다. 이는 등가선행해석 및 간편예측법은 액상화 발생 가능성을 과대평가하기 때문이다. 또한 유효응력해석에 필요한 지반해석 모델은 여러 가지 액상화 저항인자를 고려할 수 있고, 2 ‧ 3차원 모델링이 가능하며, 액상화 발생 후 과잉간극수압을 예측할 수 있는 UBC3D-PLM 모델이 가장 적합한 모델로 분석되었다.

예비평가 및 본 평가를 통한 액상화 발생 가능성 지층 평가

연구대상위치 및 지층특성분석

Fig. 2는 본 연구대상지의 시추위치와 시추주상도를 보여주고 있다. Fig. 2에서 LNG기초지반을 형성하기 위해 인공적인 매립이 이루어졌음을 알 수 있고, 좌측 및 우측 LNG 기초지반의 지층구조와 공학적 특성을 파악하기 위하여 10개의 시추조사와 다양한 실내 및 현장시험이 실시되었다. 10개의 시추조사를 통해 확인된 지층구조는 매립층(자갈/모래), 풍화토, 풍화암, 연암 및 보통암 순으로 분포하고 있었며, 매립층(자갈/모래)는 대부분 느슨한 층(N < 25)으로 확인되었다.

Fig. 2. Geological characteristics by borehole investigation.

액상화 예비평가를 통한 액상화 가능성 평가

액상화 예비평가는 지반 조건을 고려하여 액상화 평가 생략여부를 결정하는 것이다. 즉, 지하수위, 심도, 상대밀도, (N₁)₆₀, 세립토함유량 및 입도분포으로 평가하는 것이다. 액상화 예비평가 해설(KISTEC, 2020)의 각 항목에 해당하는 경우에는 “액상화에 대해 안전한 것으로 보고 액상화 평가를 생략한다”라고 제시되어 있다. 본 연구에는 이를 참조하여 액상화 예비평가를 실시하였고, 예비평가 결과를 Table 1에 나타내었다.

Table 1. Liquefaction preliminary evaluation result

(N₁)₆₀ = C_N × N₆₀

where, C_N: Overburden correction factor; N₆₀: Corrected SPT blow counts to account for the characteristics of the test.

액상화 본 평가를 통한 액상화 가능성 평가

액상화 예비평가결과 6개의 시추공에서 액상화 발생 가능성이 있는 것으로 평가되어, 액상화 본 평가가 필요한 것으로 검토되었다. 액상화 본 평가에서 액상화 발생 가능성은 대상 현장에서 액상화를 유발시키는 액상화 저항전단응력비(CRR_7.5)를 지진에 의해 발생되는 진동전단응력비(CSR)로 나눈 안전율로 평가하며, 아래 식으로 표현할 수 있다.

\(\begin{aligned}F S_{(M=6.5)}=\frac{C R R_{7.5}}{C S R} \times M S F>1.0\end{aligned}\)        (1)

여기서, CRR_7.5 : 지진규모 7.5에 해당하는 액상화 저항전단응력비

CSR : 진동전단응력비

MSF : 지진규모 보정계수(국내 지진규모 보정계수 = 1.5)

FS_(M=6.5) : 본 평가에 의한 액상화 안전율

본 연구에서는 지진구역계수(I구역), 평균재현주기(4,800년), 위험도계수(2.6) 및 지반의 분류(S4)을 고려하여 국내내진설계기준(KCSC, 2018)이 반영된 SeismoArtif 2020 프로그램으로 인공지진파를 추출하였다. 추출된 인공지진파를 ProShake 2.0 프로그램을 이용하여 지진응답해석을 수행하였고, 지진응답해석을 통해 지층별 최대전단응력을 산정하였다.

산정된 지층별최대전단응력, 현장조사 시 평가된 지층별SPT N값, 그리고 상재하중, 에너지 효율 및 세립분 함유량이 보정된 (N₁)_60cs을 이용하여 액상화 본 평가를 수행하였고, 그 결과 3개의 시추공(BT-10, BT-15 및 BT-16)에서 액상화 발생 가능성이 높은 것으로 나타났다(Table 2 참조).

(N₁)_60cs = (N₁)₆₀ + Δ(N₁)₆₀       (2)

여기서, (N₁)₆₀ : 상재하중과 에너지효율(60%)이 보정된 N값

\(\begin{aligned}\Delta\left(N_{1}\right)_{60}=\exp \left(1.63+\frac{9.7}{F C+0.01}+\left(\frac{15.7}{F C+0.01}\right)^{2}\right)\end{aligned}\)

FC : 세립토 함유량(%), 200번체 통과율

Table 2. Liquefaction main evaluation result

유효응력해석을 위한 지반모델

LNG Tank가 설치될 지반 상에서의 10개 시추공에 대해 액상화 예비평가 및 본 평가를 수행한 결과, 3개의 시추공(BT-10, BT-15, BT-16)에서 액상화 발생 가능성이 예측되었다. 이에 보다 더 신뢰성 있고 상세한 액상화 발생 가능성을 예측하기 위해 유효응력해석을 통한 액상화 평가를 수행하였다. 이에 대표적인 지반모델로는 Flac의 Finn 모델과 Plaxis의 UBC3D-PLM 모델이 있으나, Daftari(2015), Tung et al.(2020) 및 Park and Kim(2022)은 액상화 후 발생되는 과잉간 극수압을 예측할 수 있는 UBC3D-PLM 모델이 액상화 평가에 더 부합한 모델이라고 언급하였다. 따라서 본 연구에서 이용된 프로그램은 UBC3D-PLM 모델이 내장된 Plaxis2D를 선정하였고, 액상화를 평가하고자 하는 지층의 지반모델은 UBC3D-PLM 모델을, 연암의 지반모델은 Linear Elastic 모델을, 그 이외의 지층에는 HS Small Strain 모델을 적용하였다

UBC3D-PLM 모델

본 연구에서 액상화해석에 이용된 지반모델은 UBC3D-PLM 모델이다. 최초 Puebla et al.(1997)이 제안한 모델은 2차원 모델(UBCSAND)이었다. 그리고 Beaty and Byrne(1998, 2011)는 다양한 프로그램에 이를 적용하여 모델을 검증하였다. UBC3D-PLM 모델은 UBCSAND 모델을 3차원 모델화한 것이다. 초기 UBC3D-PLM 모델은 Tsegaye(2010)에 의해 개발되고 발표된 후, 유한요소 프로그램 Plaxis2D에서 User-defined Soil model로 적용되었고, 최근에는 정식모델로 이용되고 있다. Plaxis2D내에서 UBC3D-PLM 모델 Parameter는 표준관입시험((N₁)₆₀) 결과를 이용하여 산정할 수 있으며 (Bentley Systems International Ltd., 2021), 이는 아래와 같다.

① Stiffness parameter

K^e_G = 21.7 × 20.0 × (N₁)₆₀)^0.333       (3)

K^e_B = K^e_G x 0.7       (4)

k^p_G = k^e_G x ((N₁)₆₀)² x 0.003 + 100.0       (5)

여기서, K^e_G : 탄성전단계수

K^e_B : 체적전단계수

k^p_G : 소성전단계수

② Strength parameter

\(\begin{aligned}\phi_{p}=\phi_{p i}+\max \left(0.0, \frac{\left(N_{1}\right)_{60}-15}{5}\right)\end{aligned}\)       (6)

\(\begin{aligned}\phi_{p i}=\phi_{c v}+\frac{\left(N_{1}\right)_{60}}{10}\end{aligned}\)       (7)

R_f = 1.1 x (N₁)₆₀)^-.015 < 0.99       (8)

여기서, ϕ_cv : 내부마찰각

ϕ_p : 최대 내부마찰각

ϕ_pi : (N₁)₆₀이 보정된 내부마찰각

R_f : 파괴비(일반적으로 0.9)

③ Advanced parameter

(N₁)₆₀ = C_N⦁N₆₀       (9)

여기서, f_dens : 조밀화 요소(1.0)

_Epost : 액상화 후 요소(1.0)

HS Small 모델

Plaxis2D내에서 HS Small 모델은 변형경화모델, 즉 HS 모델(Hardening Soil model)에 미소변형률 강성이 고려된 확장모델이다. HS 모델의 기본이론은 다음과 같이 축변형률 ε₁과 축차응력 q의 쌍곡선 관계를 이용하는 것이다.

\(\begin{aligned}\varepsilon_{1}=\frac{1}{2 E_{50}} \frac{q}{1-\frac{q}{q_{a}}}\end{aligned}\)       (10)

여기서, q_a : 수렴축차응력

하중증가(loading)에 따른 E₅₀은 구속응력에 의존하는 탄성계수이며, 다음과 같이 표현할 수 있다.

\(\begin{aligned}E_{50}=E_{50}^{r e f}\left(\frac{c \cot \phi-\sigma_{3}^{\prime} \sin \phi}{c \cot \phi+p^{r e f}}\right)^{m}\end{aligned}\)       (11)

여기서, E^ref₅₀ : 하중증가 시 기준구속압(p^ref)에 상응하는 기준 탄성계수

m : 구속압에 대한 멱승계수(사질토 m = 0.5, 점성토 m = 1.0)

그리고 재하중(reloading) 또는 제하중(unloading)시 필요한 탄성계수 E_ur는 식 (12)로 표현할 수 있으며, 일반적으로 E^ref_ur ≈ 3E^ref₅₀값을 사용기도 한다.

\(\begin{aligned}E_{u r}=E_{u r}^{r e f}\left(\frac{c \cot \phi-\sigma_{3}^{\prime} \sin \phi}{c \cot \phi+p^{r e f}}\right)^{m}\end{aligned}\)       (12)

여기서, E^ref_ur : 재하 또는 제하 시 기준구속압(p^ref)에 상응하는 기준 탄성계수

극한 축차응력(q_f)과 수렴 축차응력(q_a)는 다음과 같이 정의할 수 있다.

\(\begin{aligned}q_{f}=\left(c \cos \phi-\sigma_{3}^{\prime}\right) \frac{2 \sin \phi}{1-\sin \phi}, q_{a}=\frac{q_{f}}{R_{f}}\end{aligned}\)       (13)

여기서, R_f : 파괴비(일반적으로 0.9으로 가정)

HS Small 모델에서는 미소변형률 강성을 고려하기 위한 입력변수를 2개로 단순화 시킨다. 일반적으로 최대 전단탄성 계수는 매우 낮은 변형률(<10^-5)에서 발생한다. 이러한 미소변형률 전단강성은 전단변형률 증가와 함께 전단탄성계수는 급격하게 감소하게 되며, 이러한 변형률 증가에 따른 전단강성변화를 고려하기 위하여 식(14)를 이용하여 전단변형률에 의존하는 비선형 전단탄성계수를 고려할 수 있다. Fig. 3은 전단변형률에 따른 비선형 전단탄성계수를 보여주고 있다.

\(\begin{aligned}G=\frac{G_{o}}{1+0.385 \gamma / \gamma_{0.7}}\end{aligned}\)       (14)

Fig. 3. Relationship between shear modulus and shear strain (Atkinson and Sallfors, 1991).

여기서, G_o : 최대전단탄성계수

γ_0.7 : G = 0.7G_o에 해당되는 전단변형률

유효응력해석을 통한 액상화 평가

1차원 유한요소망 및 경계조건

10개의 시추공에 대해 액상화 예비평가, 액상화 본 평가를 수행한 결과 BT-10, BT-15, BT-16 시추공 일부 지층에서 액상화 발생 가능성이 예측되었고, 이에 보다 더 신뢰성 있고 상세한 액상화 발생 가능성을 예측하기 위해 유효응력해석을 통한 액상화 평가를 수행하였다.

Fig. 4는 BT-10, BT-15 및 BT-16 시추주상도별 1차원 유한요소망을 보여주고 있다. 액상화 예비평가 및 본평가를 통해 액상화 발생 가능성이 있는 지층은 1 m간격으로 모델링 하였으며, 최하단 연암층도 1 m로 모델링하였다. 유한요소해석단계는 초기응력설정, 지진하중적용으로 총 2단계로 설정하였다. 해석단계별 변위 및 동적경계는 Table 3에 나타내었다. 특히 지진하중 적용 시 좌우측면의 변위경계는 Free로 동적경계는 Tied Degrees of Freedom을 적용함으로써 1차원 유효응력 액상화모델링이 가능하도록 설정하였다.

Fig. 4. Evaluation depth of liquefaction by borehole and one dimensional finite element.

Table 3. Displacement and dynamic boundary conditions for each analysis step

본 연구에서는 지진구역계수(I구역), 평균재현주기(4,800년), 지반의 분류(S4) 등을 고려하여 국내내진설계기준(KCSC, 2018)이 반영된 SeismoArtif 2020 프로그램을 이용하여 인공지진파를 추출하였다(Fig. 5 참조).

Fig. 5. Standard response spectrum and artificial seismic wave used for liquefaction evaluation.

지층별 지반모델별 해석 Parameter

Tables 4~6은 본 연구에서 적용된 지반모델별 해석 Parameter를 보여주고 있다. 대부분의 Parameter는 현장 및 실내 시험을 통해 산정하였으나, me, ne, np, f_dens, f_Epost, R_f 값은 프로그램에 내장된 기본값(Default Value)을 적용하였다. 그리고 투수계수는 문헌값에 제시된 값을 이용하였고, 매립자갈은 1.0 cm/s를, 매립모래는 0.01 cm/s를 적용하였다.

Table 4. UBC3D-PLM model parameter

Table 5. HS Small model parameter

Table 6. Linear elastic model parameter for soft rock

수치해석 결과 분석

Fig. 6은 BT-10, BT-15 및 BT-16 시추공 내 액상화 발생 가능성이 있는 지층에서의 지진하중으로 인해 발생한 과잉간 극수압을 보여주고 있다. 지진하중으로 인한 과잉간극수압의 증가량은 BT-10 > BT-15 > BT-16 시추공 순으로 나타나고 있으며, 과잉간극수압의 증가량이 큰 BT-10 시추공이 액상화에 가장 취약한 것으로 평가할 수 있다.

Fig. 6. Characteristics of time-excess pore pressure during earthquakes.

Fig. 7은 지진발생시 심도-과잉간극수압을 보여주고 있다. BT-10 시추공에서의 지진하중으로 인한 과잉간극수압 증가량(69.13 kN/m²)이 가장 크게 나고 있으며, BT-15 시추공에서 과잉간극수압 증가량은 61.66 kN/m²이며, BT-16 시추공에서 과잉간극수압 증가량(17.01 kN/m²)이 가장 작게 나타나고 있음을 알 수 있다. 이는 BT-10(N_ave = 9.8)의 N값이 BT-15(N_ave = 10.2) 및 BT-16(N_ave = 17.2)에 비해 작으며, 액상화 발생 가능성 지층(매립자갈)의 심도가 낮고, 지하수위가 높아(GL(-)2.0 m) 상대적으로 과잉간극수압이 크게 나타나고 있는 것으로 판단된다. 또한, BT-16에서의 과잉간극수압 증가량이 매우 작은 이유는 N값이 상대적으로 매우 크기 때문이다.

Fig. 7. Characteristics of depth-excess pore pressure during earthquakes.

Fig. 8은 지진발생시 심도-간극수압을 보여주고 있다. 여기서 간극수압이라 함은 “정수압+지진발생시 과잉간극수압”을 의미한다. BT-10 시추공에서의 과잉간극수압 증가량(69.13 kN/m²)으로 인해 간극수압 증가량이 매우 크게 나타나고 있다. BT-15 시추공에서의 과잉간극수압의 증가량(61.66 kN/m²)은 BT-10 시추공과 비슷하나, 심도가 깊어 간극수압 증가량은 BT-10보다 작음을 알 수 있다. BT-16에서의 간극수압 증가량은 가장 작음을 알 수 있다.

Fig. 8. Characteristics of depth-pore pressure during earthquakes.

Fig. 9는 심도-연직유효응력 관계도를 보여주고 있다. 지진하중 발생 전에는 심도가 깊이짐에 따라 연직유효응력은 증가하고 있었지만, 지진하중 발생 후에는 증가된 과잉간극수압으로 인해 액상화 발생 가능성 지층에서 연직유효응력이 감소하고 있었다. BT-15 및 BT-16에서는 지진하중 발생 후에도 연직유효응력이 (+)으로 유지하고 있으나, BT-10에서의 연직유효응력은 거의 0에 가깝거나, 일부에는 (-)으로 나타나고 있다. BT-10은 액상화 발생 가능성이 매우 높은 지층으로 추정할 수 있다.

Fig. 9. Characteristics of depth-Effective vertical stress during earthquakes.

Fig. 10은 심도-간극수압비(Ru) 관계도를 보여주고 있다. 여기서 간극수압비 = 과잉간극수압/연직유효응력을 말한다. 간극수압비(Ru)가 0.95 이상될 경우에는 액상화 발생 지층으로 판단할 수 있다. 시추공 BT-15 및 BT-16에서 최대 간극수압비는 0.46, 0.13으로 나타나 액상화에 안정한 지층으로 판단할 수 있으며, N값이 가장 큰 BT-16(N_ave = 17.2)가 액상화에 매우 안정한 것으로 나타났다. 하지만, 시추공 BT-10에서 최대 간극수압비(Ru)는 1.06으로 나타나 액상화에 불안정한 것으로 나타나 액상화 대책공법이 필요한 것으로 평가되었다.

Fig. 10. Characteristics of depth-pore pressure ratio during earthquakes.

결론

본 연구에서는 개정된 액상화 평가법(KISTEC, 2020)과 UBC3D-PLM 모델이 적용된 유효응력해석법을 이용해 지진 발생시 액상화 발생 가능성 평가를 수행한 후, 다음과 같은 결론은 얻었을 수 있다.

(1) N < 25인 느슨한 포화 모래(SC)및 자갈층(GC)이 두텁게 분포된 지층에 대해 액상화 예비평가((N₁)₆₀, 세립질함유량, 입도분포), 액상화 본 평가(표준관입시험((N₁)_60cs))를 수행한 결과 BT-10(GL(-)3.0 m~(-)7.0 m), BT-15(GL(-) 11.0 m~(-)14.0 m), BT-16(GL(-)8.0 m~(-)13.0 m) 시추공 심도에서 지진발생시 액상화 발생 가능성이 예측되었다.

(2) 보다 상세한 액상화 발생 가능성을 평가하기 위해 UBC3D-PLM 모델을 이용한 “유효응력해석을 통한 액상화 평가”를 수행하였다. 액상화 발생 가능성이 있는 지층의 지반모델은 UBC3D-PLM 모델을, 연암층은 Linear Elastic 모델을, 연암층을 제외한 지층의 지반모델은 HS Small 모델을 적용하였다.

(3) 지진하중 증가와 동시에 액상화 발생 가능성 있는 지층(UBC3D-PLM 모델이 적용된 지층)에서 지진하중에 의한 과잉간극수압이 급증하였고, BT-10 시추공에서의 과잉간극수압 증가량(69.13 kN/m²)이 BT-15(61.66 kN/m²) 및 BT-16(17.01 kN/m²) 시추공에서의 과잉간극수압 증가량보다 크게 평가되었다. 이는 BT-10 시추공에서 액상화 발생 가능성 지층(매립자갈)의 심도(GL(-)3.0 m~(-)7.0 m)가 낮고, 지하수위가 높아 상대적으로 과잉간극수압이 크게 나타나고 있는 것으로 판단된다.

(4) 지진하중에 의한 심도별 과잉간극수압, 심도별 간극수압 및 심도별 연직유효응력 분석을 통해 지진발생시 지반거동을 예측할 수 있었으며, 특히 심도별 간극수압비(Ru) 검토를 통해 지하수위 이하의 느슨한 매립층(자갈/모래)의 액상화 발생 가능성을 평가할 수 있었다.

(5) 개정된 액상화 평가법(KISTEC, 2020)은 BT-10, BT-15 및 BT-16 시추공에 액상화 보강이 필요한 것으로 평가되었으나, 유효응력해석법을 이용한 액상화 평가법은 BT-15 및 BT-16 시추공에서는 액상화가 발생하지 않는 것으로 평가되어 액상화 보강영역을 최소화할 수 있는 것으로 나타났다. 이는 개정된 액상화 평가법은 N값 및 세립분함유량으로만 평가하는 반면, 유효응력해석법은 N값, 과잉간극수압, 간극수압, 연직유효응력, 지반의 투수성 및 전단강도(점착력 및 내부마찰각) 등 보다 많은 액상화 저항요소를 고려할 수 있기 때문이다.

(6) 따라서 향후 지하수위이하의 느슨한 매립층(모래/자갈)에서 액상화 발생 가능성을 평가하기 위해서는 N값 및 세립분함유량을 이용한 액상화 평가는 가급적 지양하고, 다양한 액상화 저항인자를 고려할 수 있는 유효응력해석법을 이용한 액상화 평가를 실시하여 안전하고 보다 경제적인 기초설계를 해야 할 것이다.