Signal Analysis, Signal Demodulation and Numerical Simulation of a Quasi-Distributed Optical Fiber Sensor Based on FDM/WDM Techniques and Fabry-Pérot Interferometers

In civil engineering quasi-distributed optical fiber sensors are used for reinforced concrete monitoring, precast concrete monitoring, temperature monitoring, strain monitoring and temperature/strain monitoring. These quasi-distributed sensors necessarily apply some multiplexing technique. However, on many occasions, two or more multiplexing techniques are combined to increase the number of local sensors and then the cost of each sensing point is reduced. In this work, a signal analysis and a new signal demodulation algorithm are reported for a quasi-distributed optic fiber sensor system based on Frequency Division Multiplexing/Wavelength Division Multiplexing (FDM/WDM) and low-precision Fabry-Pérot interferometers. The mathematical analysis and the new algorithm optimize its design, its implementation, improve its functionality and reduce the cost per sensing point. The analysis was corroborated by simulating a quasi-distributed sensor in operation. Theoretical analysis and numerical simulation are in concordance. The optimization considers multiplexing techniques, signal demodulation, physical parameters, system noise, instrumentation, and detection technique. Based on our analysis and previous results reported, the optical sensing system can have more than 4000 local sensors and it has practical applications in civil engineering.


Introduction
Many practical applications for fiber optic sensors require measuring the same parameter in different spatial locations. For such applications, quasi-distributed sensors can be applied. This sensor type has received significant research attention and it necessarily applies some multiplexing technique. Several basic multiplexing techniques such as time-division multiplexing, wavelength-, coherence-, frequency-, spatial-, and code-division multiplexing have been proposed, demonstrated and applied over the last 10 years [1][2][3][4][5][6][7][8][9][10]. However, the multiplexing capacities of many of these techniques are limited to only a few sensors due to various factors including cross-talk, limited power budget and wavelength bandwidth. To increase the multiplexing capacity, several hybrid systems using two or more basic techniques have been developed [11][12][13][14][15]. These hybrid systems offer us some benefits such as cost savings per sensing point and enhancement of the competitiveness of fiber-optic sensors in their rivalry with conventional sensor technologies.
In reference [12], a quasi-distributed sensor based on wavelength-frequency division multiplexing was experimentally proposed. The sensing system had two wavelength channels and each wavelength channel had four frequency channels. Thus, the quasi-distributed sensor had eight Fabry-Pérot interferometers where each interferometer acted as a local sensor. Later, the same sensing system was applied for temperature measurement. In this case, the optical sensing system had three wavelength channels and each channel had three frequency channels. There were then nine Fabry-Pérot fiber sensors along the single array of the optical fiber. The authors reported only a resolution of 0.125 • C and a sensitivity of 10 pm • C [13]. However, each Fabry-Pérot sensor has its own resolution and its own sensitivity since each sensor has its own physical parameters as reported in references [16,17].
The optical system reported in [12,13], can be theoretically analyzed after stating three principal problems: (a) a sensing system based on frequency division multiplexing, (b) a sensing system based on wavelength division multiplexing and (c) a sensing system based on frequency-wavelength division multiplexing. Problems (a) and (b) were reported in [16,17]. Both theoretical analyses considered the capacity of the multiplexing technique, the instrumentation, the system noise, the signal demodulation, the detection technique and the local sensor properties. Both theoretical analyses had their design and functionality optimized, and both were corroborated by simulating the sensing system in operation. However, to our knowledge, there are no analytic analyses of the third problem (c), so the limitations of quasi-distributed sensors are not known. In this work, the third problem (c) is considered and therefore a quasi-distributed fiber sensor based on frequency-wavelength (FDM/WDM) division multiplexing is theoretically analyzed and numerically simulated. This analysis optimizes its design and its implementation. The analysis considers the capacity of multiplexing technique, local sensor properties, signal demodulation, detection technique, instrumentation and noise. To verify the theoretical analysis and a new demodulation algorithm, we performed a numerical simulation of a quasi-distributed sensor which has four wavelength channels and each channel has three frequency channels. Thus, the simulated quasi-distributed sensor has twelve Fabry-Pérot sensors. The local sensors are low-precision Fabry-Pérot interferometers which are formed by two identical Bragg gratings. All Bragg gratings must have low reflectivity (~1%) and the same length. From the numerical results, each Fabry-Pérot sensor has two resolutions due to the Fourier Domain Phase Analysis algorithm [18][19][20] as reported in [16,17], confirming the optimal functionality of our optical sensing system under study. Theory and simulation results are in concordance. Figure 1 shows the optical system under study. The system consists of a light source, a 50/50 optical circulator, an optical spectrum analyzer (OSA spectrometer), a personal computer used for signal processing and the quasi-distributed sensor. Particularly, the quasi-distributed sensor consists of a serial arrangement of local sensors where each sensor is an interferometer formed by two identical Bragg gratings imprinted in a single-mode fiber. All Bragg gratings have low reflectivity, ≈1% and their lengths are similar. Each interferometric sensor acts as a low-precision Fabry-Perot interferometer [13,16,17]. In the quasi-distributed sensor, Wavelength Division Multiplexing (WDM) and Frequency Division Multiplexing (FDM) techniques were combined. The spatial resolution L SR eliminates ghost interferometers.  [13,16,17]. In the quasi-distributed sensor, Wavelength Division Multiplexing (WDM) and Frequency Division Multiplexing (FDM) techniques were combined. The spatial resolution eliminates ghost interferometers. If the quasi-distributed fiber optic sensor does not have external perturbation due to temperature or strain, the reflection spectrum detected by the OSA spectrometer can be expressed as:

