Aims & Scope

이번 연구에서는 양산단층대의 남쪽에 위치하는 경남 양산시 울주군 언양읍 상천리 지역(이하 SCL)과 북쪽에 위치하는 경북 포항시 북구 송라면 보경사 인근(이하 BGS)에서 채취한 단층비지를 연구 대상으로 선정하였다. 야외노두에서 단층의 자세, 변형구조, 모암 등의 주요 특징을 파악하고, 단층핵과 주변 부차 단층(subsidiary fault)의 단층비지를 각각 채취하였다.

상천리 지역은 하양층군의 퇴적암과 유천층군의 화산암류 그리고 불국사 화강암류가 주로 분포한다. 양산단층을 기준으로 동쪽지괴에는 주로 사암과 셰일로 이루어진 하양층군 퇴적암류가 분포하고, 서쪽지괴에는 유천층군의 화산암류와 불국사화강암류가 분포한다(Choi et al., 2009). 시료를 채취한 곳은 양산단층을 가로지르는 동서방향의 가천천을 따라 노출된 노두로 단층핵과 혼합대(mixed zone) 그리고 손상대가 연속적으로 잘 노출되는 지역이다(그림 2a, 2b; Choi et al., 2009). 주단층면의 자세는 N06°E/ 73°SE이며, 단층핵의 모암은 대부분 화성기원의 안산암질 응회암으로 우수향 전단을 지시하는 S-C fab-ric이 흔히 인지된다. 실험을 위해 단층핵에서 두 개의 시료(SCL 1-1, 1-2)와 부차 단층에서 하나의 시료(SCL 1-3)를 채취하였다.

Fig. 2. 
Photograph of outcrops where fault gouges, (a) SCL 1-1, (b) SCL 1-3, (c) BGS 3-1, and (d) BGS 3-2, were collected. SCL 1-1 and BGS 3-1 are fault gouges collected from the fault core zone, and SCL 1-3 and BGS 3-2 are from the subsidiary fault within the damage zone. The white circle on each of the photographs indicates scale.

보경사 지역에서 확인되는 양산단층대는 서쪽지괴의 유천층군인 안산암질 응회암과 동쪽지괴의 하양층군 퇴적암을 가로지르고 있다(Kim et al., 1968; Yun et al., 2000). 주단층면의 자세는 N25°E/70°NW이며, 양산단층의 주향방향과 유사하다. 주단층면에서 확인되는 단층비지는 주로 경상분지 백악기 퇴적암 기원의 단층암이다(Kim et al., 2015). 시료는 북북동방향의 주단층면을 따라 확인되는 단층비지(BGS 3-1)와 퇴적암 기원의 단층비지를 가지는 북동방향의 부차단층(BGS 3-2)에서 채취하였다(그림 2c, 2d).

야외 단층 노두에서 채취한 단층비지를 실내 전단실험에 사용하기 위해 전처리 과정을 거쳐 분말형태로 만들었다. 먼저 채취한 단층비지를 #60체(250 μm)에 놓고 증류수를 부으면서 거른다. #60체를 통과한 시료들을 60℃ 오븐 내에서 건조한 후, 뭉쳐있는 단층비지 분말 덩어리를 손으로 분쇄하여 다시 #60체로 걸러 최종 단층비지 분말을 제작한다. 최종단층비지 분말은 XRD 정량분석(표 1), 입도분석(그림 3), 그리고 전단실험에 이용하였다.

Fig. 3. 
Grain size distribution of the fault gouges. The SCL 1-1 has a similar grain size distribution to SCL 1-2, and it appears they consist of relatively larger grains. Both SCL 1-3 and BGS 3-2 have a similar grain size distribution and have relatively smaller grains. The BGS 3-1 has a relatively narrower grain size distribution compared to that of other samples.

