Introducing Coherent MIMO Sensing, a fading-resilient, polarization-independent approach to φ-OTDR

Nowadays, long distance optical fibre transmission systems use polarization diversity multiplexed signals to enhance transmission performance. Distributed acoustic sensors (DAS) use the same propagation medium ie. single mode optical fibre, and aims at comparable targets such as covering the highest distance with the best signal quality. In the case of sensors, a noiseless transmission enables to monitor a large quantity of mechanical events along the fibre. This paper aims at extending the perspectives of DAS systems with regard to technology breakthroughs introduced in long haul transmission systems over the last decade. We recently developed a sensor interrogation method based on coherent phase-sensitive optical time domain reflectometry (φ-OTDR), with dual polarization multiplexing at the transmitter and polarization diversity at the receiver. We name this technique Coherent-MIMO sensing. A study is performed from a dual-polarization numerical model to compare several sensor interrogation techniques, including Coherent-MIMO. We demonstrate that dual-polarization probing of a fibre sensor makes it insensitive to polarization effects, decreases the risks of false alarms and thus strongly enhances its sensitivity. The simulations results are validated with an experiment, and finally quantitative data are given on the performance increase enabled by Coherent-MIMO sensing. © 2020 Optical Society of America under the terms of the OSA Open Access Publishing Agreement


Introduction
Optical fibres, the ones dedicated for sensing applications as much as the ones in already deployed cables in telecom infrastructure, can be employed as distributed acoustic and strain sensors. These sensors can serve in many fields, from marine environment monitoring [1] to medical applications [2]. Often chosen because they are independent of most external physical fields (electromagnetic field impacts copper cables, temperature of the fibre strongly modifies the intensity in Raman fibre sensors) [3], Rayleigh scattering-based sensors exploit the intensity, polarization and phase of light. When a disturbance occurs in the neighbourhood of the fibre, the induced intensity variations (intensity-based optical time domain reflectometry -OTDR) or phase variations (differential phase OTDR) of the backscattered light are captured to localize the disturbance source. The last years have seen the development of coherent φ-OTDR, capable of distributed strain and temperature sensing [4] over long distances and with tunable, up to sub-meter spatial resolution [5]. Coherent φ-OTDR exploits the retro-propagated optical field induced by the Rayleigh backscattering effect in optical fibres, by measuring the backscattered phase of the fibre sensor. The main advantage of phase over intensity is that it offers a linear response to strain [6,7]. It can also capture rapid changes in the disturbances, over high bandwidth. Rayleigh scattering however, is known to be polarization-dependent [8]. Besides, coherent setups are vulnerable to coherent noise also known as interference fading [9,10], not to mention laser-related noise [11] which will be out of the scope of this paper. Coherent fading noise can be suppressed using specific pulse-probing methods [12]. To overcome polarization fading effect, recent work focus on the receiver design and uses a dual-polarization coherent receiver [13,14] so that the overall received optical power is kept constant independently from the state of polarization of the incident signal at the receiver side. A wide spread assumption is that polarization-diversity receivers are the key to mitigate polarization fading [15]. Still, polarization-diversity reception can be improved [16]. What's more, the dual-polarization receiver setup is not sufficient to provide a perfect channel estimation. The transmitter setup must as well be taken into account for fully achieving a polarization insensitive sensor. Actually, polarization induced phase noise is transmitter-dependent, which was observed in [17] and further analysed in [18]. This effect depends on the input state of polarization (SOP). In the following, we compare single and dual polarization interrogation schemes, considering Single polarization Input -Single polarization Output (SISO), Single Input -Multiple Output (1 × 2 SIMO) and Multiple Input -Multiple Output (2 × 2 MIMO) probing techniques, and investigate the issues related to each of them. A numerical model is developed to investigate the separate and joint effects of physical aspects encountered along an optical fibre. The transmitter, sensed fibre, and receiver are modelled as in [19]. Finally, we demonstrate improved sensitivity of 2 × 2 MIMO, and call this sensing technique Coherent-MIMO Sensing.
The paper is organized as follows. Section 2 details the numerical model used for the simulations, under static conditions and with the inclusion of noise contributions. Section 3 investigates the different manners to recover and compute the phase from the received backscattered field, and summarizes the main fading effects that can be identified and mitigated with modern sensing schemes. Section 4 compares the performance of interrogation methods when the number of used polarization channels varies using simulations; in this section, polarization diversity is used at the receiver side, simulations are run to identify and discuss the remaining sensitivity issues. Finally, section 5 details the experimental setup that validates simulations results and assesses the sensitivity of the developed MIMO sensor.

