몬테카를로 기반의 가축 질병 확산 네트워크 시뮬레이션 모델 개발 연구

Ⅰ. 서 론

전 세계 사육 가축 두수는 2015년 41억 마리에서 2050년 58억 마리로 증가할 것으로 예상된다[1]. 국내의 경우 농업에서 차지하는 축산업 생산액은 19조 8천억원 39.4%으로 지난 54년간 65~20년 연평균 12.2%가 증가하였다. 축산업의 양적 성장은 가축 사육 밀집도를 높이는 방향으로 이루어져 가축질병 발생위험도가 증가되었고 이는 막대한 경제적 피해로 이어지고 있다. 2000년 이전까지만 해도 우리나라에서 발생하지 않았던 구제역, 고병원성 조류인플루엔자 등 해외 악성 가축전염병은 최근 매년 발생하면서 축산업에 막대한 피해를 주고 있는 실정이다. 또한 우리나라에서 발생한 적이 없는 아프리카돼지열병은 2019년에 돼지 농가에서 처음으로 확인되었고, 야생 멧돼지에서 지금까지도 지속적으로 검출되고 있어 축산농가에 큰 위협이 되고 있다. 조류인플루엔자 등 인수공통전염병을 포함한 가축전염병은 축산업뿐만 아니라 인류를 계속해서 위협할 전망이다[2]. 실제로 국내 사회재해로 발생한 손실액 중 가축전염병으로 인한 피해액은 76%로 가장 높은 비율을 차지하고 있다. 이와 같은 질병은 위험도가 높은 관리대상 질병으로 전염성이 매우 강하며 치사율 또한 상당히 높으며 아프리카돼지열병의 경우 현재까지 백신이나 치료제가 없어 더욱이 치명적이다[3]. 축산 농가뿐아니라 사회경제적 손실비용을 줄이기 위해서는 가축 질병 발생 시 감염병의 경제적·환경적 피해를 최소화하기 위한 즉각적인 방역체계 구축, 위험요소 방역관리에 대한 대책은 필수적이다. 특히 전염병의 주요 확산 경로를 파악하고 이동 경로를 예측하여 질병 확산을 차단하는 기술은 반드시 필요하다. 이에 국내에서는 고병원성 가축질병 발생빈도와 높은 상관성을 보이는 요인들의 지역별 영향력 수준을 가늠하기 위한 연구, ICT 기반 가축 감염병 대응을 위한 커뮤니티 매핑 연구, 가축 관련 운송차랑 통행 데이터를 이용한 가축전염병 확산 예측 모형 개발 연구, 현장 조사를 통한 가축전염병 대응 정책개선 방안 연구, GIS를 이용한 가축전염병 확산 경로 분석 연구 등 가축 전염병 확산에 관련된 연구가 수행되었다[4][5][6][7][8]. 최근에는 AI, 네트워크 기술을 접목한 오픈소스 기반의 기술을 적용하여 가축 전염병의 여러 형태의 데이터를 다각적인 방법을 분석하는 연구, 빅데이터와 딥러닝을 통한 동물전염병 예방 연구, 축산시설간 차량이동망의 허브 시설 연구, 양돈 농장 및 도축장 간 네트워크 분석 연구 등도 진행되고 있는 추세이다[9][10][11][12]. 가축질병, 방역관리, 위험요소에 관련된 연구로는 축산차량 방문 정보를 이용하여 축산농가 간 차량 이동에 대해 접촉 네트워크 구축 및 분석 연구. 가축질병의 발생원리와 잠복감염 종오리 및 종란에 대한 방역관리, 병원체와 접촉한 오염된 인적자원의 방역관리 연구 등이 있다[13][14]. 가축질병 요인 및 확산에 관한 기존의 연구들은 전파확산을 억제하기 위한 방역활동 대안 및 우선순위에 대해 제시하는 방식을 취하고 있으나 실제로 농가 활용 가능성 및 구체적인 시스템은 미비한 실정이다. 또한 발생농가의 특성을 바탕으로 차후 바이러스 감염 가능농장의 지리적 위치와 밀집도를 예측하였으나, 연구 과정에서 역학조사 분석을 통해 도출된 요인들의 누락 가능성과 가축질병 예찰 및 방역에 투입 가능한 행정 인력 등의 한계로 어려움을 겪고 있다[15].

