The present invention relates to an optical fiber-based distributed sensor and, more particularly, to a multicore fiber-based sensor that is able to detect the presence of microbends along the extent of a given fiber, providing three-dimensional information on the location and size of various deformations within a space surrounding the distributed sensor.
Optical fiber-based distributed sensors have emerged as an invaluable tool in performing characterizations of arbitrary deformations in three-dimensional space. Potential applications include 3D printing, surgical catheters, smart wearables, monitoring systems for fuel tanks, composite structures, and the like. The use of optical fibers for “shape sensing” offers high precision and high-speed operation, and may be particularly applicable for characterizing difficult-to-access surfaces and environments, as a result of the built-in shielding of the light beam that is used as the sensing probe.
To date, fiber-based distributed sensors have been able to reconstruct arbitrary paths and shapes only at the “macro” level (i.e., on a scale of centimeter/meter in terms of measurement). Going forward, the ability to perform distributed sensing at a smaller scale (i.e., sub-millimeter changes/bends) will become more important. For example, the influence of microbends on the attenuation of an optical communication signal propagating along a transmission fiber has been of interest for decades. As the transmission loss in optical fibers approaches the fundamental limits dictated by the intrinsic absorption and scattering in glass, the losses induced by the microscopic physical bends in the optical fibers (and cables) are becoming increasingly relevant. However, such microbends cannot be measured directly with the currently-available sensors.
The needs remaining in the prior art are addressed by the present invention, which relates to a multicore fiber-based sensor that is able to detect the presence of microbends along the extent of a given fiber, providing three-dimensional information on the location and size of various deformations within a space surrounding the distributed sensor.
In accordance with the principles of the present invention, the ability to “reconstruct” micro-deformations that are distributed along the length of an optical fiber is provided by a system that is based on the use of a twisted multicore optical fiber to probe the distributed reflection of light within the multiple waveguiding cores. The cores are formed to include continuous fiber Bragg gratings (FBGs) that all exhibit the same Bragg wavelength. A micro-scale local deformation of the sensing fiber produces a local shift in the Bragg wavelength, where the use of multiple cores allows for a complete modeling of a bend at a specific location.
In one exemplary embodiment, the present invention takes the form of a distributed system for sensing and measuring microbends and micro-deformations in a three-dimensional (3D) space that utilizes a multicore sensing fiber in combination with an optical backscatter reflectometer. In particular, the multicore sensing fiber is formed to including a plurality of offset cores that are radially spaced from a center of the multicore sensing fiber by an amount Ro and a plurality of continuous fiber Bragg gratings (FBGs) inscribed in the plurality of offset cores in a one-to-one relationship. The set of FBGs are formed to reflect light at a common Bragg wavelength λBragg. The optical backscatter reflectometer includes a tunable laser source for generating a swept wavelength output beam spanning a wavelength range surrounding λBragg, an optical beam splitter/combiner, an optical detector, and a Fourier transform analyzer for performing optical frequency domain reflectometry (OFDR). The optical beam splitter/combiner functions to split the swept wavelength output beam from the tunable laser source into a swept wavelength “probe” beam that is directed into the multicore sensing fiber and a swept wavelength reference beam directed into a reflector. The optical beam splitter/combiner is also used for combining a swept wavelength return beam from the multicore sensing fiber and a reflected swept wavelength reference beam to create an interfering FBG sensing beam. The optical detector is responsive to the interfering FBG sensing beam for creating an electronic version of the interference beam, with the Fourier transform analyzer thereafter used to perform a Fourier transform on the electronic version of the interfering FBG sensing beam to generate measurements of local changes in Bragg wavelength along the length of the multicore sensing fiber and reconstruct therefrom the shape of the three-dimensional space.
While the sensor fiber may be formed of a conventional glass material, other embodiments may utilize a sensor fiber formed of a material that is less elastic, with a smaller Young's modulus that allows for an even finer degree of measurement resolution.
Other and further embodiments and features of the present invention will become apparent during the course of the following discussion and by reference to the accompanying drawings.
Referring now to the drawings, where like numerals represent like parts in several views:
As best shown in
Each offset core 12i is formed to include a continuous FBG 14i, which may be written in the cores during the process of drawing down the preform into the final fiber. Each FBG 14i is created to exhibit the same Bragg wavelength λBragg so that they all reflect light of the same wavelength in the absence of any local bends or deformations that otherwise creates a shift in the Bragg wavelength value.
OBR 20 itself includes a tunable laser source 22 that is configured as a swept wavelength (frequency) source, centered on the Bragg wavelength (λBragg) of FBGs 14. In one exemplary embodiment, tunable laser source 22 may be configured to provide a narrow linewidth output that is swept through a wavelength range±10 nm on either side of λBragg (i.e., a wavelength range of 20 nm). For example, if λBragg=1541 nm, tunable laser source may be configured to provide an output beam that is scanned across the wavelength range of 1531 nm to 1551 nm. The tunable bandwidth of 20 nm is exemplary only, and there are instances where a larger bandwidth may be desirable, as discussed below.
The output beam from tunable laser 22 thereafter passes through a beam splitter 24 of OBR 20, which directs a majority of the beam (referred to at times as a “probe beam” or “probe signal”) out of OBR 20 and into a 1×N optical switch 30. Switch 30 is controlled to direct the probe beam into a selected offset core 12i of multicore sensing fiber 10 in a manner that will be described in detail below.
Returning to the description of OBR 20, the remaining output from beam splitter 24 (referred to at times as a “reference beam”) is directed along a reflective signal path 26. The reflected reference beam and backscattered reflections from multicore sensing fiber 10 are combined within beam splitter 24 (operating as a combiner in this direction), directing the interfering combination of these signals into an optical detector 28 included within OBR 20. The output from optical detector 28 is thereafter applied as an input to Fourier analyzer 29, which performs frequency domain analysis, converting the frequency domain measurement from optical detector 28 into a space-domain measurement of phase and amplitude as a function of length along multicore sensing fiber 10.
If multicore sensing fiber 10 is flat and straight with no microbends (or other types of micro-deformations), FBGs 14 will all maintain the reference Bragg wavelength λBragg and will consistently reflect probe beam light at only this wavelength, allowing the remaining wavelengths to continue to propagating along multicore sensing fiber 10. Therefore, inasmuch as there is no change in λBragg, there is also no change in the frequency component of the output from optical detector 28. Thus, Fourier analyzer 29 provides a constant, linear output signal indicative of an “unperturbed” multicore sensing fiber 10. Once any microbend/deformation is present within fiber 10, the Bragg wavelength of one or more offset cores 12 will change (see
Therefore, in accordance with the principles of OFDR, by illuminating each offset core 12i with a probe beam that is scanned across the defined wavelength range, deformations/microbends along multicore sensing fiber 10 will be identified within the interference signal processed by Fourier analyzer 29. That is, by performing a Fourier transformation of the interfering beams, the spectral information can be used to detect and measure micro-deformations along multicore sensing fiber 10. The Fourier transform converts the spectral information in the received interference signal into spatial (temporal) information, depicted in the form of distributed Bragg wavelength changes at locations where microbends/deformations are present.
The Fourier relation inversely relates the wavelength scanning range of tunable laser source 22 to the longitudinal spatial domain measurements of strain (and, therefore, of curvature and shape). For example, a 20 nm wavelength scanning range translates into a 40 μm measurement resolution. Increasing the wavelength scanning range to 80 nm (still centered on the defined Bragg wavelength) translates into a 10 μm resolution in the measurement of local microbends, albeit at the cost of requiring a tunable laser source 22 that is capable of generating such a large swept wavelength range.
Continuing with the description of the components of system 100, the tunable probe beam exiting OBR 20 is provided as an input to 1×N optical switch 30, as mentioned above. Optical switch 30 includes a single input/output port 32 and a plurality of N connecting ports 341-34N, with each connecting port 34i associated with a unique offset core 12i. The plurality of N outputs from optical switch 30 are coupled into a plurality of separate optical fibers 381-38N, which are associated with offset cores 121-12N in a one-to-one relationship. The far end of multicore sensing fiber 10 is immersed in an index-matching gel 50 to suppress unwanted Fresnel reflections at the far endface of fiber 10 from re-entering one or more of the multiple offset cores 12.
System 100 as shown in
By using optical switch 30 to sequentially illuminate each individual offset core 12i with the swept wavelength probe beam, the changes in Bragg wavelength associated with specific cores at a given transverse location thus allows for the type of shape deformation to be re-created. That is, the inclusion of the switching capability within system 100 allows for the collection of data from multiple offset cores 12 on a one-by-one basis so as to obtain cross-sectional deformations at selected locations along fiber 10. Repeating this process along the extent of multicore sensing fiber 10 allows for a complete reconstruction of various microbends (and other types of deformations) that occur along its span.
Once all of the measurements have been completed by Fourier analyzer 29, the distributed curvature and shape of the associated three-dimensional space may be created by a reconstruction module 27 that is coupled to the output of Fourier analyzer 29, as shown in
where Ro is the radial offset between the center of fiber 10 and the center of an offset core 12u, u defines the individual cores, εu(z) is the unit vector of the respective core 12u, and εu(z) is the strain induced in the corresponding core 12u. Using the strain-optic coefficient η of silica glass (˜0.78) and measured local changes in Bragg wavelength ΔλBragg recorded by Fourier analyzer 29, the corresponding local strain εu(z) experienced by core u can be defined by reconstruction module 29 as:
By utilizing optical switch 30 to sequentially illuminate each offset core 121-12N, the strain information from the plurality of N (for example, N=6) offset cores 12 is summed within reconstruction module 29 in the manner shown in the definition of spatially-dependent curvature κ(z), developing both the magnitude and phase of the distributed fiber curvature. The bend orientation of each curve along multicore sensing fiber 10 is represented by the phase portion of the curvature. It is to be understood that the number of individual offset cores 12 included within multicore sensing fiber 10 directly impacts the accuracy of the calculated distributed curvature, where increasing the number of offset cores will increase the amount of data that is captured and recorded by Fourier analyzer 29.
Finally, the distributed shape S of the deformed fiber may also be provided as an output from reconstruction module 27. In particular, the distributed shape is reconstructed from the calculated spatially-dependent curvature κ(z), using the Frenet-Serret formulas (which are a set of differential equations describing a three-dimensional (3D) curve) to provide the distributed shape output from module 27. Specifically, the Frenet-Serret equations relate the local shape parameters, including the tangent T(x,y,z), normal N(x,y,z) and binomial B(x,y,z) vectors, with the fiber curvature and torsion measured at the closely spaced locations. Mathematically, this is expressed as:
where, S≡[T(x,y,z); N(x,y,z); B(x,y,z)], {dot over (S)}=dS/dz, and the torsion τ(z) quantifies how rapidly the bend direction changes along the length of the curved fiber. In practice, a spatial derivative of the phase component of the distributed curvature vector leads to the amount of torsion τ(z) (=dθb(z)/dz) produced along the length of multicore sensing fiber 10. By repeatedly solving this expression for the eigenvalues and eigenvectors of the set S along the length (z-axis) of multicore sensing fiber 10, the distributed shape of the fiber may be estimated.
It is important to note that the initial conditions for solving the above expression assume an absence of curvature and torsion at the location z=0 i.e., κ(0)=τ(0)=0 at the input to multicore sensing fiber 10. Furthermore, the tangent T(x,y,z), normal N(x,y,z) and binormal B(x,y,z) vectors, at z=0, are defined as the three orthonormal unit vectors in an arbitrarily chosen three-dimensional spatial frame-of-reference. The tangent vector at any position is assumed to be “pointing” in the direction of the increasing fiber length and indicates the local fiber direction. Therefore, a concatenation of the tangent vectors at closely spaced locations along the length of multicore sensing fiber 10 represents the distributed shape of the fiber.
The measurement sensitivity of the inventive system may be increased by increasing the signal-to-noise ratio (SNR) of OBR 20, or broadening the tuning wavelength range of tunable laser 22 to increase measurement resolution (as mentioned above). The SNR depends on the spectral beating signal generated by interfering the reference beam ({right arrow over (E)}r) with the backscattered signal ({right arrow over (E)}s) in OBR 20. That is, SNR∝|{right arrow over (E)}r|×|{right arrow over (E)}s|. Therefore, the SNR can be increased two-fold (for example) by increasing the intensity of tunable laser source 22, or by simply increasing by two-fold the amplitude of refractive index modulation Δnac of Bragg gratings 14, since|{right arrow over (E)}s|∝nac. Reducing background noise present within the instrumentation of OBR 20 itself (e.g., shot-noise, dark current noise, frequency measurement noise, and the like) also increases the SNR of OBR 20 and, as a result, the measurement sensitivity of the system.
Increasing the sensitivity of measurements in the transverse plane of multicore sensing fiber 10 may also be provided by increasing the radial offset Ro between cores 12 and the central axis of fiber 10, while maintaining the same outer diameter of the fiber. The amount of Bragg wavelength shift (ΔλBragg) in the presence of bend-induced fiber strain is directly related to the value of Ro, as shown by the following relation:
where η is a fixed quantity representing the strain-optic coefficient of the silica glass, y0 is the amount of fiber displacement in the transverse plane with respect to a straight (flat) neutral plane, and kd is the period of the deformation imposed along the length of the fiber. Clearly, the amount of detected wavelength shift can be increased by proportionately increasing the radial offset (i.e., Ro) of cores 12. This leads to a linear increase in the SNR of the system, which ultimately improves the sensitivity of the measurements.
Another alternative approach to increasing measurement sensitivity is to reduce the overall diameter of multicore sensing fiber 10. Since the fiber is cylindrical in form, reducing the diameter serves to lower the moment of inertia I (I=π/4*R4), where R is the radius of fiber 10. It follows that by lowering the moment of inertia, the flexibility (and thus bending) of multicore sensing fiber 10 itself is increased, providing a larger shift in Bragg wavelength in FBGs 14. The increase in I may improve the sensitivity of the local strain, the local curvature and, ultimately, the distributed shape measurements. Specifically, by reducing the fiber diameter by 50%, the value of I50% is reduced to about 0.0625I and a factor of sixteen increase in both the resulting bend amplitude y0 and associated Bragg wavelength shift ΔλBragg.
Fabricating multicore sensing fiber 10 from an optical material with a Young's modulus (E) less than that of conventional silica glass (E=˜70 GPa) also leads to improvements in SNR. Soft glasses, such as chalcogenide and fluoride glasses present suitable platforms for the reduced Young's modulus of multicore sensing fiber 10. On the other hand, the longitudinal sensitivity of the shape-sensing measurements can be increased proportionally by increasing the precision of the estimated group delay for the distributed backscatter signal. Using a set of distributed measurements of the refractive index of fiber 10 may be used to determine the estimated group delay.
It has also been found that performing repeated measurements of each offset core allows for the noise present in the averaged values to be reduced. For example, switch 30 may be controlled to perform multiple switchings from port 341 through port 34N, forming multiple measurement scans of multicore sensing fiber 10. That is, by performing multiple scans of each offset 12, the noise contribution associated with a single scan is reduced by averaging out over multiple scans. That is, the repeated measurements result in suppressing noise present in the measurement data by averaging the data over a number of scans, thereby effectively enhancing the SNR and improving the accuracy of the fiber shape measurement.
It will be apparent to those skilled in the art that various modifications and variations can be made to the present invention without departing from the spirit or scope thereof. Thus, it is intended that the present invention cover the modifications and variations of the above-described embodiments, all of which are considered to fall within the spirit and scope of the invention as the defined by the claims appended hereto.
This application claims the benefit of U.S. Provisional Application No. 62/989,117, filed Mar. 13, 2020 and hereby incorporated by reference.
Filing Document | Filing Date | Country | Kind |
---|---|---|---|
PCT/US21/22051 | 3/15/2021 | WO |
Number | Date | Country | |
---|---|---|---|
62989117 | Mar 2020 | US |