Dual-polarization fibre model
Several existing models investigated the backscattered phase and amplitude from optical fibres. The first step is to consider the fibre under static conditions [20], the second is to add perturbations to the simulated fibre [21]. We restrict ourselves to the static approach and introduce the polarization dimension to the model: the channel is seen as N successive fibre segments of length L s , and each of these segments will be associated not only to the returned phase and amplitude (complex scalar), but to its complete Jones 2 × 2 matrix [19], denoted H, to also cover the polarization state dimension. Our scope is long SSMF (standard single mode fibres), with a highly coherent source (fibre length L f very small compared to source coherence length L coh ). A global phasor p i is derived for each segment i ∈ [1, N], which includes the local Rayleigh backscattering amplitude. In addition, A i the dual-pass attenuation in the fibre applied to segment i, is defined as e 2L i a where L i and a are fibre distance from interrogator to the segment i and linear attenuation in the fibre respectively, such that the received backscattered amplitude and phase from fibre segment i will be A i p i . Extending our scope to dual-polarization, the phase and amplitude are distributed over two orthogonal polarization axes in the fibre. To each segment i ∈ [1, N], a single-pass, unitary Jones matrix U i is associated. It comprises the random polarization effects encountered by the light travelling through the fibre segment (forward transmission). To model these polarization effects, three random unitary matrices are generated: D β i = diag(e jβ i ; e −jβ i ) and D γ i = diag(e jγ i ; e −jγ i ), diagonal phase retarders with γ i and β i uniformly drawn in [−π, π] (x → e x is the exponential function), and and ξ i uniformly drawn in [0, 1] [22]. U i is the following product: where φ is the common phase term to cover all possible unitary matrices. It results from Eq. (1) that U i is a random unitary matrix, which is accurate to describe the polarization state of any static fibre segment [23]. In the following, as the interest is on the differential phase between segments, the common phase φ won't be considered in the expression of U. Generated this way, Jones matrices cover all possible states of polarization (SOP). For verification, it is possible to switch from Jones space to Stokes vectors and project the simulated SOP on the Poincaré sphere: for a significant number of random drawings, the full unit-radius sphere is covered.
As we are interested in the backscattered light, the forward and reflection paths are modelled for each fibre segment. The reflection is placed at the end of the fibre segment. It can be modelled either as a perfect reflection M = [ 1 0 0 1 ] , or considering polarization transfers at the reflection, introducing a transfer coefficient α ∈ [0, 0.05]: . The value of this polarization transfer coefficient has to be measured and verified. The existence of such a transfer coefficient was considered for physical properties of silica [24], independently of optical fibre transmission. It could result in slight polarization changes of the Rayleigh backscattered light [25], which we estimate below 5% of light transferred from one polarization state to its orthogonal state. As it is not considered in most studies on Rayleigh backscatter, α = 0 is set in the following, if not specified otherwise. The optical fibre is reciprocal from a propagation point of view, which translates into U backwar d = U T f or war d , with T standing for the transpose operator [26]1. We assume that the forward SOP is expressed in an (x,y,z) coordinate system where the z axis points towards the propagation axis and the backward SOP is expressed in the (x,y,-z) coordinate system. Assuming that U i can be the forward matrix of segment i as well as the product of all forward matrices from segment 1 to segment i (since a product of random unitary matrices is a random unitary matrix), the full dual-pass Jones matrix from the initial segment up to fibre segment i is then: where A i and p i stand for the dual-pass attenuation and phasor from the fibre start, and where |p i | follows a Rayleigh distribution [20]. Polarization dependent loss (PDL) and polarization mode dispersion (PMD) are considered negligible in the fibre sensor. The general shape of the Jones backscatter matrix is computed as follows for all spatial segments i: 1We change the formalism for the round-trip Jones matrix with respect to the one we used in [19,27] to align with [26] which better describes the SOP evolution after a round-trip in an optical fibre for each particular realization of U forward (β, Θ, γ) (Indeed, if we are only interested in phase noise statistics averaged over several realizations of backscattered SOP, both models are statistically equivalent. However, for each particular realization, the new model better depicts the SOP evolution).
As the SOP of the emitted signal at the transmitter and the incident signal at the receiver side aren't perfectly aligned, a misalignment parameter θ will be considered. It is reported as an additional rotation, for example at the entrance of the system No phase retarder is added as the considered fibre portions at the transmitter and receiver are small compared to the fibre sensor length L f .
Finally, the polarization beat length L pb of the simulated fibre is considered. This way, the speed of polarization variations from one fibre segment to the other is controlled. For single mode fibres, L pb is of the order of magnitude of 10m [28]. By definition if L s ≥ L pb , then the SOP in two successive segments is uncorrelated: polarization parameters are randomly drawn from one segment to the other, the SOP of light returning from segment i + 1 is independent from the SOP of light returning from segment i. Else if L s /L pb → 0, then there is almost no SOP change from one fibre segment to the next. For values of L s /L pb comprised in ]0, 1[, the change in the SOP will be proportional to the ratio L s /L pb , with λ i+1 = λ i + L s L p b × λ N ew , with i the segment index, λ any polarization parameter (referring to β, γ or Θ), λ N ew randomly drawn as described in the beginning of this section and independent from λ i .
Several simulations for different L s /L pb ratios are shown in Figure 2, where it appears that the larger L s /L pb , the more random the trajectory of the SOP of light in the fibre (including attenuation A i , thus not of constant modulus on the sphere).   Figure 2(c) that from 10L s > L pb , the distribution of the SOP over the Poincaré sphere is already random. When targeting sensing applications over telecom fibres, the spatial resolution L s is typically in the range of 1m to 10m (specific fields such as medicine are out of the scope of our study). Therefore, and for sake of simplicity, the simulation study below is conducted with uncorrelated SOP from one segment to the next (L s ≥ L pb case). Notice that additional simulations performed with L s 10 −1 L pb and L s ≥ L pb in the sensor, which are not displayed here, lead to similar results.

DAS system model
Similar to [21], our sensing model is considered as an ensemble of building blocks, so that elements can be added and removed at will to test their respective influence on measurement accuracy, noise, or phase retrieval. The transmitter encloses a narrow linewidth laser, modulated by a specific pattern. Laser phase noise here is modelled as a Wiener process of variance σ 2 = 2π∆νT S where ∆ν is the laser linewidth, T S is the sampling period. Although this does not perfectly match the physical reality, for which laser frequency noise increases in 1/ f d , d ≥ 1 [11], this enables to predict the general evolution of the error on phase estimation over distance [27]. The propagation medium is the fibre, the model of which is described in the previous section.
The receiver is made of a coherent mixer, subject to insertion losses and additive white noise caused by thermal and quantum noise from the photodiodes. The white noise of the receiver is considered as σ 2 AW G N = σ 2 th + σ 2 s + σ 2 RI N , where the contributions of the shot noise σ 2 s and thermal noise σ 2 th are dominant compared to the relative intensity noise (RIN) in our case. Altogether, we previously showed [19] that the laser phase noise is an important limitation in such a system.

Phase estimation in polarization sensitive sensing
This section focuses on the phase estimation methods to be used according to the number of used polarization states at the transmitter and receiver side.
The following notations are adopted: SISO stands for Single polarization Input -Single polarization Output sensing, which means that the light sent to the fibre sensor is linearly polarized in one direction and that the backscattered light is recovered on one polarization at the receiver side. SIMO stands for Single polarization Input -Multiple polarization Output sensing, which is often referred to as "polarization diversity receiving in φ-OTDR": linearly polarized light is sent to the sensor, typically with a narrow linewidth laser, backscattered light is recovered on both polarizations using a polarization diversity receiver. Finally, MIMO stands for Multiple polarization Input -Multiple polarization Output sensing, where light is sent on two orthogonal polarization states at the transmitter side, and polarization diversity detection is used at the receiver. MISO setup is also considered: Multiple polarization Input -Single polarization Output sensing, where polarization diversity is used at the transmitter but not at the receiver.
The Jones formalism is used. The dual-pass fibre channel is noted H, or H i considering the particular segment i.
h y x h y y ], with complex scalar coefficients.