따라서 본 연구에서는 축종별 질병 발생 주요 요인을 수집하고 이를 바탕으로 가축 질병으로 인한 피해를 최소화할 수 있는 시스템을 제시하고자 하였다. 특히 본 연구에서는 요인들의 불확실성 및 한계를 보완하기 위해 몬테카를로 시뮬레이션을 활용하여 질병 확산 가능성을 확률적 변수로 지정하는 방안을 제시하였다. 질병 확산 경로를 파악하기 위해 다수의 시뮬레이션을 실행하여 가중치 및 설정값을 재추정한다면 확률적으로 정확도가 높은 결과예측이 가능할 것으로 판단된다. 몬테카를로 시뮬레이션은 무작위 추출된 난수를 이용하여 원하는 값의 확률적 분포를 계산하는 알고리즘으로 수집자료 및 통계자료가 많고 입력 값의 경향이 일치할수록 정밀한 시뮬레이션이 가능하다. 이를 바탕으로 가축전염병 확산 예측 모델을 개발하고 도출된 결과를 바탕으로 예측 모델을 평가하고자 한다. 먼저 가축질병 역학조사를 바탕으로 질병 확산에 영향을 끼치는 요인들을 도출한다. 수집된 자료를 바탕으로 질병 확산 위험도에 따라 각 요인별 가중치를 부여하였으며 몬테카를로 시뮬레이션의 적용성을 평가한다. 이후 질병 확산 예측 모델을 설계하고 질병 역학조사결과 및 실제사례를 적용하여 검증을 실시한다. 마지막으로 도출된 결과를 이용하여 질병 확산 양상을 예측할 수 있는 최종 시스템을 설계하고자 한다.

Ⅱ. 재료 및 방법

2-2 대상 데이터 수집 및 범위

본 연구에서는 구제역, 아프리카돼지열병, 조류인플루엔자와 같은 1종가축전염병을 대상으로 농림축산검역본부로부터 수집된 역학조사 분석 데이터를 활용하여 가축 질병 예찰 및 효율적인 방역대 설정 방법을 제안하였다. 제안하는 데이터는 국내에서 발병한 질병을 대상으로 발병 원인, 농장방역사항 발생농장 및 농장소재지, 축산계열화사업자명, 가축전염병 발생일, 가축의 종류 및 사육규모, 가축전염병 예방·확산 방지를 확보하였다. 특히 사육농장내로 유입된 위험요인과 초기 전파확산 차단을 위해 축산관련차량 및 사람의 이동에 의한 병원체의 전파, 야생동물의 이동 등 주요 위험요인을 효과적으로 차단할 수 있는 프로그램을 개발하고자 하였다.

본 연구에서는 시뮬레이션 모델 설계 및 네트워크 구축을 위한 자료로 가축전염병 역학조사 분석 보고서의 주요 발병자료만을 사용하여 필요한 자료만 추출하여 새롭게 구성하였다 (그림1). 하지만 고려한 역학조사 분석보고서의 데이터는 도출된 요인들의 누락 가능성 등으로 인한 수집 자료의 한계가 존재한다. 불완전한 데이터로 인해 생길 수 있는 문제점을 최소화하기 위해 수집한 요인의 영향도 및 유효성을 판단하였다 (그림2). 또한 각 요인별로 확산 위험도에 따른 가중치를 적용시켜 예측모델을 구현하였다.

Fig. 1.

visiting history report

Fig. 2.

Major variable impact levels

2-3 확산 알고리즘 모델링 (Modeling)

가축 전염병 네트워크는 방향성 (연관성, 가중치)이 존재하며 한 절점에 많은 연결선이 연결되어 있다면 질병 확산에 얼마나 큰 영향을 끼칠 수 있는지 측정할 수 있다. 즉, 최초전파자와 특정지역에서 영향력이 높은 절점을 찾아 질병의 확산을 막는 것은 중요하다. 본 연구에서는 네트워크의 절점은 농가, 연결선은 농가간을 연결하는 연결성으로 정의하였다 (그림3). 연결선은 이동 정보 (사람, 차량), 지형 정보, 거리 정보, 산 정보, 내부 정보를 고려하였다. 예를 들어 농가 간에 산 거리 요소에 대한 연관성이 없으면 연결선을 배제하였다.

Fig. 3.

Farm layout and distance information

투입된 절점 및 연결선에는 확률 통계적 분석방법인 몬테카를로 시뮬레이션 (Monte Carlo Simulations)이 적용된다. 이는 하나의 방향성에 의거하여 목표 값에 도달하는 것이 아니라, 다양한 시나리오를 고려하여 확률적 불확실성을 감안하여 예측양상을 띄기 때문에 예기치 못한 변동에 따른 위험성을 체계적으로 관리할 수 있다.

1) 1차 가중치

