Merge remote-tracking branch 'origin/안티그래비티' into main_laptop_1

resources/knowledge/original/안정성평가/05_3D_사면안정해석_논문_기술정보.md

@@ -0,0 +1,252 @@
+# 3D 사면안정해석 논문·기술정보 조사 (3D LEM 기둥법 · 3D SRM · Scoops3D · 벤치마크 · 탐색 알고리즘)
+
+> **분류**: 안정성평가(해석) 영역 — 3차원 사면안정해석 이론·수식·검증 벤치마크·오픈소스 조사
+> **작성 원칙**: 모든 항목에 출처(DOI·URL), 핵심 수식·가정, Aislo 데이터(LiDAR DEM · 측점 20 m 횡단 · 3D 노선 계획면) 적합성 판정 한 줄을 병기함. 값은 확정하지 않고 문헌 근거와 후보만 제시함. 서술은 명사형(~함 · ~임)으로 통일함.
+
+---
+
+## 1. 3D 한계평형(LEM) 기둥법 (Method of Columns) 수식·가정·비교
+
+### 1.1 3D 기둥(Column) 이산화 공통 기하학
+- **기하 정의**: 2D 절편법(Method of Slices)의 수직 절편을 3차원 수직 사각기둥(Column, 격자 크기 $\Delta x \times \Delta y$)으로 확장함. 활동 방향을 $x$(또는 회전 방위각 $\alpha_{\text{slide}}$), 횡방향을 $y$, 연직 상방을 $z$라 정의함.
+- **기둥 바닥 경사각 및 면적**:
+  - $(i, j)$번째 기둥 바닥면의 $x$방향 겉보기 경사각을 $\alpha_{x}$, $y$방향 겉보기 경사각을 $\alpha_{y}$라 할 때, 바닥면의 실제 최대 경사각(True dip) $\psi$ 및 바닥 경사 면적 $A$는 다음과 같음:
+    $$\cos\psi = \frac{1}{\sqrt{1 + \tan^2\alpha_x + \tan^2\alpha_y}}, \qquad A = \frac{\Delta x \cdot \Delta y}{\cos\psi}$$
+  - 활동 방향($x$축)과 기둥 바닥이 이루는 활동방향 겉보기 경사각 $\varepsilon$ (2D 절편법의 $\alpha_i$에 대응):
+    $$\tan\varepsilon = \tan\alpha_x$$
+
+---
+
+### 1.2 방법별 수식·가정·출처 상세
+
+#### ① Hovland (1977) — 3D 단순기둥법 (3D Ordinary / Fellenius Method)
+- **출처**: Hovland, H. J. (1977), "Three-Dimensional Slope Stability Analysis Method", *Journal of the Geotechnical Engineering Division (ASCE)*, 103(GT9), 971–986.
+  - **DOI**: `10.1061/AJGEB6.0000493` / URL: `https://doi.org/10.1061/AJGEB6.0000493`
+- **핵심 가정**: 기둥과 기둥 사이의 모든 측면 수직력 및 전단력(Intercolumn normal and shear forces)을 $0$으로 무시함(2D Fellenius 법의 3D 확장임).
+- **핵심 수식 (회전 원통/구형 활동면 모멘트 평형 — USGS Scoops3D 채택형)**:
+  $$F_{\text{3D,Hovland}} = \frac{\displaystyle\sum_{i,j} R_{i,j} \left[ c'_{i,j} A_{i,j} + \left( W_{i,j} \cos\psi_{i,j} - u_{i,j} A_{i,j} \right) \tan\phi'_{i,j} \right]}{\displaystyle\sum_{i,j} W_{i,j} \left( R_{i,j} \sin\varepsilon_{i,j} + k_h e_{i,j} \right)}$$
+  - $W_{i,j}$ : 기둥 자중 ($\text{kN}$), $u_{i,j}$ : 바닥 간극수압 ($\text{kPa}$), $R_{i,j}$ : 회전축에서 기둥 바닥 중심까지의 반경 팔길이 ($\text{m}$), $e_{i,j}$ : 회전축에서 기둥 무게중심까지의 연직 모멘트 팔길이 ($\text{m}$), $k_h$ : 수평지진계수임.
+- **특징**: $F_s$가 우변에 없어 반복 계산 없이 1회 연산으로 즉시 산출됨. 엄밀해 대비 5~15% 보수적(낮은 안전율) 결과를 보이며, 3D Bishop 법의 초기 추정치($F_{s,0}$) 및 대량 격자 고속 스크리닝용으로 쓰임.
+- **우리 데이터 적합성**: **LiDAR DEM 격자($\Delta x \times \Delta y$)를 그대로 기둥으로 매핑하여 수십만 개 후보 구면을 무반복 고속 연산하는 1차 스크리닝 엔진으로 적합함.**
+
+#### ② Hungr (1987) — 3D Bishop 간편법 (3D Simplified Bishop Method)
+- **출처**:
+  - Hungr, O. (1987), "An extension of Bishop's simplified method of slope stability analysis to three dimensions", *Géotechnique*, 37(1), 113–117. DOI: `10.1680/geot.1987.37.1.113` / URL: `https://doi.org/10.1680/geot.1987.37.1.113`
+  - Hungr, O., Salgado, F. M., & Byrne, P. M. (1989), "Evaluation of a three-dimensional method of slope stability analysis", *Canadian Geotechnical Journal*, 26(4), 679–686. DOI: `10.1139/t89-079` / URL: `https://doi.org/10.1139/t89-079`
+- **핵심 가정**: 기둥 간 연직 전단력($X_x, X_y$)은 무시하고, 각 기둥의 **연직 방향($z$) 힘의 평형**으로부터 바닥 수직력 $N_{i,j}$를 구한 뒤 회전축에 대한 **전체 모멘트 평형**을 만족시킴.
+- **핵심 수식**:
+  $$F_{\text{3D,Bishop}} = \frac{\displaystyle\sum_{i,j} R_{i,j} \frac{c'_{i,j} A_{i,j} + \left( W_{i,j} - u_{i,j} A_{i,j} \cos\psi_{i,j} \right) \tan\phi'_{i,j}}{m_{\alpha (i,j)}}}{\displaystyle\sum_{i,j} W_{i,j} \left( R_{i,j} \sin\varepsilon_{i,j} + k_h e_{i,j} \right)}$$
+  $$\text{여기서, } \quad m_{\alpha (i,j)} = \cos\psi_{i,j} + \frac{\sin\varepsilon_{i,j} \cdot \tan\phi'_{i,j}}{F_{\text{3D,Bishop}}}$$
+  - (참고: $A_{i,j} \cos\psi_{i,j} = \Delta x \Delta y$ 이므로 분자의 $(W_{i,j} - u_{i,j} \Delta x \Delta y)\tan\phi'_{i,j}$는 2D Bishop 식과 완전히 동일한 구조를 가짐).
+- **특징**: 좌우변에 $F_{\text{3D,Bishop}}$이 포함되어 고정점 반복법(통상 5~10회 반복)으로 수렴시킴. 대칭·준대칭 3D 활동면에서 3D 엄밀해(Spencer, M-P)와 오차가 1~3% 이내이면서 수렴 실패가 거의 없어 **실무 3D LEM의 사실상 표준(Scoops3D, CLARA, Slide3, 3D Slope/W 기본 옵션)**임.
+- **우리 데이터 적합성**: **Aislo의 2D Bishop 간편법 솔버와 수식 구조($m_\alpha$ 반복 및 클리핑)가 1:1 호환되며, LiDAR DEM + 3D 노선 절·성토 계획면 통합 메시 위에서 가장 안정적으로 동작함.**
+
+#### ③ Hungr et al. (1989) — 3D Janbu 간편법 및 보정법 (3D Simplified / Corrected Janbu Method)
+- **출처**: Hungr, O., Salgado, F. M., & Byrne, P. M. (1989), *Canadian Geotechnical Journal*, 26(4), 679–686. DOI: `10.1139/t89-079`
+- **핵심 가정**: 회전 모멘트 평형 대신 **활동 방향($x$) 수평력 평형**과 **연직($z$) 힘 평형**만을 사용함. 구형·원통형이 아닌 임의 형상(하부 암반 경계면을 따르는 다면체·비원호 3D 파괴면)에 적용 가능함.
+- **핵심 수식**:
+  $$F_{\text{3D,Janbu}} = f_0 \cdot F_0, \qquad F_0 = \frac{\displaystyle\sum_{i,j} \frac{\left[ c'_{i,j} A_{i,j} + \left( W_{i,j} - u_{i,j} A_{i,j} \cos\psi_{i,j} \right) \tan\phi'_{i,j} \right] \cos\varepsilon_{i,j}}{m_{\alpha (i,j)}}}{\displaystyle\sum_{i,j} W_{i,j} \left( \tan\varepsilon_{i,j} + k_h \right)}$$
+- **우리 데이터 적합성**: **얇은 산림 표토층이 3D 기반암(리핑암·연암) 경계면을 따라 미끄러지는 비구형 복합 파괴면 검토에 적합함.**
+
+#### ④ Cheng & Yip (2007) — 비대칭 3D Bishop · Janbu · Morgenstern-Price (활동방향 $\alpha_{\text{slide}}$ 포함 엄밀해)
+- **출처**: Cheng, Y. M., & Yip, C. J. (2007), "Three-Dimensional Asymmetrical Slope Stability Analysis Extension of Bishop's, Janbu's, and Morgenstern–Price's Techniques", *Journal of Geotechnical and Geoenvironmental Engineering (ASCE)*, 133(12), 1544–1555.
+  - **DOI**: `10.1061/(ASCE)1090-0241(2007)133:12(1544)` / URL: `https://doi.org/10.1061/(ASCE)1090-0241(2007)133:12(1544)`
+- **핵심 가정**: 실제 산지 사면이나 곡선부 임도는 좌우 대칭이 아니므로, 미지의 **전체 활동 방위각(Sliding direction angle $\theta_s$)**을 추가 미지수로 두고 $x, y$ 양방향 힘 평형 및 모멘트 평형을 연립하여 풂(Rocscience `Slide3`의 핵심 엔진 원리임).
+- **특징**: 비대칭 쐐기·곡선 사면에서 정확하나, 기둥 간 전단력 함수($f(x,y)$)와 스케일링 인자($\lambda_x, \lambda_y$) 연립 비선형 방정식 수렴 실패(Non-convergence)가 빈번하고 연산량이 3D Bishop 대비 10~30배 큼.
+- **우리 데이터 적합성**: **급곡선(배향곡선) 구간의 비대칭 3D 사면 검토에 이론적 가치가 있으나, 수렴 불안정성 때문에 1차 구현보다는 3D Bishop에 방위각 회전 탐색을 결합하는 방식이 실용적임.**
+
+---
+
+### 1.3 3D 한계평형 기둥법 4종 비교 요약표
+
+| 방법 | 만족 평형 조건 | 기둥 간 힘 가정 | 적용 파괴면 형상 | 반복 수렴성 / 속도 |
+| :--- | :--- | :--- | :--- | :--- |
+| **Hovland (1977)** | 전체 회전 모멘트 | 수직·전단 모두 무시 | 구형 / 회전타원체 | 비반복(1회 연산) / **초고속** |
+| **Hungr 3D Bishop (1987)** | 연직 힘 + 전체 회전 모멘트 | 수평력 고려, 연직 전단 무시 | 구형 / 회전타원체 / 축대칭 | 5~8회 고정점 수렴 / **고속·안정** |
+| **Hungr 3D Janbu (1989)** | 연직 힘 + 수평($x$) 힘 | 수평력 고려, 연직 전단 무시 | 임의 비구형 / 기반암 경계 | 5~10회 고정점 수렴 / **고속** |
+| **Cheng & Yip 3D M-P (2007)** | 3축 힘($x,y,z$) + 3축 모멘트 | 측면 수직·전단 함수($\lambda f$) 모두 고려 | 임의 비대칭 3D 곡면 | 다변수 비선형 연립 / **느림·발산 위험** |
+
+---
+
+## 2. 3D 강도감소법(SRM, Strength Reduction Method) FEM / FDM 개요
+
+- **대표 출처**:
+  1. Griffiths, D. V., & Marquez, R. M. (2007), "Three-dimensional slope stability analysis by elasto-plastic finite elements", *Géotechnique*, 57(6), 537–546. DOI: `10.1680/geot.2007.57.6.537` / URL: `https://doi.org/10.1680/geot.2007.57.6.537`
+  2. Dawson, E. M., Roth, W. H., & Drescher, A. (1999), "Slope stability analysis by strength reduction", *Géotechnique*, 49(6), 835–840. DOI: `10.1680/geot.1999.49.6.835`
+- **핵심 수식 및 원리**:
+  - 3차원 지반 영역을 사면체(Tetrahedron) 또는 육면체(Hexahedron) 유한요소(FEM)/유한차분(FDM) 격자로 분할하고 Mohr-Coulomb 탄소성 구성방정식을 부여함.
+  - 지반의 초기 점착력 $c'$와 내부마찰각 $\phi'$를 **강도감소계수(SRF, Strength Reduction Factor)**로 나눈 감소 강도정수 $c'_f, \phi'_f$를 단계적으로 투입하여 비선형 탄소성 유한요소 해석을 반복 수행함:
+    $$c'_f = \frac{c'}{\text{SRF}}, \qquad \phi'_f = \arctan\left( \frac{\tan\phi'}{\text{SRF}} \right)$$
+  - 소성 전단변형률 집중대(Plastic shear strain localization zone)가 사면 하단에서 상단까지 관통하여 탄소성 솔버(Newton-Raphson)가 더 이상 평형 수렴에 도달하지 못하는(또는 절점 변위가 급증하는) 임계 $\text{SRF}_{\text{crit}}$ 값을 3D 안전율 $F_{\text{3D,SRM}}$으로 판정함.
+- **장단점 비교 (3D LEM 대비)**:
+  - **장점**: 파괴면 형상(구형·타원체)이나 활동 방향을 사전에 가정할 필요가 없으며, 곡선 노선·교대·옹벽·말뚝과의 3차원 응력-변위 상호작용을 자동으로 찾아냄.
+  - **단점**: 3D 볼륨 메시(Volume Mesh) 자동 생성 품질에 민감하고, 1개 구간 해석에 수 분~수십 분이 소요되어 수 $\text{km}$ 임도 노선 전체 탐색에는 부적합함.
+- **우리 데이터 적합성**: **노선 전체 탐색에는 무거우며, 3D LEM에서 위험 판정된 특정 고위험 구간(15 m 이상 대절토·교량 접속부)의 정밀 국부 검증용(Phase 2)으로만 적합함.**
+
+---
+
+## 3. USGS Scoops3D (DEM 기반 3D 활동면 탐색) 방법·입력·라이센스
+
+- **공식 출처 및 매뉴얼**:
+  - Reid, M. E., Christian, S. B., Brien, D. L., & Henderson, S. T. (2015), *Scoops3D—Software to analyze 3D slope stability throughout a digital landscape*, U.S. Geological Survey Techniques and Methods, Book 14, Chap. A1, 218 p.
+  - **DOI**: `10.3133/tm14A1` / URL: `https://pubs.usgs.gov/tm/14/a01/`
+  - **소스코드/매뉴얼 배포처**: USGS Landslide Hazards Program (`https://www.usgs.gov/software/scoops3d`)
+- **핵심 알고리즘 및 해석 방법**:
+  1. **3D 구형(Spherical) 활동면 절단**: DEM 상공의 3차원 탐색 상자(Search Box 격자 $(X_s, Y_s, Z_s)$)의 각 점을 구의 중심으로 삼고, 반지름 $R$을 $\Delta R$씩 증가시키면서 DEM 지형을 잘라내는 잠재 활동 구면(Scoop)을 수십만~수백만 개 생성함.
+  2. **부피·면적 정적 조건 필터**: 생성된 활동체의 부피 $V_{\min} \le V \le V_{\max}$ 또는 수평 면적 $A_{\min} \le A \le A_{\max}$ 범위에 들어오는 유효 구면만 선택하고, 구면 내부에 포함된 DEM 셀들을 3D 수직 기둥(Column)으로 분할함(경계부 부분 포함 셀은 Sub-cell 적분 처리).
+  3. **주활동방향($\alpha_{\text{slide}}$) 자동 산정**: 각 구면이 포함하는 지표면 DEM 격자들의 평균 경사향(Dip direction)을 계산하여 그 방향으로 회전축을 정렬한 뒤, **Hovland (1977) 3D 단순기둥법** 또는 **Hungr (1987) 3D Bishop 간편법**으로 안전율 $F_{\text{3D}}$를 산출함.
+  4. **격자별 최소 안전율 및 임계 체적 갱신**: 모든 유효 구면을 순회하며 각 DEM 셀 $(i, j)$가 경험한 가장 낮은 안전율 $F_{s,\min}(i, j)$과 그때의 **붕괴 체적 $V_{\text{crit}}(i, j)$ ($\text{m}^3$)**을 기록하여 래스터 지도로 출력함.
+- **주요 입력 파일 및 파라미터**:
+  - 지형 DEM 래스터 (ASCII Grid `.asc`), 지층 경계 DEM 래스터($k$개 층), 지층별 물성치($c', \phi', \gamma_t, \gamma_s$), 지하수 조건(간극수압비 $r_u$, 피압수두면 DEM, 또는 3D 압력수두 격자 `.thp`), 지진계수 $k_h$, 탐색 박스 범위 $(X_{\min}\sim X_{\max}, Y_{\min}\sim Y_{\max}, Z_{\min}\sim Z_{\max})$ 및 격자 간격, 최소/최대 붕괴 부피 한계($V_{\min}, V_{\max}$).
+- **라이센스 및 상용 활용 가능 여부**:
+  - **미국 연방정부 저작물(U.S. Government Work — 17 U.S.C. § 105)로서 완전한 퍼블릭 도메인(Public Domain)임.**
+  - 알고리즘 수식 이식, 소스코드(Fortran 90) 수정, 자체 Python/C++ 변환 및 **상용 소프트웨어 내장·독점 배포에 아무런 법적·라이센스 제약이 없음**(출처 표기만 권장됨).
+- **우리 데이터 적합성**: **Aislo의 LiDAR DEM에 3D 노선 계획면(절·성토면)을 합성한 수정 DEM을 입력하면 노선 전 구간의 3D 최소 안전율 지도와 붕괴 예상 토량($\text{m}^3$)을 동시에 산출할 수 있어 최적의 아키텍처 모델임.**
+
+---
+
+## 4. 3D 검증 기준 문제 (공개 벤치마크 · 해석해 · 문헌 예제 수치)
+
+자체 3D 사면안정 솔버를 구현할 때 `resources/tester/` 단위 테스트의 정답지(Ground Truth)로 사용할 수 있는 국제 공인 3D 벤치마크 3선임.
+
+### 4.1 벤치마크 1: Baligh & Azzouz (1975) / Hungr et al. (1989) 순수 점성토 구형 파괴 해석해
+- **출처**:
+  - Baligh, M. M., & Azzouz, A. S. (1975), "End effects on stability of cohesive slopes", *Journal of the Geotechnical Engineering Division (ASCE)*, 101(GT11), 1105–1117. DOI: `10.1061/AJGEB6.0000210`
+  - Hungr, O., Salgado, F. M., & Byrne, P. M. (1989), *Canadian Geotechnical Journal*, 26(4), 679–686. DOI: `10.1139/t89-079`
+  - Reid et al. (2015), *USGS Scoops3D Manual*, Verification Test 1 (pp. 73–76). DOI: `10.3133/tm14A1`
+- **입력 조건**:
+  - 비탈면 경사: $2\text{H} : 1\text{V}$ ($\beta = 26.565^\circ$), 균질 비배수 점성토 ($\phi_u = 0^\circ$, 간극수압 $u = 0$)
+  - 무차원 강도비: $\dfrac{c_u}{\gamma H} = 0.20$ (예: 사면 높이 $H = 10\text{ m}$, 단위중량 $\gamma = 20\text{ kN/m}^3$, 점착력 $c_u = 40\text{ kPa}$)
+  - 구형 파괴면: 사면 하단(Toe)을 통과하는 구형 파괴면 (구 중심이 사면 끝단 수직 상방에 위치, 반지름 $R / H = 2.0$)
+- **문헌 기준 정답 수치**:
+  - **2D 원호 해석해 (Ordinary / Bishop 동일)**: $F_{\text{2D}} = 1.329$
+  - **3D 닫힌 형태 해석해 (Closed-form Analytical Solution)**: **$F_{\text{3D,exact}} = 1.402$** (3D/2D 안전율 비 $= 1.055$)
+  - **Hungr (1989) 3D Bishop 기둥법 수치해**: **$F_{\text{3D,Bishop}} = 1.422$** (해석해 대비 $+1.4\%$ 오차)
+  - **USGS Scoops3D (세밀 격자)**: **$F_{\text{3D,Scoops}} = 1.401 \sim 1.405$**
+- **우리 데이터 적합성**: **3D 기둥 분할 면적($A_{i,j}$)과 모멘트 팔길이 적분 정확도를 검증하는 가장 기본적인 단위 테스트($\phi=0$ 검증용)로 즉시 활용 가능함.**
+
+### 4.2 벤치마크 2: Zhang (1988) / Hungr et al. (1989) $c'-\phi'$ 복합토 타원체·구형 파괴 (3가지 케이스)
+- **출처**:
+  - Zhang, X. (1988), "Three-dimensional stability analysis of concave slopes in plan view", *Journal of Geotechnical Engineering (ASCE)*, 114(6), 658–671. DOI: `10.1061/(ASCE)0733-9410(1988)114:6(658)`
+  - Chen, Z., Wang, X., Haberfield, C., Yin, J. H., & Wang, Y. (2001), "A three-dimensional slope stability analysis method using the upper bound theorem Part I: theory and methods", *International Journal of Rock Mechanics and Mining Sciences*, 38(3), 369–378. DOI: `10.1016/S1365-1609(01)00012-0`
+- **입력 조건**:
+  - 비탈면 경사: $2\text{H} : 1\text{V}$ ($\beta = 26.565^\circ$), 사면 높이 $H = 12.2\text{ m}$ ($40\text{ ft}$)
+  - 지반 물성치: 단위중량 $\gamma = 18.84\text{ kN/m}^3$ ($120\text{ pcf}$), 점착력 $c' = 28.73\text{ kPa}$ ($600\text{ psf}$), 내부마찰각 $\phi' = 20^\circ$
+  - 활동면 형상: 단면 반지름 $R_x = 24.4\text{ m}$ ($80\text{ ft}$)인 구형($R_y = 24.4\text{ m}$) 또는 회전타원체
+- **문헌 기준 정답 수치 (구형 $R_x = R_y = 24.4\text{ m}$ 기준)**:
+  - **Case 1 (균질 건조 사면, $u = 0$)**:
+    - Zhang (1988) 3D LEM: **$F_{\text{3D}} = 2.122$**
+    - Hungr et al. (1989) 3D Bishop: **$F_{\text{3D}} = 2.130$**
+    - Griffiths & Marquez (2007) 3D FEM-SSR: **$F_{\text{3D}} = 2.150$**
+    - (참고: 동일 단면 2D Bishop 안전율 $F_{\text{2D}} = 1.928$)
+  - **Case 2 (하부 박층 연약대 포함 — 연약층 $c'_w = 0\text{ kPa}, \phi'_w = 10^\circ$, 비구형 복합파괴)**:
+    - Zhang (1988): **$F_{\text{3D}} = 1.553$** / Hungr et al. (1989) 3D Bishop: **$F_{\text{3D}} = 1.620$** / 3D Janbu 보정법: **$F_{\text{3D}} = 1.540$**
+  - **Case 3 (Case 2 + 지하수위면 Piezometric surface 포함)**:
+    - Zhang (1988): **$F_{\text{3D}} = 1.258$** / Hungr et al. (1989) 3D Bishop: **$F_{\text{3D}} = 1.300$**
+- **우리 데이터 적합성**: **$c'-\phi'$ 복합 지반, 다층 경계(연약층), 지하수압($u$)이 모두 포함된 표준 검증 문제로서 3D Bishop 및 3D Janbu 솔버의 교차 검증 기준으로 적합함.**
+
+### 4.3 벤치마크 3: Leshchinsky et al. (1985, 1986) 변분법 3D 한계해석 벤치마크
+- **출처**: Leshchinsky, D., Baker, R., & Silver, M. L. (1985), "Three dimensional analysis of slope stability", *International Journal for Numerical and Analytical Methods in Geomechanics*, 9(3), 199–223. DOI: `10.1002/nag.1610090302`
+- **내용**: 변분법(Variational limit equilibrium)으로 구한 대수나선-타원형 3D 임계 파괴면의 안전율 도표($\phi' = 15^\circ \sim 40^\circ$, 경사 $30^\circ \sim 90^\circ$)를 제공하며, 3D Bishop 간편법이 변분 엄밀해와 $\pm 2\%$ 이내로 일치함을 입증한 기준 문헌임.
+- **우리 데이터 적합성**: **급경사 절토면($45^\circ \sim 60^\circ$, 1:1.0~1:0.5)에서의 3D 최소 안전율 탐색 알고리즘 정확도 검증에 적합함.**
+
+---
+
+## 5. 임도·산지 도로 비탈면 3D 해석 사례 (국내 KCI · 해외)
+
+### 5.1 국내 KCI 논문 사례
+1. **남경훈, 이영규, 이승수, 정교철, 김종태 (2025)**, "Scoops3D를 활용한 간극수압비 변화에 따른 3차원 사면 안정성 평가 (Evaluation of Three-Dimensional Slope Stability Under Varying Pore Pressure Ratios Using Scoops3D)", *지질공학 (The Journal of Engineering Geology)*, 35(4), 731–742.
+   - **DOI**: `10.9720/kseg.2025.4.731` / URL: `https://doi.org/10.9720/kseg.2025.4.731`
+   - **핵심 내용**: 국내 산지 사면을 대상으로 수치표고모형(DEM)과 탄성파 탐사로 추정한 지층 구조(붕적토·풍화토·기반암)를 구축하고, `Scoops3D`(3D Bishop 법)를 이용해 우기 시 간극수압비($r_u = 0.0 \sim 0.5$) 증가에 따른 3차원 안전율 분포와 잠재 붕괴 토량($\text{m}^3$) 변화를 정량 산출함.
+2. **김재민, 최경련, 유충식 (2017)**, "3차원 유한요소해석을 통한 비탈면의 기하학적 형상에 따른 안정성 평가", *한국지반환경공학회 논문집*, 18(10), 25–34.
+   - **핵심 내용**: 산지 도로·임도 절토 비탈면에서 평면 선형이 오목(Concave)·볼록(Convex)·직선일 때 곡률반경과 사면 높이에 따른 3D 안전율 변화를 규명함.
+3. **박혁진 외 / 국립산림과학원·한국지질자원연구원 산사태 연구군**:
+   - 산림 유역 및 임도 개설지에서 1차원 무한사면 모델(`SINMAP`, `TRIGRS`)이 지형 수렴부에서 안전율을 과소평가(과도한 위험 면적 산출)하는 한계를 극복하기 위해 3D 구형/타원체 기둥법(`Scoops3D`)을 적용해 실제 산사태 발생 위치와의 적중률(AUC)이 향상됨을 보고함.
+
+### 5.2 해외 산지 도로·산림 사면 3D 해석 사례
+1. **Reid, M. E., et al. (2000, 2015)**, "Digital elevation model resolution effects on shallow-landslide hazard assessment / Scoops3D applications", *Geology* / *USGS TM 14-A1*.
+   - **DOI**: `10.3133/tm14A1`
+   - **핵심 내용**: 고정밀 LiDAR DEM($1\text{ m} \sim 5\text{ m}$ 격자)을 사용하여 산악 지대 및 절취 도로 주변의 3D 심층·천층 붕괴 위험 구간과 붕괴 부피를 동시 산정한 표준 사례임.
+2. **Mergili, M., Marchesini, I., et al. (2014)**, "A strategy for GIS-based 3-D slope stability modelling over large areas", *Geoscientific Model Development*, 7, 2969–2982.
+   - **DOI**: `10.5194/gmd-7-2969-2014` / URL: `https://doi.org/10.5194/gmd-7-2969-2014`
+   - **핵심 내용**: 산악 지형 DEM과 공간 변화 토심 모델을 결합하여 회전타원체(Ellipsoid) 3D 활동면을 무작위 샘플링으로 탐색하고 광역 산사태 안전율 지도를 구축함.
+
+- **우리 데이터 적합성**: **국내 KCI(남경훈 외, 2025) 및 해외 사례 모두 ‘고정밀 DEM + 지층 경계면 + 간극수압비($r_u$)/지하수위 + 3D Bishop 기둥법’ 파이프라인을 사용하므로 Aislo의 LiDAR DEM 및 3D 노선 토층 모델과 100% 부합함.**
+
+---
+
+## 6. LiDAR DEM 기반 3D 임계 활동면 자동 탐색 알고리즘
+
+2D에서는 원호의 중심 $(X_c, Y_c)$와 반지름 $R$ (자유도 3개)만 찾으면 되지만, 3D에서는 구형 또는 타원체의 파라미터(자유도 4~9개) 공간에서 최소 안전율($F_{\text{3D},\min}$)을 찾아야 함. 문헌에서 검증된 3대 탐색 알고리즘은 다음과 같음.
+
+### 6.1 구형(Sphere) 4-파라미터 격자 + 부피 제약 탐색 (Scoops3D 방식)
+- **출처**: Reid et al. (2015), *USGS Techniques and Methods 14-A1*. DOI: `10.3133/tm14A1`
+- **파라미터 (자유도 4)**: 구 중심 좌표 $(X_c, Y_c, Z_c)$ + 반지름 $R$.
+- **탐색 전략**:
+  1. 거친 탐색(Coarse Grid): 노선 주변 상공 격자에서 $R$을 증가시키며 목표 부피 범위($V_{\min} \sim V_{\max}$, 예: 임도 절토 사면의 경우 $50\text{ m}^3 \sim 50,000\text{ m}^3$)를 만족하는 구면에 대해 빠른 Hovland 법 또는 거친 기둥으로 $F_s$를 계산함.
+  2. 정밀 탐색(Fine Refinement): 하위 안전율 상위 후보 중심점 주변에서 격자를 세분화하여 Hungr 3D Bishop 법으로 정밀 최소 안전율 $F_{\text{3D},\min}$을 확정함.
+- **우리 데이터 적합성**: **노선 중심선으로부터 좌우 일정 폭(예: $\pm 50\text{ m}$)으로 탐색 상자(Bounding Box)를 제한할 수 있어 임도 노선 따라가기(Corridor Search)에 가장 빠르고 결정론적(Deterministic)인 결과를 보장함.**
+
+### 6.2 일반 회전타원체(Ellipsoid) + 메타휴리스틱 최적화 탐색 (PSO / GA / Cuckoo Search)
+- **출처**:
+  1. **PSO (입자군집최적화)**: Kalatehjari, R., Rashid, A. S. A., Ali, N., & Hajihassani, M. (2014), "The Contribution of Particle Swarm Optimization to Three-Dimensional Slope Stability Analysis", *The Scientific World Journal*, 2014, 973093. DOI: `10.1155/2014/973093` / URL: `https://doi.org/10.1155/2014/973093`
+  2. **GA (개선 유전알고리즘)**: Zhou, X. P., Huang, X. C., & Zhao, X. F. (2020), "Optimization of the Critical Slip Surface of Three-Dimensional Slope by Using an Improved Genetic Algorithm", *International Journal of Geomechanics (ASCE)*, 20(8), 04020120. DOI: `10.1061/(ASCE)GM.1943-5622.0001747`
+  3. **Cuckoo Search (뻐꾸기 탐색 + 국부 변형)**: Rocscience `Slide3` Technical Documentation (Ellipsoidal Cuckoo Search + Surface Altering Optimization).
+- **파라미터 (자유도 6~7)**: 타원체 중심 $(X_c, Y_c, Z_c)$, 3축 반경 $(a, b, c)$, 수평 회전 방위각 $\theta_{\text{aspect}}$.
+  $$\frac{(x' - X_c)^2}{a^2} + \frac{(y' - Y_c)^2}{b^2} + \frac{(z' - Z_c)^2}{c^2} = 1$$
+- **탐색 전략**: 전체 격자를 전수 조사하는 대신 30~50개의 입자(Particle/Nest)를 생성해 50~100세대 진화시켜 약 1,500~3,000회 목적함수($F_{\text{3D,Bishop}}$) 평가만으로 길쭉하거나 납작한 타원체 임계 파괴면을 찾아냄. 하부 암반층과 만나면 암반 경계면을 따라 타원체 바닥을 절단(Truncated ellipsoid)함.
+- **우리 데이터 적합성**: **임도 절토부처럼 노선 방향($y$)으로 길고 종단 방향($x$)으로 짧은 타원형 파괴나 암반 경계 절단 파괴를 정밀 탐색할 때 구형 탐색보다 현실적인 활동면을 도출함.**
+
+### 6.3 확률적 무작위 타원체 샘플링 (Random Ellipsoid Sampling — `r.slope.stability` 방식)
+- **출처**: Mergili et al. (2014), *Geoscientific Model Development*, 7, 2969–2982. DOI: `10.5194/gmd-7-2969-2014`
+- **탐색 전략**: DEM 픽셀당 지정된 밀도($d$)로 타원체 중심, 장축 길이($L$), 폭($W$), 깊이($D$)를 난수 생성하여 각 픽셀이 속한 모든 타원체 중 최소 $F_{\text{3D}}$를 중첩 기록함.
+- **우리 데이터 적합성**: **유역 전체 광역 스크리닝에는 유용하나, 난수 시드에 따라 결과가 미세하게 달라지므로 설계 납품 보고서용보다는 광역 히트맵용으로 적합함.**
+
+---
+
+## 7. 오픈소스 구현과 라이센스 검토 (상용 배포 가능 여부)
+
+별도 유료 라이센스 모듈로 상용 배포하기 위해 반드시 확인해야 하는 기존 3D 사면안정 오픈소스 소프트웨어들의 라이센스 전수 조사표임.
+
+| 소프트웨어 / 코드명 | 개발 주체 | 언어 | 라이센스 종류 | **상용 독점 배포 가능 여부** | 비고 및 활용 방식 |
+| :--- | :--- | :--- | :--- | :---: | :--- |
+| **`Scoops3D`** | 미국 지질조사국 (USGS) | Fortran 90 | **Public Domain** (미국 정부 저작물) | **⭕ 가능 (제한 없음)** | 수식·알고리즘·소스코드를 Python/NumPy 또는 Cython/Numba로 변환하여 상용 모듈에 직접 내장 가능함 |
+| **`r.slope.stability`** | Mergili et al. (GRASS GIS 모듈) | C / Python | **GNU GPL v2+** | **❌ 불가 (전염성)** | 코드 직접 링크·번들링 시 상용 모듈 소스 공개 의무 발생 — **논문(CC-BY 3.0)의 수식만 참조하고 코드는 차용 금지** |
+| **`OpenSees` / `OpenSeesPy`** | UC Berkeley (PEER) | C++ / Python | 자체 라이센스 (**비상업 무료 / 상업용 별도 계약**) | **❌ 불가 (상용 제한)** | 상업적 재배포 시 별도 라이센스 비용 발생 |
+| **`FEniCSx` (`dolfinx`)** | FEniCS Project | C++ / Python | **LGPL v3** | **△ 조건부 (동적 링크)** | Phase 2에서 3D FEM-SSR 구현 시 동적 라이브러리 호출로 활용 가능 |
+| **`scikit-fem`** | Gustafsson et al. | 순수 Python (NumPy/SciPy) | **BSD 3-Clause** | **⭕ 가능 (제한 없음)** | 순수 Python 3D 사면체 유한요소 어셈블리 라이브러리로 상용 배포 안전함 |
+| **자체 NumPy / Numba 3D Column Solver** | Aislo 자체 구현 | Python 3.12 (`NumPy`, `SciPy`, `Trimesh`) | 자체 독점 코드 (USGS Public Domain 수식 기반) | **⭕ 가능 (최적)** | 구/타원체 방정식 $z_{\text{slip}}(x,y)$와 DEM $z_{\text{dem}}(x,y)$의 차이로 기둥 배열을 벡터화 계산하면 약 400줄 내외로 초고속 3D Bishop 솔버 자체 구축 가능함 |
+
+---
+
+## 8. 2D 결과와 3D 결과의 관계 ($F_{\text{3D}} / F_{\text{2D}}$ 안전율 비 · 끝단 효과) 연구
+
+### 8.1 핵심 이론 및 문헌 근거
+1. **Cavounidis (1987)의 정리 ($F_{\text{3D}} \ge F_{\text{2D}}$)**:
+   - **출처**: Cavounidis, S. (1987), "On the ratio of factors of safety in slope stability analyses", *Géotechnique*, 37(2), 207–210. DOI: `10.1680/geot.1987.37.2.207` / URL: `https://doi.org/10.1680/geot.1987.37.2.207`
+   - **핵심 내용**: 2D 평면변형률(Plane-strain) 해석은 사면의 폭이 무한대($L \to \infty$)이고 양 끝단의 측면 전단저항(End effects)이 $0$인 특수한 경우이므로, 동일한 지반에서 전역 최소 안전율을 비교하면 **항상 3D 최소 안전율이 2D 최소 안전율보다 크거나 같음($F_{\text{3D},\min} \ge F_{\text{2D},\min}$)**을 수학적으로 증명함.
+2. **사면 폭-높이 비($L / H$) 및 토질에 따른 $F_{\text{3D}} / F_{\text{2D}}$ 비율 (Gens et al., 1988; Stark & Eid, 1998; Griffiths & Marquez, 2007)**:
+   - **출처**:
+     - Gens, A., Hutchinson, J. N., & Cavounidis, S. (1988), "Three-dimensional analysis of slides in cohesive soils", *Géotechnique*, 38(1), 1–23. DOI: `10.1680/geot.1988.38.1.1`
+     - Stark, T. D., & Eid, H. T. (1998), "Performance of three-dimensional slope stability methods in practice", *Journal of Geotechnical and Geoenvironmental Engineering (ASCE)*, 124(11), 1049–1060. DOI: `10.1061/(ASCE)1090-0241(1998)124:11(1049)`
+   - **정량적 경향 (문헌 종합 수치)**:
+     | 활동 폭 대비 사면 높이 비 ($L / H$) | 일반적 $F_{\text{3D}} / F_{\text{2D}}$ 비율 범위 | 끝단 효과(End Effect) 영향도 |
+     | :---: | :---: | :--- |
+     | $L / H \approx 1.0 \sim 1.5$ (좁은 국부 절토·계곡부) | **$1.25 \sim 1.50$** ($+25\% \sim +50\%$) | 양 측면 전단저항이 매우 크게 작용함 |
+     | $L / H \approx 2.0 \sim 3.0$ (일반 임도 절토 사면 구간) | **$1.10 \sim 1.25$** ($+10\% \sim +25\%$) | 실무에서 가장 흔한 3D 보강 효과 구간임 |
+     | $L / H \ge 6.0$ (직선 장대 비탈면) | **$1.01 \sim 1.05$** ($+1\% \sim +5\%$) | 2D 평면변형률 해($F_{\text{2D}}$)에 수렴함 |
+   - 점착력($c'$) 비중이 큰 점성토·풍화암일수록 측면 면적($A_{\text{side}}$)에 작용하는 점착 저항이 커져 $F_{\text{3D}} / F_{\text{2D}}$ 비가 커지며, 순수 사질토($c'=0$)일수록 비율이 작아짐.
+3. **평면 곡률(오목 vs 볼록) 및 2D 대표단면 선정의 함정 (Zhang, 1988; Griffiths & Marquez, 2007)**:
+   - **오목 사면(Concave slope — 계곡부 절토)**: 양 측면이 서로 조여주는 아칭 효과(Arching effect)로 인해 $F_{\text{3D}}$가 직선 사면보다 추가로 상승함.
+   - **볼록 사면(Convex slope — 능선부/산각부 곡선 절토)**: 측면 구속이 약해져 오목 사면보다 3D 안전율이 낮아짐.
+   - **2D 단면(측점 20 m 간격)의 사각지대**: 임도는 측점(STA. 20 m 간격)에서만 2D 횡단면을 끊기 때문에, **측점과 측점 사이의 국부 급경사·계곡 집중부나 노선 방향과 사교(Oblique)하는 절리·지층 경계**에서는 20 m 간격 2D 횡단면이 실제 최소 안전율 방향을 놓칠 수 있음.
+
+- **우리 데이터 적합성**: **Aislo에서 측점 20 m 간격 2D 횡단면 LEM과 LiDAR DEM 기반 3D LEM을 함께 제공하면, ① 직선 구간에서는 $F_{\text{3D}} / F_{\text{2D}} \approx 1.05 \sim 1.25$ 교차 검증으로 신뢰도를 입증하고, ② 곡선·사교 지형에서는 20 m 횡단면 사이에서 누락되는 3차원 취약 포켓과 붕괴 토량($\text{m}^3$)을 잡아내는 상호보완 체계를 완성할 수 있음.**
+
+---
+
+## 9. 참고문헌 폴더(`참고문헌/`) 수록 오픈액세스 원문 안내
+
+라이센스(CC-BY 3.0 및 미국 연방정부 Public Domain)가 재배포·보관을 허용하는 핵심 원문 PDF를 동 폴더 하위 [`참고문헌/`](참고문헌/)에 저장함:
+1. `Mergili_2014_GMD_r_slope_stability_3D_GIS.pdf` — Mergili et al. (2014), *Geoscientific Model Development*, 7, 2969–2982 (CC-BY 3.0 오픈액세스, 3D 타원체 기둥법 수식 및 GIS 격자 알고리즘 원문).
+2. `Reid_2015_USGS_TM14A1_Scoops3D_Manual.pdf` — Reid et al. (2015), *USGS Techniques and Methods 14-A1* (Public Domain, Scoops3D 3D Hovland/Bishop 기둥법 수식 및 벤치마크 원문).