Phase extraction
MIMO allows to recover the full Jones matrix H of the channel. Assuming no attenuation A i and assuming H has a common backscattered phase φ such that With H being the dual-pass matrix of the channel, φ is the dual-pass phase. Hence the use of the following phase estimator:φ where det H is the determinant of H. SIMO probes a single polarization. For a recovered signal on both polarizations, if the first row of H is sent: Note that in Eq. (4) we also haveφ SI MO = arctan( ) with , the real and imaginary parts of h xx and h yx data [13]. SISO recovers information on one polarization only (h xx for example), such that the estimated phase isφ SI SO = ∠(h xx ). MISO proceeds the same way, and recovers h xx + h xy with no additional information to differentiate both components. Estimated MISO phase is then Existing SISO and MISO configurations do not necessarily include a coherent receiver with polarization diversity but rather a single polarization or single 90 • hybrid direct detection receiver, and use all-optical techniques to recover the phase such as delay lines, dual pulses, chirped pulses... Using direct detection, there is no sensitivity to SOP variations of the incoming light: the full intensity of the signal is received [29]. Moreover, in direct detection and coherent source scheme, it is possible to have a rather stable φ-OTDR trace from comparisons of successive power-versus-time traces. Further trace-averaging can be performed if a speckle pattern is still noticeable [30]. Also, it was noticed that considering polarizations separately could bring a better accuracy to the system [31].
Here, coherent detection -coherent φ-OTDR is used. Coherent receivers are sensitive to the SOP of incoming light. To recover the full signal, there is a need to use a polarization diversity scheme, where each polarization component is detected by projecting the incident signal over two orthogonal polarization states and beating each projection with the same local oscillator (LO) [29]. With only one polarization projected onto the LO laser, a part of incoming light can be lost. Our SISO case refers to the latter situation.