이번 연구에서는 단층비지 내의 주된 점토광물인 일라이트(illite), 녹니석(chlorite) 이외에 몬모릴로나이트(montmorillonite)를 식별하기 위해 별도의 전처리(즉, 에틸렌글리콜 처리) 과정을 실시하였고, 절차는 다음과 같다. 먼저 전체 시료에 대해 XRD 정량분석 후, 몬모릴로나이트의 피크가 인지되는 시료들을 따로 선별하였다. 이후 선별한 시료에서 점토광물의 층간을 최대로 팽창시키는 에틸렌글리콜 처리 후, 정성분석을 통해 몬모릴로나이트의 존재 여부를 식별하였다. 몬모릴로나이트의 존재가 확인된 시료에 대해서만 XRD 정량분석 시, 몬모릴로나이트를 추가로 입력하여 산출하였다.

에틸렌글리콜 처리 방법은 앞서 설명한 처리과정으로 제작한 분말시료를 물과 섞어 현탁액을 만든 후, 원심분리를 통해 비커의 맨 위에 떠다니는 점토광물을 따로 건져내었다. 이후, 슬라이드글라스에 올려서 건조시킨 후, 에틸렌글리콜을 한방울 떨어뜨리고 다시 5시간 건조시켰다.

XRD 정량분석결과는 표 1에 제시하였다. SCL 1-1는 석영 함량이 51.4 wt%로 다른 시료들에 비해 가장 높은 함량으로 산출되었고, 장석 함량과 합하여 가장 높은 비점토광물 함량을 보였다. 점토광물의 총 함량은 35.6 wt%로 가장 낮게 산출되었다. SCL 1-2는 석영 함량이 SCL 1-1보다 적지만, 장석을 많이 함유하여, 총 비점토광물과 점토광물 함량은 SCL 1-1과 유사하게 각각 62.5 wt%와 37.5 wt%로 산출되었다. SCL 1-3의 경우, 비점토광물 함량은 51.2 wt%로 앞선 시료들보다 적지만, 총 점토광물 함량은 48.8 wt%로 상대적으로 많이 산출되었다. BGS 3-1의 경우, 장석 함량은 3 wt%미만으로 거의 산출되지 않았지만 석영이 43.0 wt%로 산출되었다. 총 비점토광물 함량은 52.6 wt%이고, 점토광물의 총 함량은 47.4 wt%로 산출되었으며, 특징적으로 몬모릴로나이트가 9.5 wt%로 산출되었다. BGS 3-2는 BGS 3-1과 비점토광물과 점토광물의 총 함량은 비슷하지만, 몬모릴로나이트는 존재하지 않았다. 따라서 XRD 정량분석에 근거할 때, 광물조성이 유사한 시료들은 총 비점토광물 함량과 점토광물 함량이 유사한 상천리 주 단층비지대(SCL 1-1, 1-2)에서 채취한 두 시료와 상천리와 보경사 인근의 부차 단층(SCL 1-3과 BGS 3-2)에서 채취한 두 시료로 구분 가능하다.

입도분석은 Malvern 사의 Mastersizer 2000 장비를 이용하여 실시하였고, 그 결과를 입도크기와 부피 비율의 함수로 그림 3에 제시하였다. 그래프에서 정점(peak)의 X축 위치는 가장 높은 비율을 차지하는 평균 입자의 크기를 의미하고, Y축 위치는 평균 입자들이 차지하고 있는 비율을 나타낸다.

상천리에서 채취한 시료들(SCL 1-1, 1-2, 1-3)은 각각 0.025, 0.018 그리고 0.004 mm의 평균 입자크기를 가지며, 시료마다 정점의 양 옆으로 분산되어 있는 것으로 다양한 크기의 입자들이 분포함을 보여준다. 이는 SCL 1-3에서 가장 현저하게 나타나며, 분급이 좋지 않음을 의미한다. 보경사 인근에서 채취한 시료(BGS 3-1, 3-2)는 각각 0.004, 0.003 mm의 평균 입자크기 갖는 것으로 나타났다. BGS 3-1의 경우 다른 시료들과 비교했을 때, 정점 주변으로 집중되는 경향이 가장 현저하게 나타나며, 이는 평균입자크기로 입자들이 균질하게 분포하여 분급이 좋음을 의미한다. 따라서 다른 시료들에 비해 BGS 3-1은 상대적으로 세립의 입자들이 높은 비율을 차지하고, 분급도 우세하다. 반면, BGS 3-2의 평균 입자크기는 BGS 3-1과 유사하지만, 정점에 집중되는 경향을 보이지 않기 때문에 불량한 분급을 갖고, SCL 1-3과 유사한 입도분포를 보인다.

