A Novel Application of Strain Energy to Estimate Failure Pattern of Gneiss Rock using Genetic Algorithm and Triaxial System

기술 요약: 유전 알고리즘과 삼축 압축 시스템을 이용한 편마암의 파괴 패턴 추정을 위한 변형 에너지의 새로운 적용

문제 정의
인도 우타라칸드주 찰몰리 지역의 편마암과 같은 변성암은 광물 엽리(foliation)로 인해 뚜렷한 구조적 이방성을 나타낸다. 이러한 이방성은 결정론적 접근 방식에 의존하는 기존의 강도 예측을 복잡하게 만드는데, 기존 방식은 흔히 응력-변형률 관계에만 의존한다. 이러한 전통적인 방법들은 비선형적이고 규모 의존적인 암석 거동과 엽리 방향, 응력 조건, 그리고 에너지 진화 사이의 복잡한 결합을 포착하는 데 자주 실패한다. 결과적으로, 재해 발생 가능성이 높은 엽리 발달 암반의 안정성을 평가하기 위해 에너지 기반 파괴 이론과 진화적 최적화를 통합하는 예측 프레임워크가 필요하다.

연구 방법론
본 연구는 통합된 실험 및 계산적 접근 방식을 채택하였다:

• 시료 채취 및 준비: 인도 찰몰리의 파갈 날라(Pagal Nala) 산사태 지점에서 편마암 암괴를 채취하였다. 원통형 코어 시료는 길이 대 직경 비율이 2가 되도록 추출되었으며, 7개의 특정 엽리 각도(0°, 15°, 30°, 45°, 60°, 75°, 90°)에 맞춰 배향되었다.
• 물리적 및 암석학적 특성 분석: 물리적 성질(단위 중량, 공극률, 흡수율)을 측정하였다. 광물 조성은 X선 회절법(XRD)을 통해 분석되었으며, 방해석, 석영, 장석가 주요 광물임을 확인하였다. 지열-압력계(Geothermobarometry)를 사용하여 암석 형성의 최고 변성 조건을 추정하였다.
• 역학적 시험:
  • 초음파 펄스 속도(UPV): P파 및 S파 속도를 측정하여 동적 탄성 계수(영률, 전단 계수, 부피 탄성 계수)를 결정하고 방향성 강성을 평가하였다.
  • 일축 압축 시험: 표준 일축 압축 강도(UCS) 시험을 수행하였다.
  • 단일 사이클 하중-제거(SC-UCS): 시편을 최대 강도의 95%까지 하중을 가한 후 제거하였다. 이를 통해 총 변형 에너지($U$), 탄성 변형 에너지($U_e$), 소산 변형 에너지($U_d$)를 정량화할 수 있었다. 에너지 저장 지수($W_{et} = U_e / U_d$)를 계산하였다.
  • 삼축 압축 시험: 1, 2, 3 MPa의 구속 압력 하에서 시험을 수행하여 전단 강도 파라미터(점착력 $c$ 및 내부 마찰각 $\phi$)를 결정하고 파괴 모드를 관찰하였다.
• 유전 알고리즘(GA) 모델링: 암석 강도와 전단 파라미터를 최적화하고 예측하기 위해 MATLAB을 사용하여 GA 기반 모델을 개발하였다. 이 모델은 다양한 엽리 각도에 따른 지표 특성, 동적 계수, 변형 에너지 상수를 포함한 입력 파라미터를 활용하였다. 모델은 실험 데이터에 훈련되어 예측 오차를 최소와하고 지배적인 파괴 패턴을 식별하였다.

주요 결과

• 이방성 및 탄성 특성: 찰몰리 편마암은 상당한 탄성 이방성을 보였다. P파 속도는 90°일 때와 비교하여 30° 엽리 각도에서 27.83% 감소하였는데, 이는 약한 광물 경계를 따른 파동 산란에 기인한다. 동적 영률은 방향에 따라 약 39 GPa에서 66 GPa까지 변화하였다.
• 전단 강도 및 파괴 모드: 전단 강도는 엽리 방향에 크게 의존하였다. 30° 방향은 가장 낮은 전단 강도를 보였으며 90° 방향 대비 가장 큰 강도 감소를 나타내어, 운모가 풍부한 엽리면을 따라 미끄러짐이 지배적인 임계 약면임을 확인하였다. 파괴 모드는 90°에서의 축 방향 분열 및 Y자형 파괴에서 중간 각도에서의 단일 평면 전단 또는 이중 전단으로 변화하였다.
• 에너지 진화: 엽리에 수직으로 하중을 받은 시편(90°)은 최대 총 변형 에너지와 탄성 변형 에너지를 저장한 반면, 30° 시편은 가장 적게 저장하였다. 90° 방향은 강한 석영-장석 밴드를 통한 인장 균열 전파를 통해 에너지 저장을 촉진한 반면, 중간 각도에서는 더 높은 에너지 소산을 동반한 전단 지배적 파괴를 촉진하였다. $W_{et}$ 비율은 주로 취성 파괴 반응($W_{et} > 1.3$)을 나타냈으나, 격렬한 암반 폭발(rockburst) 임계값($W_{et} < 2$) 미만으로 유지되었다.
• 모델 성능: GA 기반 모델은 응력에 대해 75.29%, 전단 파라미터에 대해 66.89%의 예측 일치도를 달了성하였다. 모델은 일반적인 이방성 경향(30°에서 강도 최소, 90°에서 강도 최대)을 성공적으로 재현하였으나, 30° 방향에 대한 일치도는 낮았는데, 이는 해당 임계 각도에서 발생하는 복잡한 불연속면 제어 미끄러짐 메커니즘 때문인 것으로 보인다.

의의 및 주장
저자들은 본 연구가 재해 발생 가능성이 높은 지형에서 암반 안정성 평가를 위한 새롭고 빠르며 예측 가능한 의사결정 지원 도구를 제공한다고 주장한다. 광물학적 특성 분석, 에너지 기반 파괴 분석, 유전 알고리즘 최적화를 통합함으로써, 본 프레임워크는 전통적인 응력-변형률 분석을 넘어 이방성 강도와 에너지 진화를 정량화하는 방법을 제시한다.

본 연구는 데이터 기반 지반 설계의 기초로서 엽리 각도, 강도 이방성, 그리고 에너지 소산을 연결하여 파괴 패턴을 예측하는 위치를 점한다. 저자들은 이 모델이 예비 스크리닝 도구 역할을 한다는 점을 겸허히 언급하며, 향후 실험 데이터셋(특히 임계 엽리 각도 주변)을 확장하고 하이브리드 머신 러닝 접근 방식이나 방향 의존적 가중치를 통합함으로써 정확도를 높일 수 있다고 명시하였다. 본 연구는 터널 공학 및 사면 안정성과 같은 공학적 응용 분야에서 변성암 암반 평가를 위한 전통적인 암석 시험과 지능형 모델링 사이의 간극을 메우는 것을 목표로 한다.