Common fading phenomenons
Several fading types were mentioned in the introduction. Now that the phase estimators are defined, this section aims at linking phase estimators to their respective fading risks.
Interference fading, or coherent fading [10], also referred to as Rayleigh fading and fading noise [32], is caused by coherent interference within the same pulse due to the presence of scatterers in the fibre and to the use of a highly coherent source, resulting in a jagged intensity trace at the reception [9]. If a vibration occurs, the optical path changes and so does the interference figure. It is the reason why intensity variations due to a perturbation have no linear relationship with that perturbation. Such variations are reported as an issue for OTDR techniques that rely on intensity, or that need high spatial resolution [12,30], since intensity variations affect locally the sensor response. However, coherent φ-OTDR is based on that effect, allowing for precise phase measurement within a short time and thus over a high bandwidth, hence it shouldn't be mitigated using low-coherence sources [3,30] but rather avoided. A way for avoiding such an effect would be to use multiple channels [33] to probe the line.
Polarization fading was identified in interferometric sensors [34,35], and remains an issue in coherent detection DAS. In more recent fibre sensors with coherent detection, polarization effects can prevent detection of the strength of perturbations [36]. Polarization fading occurs when the incident SOP and the SOP of the LO become orthogonal, which prevents mixing on the photo-detector in interferometers, and which fades the measure in heterodyne or homodyne coherent systems.
Referring to the notations in subsection 3.1, polarization fading strongly impacts SISO configuration: when detecting h xx , an intensity fading will occur for periodic (Θ, β) pairs corresponding to orthogonality between the emitted x-polarization and the received x-polarization, resulting in no reliable phase estimation if the detected intensity is too low. Conversely, SIMO and MIMO modes are supposed to be resilient to polarization fading, meaning there will always be a phase estimation that is stable in time: φ SI MO = ∠(p i × (e j2γ i cos 2β i − j sin 2β i (e j2γ i cos 2Θ i + sin 2Θ i ))) (here Eq. (4) is applied to H with α = 0) and φ M I MO = ∠(p i ). MISO also allows to recover an estimation as its expression is similar to that of SIMO using Eq. (4): even if the receiver SOP is orthogonal to polarization Y, polarization X will be recovered, and vice-versa. Of course, MISO and SIMO-estimated absolute phase differ from the real common phase due to polarization effects (Θ, γ, β appear in their expression) but when sensing dynamic events, we are only interested in phase variations compared to a stable reference.
Recent and older results, together with our definitions of SISO and SIMO, agree on the need for multiple receivers, should it be in direct detection [37] or by employing polarization diversity system in coherent OTDR [15,38]. In fibre interferometers, using depolarized light at the transmitter limits the polarization fading effects [18]. MISO case is not considered in the literature. The polarization beat length was also mentioned, as it could play a role in the quality of the received signal when confronted to polarization fading [21]. This parameter has to be considered in future works. The interest of the mentioned solutions against fading is discussed below, and comments are given on the remaining fading types along the sensor.