1차 가중치에 고려되는 인자는 축종별 가축 전염병 역학조사 분석 보고서에 따라 농가와 농가 간 질병 주요 전파 요인으로 산정했으며 이는 사람 이동 정보, 차량 이동 정보, 농장 거리 정보, 주변 지형 정보, 산 정보, 내부 정보에 해당한다. 1차 가중치는 농가 간 실제적인 확산 요인들을 중요도에 가중치를 고려하기 때문에 기존의 전염병을 예측하는 방법들과 달리 실제 영향을 끼치는 요인 및 중요도를 반영할 수 있다는 강점을 갖는다.

Fig. 4.

Weighting factors by risk sources

1차 가중치 값은 각 요소별 중요도에 따라 0~10까지 설정할 수 있으며 영향도가 없는 요소의 가중치는 0에 해당한다. 각 요소별 최대값은 1이며, 연결선을 가지는 요소는 총 6개에 해당한다. 예를 들어 모든 요소의 가중치가 최대값으로 산정되면 적용되는 1차 가중치는 6에 달한다. 가중치에 따른 1차 가중치를 산정하는 식은 아래와 같다.

$\[ \mathrm{실}\mathrm{제}\mathrm{적}\mathrm{용}\mathrm{되}\mathrm{는}1\mathrm{차}\mathrm{가}\mathrm{중}\mathrm{치}=\frac{\mathrm{각}\mathrm{요}\mathrm{소}\mathrm{가}\mathrm{중}\mathrm{치}\mathrm{값}}{\left(10\right)} \]$

(1)

2) 2차 가중치

2차 가중치란 농장과 농장 사이 간 중요도를 의미하며, 실제 적용되는 2차 가중치는 1차 가중치를 산정 후 요소별 중요도에 따라 결정된다 (그림5).

Fig. 5.

Weighting factors for the detailed options

2차 가중치 값은 각 요소별 중요도에 따라 사용자가 0부터 10까지 설정할 수 있다. 여기서, 사용자가 특정 요소에 대해 연결도가 없다고 판단하면 2차 가중치는 고려하지 않으며 가중치 값은 0에 해당한다.

3) 2차 가중치 (농장 거리 정보)

구제역 발생 원인과 바이러스가 전파되는 과정에 대해서 다양한 추측설들 (자연발생, 공기에 의한 전파, 사료차 및 이동수단 등)이 있지만 아직까지는 명확한 원인과 근거를 찾을 수 없는 실정이다. 현재 가축질병 발생 시 방역체계를 구축하기 위해 구제역 발생농가를 중심으로 한 오염지역 반경 500m, 위험지역 반경 3Km을 각각 설정하고 오염지역 내 가축은 모든 살처분 대상이 된다. 이에 발생농가의 지형적인 특정을 고려하여 지형분류도 작성 및 지형분리 체계를 수립함으로서 가축전염병 발생 시 불명확한 경계선으로 인한 피해를 최소화하였다[16]. 본 논문에서는 농장과 농장 간 거리 정보에 따른 연결성을 판별하기 위해 농장 위치에 대한 정보를 수집하였다. 위치 정보는 위도와 경도를 기준으로 측정하였으며, 각 농장에 ID를 부여하여 분석하였다 (그림6). 특정 농장을 기준으로 반경 10km 이하에 들어가는 농장들은 그림7과 같이 서로 네트워크로 연결하였으며 농장 거리를 기준으로 가중치를 달리 하였다. 농장 거리에 대한 연결도를 산정하는 식은 아래와 같다.

$\[ \mathrm{농}\mathrm{장}\mathrm{거}\mathrm{리}\mathrm{정}\mathrm{보}\mathrm{연}\mathrm{결}\mathrm{도}=1\mathrm{차}\mathrm{가}\mathrm{중}\mathrm{치}\times \mathrm{농}\mathrm{장}\mathrm{거}\mathrm{리}\mathrm{정}\mathrm{보}\mathrm{점}\mathrm{수} \]$

(2)

Table 1.

Farm distance score by the distance between each two farms

Fig. 6.

Farm location Information

Fig. 7.

Node connection by distance among farms

4) 2차 가중치 (산 정보)

가축 질병 발생 시, 질병 특성 파악을 통해 바이러스 전파와 지형적인 요소의 상관여부를 중점적으로 분석하여 피해를 최소화하는 것은 중요하다. 사람이나 차량이 이동하는 경로를 소독하고 통제해도 전염병에 감염된 야생동물이 산을 이용해 장거리를 이동할 경우 확산 양상을 파악하기 어렵다. 이들의 행동반경은 하룻밤 사이 50km 내외를 오갈 수 있으며 축산농가에 가축들에게 전파될 가능성이 매우 높다. 구제역의 경우 발생농가와 비발생농가의 고도를 분석한 결과, 구제역 발생 농가가 비발생 농가보다 평균 24m 높게 분포되어 있었다[17]. 이에 본 연구에서는 산의 위치 정보를 고려하여 질병 예측 모델을 시뮬레이션 하였다. 산 정보는 위도, 경도, 높이를 포함한다 (그림8).