Optical System and its Reflection Spectrum
or: where is the total reflection spectrum, = 1,2, … , = 1,2, … , are the interference patterns produced by the Fabry-Pérot sensors. The light spectrum consists of a series of wavelength channels and each wavelength channel contains a series of frequency channels. Therefore, based on our previous work [16,17], the maximum number of interference patterns is: K × M indicates the maximum number of local optical sensors and the parameters are: Δ is the dynamic range for all Fabry-Pérot sensors, is the k-th Bragg wavelength (wavelength channel), is the effective refraction index of the core, is the length of the gratings, Δλ is the spectrometer resolution and = − is the working interval: is the maximum wavelength and is the minimum wavelength. In this sensing system, all interference patterns have the same form: If the quasi-distributed fiber optic sensor does not have external perturbation due to temperature or strain, the reflection spectrum detected by the OSA spectrometer can be expressed as: or: where R T (λ) is the total reflection spectrum, R km (λ)(k = 1, 2, . . . , K y m = 1, 2, . . . , M) are the interference patterns produced by the Fabry-Pérot sensors. The light spectrum R T (λ) consists of a series of wavelength channels and each wavelength channel contains a series of frequency channels. Therefore, based on our previous work [16,17], the maximum number of interference patterns is: K × M indicates the maximum number of local optical sensors and the parameters are: ∆λ op is the dynamic range for all Fabry-Pérot sensors, λ BGk is the k-th Bragg wavelength (wavelength channel), n is the effective refraction index of the core, L BG is the length of the gratings, ∆λ is the spectrometer resolution and λ w = λ max − λ min is the working interval: λ max is the maximum wavelength and λ min is the minimum wavelength. In this sensing system, all interference patterns have the same form: The physical parameters are: a nm , the amplitude factors, π the constant 3.1415, n 1 the amplitude of the effective refractive index modulation of the gratings, λ is the wavelength and L FPkm is km-th cavity length. Analysis of the interference pattern (4); each interference pattern consists of two functions: enveloping and modulating. The enveloped is a sinc function which is the reflection spectrum of the gratings, the width ∆ BGk is the spectral distance between its +1 and −1 zeros: The modulating term of the cosine function whose frequency ν FPkm is: The cosine function was produced due to the interference between two beams generated from the interferometry sensor. On the other hand, if the quasi-distributed sensor experiences an external perturbation due to temperature or strain, the physical variable affects the Fabry-Pérot sensor [18]. In contrast, the interference patterns have a small shift which is proportional to the magnitude of the physical parameter. In this case, the optical spectrum detected by the OSA spectrometer can be expressed as: or: R T (λ, δλ) is the reflection spectrum due to the external disturbances and δλ km is the km-th Bragg wavelength shift due to the measured change.