Impact of polarization diversity on phase estimation
With phase OTDR, a local phase information is processed as a differentiation between the phase estimated at one fibre segment location and that estimated at one previous location. The first fibre segment serves as a reference. Repeating this phase estimation process in time allows to capture and localize mechanical perturbations that may affect the fibre. The phase retrieval is detailed below.
In static mode, for two segments i and j with j > i, the measured differential phase is constant over time and given by φ st atic = ϕ j − ϕ i . In dynamic mode, assuming a strain is applied locally at segment position i, the phase is linearly disturbed relatively to the perturbation, following: where n is the refraction index, ξ photo-elasticity coefficient (= 0.79 for silica), δl the fibre strain. So if i is disturbed, the measured phase is φ dyn = ϕ j −ϕ i +2∆ϕ single−pass = φ st atic +∆ϕ dual−pass . The phase delay is applied twice since the light does a round-trip in the segment (way faster than the mechanical disturbance action, so the measurement of mechanical perturbations is considered instantaneous). For the detection of an event, we chose to monitor the standard deviation of the differential phase over time, along the fibre distance (denoted StDv). When an event occurs along the sensor, phase variations appear in time, thus increasing the phase standard deviation locally. The displayed traces in this section are simulations of a coherent detection -coherent φ-OTDR, where the phase is recovered according to the estimators defined in subsection 3.1. All simulations are static, e.g with no strain applied to the sensor. Therefore, this study focuses on the false alarm issue, not onto the non-detection case: the main requirement for phase sensors is that no variation or "false alarm" occurs when the fibre is not disturbed, which is also the condition to achieve the highest sensitivity. A simulation study of non-detection case would require a dynamic model, which is out of the scope of this paper.