전단실험에 사용되는 단층비지 분말은 그림 4에서 보는 바와 같이 암석 블록(측면블록, 4×5×2 cm, 2개, 중앙블록, 4×8×2 cm, 1개) 사이에 2 mm 두께의 샌드위치 형태로 채워 넣어 샘플 조립체(sample as-sembly)를 만들었다. 암석 블록으로 제작한 재료는 암석역학 실내시험에서 보편적으로 이용되는 표준시료인 Berea 사암 또는 전북 황등 지역에서 산출되는 석영, 사장석, K-장석 및 흑운모가 주 구성광물인 세립질의 황등화강암을 이용하였다. 암석 블록의 슬립면(slip surface)은 샘플결합에 앞서 실리콘 카바이드(# 80)로 연마하였다. 전단면의 연마는 암석 블록의 표면을 거칠게 만듦으로써, 단층의 전단면에 형성되는 굴곡면(asperity)을 재현하여, 단층비지 내부에서만의 미끌림을 유도하기 위함이다(Marone et al., 1990). 또한 단층비지의 손실을 막기 위해, 암석블록의 접합부 주변에 테이핑을 하여, 24시간 동안 증류수에 침수시켜 단층비지 분말을 포화 상태로 만들었다.

Fig. 4. 
Schematic diagram of biaxial loading apparatus consisting of a vertical loading frame for applying shear force and a horizontal frame for normal force. The inset shows a magnified sample assembly which is made up by sandwiching a gouge layer of 2 mm thickness in between the center forcing block (4 cm × 8 cm × 2 cm) and the side forcing block (4 cm × 5 cm × 2 cm). For running a double direct shear test, the sample assembly is positioned just below the loading piston equipped under the shear load cell.

그림 4에 삽입된 모식도와 같이 전단실험 동안 옆으로 단층비지가 새어 나오는 것을 최소화하기 위해, 마찰력에 영향을 주지 않는 알루미늄으로 제작한 덮개판을 부착하였다. 이후 전단면에 수직응력을 가하는데 이용된 수평으로 놓인 로딩프레임(loading frame)에 장착하였다. 또한 전단실험 동안 수직 및 전단변위(normal and shear displacement)를 측정하기 위해 LVDT (linear variable differential trans-former)를 장착하였다.

이후 일정한 가력 속도(loading rate, 0.01 tonf/s)로 전단면에 수직응력(normal stress)을 증가시키고, 정해진 수직응력에 도달한 후, 서보-컨트롤러(servo-con-troller)의 제어를 통해 전단실험을 마칠 때까지 일정하게 유지하였다. 실험에 사용된 수직응력은 5, 10, 20 MPa이다. 이후 수직방향으로 놓인 로딩프레임을 이용하여 전단변위속도를 제어, 변화시키면서 전단면에 전단력(shear force)을 가하였다(velocity stepping test). 수직응력 5, 10 MPa에서는 전단변위속도(1-10-100-10-1 μm/s)를 이용하였고, 수직응력 20 MPa에서는 전단변위속도(1-3-10-30-60-100 μm/s)를 이용하였다. 각 실험에서 정해진 전단변위속도를 총 4회 반복하면서 전단실험을 실시하였다. 전단실험 동안의 수직, 수평 방향의 응력과 변위를 80-100 Hz로 기록하였다.

실내 전단실험을 통해 얻은 자료로부터, 마찰거동 특성을 지시하는 파라미터(constitutive parameter)를 산출하기 위해서 Dieterich (1978, 1979)와 Ruina (1983)가 제안한 Rate and state friction law (수식 1)를 적용하였다. Rate and state friction law는 (1) 전단실험 동안 전단 속도의 변화에 따라 마찰계수도 동시에 변화하는 효과인 direct effect (rate)와 (2) direct effect 이후, 변화하는 마찰계수를 나타내는 효과인 evolution effect (state)의 조합으로 표현된다(그림 5a; 수식 1).

Fig. 5. 
(a) Behavior of sliding friction coefficient in response to sliding velocity increase. “a” and “b” are experimental constants. “a” represents the magnitude of the instantaneous increase in the sliding friction coefficient. “b” represents the magnitude of the subsequent decay in the sliding friction coefficient. D_c is the critical distance over which sliding friction coefficient evolves from one steady state value to another. Note that the discrepancy between velocity weakening (indicated by solid line) and velocity strengthening (indicated by dotted line). Velocity weakening (a-b<0) indicates unstable sliding, and velocity strengthening (a-b>0) indicates stable sliding (modified from Samuelson, 2010). (b) Enlargement of a velocity step test (SCL 1-3) showing the variation of sliding friction coefficient in response to an upstep velocity and the accompanying dilation of the gouge layer (ΔH).

수식 1에서 μ₀ 는 전단 속도가 변화하기 전의 마찰계수이고, V₀는 전단 속도가 변화하기 전의 속도, V는 전단 속도가 변화한 후의 속도이다. D_c는 임계전단거리(critical slip distance)로 전단속도 변화 후, 마찰계수가 새로운 평형상태에 도달하는데 필요한 거리를 의미한다. θ는 상태변수(state valuable)로써 임계전단거리(D_c)와 속도의 함수로 표현된다. 따라서 θ는 임계전단거리 동안 어떤 속도로 움직인 평균 이동 시간(average time of asperity contact, Dieterich, 1979)을 의미한다. “a” 와 “b” 는 전단실험을 통해 산출되는 상수(experimental constants)이다. 전단속도의 제어를 통한 속도단계실험(velocity stepping test)에서 전단속도의 증가 혹은 감소와 동시에 마찰계수가 증가 혹은 감소하는 direct effect의 정도를 나타내는 실험상수는 “a” 이고, D_c를 이동하는 동안 마찰계수가 감소 혹은 증가 하는 evolution effect의 정도를 나타내는 실험상수는 “b” 이다(그림 5a).

또한 “a” 와 “b” 는 단층의 안정성 판별에 이용된다. “velocity weakening”, (a-b<0)인 경우 파괴 후 단층은 불안정한 지진성 전단거동을 유발하게 된다. “velocity strengthening”, (a-b>0)인 경우, 단층은 안정한 비지진성(aseismic) 전단거동을 보인다. 위에서 설명한 파라미터들과 더불어 전단실험 동안의 단층비지 두께변화를 측정하여, 산출할 수 있는 ΔH가 있다(그림 5b). ΔH는 전단속도 변화에 따라 샘플 조립체 사이에 끼운 단층비지 체적의 증가(dilatancy) 또는 감소(compaction)현상을 표현한 것이다. 이는 전단실험 동안, 단층비지 두께변화로 단층비지 내부에서의 변형 정도를 지시한다. 이러한 파라미터를 구하기 위해 전단실험을 마친 데이터에서 속도변화가 일어난 구간을 따로 선별하여, 곡선 맞춤(curve fitting)을 하여 계수를 산출하였다.

전단실험 중 변형된 단층비지의 미구조 분석을 위해 실험 후 각 샘플조립체를 회수하고 이를 슬라이드글라스 크기로 절단하여 연마 박편(polished thin section)을 제작하였다. 이후 전계방출형 주사전자현미경(field-emission scanning electron mi-croscope) 하에서 단층비지를 관찰하였다. 관찰한 연마 박편 시료는 모두 수직응력 5 MPa 조건에서 최종 전단변위가 20-25 mm에 도달할 때까지 변형된 것들이다.

자연 단층비지를 이용한 실내 전단실험을 통해 얻어진 마찰 거동은 그림 6과 그림 7에 제시되어 있고, 각 실험에서 얻어진 최대마찰계수(peak friction co-efficient)는 표 2에 정리되어 있다. 실험 결과, 수직응력(5, 10, 20 MPa)에 따라 다른 마찰거동이 관찰되었다. 상천리 인근에서 채취한 단층비지 시료의 경우, 대체적으로 수직응력이 5, 10 MPa 조건에서 모든 시료(SCL 1-1, 1-2, 1-3)는 “slip weakening” 거동을 보였다. 즉 마찰계수가 최대마찰강도까지 급격하게 증가한 뒤, 점진적으로 마찰계수가 감소하여 새롭게 안정된 정상마찰계수에 도달하게 된다. 한편 수직응력이 20 MPa 조건에서는 마찰계수가 최대마찰강도에 도달한 이후, 전단변위가 증가함에 따라 계속 증가하는 “slip hardening” 거동을 보였다(그림 6).

Fig. 6. 
Coefficient of sliding friction as a function of shear displacement for Sangcheon-ri fault gouges, (a) SCL 1-1, and (b) SCL 1-2, and (c) SCL 1-3, tested at σ_n=5, 10, and 20 MPa and with varying shear velocity.

Fig. 7. 
Coefficient of sliding friction as a function of shear displacement for Bogyeongsa fault gouges, (a) BGS 3-1 and (b) BGS 3-2, tested at σ_n=5, 10, and 20 MPa and with varying shear velocity.

한편, 보경사 인근 단층노두에서 채취한 BGS 3-1의 경우 수직응력 5, 10, 20 MPa 조건에서 모두 “slip hardening” 거동을 보였다(그림 7a). 하지만, BGS 3-2의 경우 상천리 시료들과 유사하게 수직응력 5, 10 MPa 조건에서는 “slip weakening” 거동을 보였고, 수직응력 20 MPa 조건에서는 “slip hard-ening” 거동을 보였다(그림 7b).

최대마찰계수(peak friction coefficient)는 그림 8과 표 2에 제시하였다. 그림 8에서 보는 바와 같이 전반적으로 수직응력이 클수록 최대마찰계수가 작아지는 경향을 보인다. 단층핵으로부터 채취한 상천리 시료(SCL 1-1, 1-2)의 경우, 실험에 이용된 모든 수직응력(5, 10, 20 MPa) 하에서 최대마찰계수(μpeak)는 0.284-0.364의 범위로 측정되었고, 두 시료의 마찰거동과 마찰계수는 유사한 것으로 나타났다(표 2). 반면에, 단층손상대내 부차단층으로부터 채취한 SCL 1-3은 수직응력 5, 10 MPa에서 진행된 전단실험에서 최대마찰계수가 가장 높게 측정되었다. 같은 수직응력 조건(5-20 MPa)에서 보경사 인근 단층노두의 단층핵으로부터 채취한 시료 BGS 3-1은 모든 시료 중에서 가장 낮은 최대마찰계수가 측정되었고, 단층손상대내 부차단층으로부터 채취한 BGS 3-2는 SCL 1-1, 1-2와 유사한 범위로 측정되었다(그림 8; 표 2).

Fig. 8. 
Peak friction coefficient versus normal stress. Note that the larger normal stress is applied, the lower peak friction coefficient is measured. The lowest peak friction coefficient of 0.2 was observed in BGS 3-1. coefficients obtained from each of the tests (SCL: Sangcheon-ri, BGS: Bogyeongsa).

일반적으로 단층비지 내에서 총 입자(total clast)의 함량이 많아지고, 총 점토광물(total clay)의 함량이 적어지면 마찰계수는 증가한다(Brown et al., 2003; Tembe, 2007). XRD 정량분석결과와 전단실험결과를 비교해보면, SCL 1-1과 SCL 1-2의 경우, 총 비점토광물의 함량과 점토광물의 함량이 서로 유사하게 산출되었다(표 1). 따라서 SCL 1-1과 SCL 1-2는 유사한 마찰계수를 가질 것으로 예상되며, 전단실험결과에서도 마찰계수가 큰 차이를 보이지 않았다. SCL 1-3의 경우, 총 비점토광물의 함량이 SCL 1-1과 SCL 1-2 보다 약 10 wt% 적고, 총 점토광물의 함량은 약 10 wt% 많게 산출되었다. 따라서 단층암 구성광물에 따른 마찰계수의 일반적인 경향으로 볼 때, SCL 1-3에서 두 시료(SCL 1-1, 1-2)보다 더 낮은 마찰계수가 측정되어야 한다. 하지만, 전단실험결과 SCL 1-3의 마찰계수는 모든 시료 중 가장 높은 값을 보였다(표 2). 이를 규명하기 위해서는 추가적인 연구가 필요하다.

BGS 3-1과 BGS 3-2의 경우, 비점토광물의 총 함량이 SCL 1-1과 SCL 1-2 보다 약 10 wt% 낮게 산출되었고, 점토광물의 총 함량은 약 10 wt% 높게 산출되었다(표 1). 단층암 구성광물에 따른 마찰계수의 일반적인 경향으로 볼 때, BGS 3-1과 BGS 3-2은 SCL 1-1과 SCL 1-2보다 마찰계수가 낮게 측정되어야 한다. 하지만 BGS 3-1은 낮게 측정된 반면, BGS 3-2는 SCL 1-1과 SCL 1-2와 유사하게 측정되었다. 따라서 이번 연구에서 XRD 정량분석에서 산출된 총 비점토광물의 함량과 총 점토광물의 함량에서 10 wt% 차이는 마찰계수에 큰 영향을 미치지 않는 것으로 보인다. 반면, BGS 3-2와 유사한 총 조암 및 점토광물 함량을 보인 BGS 3-1은 모든 시료 중에서 가장 낮은 마찰계수가 측정되었다. 이는 몬모릴로나이트 함량에 기인한 것으로 보인다(표 1).

XRD 정량분석에서 BGS 3-1은 몬모릴로나이트가 9.5 wt%로 가장 많이 산출되었다. SCL 1-2에서 2.1 wt%로 산출된 것을 제외하고, 다른 모든 시료들(SCL 1-1, SCL 1-3 그리고 BGS 3-2)에는 몬모릴로나이트가 없는 것으로 나타났다.

두 지역(SCL, BGS)의 시료 모두 기질부(matrix)는 주로 판상의 점토광물(일라이트, 녹니석)로 이루어져 있고, 입자들은 대부분 XRD 정량분석결과와 일치하는 석영, 장석 등으로 관찰되었다. 미구조분석은 전단방향에 따라 발달하는 전단띠 형태로 구분하였고, 명명법(Logan et al., 1992)과 전단방향은 그림 9에 삽입된 모식도와 같다.

Fig. 9. 
BSE photograph of the samples subjected to σ_n=5 MPa and deformed up to a shear displacement of 20-25 mm. (a) SCL 1-1, (b) SCL 1-3, (c) BGS 3-1, and (d) BGS 3-2. The sense of shear is indicated by the arrows on each of the photographs. A schematic diagram (following Logan et al., 1992) shows typical shear localizations formed in a deformed gouge layer. Following distributed shear, evolution progress form R₁ (Riedel shear), to B (boundary shear), and Y (boundary-parallel shears). Note that BGS 3-1 shows a noticeable development of R₁ and Y-shear compared to other samples.

SCL 1-1의 경우, 10-100 μm 크기의 입자들이 관찰되고, 주로 20 μm 크기의 입자가 주로 분포한다. 전단띠의 발달 정도는 미미하며, R₁-전단띠가 발달한 것이 관찰되었다(그림 9a). SCL 1-3은 5-200 μm 정도의 입자들이 불량한 분급을 갖는 것으로 관찰되며, 30-50 μm 크기의 입자가 주로 관찰된다. 또한 단층비지 전면에 걸쳐 입자들이 산발적으로 분포하는 것이 현저하게 관찰된다. 전단띠의 발달은 미미하고, R₁-전단띠가 발달한 것이 관찰되었다(그림 9b).

반면, BGS 3-1의 경우, 다른 시료들에 비해 상대적으로 단층비지 내에 입자들이 차지하는 비율이 뚜렷하게 적으며, 좋은 분급을 보이는 주로 5 μm 크기의 입자들이 관찰된다. 또한 다른 시료들에 비해 현저한 전단띠의 발달이 관찰된다. 주로 R₁-전단 방향의 얇고 짧은 다수의 전단띠가 전단면 전면에 걸쳐 발달해 있는 것이 특징이다(그림 9c). 또한 R₁-전단띠의 발달 이후 계속적인 변형으로 서로 연결되어 전단방향과 거의 평행한 Y-전단띠의 발달도 확인된다. BGS 3-2의 경우 입도분석 결과에서 SCL 1-3과 유사한 입자크기분포를 보였지만, 입자크기가 10-20 μm 정도의 입자들이 주로 분포하는 것이 관찰되며, R₁-전단띠의 발달이 관찰되었다(그림 9d).

다른 시료들에 비해 상대적으로 R₁-전단띠와 Y-전단띠의 발달이 현저하게 관찰되는 BGS 3-1의 경우, 기질부를 이루고 있는 점토광물이 전단변위 증가와 함께 전단띠 형성에 용이한 방향으로 배열되었음을 의미한다. 또한 단층비지에 존재하는 포화된 몬모릴로나이트도 전단띠를 발달시키는 과정을 더욱 심화시켰을 것으로 판단된다. 몬모릴로나이트는 물에 노출될 경우, 층간(inter-layer) 사이에 물을 함유하게 됨으로써 팽창되는 광물(expendable min-eral)이다. 몬모릴로나이트의 층간에서 흡착된 물은 얇은 막을 형성하여 윤활제 역할을 한다. 이러한 물성 때문에 단층면에 전단력이 가해질 때, 단층비지 내에서 R₁-전단띠와 Y-전단띠과 같은 변형구조의 발달을 촉진시켰을 것으로 해석된다. 즉, 단층비지 내부에서 마찰력의 약화를 유발하여, 마찰계수를 낮추게 된다(Summers and Byerlee, 1977; Moore and Lockner, 2004).

그러므로 BGS 3-1의 상대적으로 낮은 마찰계수는 다음과 같이 설명 가능하다. 단층비지 내에 존재하는 물에 포화된 몬모릴로나이트는 전단변위가 증가함에 따라 전단에 유리한 방향으로 배열되는 것을 심화시켜 마찰계수를 낮추는 역할을 한 것으로 보인다. 즉, 마찰계수가 낮다는 것은 단층비지가 잘 미끄러진다는 것을 의미하고, 이는 곧 전단띠의 발달이 잘 이루어 졌음을 지시한다. 미구조 관찰결과 역시 이를 잘 반영하고 있으며, BGS 3-1은 다른 시료들에 비해 전단띠의 발달이 현저하게 관찰된다(그림 9c).

마찰율 파라미터(friction rate parameter, 이하 a-b)는 기존 단층이 재활성(즉, 파괴) 된 후에 발생하는 슬립거동의 안정성을 해석하는데 이용된다. a-b 값은 전단속도변화(velocity stepping)에 따른 마찰계수의 변화를 실험적으로 측정함으로써 산출된다(그림 5a). 만약 전단속도를 증가시켰을 때 발생하는 마찰계수가 속도변화 전의 마찰계수(μ₀)보다 높으면, 즉 a-b>0 (velocity strengthening), 그 단층은 재활성 된 후 안정한 슬립거동(stable slip or aseis-mic slip)을 보인다. 하지만 전단속도를 증가시켰을 때 발생하는 마찰계수가 속도변화 전의 마찰계수(μ₀)보다 낮아진다면, 즉 a-b<0 (velocity weakening), 그 단층은 재활성 되어 지진을 유발할 수 있는 불안정한 거동(seismic slip)을 보일 가능성이 있다.

이번 연구에서 각 단층비지의 a-b 값을 산출하기 위해 전단속도변화 실험(velocity stepping test)을 실시하였으며, 그 결과는 그림 10에 제시되어 있다. 수직응력 5 MPa 하에서 BGS 3-1에서 음의 a-b 값이 관찰된 것을 제외하고, 다른 모든 시료 및 실험조건 하에서 양의 a-b 값이 관찰되었다. 즉, 시료를 채취한 단층 분절(segment)에서 단층의 재활이 발생하더라도 지진을 유발하기보다는 안정한 마찰거동을 보일 것으로 판단된다.

Fig. 10. 
Friction rate parameter (a-b) as a function of shear displacement for the normal stresses, (a) σ_n=5 MPa, (b) σ_n=10 MPa, and (c) σ_n=20 MPa. Note that besides the case of BGS 3-1 at σ_n=5 MPa, friction rate parameters measured are positive for all the samples and test conditions, suggesting frictional slip after the fault fails exhibits an aseismic (or stable) slip. Additionally, except the case of σ_n=20 MPa, the friction rate parameter decreases with shear displacement.

또한 전단변위 증가에 따른 마찰율 파라미터의 변화와 관련하여, 수직응력이 5, 10 MPa의 경우 모든 시료에서 전단변위가 증가함에 따라 a-b 값이 감소하는 경향이 관찰(그림 10a, 10b)되었다. 한편 수직응력이 20 MPa의 경우, 전단변위가 증가함에 따라 오직 SCL 1-2와 BGS 3-1에서만 a-b 값이 감소하는 경향이 관찰된다(그림 10c). 전반적으로 전단변위가 증가하면서 a-b 값은 점진적으로 감소하는 경향이 관찰된다. 이는 기존 연구결과와 일치한다. a-b 값은 전단변위가 증가할수록 감소하는 경향을 보이는 것으로 알려져 있고, 이는 전단변위가 지속적으로 단층비지에 축적되면서 전단띠가 발달하는 것과 관련이 있다(Marone, 1998; Mair and Marone, 1999; Saffer and Marone, 2003; Numelin et al., 2007; Ikari et al., 2009, 2011). 즉 이번 연구에서 실험한 시료들은 전반적으로 단층비지 내에 전단띠의 발달이 이루어졌음을 지시하고, 이는 미구조분석에서도 R₁-과 Y-전단띠의 형태로 관찰된다.

한편 이상의 결과는 이산화탄소 지중저장과 관련하여 다음과 같은 시사점이 있다. 지하 저장층 내로 CO₂ 주입은 저장층 내 공극수압의 상승을 초래할 뿐만 아니라 저장층 주변에 존재하는 단층대 내에 공극수압의 상승을 유발한다. 공극수압이 임계압 이상으로 상승할 경우, 단층의 재활성을 초래 할 수 있으며, 이는 지진을 수반할 수 있다. 따라서 지하에서의 단층암이 지표 노두에서의 단층암과 특성이 유사하다면, 이번 연구에서 관찰된 양의 a-b 값은 단층비지를 채취한 양산단층 분절에서 CO₂ 지중주입에 따른 공극수압 상승으로 인한 단층 재활성이 발생되더라도, 지진을 유발하기 보다는 안정한 슬립거동을 보일 것이라는 것을 지시한다.