Frequency Spectrums
To calculate all frequency components of our optical spectrum R T (λ), we apply the Fourier transform: R T (ν) is the frequency spectrum for the optical signal detected by OSA spectrometer if and only if the quasi-distributed fiber optic sensor does not have external perturbations. Substituting Equation (2) and Equation (4) into Equation (9), the frequency spectrum can be calculated through: invoking the Fourier transform properties, convolution properties, series properties and using the identities: cos 2 (ϕ) = 1 2 [1 + cos(2ϕ)] and cos(ϕ) = e iϕ +e −iϕ 2 , then solving the frequency spectrum R T (ν) is: R T (ν) is two series of triangle functions and the function tri(x) has the following definition tri( otherwise , c km are amplitude factors, ν FPkm is the km-th central frequency which was defined by Equation (6) and ν BGk is the bandwidth: From Equation (11), the spectrum R T (ν) consists of 2(K × M) + 1 triangle functions (frequency components) and their central frequencies are between −ν FPkm and ν FPkm . The component ν FP00 = 0 contains information from all Fabry-Perot sensors while the lateral components contain information from each Fabry-Perot sensor. Positive semi-plane and negative semi-plane contain the same information. If the quasi-distributed sensor has external perturbation, the perturbed frequency spectrum can be calculated through: Substituting Equation (8) into Equation (13): invoking the shift property and solving the transformation, the perturbed frequency spectrum will be: Analyzing the Equation (15), the spectrum R T (ν, δλ) is the product of R T (ν) and two sets of phases. We notice that the phases contain the information about the perturbations and then the Fourier Domain Phase Analysis (FDPA) algorithm can be applied for the signal demodulation [18][19][20].

Sensor Conditions
To obtain the optimal signal demodulation is important to know the sensor limits. Therefore, we mention some requirements for the optical sensing system presented in Figure 1 [16][17][18]: First condition: The cavity length L FP should be in the interval of: This condition ensures that the OSA spectrometer is able to detect the optical signal R T (λ) and there is not overlapping between the components ν FP11 and ν FP00 .
Second condition: The Bragg gratings should have approximately the same length and their reflectivity is ≈ 1%. Thus, the cross-talk noise is eliminated. Third Condition: The spatial resolution should be given by This condition eliminates ghost interferometers. Fourth condition: The minimum distance between centers of gratings for the shortest interferometers is given by: This condition eliminates overlapping between triangle functions. Fifth condition: The number of samples is given by: where λ BG1 is the first Bragg wavelength (first wavelength channel) and L FPKM is the maximum cavity length (last frequency channel and last wavelength channel). The relation (19) satisfices the Nyquist theorem, considers the detection technique, the instrumentation and the sensor parameters.
Sixth condition: The maximum number of Fabry-Pérot sensors is given by Equation (3). This Equation considers both multiplexing techniques (WDM and FDM), sensor parameters, optical source, optical instrumentation and signal demodulation. If the six conditions are satisfied, then the sensing system would have the optimal operation. Figure 2 shows, schematically, the signal demodulation procedure for the sensing system presented in Figure 1. The signal demodulation applies the Fourier Domain Phases Analysis (FDPA) algorithm and two banks of filters. The FDPA algorithm was described and also was applied in references [16][17][18][19][20]. The first bank of K filters is defined as [16]:

Signal Demodulation
The bank is a series of rect(·) functions in the wavelength domain where the rect function is , ∆λ op is the dynamic range and λ BGk is the wavelength channel.
The second bank of M filters is given by [17]: The bank is a series of rect(·) functions in frequency domain and each rect function definition is , ν BGk is the bandwidth and ν FPkm is the centering frequency.
From Figure 2, the signal demodulation algorithm consists of two phases: calibration and measurement. In the calibration phase, the references are estimated and five steps are required: (a) the optical signal R T (λ) is acquired with the OSA spectrometer; (b) using the bank of K filters is filtered a wavelength channel, R m (λ) = F(λ)R T (λ); (c) the frequency spectrum R m (ν) is estimated through: , is a series of triangle functions; d) a frequency channel (a triangle function) , is filtered developing the operation , = F , ; e) the relative phase , is calculated using both spectra * and , ; f) the ambiguity 2 is eliminated through a linear regression and then the absolute phase , is calculated; g) an adaptive filter is applied, and h) the Bragg wavelength shift is calculated. To minimize the noise influence and provide the best estimate, the absolute phase is multiplied by a set of coefficients (step g). Those coefficients act as adaptive filters [20]. R m (ν) is a series of triangle functions [13,16]; (d) a frequency channel (a triangle function) R m (ν) is filtered using R m (ν) = F(ν)R m (ν) and finally (e) its complex conjugate is estimated R * m (ν) where the symbol * indicates the complex conjugate. R * m (ν) is the reference for one Fabry-Pérot sensor and then steps a-d are essential for each interferometer sensor.
In the measurement phase, eight steps are required: (a) the optical spectrum R T (λ, δλ) is acquired using the OSA spectrometer; (b) applying the bank of K filters is filtered the signal R m (λ, δλ) = F(λ)R T (λ, δλ); (c) the frequency spectrum R m (ν, δλ) is determined by: R m (ν, δλ) is a series of triangle functions; (d) a frequency channel (a triangle function) R m (ν, δλ) is filtered developing the operation R m (ν, δλ) = F(ν)R m (ν, δλ); (e) the relative phase ϕ m,rel is calculated using both spectra R * m (ν) and R m (ν, δλ); (f) the ambiguity 2πP is eliminated through a linear regression and then the absolute phase ϕ m,abs is calculated; (g) an adaptive filter is applied, and (h) the Bragg wavelength shift δλ m is calculated.
To minimize the noise influence and provide the best estimate, the absolute phase is multiplied by a set of coefficients (step g). Those coefficients act as adaptive filters [20].

