光纤能测冰厚吗？ | Can optical fibers measure ice thickness?

OFSCN Global Technical Center OFSCN Technical Knowledge Base
1

How does a fiber optic cable know how much ice has accumulated on it while hanging from a wire?

2

In electrical engineering, optical fibers do not directly “sense” or “feel” the physical thickness of ice accumulation. Instead, they indirectly calculate it through mechanical sensing and mechanical model inversion (catenary equation).

The core physics and mathematical logic for measuring ice thickness on overhead transmission lines (including OPGW fiber optic composite overhead ground wires containing optical fibers, or externally attached fiber optic sensors) are as follows:

I. Fundamental Measurement Physics Principles

1. Ice Load Conversion to Axial Tension and Strain

The transmission line or ground wire is suspended between two towers, and its natural shape conforms to the catenary (Catenary) equation in mechanics.

  • No-ice state: The conductor bears its own gravitational load and pre-tension.
  • Icing state: Ice layers adhere to the conductor, significantly increasing the total line load per unit length (q):
    q_\text{total} = q_\text{line} + q_\text{ice} + q_\text{wind}
  • As the weight of ice accumulation (q_\text{ice}) increases, the conductor’s sag increases, and the horizontal tension at the suspension points and the axial tensile stress/strain (\varepsilon) within the conductor increase synchronously.

2. Optical Fiber’s Perception of Mechanical Strain

Fiber Bragg Gratings (FBGs) or distributed sensing media (like strain optical fibers based on BOTDA/BOTDR/OFDR technology) inside optical fibers are extremely sensitive to minute axial strains.
When the transmission line is stretched by ice load, the optical fiber experiences mechanical strain, causing a shift in the center wavelength of its grating reflection:

\frac{\Delta \lambda_B}{\lambda_B} = (1 - p_e) \varepsilon + (\alpha_\Lambda + \xi) \Delta T

Where:

  • p_e is the effective photoelastic coefficient;
  • \varepsilon is the axial strain;
  • \alpha_\Lambda and \xi are the coefficients of thermal expansion and thermo-optic effect, respectively;
  • \Delta T is the change in ambient temperature.

3. Strict Temperature Compensation and Decoupling

In extremely cold winter environments, conductor contraction due to cold can also lead to increased mechanical tension and wavelength drift. To distinguish between strain caused by “cold contraction” and strain caused by “ice weight,” the system typically installs parallel fiber optic temperature measurement points in locations not subjected to external pulling forces. This allows for precise subtraction of wavelength drift \Delta \lambda_T caused by temperature changes, extracting the mechanical strain \varepsilon_\text{mechanical} purely from ice weight.

II. Mathematical Model for Inverting “Ice Thickness” from “Tension/Strain”

After the demodulator obtains strain and line tension data, the backend algorithm derives the equivalent ice thickness through the following three steps:

  1. Calculate Tension Increment: Based on the conductor’s cross-sectional area and elastic modulus (\sigma = E \cdot \varepsilon), calculate the additional tension \Delta T on the entire span due to ice accumulation.
  2. Invert Ice Weight per Unit Length (W_{\text{ice}}): Using the span length (L), height difference, and the initial geometric state of the conductor, the catenary mechanical equilibrium equation is used to inversely calculate the mass of ice added per meter of line m_\text{ice}\ (\text{kg/m}).
  3. Calculate Equivalent Cylindrical Ice Thickness (d_\text{ice}):
    Power engineering standards (such as Chinese power standards and IEC specifications) typically simplify irregular ice layers into a model of a uniform cylindrical sheath covering the outside of the conductor:
    m_\text{ice} = \rho_\text{ice} \cdot \pi \cdot \left[ \left( \frac{D}{2} + d_\text{ice} \right)^2 - \left( \frac{D}{2} \right)^2 \right]
    After simplification, this becomes:
    m_\text{ice} = \rho_\text{ice} \cdot \pi \cdot d_\text{ice} (D + d_\text{ice})
    Where:
    • D is the outer diameter of the bare conductor;
    • \rho_\text{ice} is the density of the ice layer (typically taken as 0.9\ \text{g/cm}^3 for standard rime ice in industrial calculations);
    • d_\text{ice} is the equivalent ice thickness.

By solving this quadratic equation, the system can output the equivalent ice thickness (mm) on the conductor in real-time.

III. Related Fiber Optic Sensing Components

In actual power online monitoring projects, the primary method relies on fiber optic tension/stress sensors installed at both ends of tension insulator strings, coupled with temperature sensors for real-time strain acquisition: