고대 제국 역학 시뮬레이션 평가와 개선 방향

초록

**
본 연구는 Turchin et al. (2013) 모델을 재구현하고, 제국 밀도, 인구 분포, 전쟁 기록 등 세 가지 역사 데이터와 비교하여 검증하였다. 전반적으로 인구와의 상관관계는 높은 편(R² ≈ 0.75)이나 개별 제국의 형태와 전쟁 발생 패턴은 충분히 재현하지 못한다. 약탈적 공격 전략을 도입한 변형 모델은 성능이 크게 저하돼, 원 모델의 무작위 공격 메커니즘이 핵심임을 확인하였다. 향후 실제 정착지 네트워크와 연계한 확률적 제국 강도 매핑이 필요하다.

**

상세 분석

**
이 논문은 클리오다이내믹스 분야에서 가장 널리 인용되는 Turchin et al. (2013) 모델을 독립적으로 구현하고, 재현 가능성을 검증한 뒤, 모델의 한계와 개선 가능성을 체계적으로 탐색한다. 구현 단계에서는 100 km × 100 km 격자를 사용해 지형(사막·초원·농경)과 고도 정보를 초기화하고, 군사 기술을 스텝 지역에서 농경지로 확산시키는 과정을 그대로 재현하였다. 각 셀은 초사회적 특성을 부여받으며, 전쟁은 인접한 농경 셀 사이에서 무작위로 발생하고, 제국 규모와 초사회적 특성에 따라 승패 확률이 조정된다.

재현 결과는 원 논문이 제시한 제국 밀도 지도와 시각적으로 거의 일치했으며, 3개 시대(1500 BC‑500 BC, 500 BC‑500 AD, 500 AD‑1500 AD)에 대해 셀별 회귀 분석을 수행한 R² 값도 원 논문과 비슷하거나 더 높았다(예: 1500 BC‑500 BC 구간 R² = 0.66 vs 0.56). 이는 모델 구현이 정확함을 의미한다.

하지만 모델이 실제 제국의 형태와 규모를 재현하는 데는 한계가 있다. 시뮬레이션에서는 소규모 혹은 단일 셀 제국이 다수 존재하고, 대규모 제국은 매우 드물다. 이는 제국 붕괴 확률이 매 단계 5 % 이상으로 설정된 결과이며, 실제 역사적 제국(예: 로마 제국)의 연속성과 영역 확장을 충분히 설명하지 못한다.

인구와의 상관관계 분석에서는 누적 제국 밀도와 인구 데이터(Reba et al., 2016) 사이에 R² ≈ 0.5‑0.75 수준의 강한 양의 상관이 관찰되었다. 특히 1400 AD‑1500 AD 구간에서 0.75에 근접했으며, 이는 인구 데이터의 품질 향상과 모델이 인구 집중 지역을 잘 포착한다는 점을 시사한다. 그러나 인구가 높은 사막 셀 등 모델에서 제국에 포함될 수 없는 지역은 실제와 차이를 만든다.

전쟁 기록과의 비교에서는 공격 횟수와 역사적 전투 데이터 간 상관이 낮았다. 모델은 제국 내부에서의 공격을 금지하고, 제국이 없는 지역에서 공격이 집중되기 때문에, 인구가 많고 제국이 발달한 유럽 지역에서는 실제 전투보다 적은 공격이 발생한다. 반면 사하라 이남 아프리카에서는 모델이 과도하게 많은 공격을 기록했으며, 이는 데이터 편향과 모델의 전쟁 메커니즘이 실제 전쟁 동기를 충분히 반영하지 못함을 보여준다.

또한, 저자는 공격 전략을 ‘약탈적(greedy)’으로 바꾸어 약한 이웃을 우선 공격하도록 수정한 변형 모델을 실험했다. 이 경우 R² 값이 전반적으로 감소하고, 제국 밀도와 인구·전쟁과의 상관이 크게 악화되었다. 이는 원 모델의 무작위 공격 가정이 제국 확산과 초사회적 특성 전파에 중요한 역할을 함을 뒷받침한다.

마지막으로 논문은 현재 모델의 구조적 한계를 짚으며, 정착지 네트워크 기반의 ‘확률적 제국 강도’를 추정하고, 인구 압력·내부 붕괴·문화적 상호작용을 명시적으로 포함하는 다중 레이어 에이전트 기반 모델로 확장할 것을 제안한다. 이러한 방향은 제국의 형성·소멸 메커니즘을 보다 정밀하게 재현하고, 역사적 사건과의 정량적 매칭을 가능하게 할 것으로 기대된다.

**