충격 방출 최적화를 통한 종 경쟁 제어 모델

초록

본 논문은 두 종이 자원 경쟁을 하는 상황에서, 경쟁 종 S₂를 충격적으로 방출하여 목표 종 S₁을 제거하고자 하는 최적 제어 전략을 제시한다. 모델의 존재·유일성·양성·유계성을 증명하고, S₁‑무존재 해의 전역 안정성을 위한 충분조건을 도출하였다. 또한 방출 횟수와 규모를 최소화하면서 S₁을 알리 효과 이하로 억제하는 최적 제어 문제를 수학적으로 정립하고, 해의 존재성을 보였다. 마지막으로 Aedes aegypti 모기의 야생형과 Wolbachia 감염형을 사례로 삼아 수치 실험을 수행, 이론 결과와 최적 방출 정책의 실효성을 확인하였다.

상세 분석

본 연구는 기존 연속형 경쟁 모델을 충격 방출 형태로 변형함으로써 실제 생물학적 방제 현장의 비연속적 특성을 반영하였다. 모델식(1a‑1b)은 두 종 S₁, S₂의 성장·사망·경쟁 항을 포함하고, S₁에만 알리 효과를 도입해 최소 생존 인구(K₀) 이하에서는 성장률이 음수가 되도록 설계하였다. 파라미터 가정 ψ₁>δ₁, ψ₂>δ₂와 ψ₂<ψ₁, δ₂>δ₁, r₂<r₁은 S₁이 내재적으로 더 강인함을 의미한다. 이러한 비대칭성은 실제 야생 모기가 감염 모기보다 번식력이 높다는 생태학적 사실과 일치한다.

수학적 분석에서는 먼저 V₀ 클래스와 상한 미분 연산자 D⁺를 도입해 충격 시스템의 연속 구간과 충격 순간 모두에서 해의 존재·유일성을 확보하였다. 정리 2(비교정리)를 이용해 해가 비음수가 되도록 보장하고, 정리 1(정리 1의 적응)에서 네 개의 평형점(노드·반발점·안정점)을 도출하였다. 특히 K_b와 K_* 사이의 초기값에 따라 S₁이 지속되거나 S₂가 지배하는 두 가지 수렴 시나리오가 존재함을 확인했다.

핵심 기여는 S₁‑무존재(0, S₂) 해의 전역 안정성을 보장하는 충분조건을 제시한 점이다. 이 조건은 충격 주기 τ와 방출량 u_k가 K₀, K₁, K₂, r_i 등 모델 파라미터와 어떻게 상호작용해야 하는지를 명시한다. 구체적으로, 매 주기마다 최소 u_min ≥ f(τ, 파라미터) 만큼의 S₂를 투입하면 S₁의 인구가 알리 임계치 이하로 떨어져 자연 소멸하게 된다. 이때 Lyapunov 함수와 비교정리를 결합해 전역 수렴을 증명하였다.

다음으로 최적 제어 문제를 설정하였다. 목적함수는 ∑_{k=0}^{N-1} u_k (방출 총량)와 최종 시점 T에서 S₁이 K₀ 이하가 되도록 하는 제약조건을 포함한다. 제어 변수는 각 충격 시점의 방출량 u_k이며, 제약집합 U =