Fig. 8.

Mountain Information

Fig. 9.

Node connection by mountain information

산 정보는 특정 산 기준 반경 10km 이하에 들어가는 농장들을 서로 네트워크 연결하였다. 농장 간 네트워크 점수는 산으로부터 두 농장을 거리를 기준으로 계산한다. 단, 여러 네트워크 연결선 형성 시 산으로부터 해당 농장 거리에 대해 합하여 적용한다. 산으로부터 두 농장 거리의 합 점수는 거리 기준으로 하며 표2와 같다.

$\[ \mathrm{농}\mathrm{장}\mathrm{산}\mathrm{정}\mathrm{보}\mathrm{연}\mathrm{결}\mathrm{도}=1\mathrm{차}\mathrm{가}\mathrm{중}\mathrm{치}\times 2\mathrm{차}\mathrm{가}\mathrm{중}\mathrm{치}\times \mathrm{산}\mathrm{정}\mathrm{보}\mathrm{점}\mathrm{수} \]$

(3)

Table 2.

Scores by mountain range

5) 2차 가중치 (사람 및 차량 이동 정보)

가축관련 운송차량은 가축전염병 확산의 주요 원인으로 제기되고 있으나 국내에서는 가축관련운송차량의 이동성과 관련한 실증적인 방역절차와 기준이 마련되지 않고 있는 실정이다. 또한 가축관련 운송차량은 가축전염병 확산의 주요 원인으로 제기되고 있으나 국내에서는 가축관련 운송차량의 이동성과 관련한 실증적인 방역절차와 기준이 마련되지 않고 있는 실정이다[18]. 이에 본 연구에서는 Type (차량, 사람, ID) 기준으로 방문 농장 (이동) 농장의 네트워크를 연결하였다 (그림10, 11).

Fig. 10.

People and vehicle movement information

Fig. 11.

Node connection by people and vehicle information

농장 간 네트워크 연결은 시간 내림차순에 따라 ID별로 순서대로 연결하게 된다. (차량, 사람) 정보 가중치의 최대값은 각 1이며, 연결도가 겹치는 경우 각 가중치를 합하여 산정된다. 이를 식으로 나타내면 다음과 같다.

$\[ \mathrm{이}\mathrm{동}\mathrm{정}\mathrm{보}\mathrm{연}\mathrm{결}\mathrm{도}=1\mathrm{차}\mathrm{가}\mathrm{중}\mathrm{치}\times 2\mathrm{차}\mathrm{가}\mathrm{중}\mathrm{치} \times \mathrm{이}\mathrm{동}\mathrm{정}\mathrm{보}\mathrm{정}\mathrm{수}\mathrm{의}\mathrm{합} \]$

(4)

6) 2차 가중치 (내부 농장 관리 정보)

농림축산식품부에 의하면 방역시설 미흡, 과밀사육 등 사육환경이 방역 취약요인으로 작용하며 소독자원, 방역관련 업체 관리 등 방역지원체계 강화시키는 것은 필수적이다. 이에 본 논문에서는 농장 내부 관리 정보에 대해 가중치를 적용시켜 2차 가중치를 산정하였다. 가중에 따른 내부정보를 산정하는 식은 아래와 같다.

$\[ \mathrm{내}\mathrm{부}\mathrm{정}\mathrm{보}=\sum _{\mathit{x}=0}^{\mathit{n}}\frac{\left\{1\mathrm{차}\mathrm{ }\mathrm{가}\mathrm{중}\mathrm{치}\times 2\mathrm{차}\mathrm{ }\mathrm{가}\mathrm{중}\mathrm{치}\times {\mathit{N}\mathit{O}}_{\mathit{i}}\right\}}{10\times \mathrm{ }\mathrm{사}\mathrm{용}\mathrm{자}\mathrm{가}\mathrm{ }\mathrm{선}\mathrm{택}\mathrm{한}\mathrm{ }\mathrm{요}\mathrm{인}\mathrm{ }\mathrm{수}} \]$

(5)

7) 2차 가중치 (주변 지형 정보)

