Numerical Model of a Multiple-Input-Multiple-Output Distributed Acoustic Sensing System with Joint Phase and Birefringence Estimation

기술 요약: 공동 위상 및 복굴절 추정을 포함한 MIMO-DAS의 수치 모델

문제 정의
분산형 광섬유 센싱(DAS)은 교통 감시, 구조물 건전성 모니터링, 지진 탐지와 같은 모니터링 응용 분야를 위해 기존의 통신 인프라를 활용합니다. 전통적인 접근 방식은 일반적으로 다음 두 가지 지표 중 하나에 집중합니다:

1. 차동 위상 OTDR ($\Delta\phi$-OTDR): 후방 산란광의 공통 광학 위상을 추정하며, 이는 종방향 변형(longitudinal strain)에 매우 민감하지만 편광 측면을 간과하는 경향이 있어, 순수 횡방향 또는 비등방성 섭동(perturbation)을 감지할 수 없습니다.
2. 편광-OTDR (P-OTDR): 편광 특성을 활용하여 횡방향 변형(굽힘, 뒤틀림, 횡방향 응력)을 감지합니다. 그러나 전통적인 P-OTDR은 종종 여러 개의 입력 상태 편광(SOP)을 필요로 하므로, 동적 모니터링을 제한하는 긴 측정 시간을 초래합니다. 또한, 이는 일반적으로 공통 위상을 활용하지 않으므로 축대칭 섭동에는 무감각합니다.

본 논문은 서로 다른 유형의 변형 이벤트(종방향 vs 횡방향/비등방성)를 높은 감도와 공간 해상도로 식별하기 위해, 공통 광학 위상과 편광 관련 양(특히 복굴절)을 동시에 추정하는 다중 파라미터 동적 센싱에 대한 필요성을 다룹니다.

방법론
저자들은 MIMO-DAS(Multiple-Input-Multiple-Output Distributed Acoustic Sensing) 시스템을 위한 수치 모델을 개발하고 실험적으로 검증하였습니다. 방법론은 세 가지 주요 구성 요소로 이루어집니다:

1. 광섬유 전파 및 산란 모델:

   • 광섬유는 선형 복굴절을 가진 웨이브플레이트(waveplate)로 취급되는 일련의 세그먼트로 모델링됩니다.
   • 정적 조건: 고유 복굴절은 스토크스 복굴절 벡터 성분이 독립적인 랑주뱅 과정(Langevin processes)을 따르는 모델을 사용하여 생성되며, 이는 섬유를 따라 복굴절 크기의 레이leigh 분포와 무작위 방향성을 갖습니다.
   • 동적 섭동: 모델은 광탄성 효과(photoelastic effect)와 광섬유 신장을 포함합니다. 모델은 두 가지 이벤트 유형을 구분합니다:
     • 종방향 변형(Longitudinal Strain): 포아송 수축($\epsilon_x = \epsilon_y = -\nu\epsilon_z$)을 동반한 광섬유 신장($\epsilon_z$)으로 모델링됩니다. 이는 광탄성 효과를 통해 공통 위상을 변화시키지만, 추가적인 복굴절을 유발하지는 않습니다 ($\epsilon_x = \epsilon_y$ 이기 때문).
     • 횡방향/비등방성 변형(Transverse/Anisotropic Strain): 불균일한 횡방향 변형($\epsilon_x \neq \epsilon_y$)을 모델링하며, 이는 추가적인 복굴절을 유발합니다. 모델은 스토크스 공간에서 고유 복굴절과 이벤트 유발 복굴절의 벡터 합을 계산하며, 이때 결과적인 크기는 상대적인 fast axis의 방향에 따라 비선형적임을 명시합니다.

2. MIMO-DAS 시스템 아키텍처:

   • 시스템은 섬유를 탐사하기 위해 편광 분할 다중화(PDM)된 BPSK 매핑 직교 골레이 시퀀스(Golay sequences)를 사용합니다. 이 코딩 방식은 피크 광출력을 낮추어 WDM 데이터 전송과 공존하면서도 신호 대 잡음비(SNR)를 유지할 수 있게 합니다.
   • 시스템은 이중 편광 코히어런트 수신기를 사용하여 후방 산란 신호를 포착합니다.
   • 신호 처리: 시스템은 각 세그먼트에 대한 왕복 존스 행렬(Jones matrices)을 추정합니다. 이 행렬들로부터 시스템은 다음을 공동으로 추출합니다:
     • 공통 위상 ($\phi$): 존스 행렬의 행렬식(determinant)으로부터 도출됩니다.
     • 유효 복굴절 ($\Delta\beta_{eff}$): 연속적인 존스 행렬 곱의 고윳값(eigenvalues)으로부터 도출되며, 게이지 길이(gauge length)에 걸쳐 누적된 지연(retardance)을 나타냅니다.