Simulations: accuracy and false alarms
Firstly, the global behaviour of the StDv along sensor distance is considered to highlight the difference between the given phase estimators. Secondly, paying attention to the StDv peaks, we show how they relate to the polarization rotation and birefringence in the fibre segments. Figure 4 gives the phase StDv in time along fibre distance for static simulations. The coded sequence duration is 5.24ms at a symbol rate f symb = 25Mbaud, with a constant power of 7dBm at the entrance of the fibre sensor and σ AW G N = 1.7mV RM S at the receiver. The native gauge length is 4m. Figure 4(a) is a zoom on a 1km section with 4 different phase estimators (SISO, SIMO, MISO, MIMO) for the same simulated fibre, with a spatial step of 20m ; Fig. 4(b) displays the phase StDv averaged after 50 independent 50km long fibre simulations for SISO, SIMO and MIMO phase estimators, and for a theoretical case with no polarization effects along the fibre with a spatial step of 200m. This theoretical case is derived by taking only coherent noise into account, i.e. U i in Eq. (1) becomes the identity matrix for all fibre segments. Figure 4 shows a clear hierarchy between no polarization diversity (DP) at all, DP at one side, and DP at both TX (transmitter) and RX (receiver) sides. The MISO case is not displayed in Figure 4(b) since it is comparable to SIMO. As depicted in Figure 4(b), the general trend of the phase StDv over long distance is to grow exponentially with distance, with a slope depending on the laser linewidth ∆ν [19,27]. The value ∆ν = 75Hz is used in the simulations since it corresponds to the laser used in the experiments in next section [19]. What differentiates the phase StDv traces of the different phase estimators is the amount and dynamic of StDv peaks, or false alarms, as shown in the single fibre simulation from Figure 4(a) (e.g SISO StDv peaks are almost two orders of magnitude higher than MIMO StDv peaks). Deriving the mean value over several simulations as displayed in Figure 4(b) shows how variable the phase estimation tends to be according to the chosen estimator. We also demonstrate in Figure 4(b) that considering no polarization effects or using MIMO probing will lead to the same phase StDv values. On average, MIMO is matching the theoretical fading-free phase StDv profile: it is polarization-independent. What comes out is that MIMO estimated phases are the most stable with no StDv value above 5.10 −2 rad over the first 5km, while SISO can vary a lot locally and quite often over fibre distance with 6 peaks above 10 −1 rad within one km in Figure 4(a) for example. SIMO and MISO show an intermediate behaviour with randomly distributed StDv peaks, here 1 to 3 StDv peaks above 10 −1 rad within 1km. It was shown in [39] that SISO performance varies with the input as well as with the output polarization state when probing the line. Figure 4(b) shows that after 50km of fibre sensor, the StDv decreases by a factor of 5 and 2 when going from SISO to MIMO and SIMO to MIMO respectively. Furthermore, several previous works discussed in subsection 3.2 have already introduced polarization diversity at the receiver side to improve the sensing quality. Therefore, for sake of clarity, the rest of the study concentrates on the comparison between SIMO and MIMO cases only.
Regardless of the phase estimation method, the evolution of the state of polarization in the fibre is described mainly by Θ, β, γ and θ defined in subsection 2.1 provided that the Rayeligh scattering is an ideal reflector (no polarization crosstalk at reflection stage). In Figure 5(a), a 3-D view involving the evolution of angles Θ and θ for fixed randomly drawn values of γ and β together at a single segment location with a fixed reflected intensity is used to track the StDv behaviour. It shows that the occurence of sudden phase variations (errors on the phase estimation, or "false alarms") in SIMO is not random, but rather concentrated at some polarization parameters. Specific (Θ, θ) pairs at fixed random γ and β trigger StDv peaks in SIMO, which is quite clear in Fig. 5(a) and (b). On the same fibre simulations, the corresponding MIMO estimation is performed in Fig. 5(c) and (d): the estimated phase is stable for all possible (Θ, θ) pairs. If α 0 where α is the polarization transfer coefficient defined in subsection 2.1 (i.e. M is not an ideal reflector and introduces some crosstalk), then the locations of StDv peaks shift as shown in Figure 5(b). The noise floor slightly rises in MIMO compared to α = 0. The laser phase noise ∆ν impacts the height of the StDv peaks and StDv values and not the general shape (position of the StDv peaks in the (Θ, θ) plane), so the case ∆ν = 0Hz is not displayed. From Figure 5, we may suspect the phase estimation to be modulated by some (Θ, θ) "pathological" pairs which degrade the backscattered information. This assumption is further investigated below.
Back to subsection 3.1, the SIMO phase estimator was defined as: φ SI MO = ∠(p i × (e j2γ i cos 2β i − j sin 2β i (e j2γ i cos 2Θ i + sin 2Θ i ))) if we dismiss α and θ. It is recalled that Θ stands for the random polarization rotation in the fibre segment, and β, γ are phase retarder coefficients, defined in subsection 2.1. This estimator is declared immune to intensity fading issues as it allows to get a phase estimation at all locations in the fibre, regardless of the state of polarization.
However, the phase ∠p is modulated by the coefficient c mult,SIMO = e j2γ cos 2β − j sin 2β × (e j2γ cos 2Θ + sin 2Θ) This coefficient is plotted in Figure 6(a) (for a given input angle θ and a randomly chosen γ). Some (Θ, β) pairs lead to extremely low values for that coefficient, which highly attenuate or "fade" the received phase value. For α 0, this phase-fading coefficient c mult,SIMO becomes a function of γ, β, Θ, α and produces a similar, still less regular pattern than in Figure 6(a). This is illustrated in Figure 6(b), where a constant γ is chosen in [−π, π], and α is set to 0.15. This configuration also has fading combinations of γ, β, Θ, shown here for Θ, β varying in [−π, π]. It is possible to predict the values of the polarization parameters for which they cause such phase-fading.
Note that there is no phase fading with MIMO phase estimation since the equivalent coefficient for MIMO is c mult,MIMO = abs(det H) = 1 since H is a unitary 2 × 2 matrix.