가축 전염병의 확산은 발병 농장 주변의 지형 및 형태에 따라 확산 양상 및 규모가 결정된다. 주변지형 정보를 고려한 2차 가중치는 농장 사이 네트워크 연결이 있는 경우 주변 지형 점수를 0부터 최대 10까지 적용하였다. 반면 농장 사이 네트워크 연결이 없는 경우 주변 지형 정보 점수를 적용하지 않고 연결선 또한 생성하지 않는다. 가중치는 특정 농가를 기준으로 주변을 8방으로 나누었다 (그림12). 주변 지형 정보 2차 가중치는 8방 셀에 포함된 지목들을 각각 적용하여 1차 가중치, 2차 가중치를 산정한 후 주변 지형 점수의 최대값을 나누어 구한다. 주변 지형 정보를 산정하는 식은 아래와 같다.

$\[ \begin{array}{c}\mathrm{주}\mathrm{변}\mathrm{지}\mathrm{형}\mathrm{정}\mathrm{보}{=}_{} \hfill \\ \sum _{x=0}^n\frac{\left\{1\mathrm{차}\mathrm{ }\mathrm{가}\mathrm{중}\mathrm{치}\times \mathrm{ }2\mathrm{차}\mathrm{ }\mathrm{가}\mathrm{중}\mathrm{치}\mathrm{ }\left(x\mathrm{의}\mathrm{ }\mathrm{정}\mathrm{수}\right)\times x\mathrm{의}\mathrm{ }\mathrm{정}\mathrm{수}\right\}}{x\times 8\times 10}\hfill \end{array} \]$

(6)

Table 3.

Internal Information Score Example

Fig. 12.

Estimation of surrounding terrain weights

8) 최종 연결선 도출

본 연구에서는 절점 (농가)간의 연결선을 통해 위험도가 높은 관계를 파악하여 가축 전염병 발생 시 확산 위험도를 분석하고 통계적 모형을 이용해 다양한 의미를 도출하고자 하였다. 네트워크 연결선을 도출하기 위해 농장과 농장 사이의 연결 여부 혹은 연결 강도를 산정하였다. 네트워크 연결선을 산정식은 식7과 같다. 농장 거리 정보, 이동 정보, 산 정보, 내부 정보에서 최종적으로 도출한 점수에 제시한 1차, 2차 가중치를 각각 곱한 뒤 더하여 최종 네트워크 연결도를 결정하였다.

$\[ \begin{array}{c}\mathrm{최}\mathrm{종}\mathrm{연}\mathrm{결}\mathrm{도}=\hfill \\ \mathrm{농}\mathrm{장}\mathrm{거}\mathrm{리}\mathrm{정}\mathrm{보}1\mathrm{차}\mathrm{가}\mathrm{중}\mathrm{치}\times \mathrm{농}\mathrm{장}\mathrm{거}\mathrm{리}\mathrm{정}\mathrm{보}\mathrm{점}\mathrm{수}\hfill \\ +\mathrm{이}\mathrm{동}\mathrm{정}\mathrm{보}1\mathrm{차}\mathrm{가}\mathrm{중}\mathrm{치}\times \mathrm{이}\mathrm{동}\mathrm{정}\mathrm{보}2\mathrm{차}\mathrm{가}\mathrm{중}\mathrm{치}\times \mathrm{이}\mathrm{동}\mathrm{정}\mathrm{보}\mathrm{점}\mathrm{수}\hfill \\ +\sum \mathrm{지}\mathrm{형}\mathrm{정}\mathrm{보}1\mathrm{차}\mathrm{가}\mathrm{중}\mathrm{치}\times \mathrm{각}\mathrm{주}\mathrm{변}\mathrm{지}\mathrm{형}\mathrm{정}\mathrm{보}2\mathrm{차}\mathrm{가}\mathrm{중}\mathrm{치}\hfill \\ \times \mathrm{각}\mathrm{주}\mathrm{변}\mathrm{지}\mathrm{형}\mathrm{정}\mathrm{보}\mathrm{점}\mathrm{수}\hfill \\ +\mathrm{산}\mathrm{정}\mathrm{보}1\mathrm{차}\mathrm{가}\mathrm{중}\mathrm{치}\times \mathrm{산}\mathrm{정}\mathrm{보}2\mathrm{차}\mathrm{가}\mathrm{중}\mathrm{치}\times \mathrm{산}\mathrm{거}\mathrm{리}\mathrm{합}\mathrm{정}\mathrm{보}\mathrm{점}\mathrm{수}\hfill \\ +\sum \mathrm{내}\mathrm{부}\mathrm{정}\mathrm{보}1\mathrm{차}\mathrm{가}\mathrm{중}\mathrm{치}\times \mathrm{각}\mathrm{내}\mathrm{부}\mathrm{정}\mathrm{보}2\mathrm{차}\mathrm{가}\mathrm{중}\mathrm{치}\times \mathrm{각}\mathrm{내}\mathrm{부}\mathrm{정}\mathrm{보}\mathrm{점}\mathrm{수}\hfill \end{array} \]$