3. 검증 전략:

   • 수치 시뮬레이션: 모델은 시뮬레이션된 종방향 및 횡방향 변형 이벤트가 적용된 3 km 광섬유에서 테스트되었습니다.
   • 실험적 검증: 두 개의 광섬유 스풀(235 m 및 1956 m)과 제어된 종방향 변형(압전 스테이지를 통해) 및 횡방향 변형(광섬유 스퀴저를 통해)을 가하는 셋업을 사용하여 실험실 실험이 수행되었습니다.

주요 기여

• 공동 추정 프레임워크: 본 논문은 코딩된 MIMO-DAS 시스템의 분산 존스 행렬로부터 공통 광학 위상과 유효 선형 복굴절 크기를 공동으로 추정하는 방법을 제안합니다.
• 이벤트 식별: 연구는 두 가지 이벤트 유형의 뚜렷한 시그니처를 기반으로 이들을 식별할 수 있음을 보여줍니다:
  • 종방향 변형: 공통 위상에는 큰 영향을 미치지만, 복굴절(retardance) 궤적에는 거의 영향을 주지 않습니다.
  • 횡방향/비등방성 변형: 공통 위상과 복굴절 궤적 모두에 감지 가능한 시그니처를 생성합니다.
• 정적 복굴절 특성 규명: 모델은 유효 복굴절의 정적 분포를 성공적으로 추정하였으며, 이는 표준 단일 모드 광섬유에 대한 이론적 기대치와 일치하는 레이leigh 분포를 보여주었습니다.
• 실험적 검증: 수치 모델은 실험실 실험을 통해 검증되었으며, 약 1.3 미터의 평균 공간 해상도로 종방향 및 횡방향 변형을 탐지, 국지화 및 식별할 수 있는 능력을 확인했습니다.

결과

• 정적 조건: 시스템은 광섬유를 따라 유효 복굴절 크기를 추정했습니다. 1956 m 광섬유 스풀에 대한 추정값의 분포는 이론적 예측과 일치하는 평균 $\approx 0.57$ rad/m의 레이leigh 분포를 따랐습니다. 또한 모델은 서로 다른 평균 복굴절 값을 가진 연결된 광섬유 스풀들을 성공적으로 구분했습니다.
• 동적 종방향 변형: 시뮬레이션(800 Hz)과 실험(100 Hz) 모두에서 사인파 형태의 종방향 변형은 차동 위상 추정에서는 명확한 시그니처를 생성했으나, 지연(retardance) 궤적에서는 감지되지 않았습니다.
• 동적 횡방향 변형: 비등방성 횡방향 변형(시뮬레이션 및 광섬유 스퀴저를 통해 적용됨)은 차동 위상과 지연 궤적 모두에서 감지 가능한 시그니처를 생성했습니다.
• 공간 해상도: 저해상도 계수(5개 그룹 중 가장 높은 반사 세그먼트를 선택)를 사용할 때, 시스템은 1.28 m(약 1.3 m)의 평균 공간 해상도를 달erm했습니다. 저자들은 민감도와 해상도 사이의 트레이드오프를 언급하며, 해상도 계수를 높이면(예: 10) 대비(contrast)는 개선되지만 공간 해상도는 감소한다고 설명했습니다.

의의 및 주장
본 논문은 코딩된 MIMO-DAS 아키텍처에서 위상과 편광 추정을 모두 활용함으로써, 단일 파라미터 시스템에 비해 이벤트 식별 능력과 횡방향 섭동에 대한 민감도를 향상시킬 수 있다고 주장합니다. 저자들은 이 접근 방식이 등방성(종방향) 및 비등방성(횡방향) 섭동을 모두 탐지할 수 있게 한다고 강조합니다.

그러나 본 논문은 시스템의 정량적 능력에 대해 신중한 입장을 유지합니다. 저자들은 추정된 유효 복굴절이 섭동의 진폭에 따라 선형적으로 변하지 않는다는 점을 명시적으로 밝힙니다. 즉, 그 반응은 고유 광섬유 복굴절과 이벤트 유발 복굴절 사이의 상대적 방향에 따라 비선형적입니다. 결과적으로, 시스템은 변형의 성격(종방향 vs 횡방향)은 식별할 수 있지만, 방향에 대한 추가적인 지식 없이는 현재 변형의 크기를 정밀하게 정량화할 수는 없습니다. 또한 저자들은 실제 환경 구축 시 케이블에서 광섬유로 전달되는 횡방향 변형의 낮은 전달율과 같은 실질적인 한계를 언급하며, 이것이 작은 횡방향 이벤트를 감지하는 데 제한이 될 수 있다고 지적했습니다. 향후 연구 과제로 노이즈 효과에 대한 이해와 현장 테스트 수행을 제안했습니다.