Expectations on sensing quality
To begin quantifying the performance increase brought by MIMO probing and its associated phase estimation compared to SIMO, a study is conducted by means of simulations. Here, the fibres are probed using complementary Golay sequences on either one or two polarizations. The use of these complementary sequences allows for a mathematically perfect estimation of the backscattered optical field, considering either a single polarization channel, or a dual-polarization one, as demonstrated in [6]. The statistical study is conducted on a high number of simulations (N = 10 6 simulated points, corresponding to 2000 fibre generations for 340 m length, 40 fibres for 25 km and 20 for 50 km), with f symb = 50Mbaud, corresponding to 2m gauge length. MIMO and SIMO probing simulations are run on randomly generated fibres, and the differential phase StDv is derived for every fibre segment. The polarization beat length L pb is selected such that L pb ≤ L s , assuming independent polarization state between consecutive segments (see subsection 2.1).
Three different fibre lengths are simulated to investigate the influence of distance on the StDv distribution according to the phase estimation method, so to compare the relative robustness of the SIMO and MIMO estimators over distance.  Figure 7(b) and (c)). The distributions shape change with distance, however the same differences are observed between SIMO and MIMO cases. We show that in Figure 7, between 35% and 50% of the SIMO estimated phases vary more than 75% of MIMO estimated phases over 50km and 340m respectively, 9% to 15% of the SIMO estimated phases vary more than 95% of MIMO estimated phases over 50km and 340m respectively.
Considering that the phase StDv is a measure of noise, we define the signal to noise ratio as SNR phase = 10 log 10 (1/σ 2 phase ). The mean value of the SNR phase for MIMO is 1dB above SIMO for all simulated lengths (50km, 25km, 340m). Moreover, the additional variance of SNR phase in SIMO compared to MIMO method is of 3dBrad 2 at 340m, 11dBrad 2 at 25km and 7dBrad 2 at 50km. We conclude that MIMO enhances the mean achieved SNR and also reduces its fluctuations with respect to SIMO.
From a high number of random simulations, it appears that SIMO phase estimations bring in more variations than MIMO, for the same simulated fibre with same scatterers distribution and birefringence, and for all distance ranges. The performance gap between SIMO and MIMO methods is striking for short distances, and remains significant at very long distances (50km here) when coherent fading and attenuation also affect the measurement. MIMO sensing appears as a strong, stable phase estimator on short fibre lengths as well as on longer distances. An experiment is performed on 340m SMF, with no perturbations around the fibre so to stay in static mode as in the simulations. The fibre is interrogated with BPSK sequences on one (SIMO) or both polarizations (MIMO). The setup is described in Figure 8 and its operating mode further detailed in [6]. SIMO and MIMO measurements could not be performed simultaneously, however they are performed successively to ensure close experimental conditions for both methods. The mechanical bandwidth is of 1.5 kHz (50 Mbaud, 8192 symbols code length), from which the first 50 Hz bandwidth is filtered out so to avoid capturing the noise of air conditioning and other low frequency background noise sources usually met in buildings. The measured phase StDv, calculated from a 1s measurement window with 2m gauge length, along with its distribution are plotted in Figure 9.