(7)

Ⅲ. 결과 및 고찰

3-1 전처리 데이터 결과

본 연구에서 데이터를 활용하기 위하여 전처리 과정을 통하여 데이터들을 분류하였다. 앞서 말했듯이 이동 정보, 거리 정보, 지형 정보, 내부 정보로 나누어 전처리를 진행하였다. 전처리 된 데이터의 예시로 지형 정보와 같은 경우에 그림13은 각 농가별 점수를 나타내는 전처리 과정 결과이다. 그림14는 농가간의 평균값을 이용하여 연결가중치를 구하는 행렬형태의 데이터 처리 방식이며 그림15는 일반적으로 행렬의 값이 대부분 0의 값을 갖고 일부만 연결이 되기 때문에 0을 제외하고 값이 있는 정보를 도출해서 정리한 데이터 결과이다. 이처럼 전처리 된 데이터들의 가중치를 결정하여 시뮬레이션 시나리오를 결정한다면 사용자의 의해 원하는 시나리오 결과를 도출해낼수 있다.

Fig. 13.

Weight data for each farm

Fig. 14.

Linked weight calculation matrix data processing

Fig. 15.

Data preprocessing results

3-3 요인별 시뮬레이션 분석 결과

본 연구에서는 가축 전염병 발생 원인과 바이러스가 전파되는 명확한 원인와 근거를 찾기 위해 몬테카를로 시뮬레이션을 이용하여 분석하였다. 분석 요인에는 거리 정보, 내부 관리 정보, 지형 정보 총 3가지 요인으로 분석을 진행하였다. 가축 전염병 발생 기준 농가는 한 곳으로 고정 설정하여 분석 요인별 몬테카를로 시뮬레이션을 진행하였으며 기준 농가 위치는 다음 그림16과 같다.

Fig. 16.

Infectious Disease Outbreak Criteria Farmhouses

1) 거리 정보에 따른 분석 결과

그림17은 거리 정보에 따른 질병 확산률을 나타낸 그림이다. 기준 농장을 기점으로 반경 10km 이하에 들어가는 농장들은 서로 네트워크로 연결하였으며, 농장 거리를 기준으로 가중치를 달리 적용하였기 때문에 기준 농장과 거리가 가까울수록 전염병 확산률은 높았으며, 거리가 멀수록 전염병 확산률은 낮게 나타났다.

Fig. 17.

Epidemic spread rate according to distance information

가축 전염병 발생 기준 농가를 기점으로 전국 농가 데이터 및 위치 정보를 이용하여 거리 정보 요인의 시간에 따른 질병 확산 양상을 분석하였으며 그 결과는 그림18과 같다. 전염병 발생 1일차에는 기준 농가를 기점으로 거리 정보상 가중치가 높은 주변 농가 3곳으로 전염병이 확산되었으며, 2일차부터 4일차까지 전염병 확산 속도가 증가하였는데 이는 추가 감염된 농장을 기점으로 거리 요인에 따른 가중치가 다시 설정되었기 때문이라고 판단된다. 전염병 발생 4일차 이후부터는 전염병 확산 속도가 점점 감소하기 시작하여 7일차때는 추가적인 전염병 확산은 나타나지 않았다. 이는 추가 감염된 농장으로부터의 주변 농장들의 거리 정보 가중치가 낮아 추가적인 감염병 전파는 없다고 판단된다.

Fig. 18.

Epidemic spread over time of distance information factor (From the top left, day 1, 2, 4, 8)

2) 내부 관리 정보에 따른 분석 결과

그림19는 내부 관리 정보에 따른 질병 확산률을 나타낸 그림이다. 분석 결과 전염병 발생 기준 농가 포함하여 모든 농가의 전염병 확산률은 0%로 나타났다. 이는 내부 관리 정보 요인에는 전염병의 주 확산 원인인 사람, 차량, 야생동물과 같은 매개체의 이동 등이 포함되지 않기 때문에 내부 관리 정보 요인의 가중치는 다른 요인들의 가중치에 비해 작게 산정되어 영향이 미비하다고 판단된다.

Fig. 19.

Epidemic spread rate according to interna management information

따라서 내부 관리 정보 요인 하나만을 가중치로 선정하지 않고 거리 정보 요인, 지형 정보 요인, 사람과 차량의 이동 정보 요인 등 다양한 전파 요인을 고려하여 시뮬레이션을 진행하면 내부 관리 정보 요인 하나만을 가중치로 선정하였을때보다 가축 전염병 확산에 영향을 미칠 것이라고 판단된다.