Parameters
To verify our signal analysis and our signal demodulation algorithm, a quasi-distributed fiber optic sensor based on the wavelength/frequency multiplexing techniques and low-precision Fabry-Pérot interferometers was numerically simulated. The simulated quasi-distributed sensor can be seen in Figure 1. The sensing system consists of four wavelength channels and each wavelength channel has three frequency channels, therefore, the sensing system has twelve Fabry-Pérot sensors. Their physical parameters are listed in Table 1.  10 0 to √ SNR = 10 4 . Typical Bragg gratings with rectangular profiles, a refractive index modulation were used. The number of samples was 4096. Given that, the Fast Fourier transform algorithm was considered. The spatial resolution was L SR = 50 cm. For each Fabry-Pérot sensor, the reference spectrum and 50 measurements were simulated. These measurements were within the interval in Table 2. In the simulation a GHIA computer with 8 GB RAM memory and CPU frequency 3.6 GHz was used.

Reflection Spectrum
Applying the parameters presented in Section 6.1, the optical signal R T (λ) was computed. The normalized spectra can be observed in Figure 3. The composite reflection spectrum of a multi-point Fabry-Perot sensor is a superposition of reflection spectra of all local sensors. From Figure 3, the first wavelength channel (λ BG1 = 1536 nm) has its bandwidth of ∆ BG1 ≈ 3.2 nm, the second wavelength channel (λ BG2 = 1542 nm) has its bandwidth of ∆ BG2 ≈ 3.24 nm, the third wavelength channel (λ BG3 = 1548 nm) has its bandwidth of ∆ BG2 ≈ 3.3 nm and the fourth wavelength channel (λ BG4 = 1554 nm) has its bandwidth of ∆ BG4 ≈ 3.308 nm. Each wavelength channel has its own value and they have a very small variation between them. The working interval is λ w = 24 nm: λ max = 1556 nm and λ min = 1532 nm. spectrum and 50 measurements were simulated. These measurements were within the interval in Table 2. In the simulation a GHIA computer with 8 GB RAM memory and CPU frequency 3.6 GHz was used.