Experiment and sensitivity
Note that there is only a short distance of fibre here, so the overall measurement shows low phase standard deviation values, well below 0.1rad. A first look at Figure 9(a) already confirms the gap of performance between SIMO and MIMO interrogation methods. The peaks of StDv are not perfectly aligned due to mere synchronization issues, but overall the highest StDv values for MIMO are located where the StDv is also high for SIMO. In these locations, the speckle pattern is suspected: coherent fading is still an issue in the fibre sensor, though it is the only one remaining with MIMO.
The SIMO estimated phase on half of the fibre segments varies more than the MIMO estimated phase on 75% of the (same) fibre segments. It is consistent with the values obtained through As for the SNR phase values, and despite the restricted number of measured points, the variance of SNR phase is 2.3dBrad 2 higher in SIMO compared to MIMO, which is consistent with the simulation value (3dBrad 2 at 340m). Although no SIMO probing was experimented on longer fibre, it is expected that the relative distribution of SIMO and MIMO is not supposed to vary much. In a recent work [40], we demonstrated a good stability of MIMO phase estimation on several tens of km. Therefore, these experimental results, in line with previous simulation ones, confirm the better sensitivity gain brought by a Coherent-MIMO approach in phase OTDR.

Conclusion
The exploitation of two orthogonal polarization channels in an optical fibre was highlighted in the case of optical fibre sensing. From a numerical model of the Rayleigh backscattering in singlemode fibres, simulations were run for different configurations of transmitter and receiver. Several phase estimation methods for coherent φ-OTDR were compared, among which Single polarization Input -Multiple polarization Output sensing (SIMO), often referred to as "polarization diversity receiving in φ-OTDR" in the literature, and Multiple polarization Input -Multiple polarization Output sensing (MIMO). It was demonstrated that some sensing configurations, in particular those with no polarization diversity, were sensitive to fading effects and could neither allow to correctly recover the intensity nor the phase information from all fibre locations. Even with polarization diversity at the receiver, we pointed out that SIMO estimated phase is subject to a modulation which biases the phase estimation for some fibre segments and triggers local phase variations. It was shown both in simulation and through experiment performed in the absence of mechanical perturbation that SIMO estimated phases vary considerably more than MIMO estimated ones, for the same fibre in the same conditions: Coherent-MIMO, by completely mitigating the effects of polarization, appears twice more sensitive and more reliable.
In essence, we were able to thoroughly compare existing interrogation methods for coherent optical fibre sensors, and so to demonstrate the superiority of Coherent-MIMO over partial polarization diversity sensing techniques in terms of sensitivity. It proves that further research on dual-polarization interrogation of fibre sensors is of great interest, and gives insights on the key concerns towards a more reliable and sensitive fibre sensing.