농장 내부 관리 정보에 대해 가중치를 적용시켜 가축 전염병 발생 기준 농가를 기점으로 2차 가중치를 산정하여 시간에 따른 내부 관리 정보 요인의 시간에 따른 질병 확산 양상을 분석하였으며 그 결과는 그림20과 같다. 전염병 발생 1일차부터 9일차까지 추가적인 전염병 확산은 없었으며 이는 위의 내용과 같이 내부 관리 정보 요인에는 전염병의 주 확산 원인인 매개체의 이동 등이 포함되지 않기 때문에 전염병 확산에 미치는 영향이 적어 추가적인 전염병 확산이 없었다고 판단된다. 또한 내부 관리 정보 요인만으로는 전염병 확산도를 파악하기에는 문제점이 있어 다른 전파 요인을 고려하여 분석해야 한다고 사료된다.

Fig. 20.

Epidemic spread over time of internal management information factors (From the top left, day 1, 2, 4, 8)

3) 지형 정보에 따른 분석 결과

그림21은 지형 정보에 따른 질병 확산률을 나타낸 그림이다. 분석 결과 전염병 발생 기준 농가를 포함하여 모든 농가의 전염병 확산률은 100%로 나타났다. 지형 정보 요인은 주변 지형 정보를 고려하여 전염병 발생 기준 농가를 기준으로 주변을 8방으로 나눈후 8방에 포함된 지목들을 각각 적용하여 가중치를 산정하는 방식을 이용하였다. 이에 따라 전염병 발생 기준 농가로부터의 거리와 상관없이 감염 농가의 지형 정보 요인의 가중치가 높게 산정되어 지형 정보에 따른 전염병 확산률이 100%로 나타났다고 판단된다. 따라서 지형 정보 요인에 따른 전염병 확산률을 판단하기위해서는 전염병 발생 기준 농가로부터 특정범위를 한정하여 분석하여야 한다고 판단된다.

Fig. 21.

Epidemic spread rate according to geographic information

지형 정보에 따른 가중치를 적용시켜 가축 전염병 발생 기준 농가를 기점으로 2차 가중치를 산정하여 시간에 따른 지형 정보 요인의 시간에 따른 질병 확산 양상을 분석하였으며 그 결과는 그림22과 같다. 분석 결과, 전염병 발생 1일차에는 기준 농가로부터 거리가 먼 농가로 전염병이 확산되었으며, 2일차부터 3일차까지 전염병 확산 속도가 증가하여 4일차에는 모든 농가가 감염되었다. 지형 정보 요인에 따른 질병 확산 시간은 은 3가지 분석 요인들중 가장 짧았으며 이는 앞선 내용과 같이 지형 정보 요인이 전염병 발생 기준 농가로부터의 거리와 상관없이 가중치가 높게 산정되어 확산 속도가 빠르다고 판단된다. 따라서 전국을 대상으로 지형 정보 요인에 따른 질병 확산 양상을 분석하기 위해서는 전염병 발생 기준 농가로부터 특정범위를 한정하여 분석하거나, 전염병 발생 기준 농가로부터의 거리 정보를 명확히 하여 추가적인 분석이 필요하다고 판단된다.

Fig. 22.

Infectious disease spread over time of geospatial information factors (From the top left, day 1, 2, 4, 8)

Ⅳ. 결 론

본 연구에서는 몬테카를로 시뮬레이션을 이용하여 축종별 질병 발생 주요 요인을 수집하고 이를 바탕으로 가축 질병으로 인한 피해를 최소화할 수 있는 시스템을 개발하여 전염병 확산 요인별 전염병 확산률과 시간에 따른 전염병 확산 양상을 분석하였다. 이동정보, 거리정보, 지형정보, 내부정보로 나누어 전처리를 진행하였으며 전처리 된 데이터들의 가중치를 산정한 후 몬테카를로 기반의 네트워크 시뮬레이션 구현을 위해 pyCharm IDE를 이용하여 python code를 작성하였다. 취합한 데이터에 각각 1차, 2차 가중치를 부여하였고 streamlit 라이브러리를 이용하여 결과를 분석을 진행하였다.

이 연구를 통해 전염병 발생 주요 요인별 전염병 확산률과 시간에 따른 전염병 확산 양상을 분석할 수 있었다. 향후, 지형 정보 요인 분석시 전염병 발생 기준 농가로부터의 거리 정보를 명확히 하여 분석이 필요하다고 판단되며 추가적으로 기존 가축 전염병 확산 네트워크 모델을 보완하여 가축 전염병에 영향을 미치는 다양한 주요 요인들이 포함된 개선된 모델을 개발하고자 하였다.

Acknowledgments