Reflection Spectrum
Applying the parameters presented in Section 6.1, the optical signal was computed. The normalized spectra can be observed in Figure 3. The composite reflection spectrum of a multi-point Fabry-Perot sensor is a superposition of reflection spectra of all local sensors. From Figure 3, the first wavelength channel = 1536 nm has its bandwidth of Δ ≈ 3.2 nm, the second wavelength channel = 1542 nm has its bandwidth of Δ ≈ 3.24 nm , the third wavelength channel = 1548 nm has its bandwidth of Δ ≈ 3.3 nm and the fourth wavelength channel = 1554 nm has its bandwidth of Δ ≈ 3.308 nm. Each wavelength channel has its own value and they have a very small variation between them. The working interval is = 24 nm : = 1556 nm and = 1532 nm.   Figure 4 shows the positive components for our frequency spectrum R T (ν). Analyzing Figure 4, the quasi-distributed fiber optic sensor has twelve Fabry-Pérot sensors and each local sensor produces one triangle function (frequency component, peaks). The red color corresponds to the frequency channels, ν FP11 − ν FP13 . Three peaks have the same bandwidth ν BG1 = 1.23 and their central frequencies are ν FP11 = 9.90, ν FP12 = 19.80, and ν FP13 = 39.60. The blue color corresponds to the frequency channels ν FP21 − ν FP23 . Again, three peaks have the same bandwidth ν BG2 = 1.22 and their central frequencies are ν FP21 = 6.14, ν FP22 = 14.73 and ν FP23 = 29.47. The black color corresponds to the channels ν FP31 − ν FP33 . Now, the bandwidth is ν BG3 = 1.218 and their central frequencies are ν FP31 = 8.52, ν FP32 = 20.71 and ν FP33 = 32.90. Finally, the cyan color corresponds to the channels ν FP41 − ν FP43 . The bandwidth is ν BG4 = 1.209 and the central frequencies are ν FP41 = 14.50, ν FP42 = 26.60 and ν FP43 = 45.94, respectively. The bandwidths ν BG1 . . . ν BG4 have small variations because each wavelength channel has its own Bragg wavelength, see Equation (12). Each interference pattern has its own frequency because each Fabry-Pérot sensor has its own cavity length and wavelength channel (see Equation (6)), see Table 1 and Figure 4.

