균열 매질을 위한 고속 GPR 모델링 해석적 파동 전파와 얇은 층 상호작용 및 쌍극자 산란

초록

본 논문은 균열이 포함된 균일·등방성 매질에서 지표면 침투 레이더(GPR) 신호를 해석적 식으로 빠르게 모델링하는 방법을 제시한다. 균열을 얇은 평면으로 가정하고, 각 면을 전기장에 선형 편극되는 무한소 쌍극자로 이산화한다. 전파는 Jonscher 모델을 이용한 주파수 의존 감쇠·분산을 포함한 균일 배경에서 해석식으로 전파하고, 얇은 층의 유효 반사계수를 이용해 정상방향 효과를 보정한다. 결과는 전통적인 FDTD와 비교해 정확도는 유지하면서 계산 속도와 기하학적 유연성이 크게 향상됨을 보인다.

상세 분석

이 연구는 GPR 신호가 균일·등방성 배경 매질을 통과한 뒤, 임의의 경사와 방위, 두께, 충전 물질을 가진 얇은 균열면에 도달했을 때 발생하는 전자기 산란을 정확히 기술하는 새로운 전방 모델링 프레임워크를 제시한다. 핵심 아이디어는 균열을 무한소 전기 쌍극자들의 집합으로 보는 미시적 맥락과, 매질 전체를 기술하는 거시적 맥락을 Ewald‑Oseen 소멸 정리를 통해 연결한다는 점이다. 구체적으로, 배경 매질의 복소 유전율 εₑ,b와 투자율 μ₀를 이용해 파수 k_b = ω√(εₑ,b μ₀)를 정의하고, 이 파수를 사용해 점쌍극자 방사식(near‑, intermediate‑, far‑field 포함)을 배경 매질에 맞게 변형한다. 균열은 일정 두께 d와 전기적 특성(εₑ,t, σ_t)을 갖는 얇은 층으로 모델링되며, 이때 TE·TM 모드에 대한 유효 반사계수 R_e,TE/TM을 광학에서 유도된 식(2‑4)을 그대로 적용한다.

균열 면을 “요소(element)”라 부르는 작은 사각형으로 이산화하고, 각 요소의 가로·세로 길이를 파장보다 충분히 작게 설정한다. 이렇게 하면 해당 요소 전체에 걸쳐 입사 전자기장이 거의 평면파라고 가정할 수 있다. 각 요소 내부에서는 전기장이 일정하므로 편극밀도 P = χ_e E 로 표현되며, 이 P를 표면 편극밀도 𝑃̂ 로 변환해 첫 번째 경계면(z_c)에서만 반사 효과를 적용한다. 결과적으로, 요소가 방출하는 전기 변위장 D_m 은

D_m = A_m 𝑃̂ G(r−r_c)

와 같이 요소 면적 A_m, 중심 좌표 r_c, 그리고 배경 매질의 그린함수 G를 곱한 형태가 된다. 여기서 G는 점쌍극자 방사식에 배경 파수를 대입한 형태이며, 복소 파수에 의해 감쇠와 위상 변이가 자동으로 반영된다.

전파와 산란을 모두 해석식으로 기술함으로써, 전통적인 유한차분시간영역(FDTD) 방식이 겪는 “계단화(staircasing)”와 시간 샘플링에 의한 수치 분산 문제를 회피한다. 또한, 요소 간 상호작용을 무시하고 외부 소스에 의한 입사장만을 고려함으로써 계산 복잡도를 O(N) 수준으로 낮춘다(여기서 N은 요소 수). 실험적으로는 동일한 매개변수를 가진 FDTD 시뮬레이션과 비교했을 때, 파형 형태와 진폭·위상 차이가 5 % 이내에 머물며, 계산 시간은 수십 배에서 수백 배까지 단축되는 것을 확인했다.

이 모델은 특히 얇은 균열(수 mm 수준)의 두께와 채워진 물질(수분, 염분 등)의 전도성 변화를 시간에 따라 동적으로 반영할 수 있다는 장점이 있다. 따라서 트레이서 실험이나 현장 교차보어홀 GPR 관측에서, 균열 내 유체 이동에 따른 전기적 특성 변화를 실시간으로 역추정하는 데 유용하게 활용될 수 있다.