본 연구는 농림축산식품부 스마트팜다부처패키지혁신개발사업(421042042SB010)의 지원에 의해 이루어진 것임.

References

• Food and Agriculture Organization of the United Nations (FAO) [Internet]. Available: https://www.fao.org/home/en
• Livestock Health Control Association (LHCA) [Internet]. Available:https://www.lhca.or.kr/business/index.do
• World Organization for Animal Health (OIE) [Internet]. Available: http://www.oie.int
• S. J. Hong, “Spatial Variations of Risk Factors Associated with the Diffusion of Highly Infectious Animal Diseases”, Journal of the Korean cartographic association, Vol. 18, No. 1, pp. 81-91, January 2018. [https://doi.org/10.16879/jkca.2018.18.1.081]
• J. H. Koo, “A Study on Community Mapping for ICT-Based Livestock Infectious Disease Response”, Journal of the Korean Surveying Society, Vol. 38, No. 3, pp.247-257, June, 2020. [https://doi.org/10.7848/ksgpc.2020.38.3.247]
• W. Kang, “Development of Predicting Model for Livestock Infectious Disease Spread Using Movement Data of Livestock Transport Vehicle”, The Journal of The Korea Institute of Intelligent Transportation Systems, Vol. 21, No. 4, pp. 78-95, August, 2022. [https://doi.org/10.12815/kits]
• H. S. Park, “Study on Policy Improvement Measures to Respond to Infectious Diseases of Livestock through Field Investigation”, Journal of The Korean Society of Disaster Information, Vol. 17, No.2, pp. 275-288, June, 2021. [https://doi.org/10.15683/kosdi.2021.6.30.275]
• S. Y. Lee, “Analysis of Animal Infected Diseases Dispersion Route by GIS”, Korean Society of Agricultural Machinery and Machinery, Vol.17, No.1, pp.314-317, February, 2017.
• M. J. Kyung, “Implementation of Open Source SOLAP Decision-Making System for Livestock Epidemic Surveillance and Prevention”, Korean Society Of Subveying Geodecy, Photogrammetry And Cartography, Vol. 30, No. 3, pp. 287-294, March 2012. [https://doi.org/10.7848/ksgpc.2012.30.3.287]
• S. H. Kim, “Animal Infectious Diseases Prevention through Big Data and Deep Learning”, Korea Institute of Intelligent Information Systems, Vol. 24, No. 4, pp. 137-154, December, 2018. [https://doi.org/10.13088/jiis.2018.24.4.137]
• G. J. Lee, “Hub Facilities in Vehicle Movement Network between Livestock Facilities ”, Journal of the Korea Academia-Industrial cooperation Society, Vol. 19, No. 6, pp. 137-146, May, 201. [https://doi.org/10.5762/KAIS.2018.19.6.137]
• M. J. Jeong, “Network Analysis of Swine Farms and Slaughters : Based on Automobile GPS Data”, Journal of the Korean Association of Cartography, Vol. 19, No. 1, pp. 69-79, April, 2019. [https://doi.org/10.16879/jkca.2019.19.1.069]
• G. J. Lee, “Hub Facilities in Vehicle Movement Network between Livestock Facilities”, Journal of Korea Academia-Industrial cooperation Society, Vol. 19, No. 6, pp. 137-146, June 2018. [https://doi.org/10.5762/KAIS.2018.19.6.137]
• U. K. Mun, “Column - Principles of the occurrence of livestock diseases and the importance of prevention and prevention”, Korea Duck Association, Vol. 102, No. 102, pp. 34-49, December 2011.
• C. H. Eom, “Analysis of Potential Infection Site by Highly PathogenicAvian Influenza Using Model Patterns of Avian Influenza Outbreak Area in Republic of Korea”, Journal of the Korean Association of Geographic Information Studies, Vol. 20, No. 2, pp. 60-74, February 2017. [https://doi.org/10.11108/kagis.2017.20.2.060]
• G.Y. Baek, “Analysis of foot and mouth disease area according to terrain classification”, Proceedings of the Korean Society for Agricultural Machinery Conference, 2012, 17.2: 362-366.
• Hyeon-Tae Kim, “Type Analysis of Gyeongbuk Region Occurred Foot and Mouth Disease According to Terrain Classification”, Journal of Agriculture & Life science, NOV. 28. 2013
• Jeong HeeHyeon, “ Assessment of Livestock Infectious Diseases Exposure by Analyzing the Livestock Transport Vehicle's Trajectory Using Big Data”, J. Korea Inst. Intell. Transp. Syst, 2020, vol.19, no.6, 92 pp. 134-143 (10 pages) [https://doi.org/10.12815/kits.2020.19.6.134]