Frequency Spectrum
Here, each frequency channel contains information from a specific Fabry-Pérot sensor and then the demodulation algorithm described in Section 5 can be applied.  Figure 4 shows the positive components for our frequency spectrum . Analyzing Figure 4, the quasi-distributed fiber optic sensor has twelve Fabry-Pérot sensors and each local sensor produces one triangle function (frequency component, peaks  (12). Each interference pattern has its own frequency because each Fabry-Pérot sensor has its own cavity length and wavelength channel (see Equation 6), see Table 1 and Figure 4.

Frequency Spectrum
Here, each frequency channel contains information from a specific Fabry-Pérot sensor and then the demodulation algorithm described in Section 5 can be applied.

Numerical Results
Applying the demodulation algorithm described in Section 5 and using the parameters presented in Section 6.1, the quasi-distributed fiber optic sensor (Figure 1) was numerically simulated. Our numerical results are shown in Figure 5. We present the behavior of demodulation errors vs. signal-to-noise rate . If demodulation errors are called the "resolution" as in our previous work [16,17], then all low-precision Fabry-Pérot sensors have two resolutions. Both resolutions are possible because the signal demodulation ( Figure 2) is based on the FDPA algorithm and this algorithm evaluates the Bragg wavelength shift twice. Observing Figure 5, each Fabry-Pérot sensor has its own high resolution; however each wavelength channel produces one low resolution because this resolution does not depend on the cavity length. Additionally, the transition from high resolution to low resolution was reported [16,17,19], again, our numerical results were presented. These results corroborate our signal analysis and our demodulation algorithm. Wavelength channel BG2 =1542nm Wavelength channel BG1 =1536nm Wavelength channel BG4 =1554nm

Numerical Results
Applying the demodulation algorithm described in Section 5 and using the parameters presented in Section 6.1, the quasi-distributed fiber optic sensor (Figure 1) was numerically simulated. Our numerical results are shown in Figure 5. We present the behavior of demodulation errors vs. signal-to-noise rate SNR 1 2 . If demodulation errors are called the "resolution" as in our previous work [16,17], then all low-precision Fabry-Pérot sensors have two resolutions. Both resolutions are possible because the signal demodulation (Figure 2) is based on the FDPA algorithm and this algorithm evaluates the Bragg wavelength shift twice. Observing Figure 5, each Fabry-Pérot sensor has its own high resolution; however each wavelength channel produces one low resolution because this resolution does not depend on the cavity length. Additionally, the transition from high resolution to low resolution was reported [16,17,19], again, our numerical results were presented. These results corroborate our signal analysis and our demodulation algorithm. For example: if the OSA spectrometer has ∆ = 10 pm (a typical value), the broadband source has = 100 nm and as ∆ = − = 1554 nm − 1548 nm = 6 nm (see Figure 3), the quasi-distributed sensor will have their limits as Table 3 illustrates.  (18) From Tables 1-3, the simulated quasi-distributed sensor satisfies the instrumentation and the signal requirements. Observing Tables 1, 3 and Figures 3 For example: if the OSA spectrometer has ∆λ = 10 pm (a typical value), the broadband source has λ w = 100 nm and as ∆λ op = λ BG4 − λ BG3 = 1554 nm − 1548 nm = 6 nm (see Figure 3), the quasi-distributed sensor will have their limits as Table 3 illustrates. Table 3. Quasi-distributed sensor limits ∆λ = 10 pm and λ w = 100 nm .

Parameter
Value Equation From Tables 1-3, the simulated quasi-distributed sensor satisfies the instrumentation and the signal requirements. Observing Tables 1 and 3 and Figures 3-5, the numerical results are in concordance with the theory. Thus, we confirm our theoretical analysis and our new demodulation algorithm. The numerical results are shown in Figure 5. The theoretical analysis and our numerical results are in concordance with experimental results reported by Shlyagin et al. [12] and Della-Tamin et al. [13]. Then, WDM/FDM techniques can be implemented on low-precision Fabry-Pérot sensors and our new algorithm demodulates its optical signal. The presented study optimizes significantly the quasi-distributed sensor implementation, its design and the sensibility of all local sensors.

Discussion
Based on our signal analysis and our numerical results, the quasi-distributed sensor would be built through WDM/FDM techniques and low-finesse Fabry-Pérot interferometers. The theoretical analysis optimizes the quasi-distributed sensor presented in Figure 1 which was experimentally proposed in [12,13]. The optimization considers the multiplexing technique, the optical instrumentation, the detection technique, the local sensor properties, the noise system, the spatial resolution and the signal demodulation. Observing Figure 5, the quasi-distributed sensor has good functionality because each Fabry-Pérot sensor is performing very well and our new algorithm can demodulate the optical signal. The Fabry-Pérot sensors have high resolution when the signal-to-noise rate (SNR) is big, but the same sensors have low resolution if and only if the signal-to-noise rate is low (the optical system has many noise). The threshold value between both resolutions (high and low resolutions) can be determined by: where σ envkm is the low resolution for the km-th Fabry-Pérot sensor. This value depends of the Bragg wavelength and the cavity length. Applying the parameters presented in Table 1, the threshold values were calculated as Table 4 shows. From Table 4, each Fabry-Pérot sensor has its own threshold value. Combining the multiplexing techniques (WDM and FDM), the number of local sensors was increased for the optical sensing system presented in Figure 1. For example: if the frequency division multiplexing has forty frequency channels (M = 40) [16] and if the wavelength division multiplexing has 100 wavelength channels (K = 100) [21,22], then our sensing system would have MXK = 40 × 100 = 4000 local sensors. This confirms that the quasi-distributed sensor has high capacity to measure some physical parameters along the optic fiber. Additionally, applying the conditions 1-6 described in Section 4, ghost interferometers, cross-talk noise and overlapping between two frequency components are eliminated, proving that the sensing system will have optimal functionality. Figure 1 showed a serial array of interferometer sensors and the topology permits the measurements along of a cable. For the measurement over any surface, a serial/parallel topology can be applied. In this case, the optical system consists of an optical broadband source, an OSA spectrometer, a 2XK-splitter and a quasi-distributed sensor. The sensor would have K-fiber optics and each fiber optic can have M-Fabry-Perot interferometers. The proposed optical sensor can be observed in Figure 6. Their practical application can be on: civil engineering, mechanical engineering, military applications, civil protection and disaster risk reduction. The quasi-distributed fiber sensor presented in this work applied the direct spectroscopic detection (DSD) technique. On the one hand, this detection method has a simple configuration and it just basically required an optical broadband source, an optical circulator, an optical spectrometer analyzer (OSA spectrometer) and a PC computer, see Figure 1. On the other hand, the authors described that the Optical Frequency Domain Reflectometry (OFDR) fell into two main categories [23]: incoherent (or direct detection) OFDR (I-OFDR) and the coherent OFDR (C-OFDR). From the first category, there are four derivations known as Network Analysis OFDR (NA-OFDR) [24], Incoherent Frequency-Modulate Continuous Wave (I-FMCW) [25], Step-frequency method [26] and Sweep Frequency method [27]. Each method was previously applied on other sensing systems. Comparing the sensing system presented in Figure 1 with the sensing systems reported by the authors of reference [23], the direct spectroscopic detection technique has some benefits, for example: simple configuration, low cost, less complexity and fewer instruments used, therefore, our sensing system has some advantages over other detection techniques.
Our future work is in the following direction: the transition between both resolutions could be determined using the physical parameters, instrumentation, the FDPA algorithm and noise. The quasi-distributed sensor application would be in this other direction.

Conclusions
In this work, a quasi-distributed fiber sensor was theoretically analyzed and also was numerically simulated. The quasi-distributed sensor was based on wavelength/frequency division multiplexing and low-finesse Fabry-Pérot interferometers. Theory and simulation were in concordance. During the analysis and the simulation, we considered the signal processing, multiplexing techniques, the optical instrumentation, the system noise, the detection technique and the local sensor parameters. The signal analysis and signal demodulation algorithm optimize the sensor implementation while the numerical simulation demonstrated its excellent functionality. From our numerical results, we confirmed that each Fabry-Pérot sensor has two resolutions since the Fourier Domain Phase Analysis algorithm makes two evaluations of the Bragg wavelength shift as previously reported by us. This quasi-distributed sensor finds many potential industrial applications due to its functionality, low cost by sensing point, high resolution and high sensitivity. Funding: The authors thank the Mexico's National Council of Science and Technology (CONACyT) and University of Guadalajara for the support granted. This investigation was carried out following the line of The quasi-distributed fiber sensor presented in this work applied the direct spectroscopic detection (DSD) technique. On the one hand, this detection method has a simple configuration and it just basically required an optical broadband source, an optical circulator, an optical spectrometer analyzer (OSA spectrometer) and a PC computer, see Figure 1. On the other hand, the authors described that the Optical Frequency Domain Reflectometry (OFDR) fell into two main categories [23]: incoherent (or direct detection) OFDR (I-OFDR) and the coherent OFDR (C-OFDR). From the first category, there are four derivations known as Network Analysis OFDR (NA-OFDR) [24], Incoherent Frequency-Modulate Continuous Wave (I-FMCW) [25], Step-frequency method [26] and Sweep Frequency method [27]. Each method was previously applied on other sensing systems. Comparing the sensing system presented in Figure 1 with the sensing systems reported by the authors of reference [23], the direct spectroscopic detection technique has some benefits, for example: simple configuration, low cost, less complexity and fewer instruments used, therefore, our sensing system has some advantages over other detection techniques.
Our future work is in the following direction: the transition between both resolutions could be determined using the physical parameters, instrumentation, the FDPA algorithm and noise. The quasi-distributed sensor application would be in this other direction.

Conclusions
In this work, a quasi-distributed fiber sensor was theoretically analyzed and also was numerically simulated. The quasi-distributed sensor was based on wavelength/frequency division multiplexing and low-finesse Fabry-Pérot interferometers. Theory and simulation were in concordance. During the analysis and the simulation, we considered the signal processing, multiplexing techniques, the optical instrumentation, the system noise, the detection technique and the local sensor parameters. The signal analysis and signal demodulation algorithm optimize the sensor implementation while the numerical simulation demonstrated its excellent functionality. From our numerical results, we confirmed that each Fabry-Pérot sensor has two resolutions since the Fourier Domain Phase Analysis algorithm makes two evaluations of the Bragg wavelength shift as previously reported by us. This quasi-distributed sensor finds many potential industrial applications due to its functionality, low cost by sensing point, high resolution and high sensitivity. Funding: The authors thank the Mexico's National Council of Science and Technology (CONACyT) and University of Guadalajara for the support granted. This investigation was carried out following the line of research "Nanostructured Semiconductor Oxides" of the academic group UDG-CA-895 "Nanostructured Semiconductors" of CUCEI, University of Guadalajara.

Conflicts of Interest:
The authors declare no conflict of interest.