WO2021246497A1 - 形状測定システム及び形状測定方法 - Google Patents

形状測定システム及び形状測定方法 Download PDF

Info

Publication number
WO2021246497A1
WO2021246497A1 PCT/JP2021/021283 JP2021021283W WO2021246497A1 WO 2021246497 A1 WO2021246497 A1 WO 2021246497A1 JP 2021021283 W JP2021021283 W JP 2021021283W WO 2021246497 A1 WO2021246497 A1 WO 2021246497A1
Authority
WO
WIPO (PCT)
Prior art keywords
core
optical fiber
core optical
outer peripheral
strain
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Ceased
Application number
PCT/JP2021/021283
Other languages
English (en)
French (fr)
Japanese (ja)
Inventor
信智 半澤
和秀 中島
隆 松井
英晶 村山
良太 和田
真輝人 小林
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
University of Tokyo NUC
NTT Inc
Original Assignee
Nippon Telegraph and Telephone Corp
University of Tokyo NUC
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Nippon Telegraph and Telephone Corp, University of Tokyo NUC filed Critical Nippon Telegraph and Telephone Corp
Priority to JP2022528897A priority Critical patent/JP7376052B2/ja
Priority to EP21818720.1A priority patent/EP4163586B1/en
Priority to US17/915,924 priority patent/US12203847B2/en
Publication of WO2021246497A1 publication Critical patent/WO2021246497A1/ja
Anticipated expiration legal-status Critical
Ceased legal-status Critical Current

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N21/47—Scattering, i.e. diffuse reflection
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B11/00—Measuring arrangements characterised by the use of optical techniques
    • G01B11/16—Measuring arrangements characterised by the use of optical techniques for measuring the deformation in a solid, e.g. optical strain gauge
    • G01B11/18—Measuring arrangements characterised by the use of optical techniques for measuring the deformation in a solid, e.g. optical strain gauge using photoelastic elements
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01L—MEASURING FORCE, STRESS, TORQUE, WORK, MECHANICAL POWER, MECHANICAL EFFICIENCY, OR FLUID PRESSURE
    • G01L1/00—Measuring force or stress, in general
    • G01L1/24—Measuring force or stress, in general by measuring variations of optical properties of material when it is stressed, e.g. by photoelastic stress analysis using infrared, visible light, ultraviolet
    • G01L1/242—Measuring force or stress, in general by measuring variations of optical properties of material when it is stressed, e.g. by photoelastic stress analysis using infrared, visible light, ultraviolet the material being an optical fibre
    • G01L1/246—Measuring force or stress, in general by measuring variations of optical properties of material when it is stressed, e.g. by photoelastic stress analysis using infrared, visible light, ultraviolet the material being an optical fibre using integrated gratings, e.g. Bragg gratings
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/62—Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light
    • G01N21/63—Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light optically excited
    • G01N21/636—Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light optically excited using an arrangement of pump beam and probe beam; using the measurement of optical non-linear properties
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/84—Systems specially adapted for particular applications
    • G01N21/88—Investigating the presence of flaws or contamination
    • G01N21/95—Investigating the presence of flaws or contamination characterised by the material or shape of the object to be examined
    • G01N21/954—Inspecting the inner surface of hollow bodies, e.g. bores
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N21/47—Scattering, i.e. diffuse reflection
    • G01N2021/4704—Angular selective
    • G01N2021/4709—Backscatter
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/62—Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light
    • G01N21/63—Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light optically excited
    • G01N21/636—Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light optically excited using an arrangement of pump beam and probe beam; using the measurement of optical non-linear properties
    • G01N2021/638—Brillouin effect, e.g. stimulated Brillouin effect
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/84—Systems specially adapted for particular applications
    • G01N21/88—Investigating the presence of flaws or contamination
    • G01N21/95—Investigating the presence of flaws or contamination characterised by the material or shape of the object to be examined
    • G01N21/954—Inspecting the inner surface of hollow bodies, e.g. bores
    • G01N2021/9542—Inspecting the inner surface of hollow bodies, e.g. bores using a probe
    • G01N2021/9546—Inspecting the inner surface of hollow bodies, e.g. bores using a probe with remote light transmitting, e.g. optical fibres

Definitions

  • the present disclosure relates to an apparatus and a method for deriving a three-dimensional shape of an object to be measured, such as a pipeline or a submarine cable, by a reflectometry technique using an optical fiber.
  • Non-Patent Document 1 A technique has been proposed in which the reflection spectrum in the frequency domain of each core of a multi-core optical fiber is measured and analyzed by OFDR (Optical Frequency Domain Reflectometry) to derive a three-dimensional shape to be measured (for example, Non-Patent Document 1). reference.). Further, a method of imparting FBG (Fiber Bragg Grating) to the entire length of the sensing medium (optical fiber) to improve the measurement resolution has also been proposed (see, for example, Non-Patent Document 2).
  • FBG Fiber Bragg Grating
  • OFDR When trying to derive the three-dimensional shape of the object to be measured by OFDR, the following difficulties occur. (1) OFDR can realize high resolution on the order of several tens of mm, but the measurement distance is limited to about several tens of meters, and it is difficult to derive a long-distance three-dimensional shape. (2) If an FBG is added to an optical fiber, it becomes difficult to improve the ease of manufacturing and the economic efficiency of the optical fiber. (3) The technique of performing shape identification by quasi-distribution measurement such as addition of FBG has a limitation on the number of measurement points, and it is difficult to identify the shape over a long distance.
  • an object of the present invention is to provide a shape measurement system and a shape measurement method capable of deriving a three-dimensional shape of a linear object to be measured over a long distance and with high resolution in order to solve the above problems. do.
  • the shape measurement system uses a multi-core optical fiber having a predetermined core arrangement as a sensing medium, and analyzes data acquired by BOTDR (Brillouin Optical Time Domain Reflectometer). ..
  • the shape measuring system is A central core arranged in the center of the cross section, and a multi-core optical fiber having three or more outer peripheral cores arranged concentrically and outside the central core at equal intervals.
  • a measuring device for measuring the rear Brillouin scattered light distribution in the propagation direction of each core of the multi-core optical fiber and a measuring device.
  • the shape measuring method according to the present invention is A central core arranged in the center of the cross section and a multi-core optical fiber having 3 or more and 6 or less outer peripheral cores arranged outside the central core and concentrically at equal intervals are arranged along the linear structure.
  • the rear Brilluan scattered light distribution in the propagation direction of each core of the multi-core optical fiber is measured, and the multi-core optical fiber and the three-dimensional shape arranged along the linear structure whose three-dimensional shape is unknown are known. From the rear Brilluan scattered light distribution of the multi-core optical fiber arranged along a certain linear structure, the position coordinates in the three-dimensional space of the linear structure whose three-dimensional shape is unknown are calculated, and the time change thereof is calculated. To identify, I do.
  • the present invention can provide a shape measurement system and a shape measurement method capable of deriving a three-dimensional shape of a linear object to be measured over a long distance and with high resolution.
  • the specific analysis method is as follows. Let z be the position in the longitudinal direction of the multi-core optical fiber.
  • the analyzer is The amount of strain at the position z of each core of the multi-core optical fiber obtained from the rear Brillouin scattered light distribution when the multi-core optical fiber is arranged along a linear structure whose three-dimensional shape is unknown, and the above-mentioned The difference from the strain amount of the position z of each core of the multi-core optical fiber obtained from the rear Brillouin scattered light distribution when the multi-core optical fiber is arranged along a linear structure whose three-dimensional shape is known.
  • the bending angle ⁇ represents an angle formed by a vector connecting the MCF center and the center of curvature and a reference direction vector, for example, a vector connecting the MCF center and a specific core in the MCF cross section.
  • ⁇ r ⁇ ⁇ ⁇ cos ( ⁇ - ⁇ )
  • r is the distance between the centers of the central core and the outer peripheral core
  • ⁇ is an angle representing the position of the outer peripheral core in the cross section of the multi-core optical fiber.
  • a multi-core optical fiber when a multi-core optical fiber is placed along a structure, an unintended twist may occur, which may cause an error in the measurement result.
  • the multi-core optical fiber is given a known torsion, and the analyzer is known to have a torsional strain occurring in the multi-core optical fiber when calculating the position coordinates.
  • the unintended torsion that occurs when the three-dimensional shape is placed along a linear structure whose three-dimensional shape is unknown is estimated based on the strain due to the torsion, and the influence of the unintended torsion is excluded. matter, Is preferable.
  • the specific analysis method in this case is as follows. Let z be the position in the longitudinal direction of the multi-core optical fiber.
  • the analyzer is The amount of strain at the position z of each core of the multi-core optical fiber obtained from the rear Brillouin scattered light distribution when the multi-core optical fiber is arranged along a linear structure whose three-dimensional shape is unknown, and the above-mentioned The difference from the strain amount of the position z of each core of the multi-core optical fiber obtained from the rear Brillouin scattered light distribution when the multi-core optical fiber is arranged along a linear structure whose three-dimensional shape is known.
  • r is the distance between the centers of the central core and the outer peripheral core
  • k 1 is the twist correction coefficient represented by the equation (3)
  • ⁇ is the Poisson's ratio of the multi-core optical fiber
  • a i is the core i.
  • the initial angle of, ⁇ i is the angle representing the position of the outer peripheral core in the cross section at the position z of the multi-core optical fiber represented by the equation (4)
  • p is the spin rate of the outer peripheral core
  • ⁇ twisting is the multi-core light.
  • the torsional strain generated at the position z of the fiber, ⁇ (z) is the specific torsion angle at the position z of the multi-core optical fiber represented by the equation (5)
  • k 2 is the equation (6). It is a correction coefficient of twist represented by.
  • the analyzer calculates the spin rate p of the outer peripheral core from the periodic fluctuations of the bending strain ⁇ bending and i of each of the outer peripheral cores.
  • the error in the number of twists can be taken into consideration due to the manufacturing error of the multi-core optical fiber, and the shape change of the object to be measured can be specified with high accuracy.
  • the multi-core optical fiber preferably has 3 or more outer peripheral cores and a clad diameter of 375 ⁇ m or more. Further, in the multi-core optical fiber, the distance between the center core and the outer peripheral core is preferably 120 ⁇ m or more. Fan-in / fan-out (FIFO) devices can be easily realized by bundling existing single-mode fibers (clad diameter 125 ⁇ m). Further, the wider the distance between the center of the central core and the center of the outer peripheral core, the larger the distortion generated in the outer peripheral core, and the higher the dynamic range can be achieved.
  • FIFO Fan-in / fan-out
  • the present invention can provide a shape measurement system and a shape measurement method capable of deriving a three-dimensional shape of a linear object to be measured over a long distance and with high resolution.
  • FIG. 1 is a diagram illustrating a shape measuring system of the present embodiment.
  • This shape measurement system is A central core 11 arranged in the center of the cross section, and a multi-core optical fiber 10 having three or more outer peripheral cores 12 arranged concentrically and outside the central core 11 at equal intervals.
  • the analysis device 30 that calculates the position coordinates in the three-dimensional space of the linear structure whose three-dimensional shape is unknown from the above and identifies the time change thereof.
  • FIG. 1 shows an example in which a specific one core of the multi-core optical fiber 10 and a measuring device 20 for measuring the rear Brillouin scattered light distribution are connected.
  • a fan-out mechanism that separates each core of the multi-core optical fiber into a single-core optical fiber, and a form in which the multi-core optical fiber 10 and the measuring device 20 are connected via an optical switch may be used.
  • This shape measurement system uses a multi-core optical fiber 10 as a sensing medium, a measuring device (BOTDR) 20 for detecting rear Brillouin scattered light in the propagation direction of each core of the multi-core optical fiber 10, and measurement data acquired by the BOTDR 12. It is provided with an analysis device 30 for analysis.
  • BOTDR measuring device
  • the multi-core optical fiber 10 is installed along the longitudinal direction of the linear structure to be measured.
  • the linear structure is, for example, a pipeline, a riser (a pipe through which a fluid can pass from the seabed to equipment on the sea surface in offshore drilling and offshore production), a submarine cable, and the like.
  • FIG. 2 is a diagram illustrating a cross section of the multi-core optical fiber 10.
  • the multi-core optical fiber 10 has a central core in the center of the cross section and an outer peripheral core arranged at substantially equal intervals on the same circumference from the center of the cross section, for a total of four cores.
  • the four cores have substantially the same refractive index distribution and optical characteristics, but the structure may be arranged so that the refractive index distributions and optical characteristics of the cores are intentionally different. ..
  • the limits on the number of cores are as follows. If the number of cores is 3 or less (2 or less outer peripheral cores), the shape cannot be identified. On the other hand, when the number of cores is 5 or more (4 or more outer peripheral cores), the measurement accuracy is improved, but the measurement time is increased by the number of cores. Also, a circle of a certain radius can be surrounded by six circles of the same radius in contact with it. In this case, the total number of circles is seven. That is, a single-core optical fiber having the same diameter can be used as a multi-core optical fiber at the core position when seven optical fibers are closely arranged, and a FIFO can be easily produced.
  • the number of outer circles is 7 or more (the total number of circles is 8 or more)
  • the clad diameter of the multi-core optical fiber becomes large, and it becomes difficult to create a fan-in / fan-out (FIFO). Therefore, it is desirable that the number of outer peripheral cores of the multi-core optical fiber is 3 or more and 6 or less.
  • the multi-core optical fiber 10 has a larger clad outer diameter (> 125 ⁇ m) than a general communication optical fiber.
  • the distance between the center core 11 arranged at the center of the optical fiber and the outer peripheral core 12 arranged on the outer periphery of the central core 11 is set to be larger than 30 ⁇ m. The reason is that in a standard multi-core optical fiber having a clad outer diameter of 125 ⁇ m, the distance between the center core and the outer peripheral core that can be arranged while suppressing the influence of leakage loss is about 30 ⁇ m.
  • FIG. 3 is a diagram for explaining the relationship between the center-to-center distance and the strain amount of the outer peripheral core 12 with respect to the curvature of the center of the multi-core optical fiber 10 in which ( ⁇ - ⁇ ) is 0 in the relational expression (1).
  • the amount of strain obtained by a standard multi-core optical fiber having a clad outer diameter of 125 ⁇ m even when the distance between the centers of the central core 11 and the outer peripheral core 12 is 120 ⁇ m or more and the curvature is ⁇ 0.5 [1 / m]. It can be seen that a strain amount of 5 times or more can be obtained. This means that by increasing the distance between the centers of the central core 11 and the outer peripheral core 12, the sensitivity to the shape change of the small curvature generated in the multi-core optical fiber 10 can be improved.
  • This shape measurement method is Arranging the multi-core optical fiber 10 along the linear structure (step S01), Measuring the rear Brilluan scattered light distribution in the propagation direction of each core of the multi-core optical fiber 10 (step S02), and the multi-core optical fiber 10 and three-dimensionally arranged along a linear structure whose three-dimensional shape is unknown. From the rear Brilluan scattered light distribution of the multi-core optical fiber 10 arranged along the linear structure whose shape is known, the position coordinates in the three-dimensional space of the linear structure whose three-dimensional shape is unknown are calculated and changed over time. (Step S03), I do.
  • the measuring device 20 performs step S02, and the analysis device 30 performs step S03. In step S02, it is desirable to correct the temperature dependence of the Brilluan frequency shift.
  • step S03 The amount of strain at the position z of each core of the multi-core optical fiber obtained from the rear Brillouin scattered light distribution when the multi-core optical fiber is arranged along a linear structure whose three-dimensional shape is unknown, and the above 3 It is a difference from the strain amount of the position z of each core of the multi-core optical fiber obtained from the rear Brillouin scattered light distribution when the multi-core optical fiber is arranged along a linear structure whose dimensional shape is known.
  • the bending strain ⁇ of each of the outer peripheral cores is calculated by subtracting the differential strain of the central core from the differential strain of each of the outer peripheral cores (step S32).
  • the curvature ⁇ and the bending angle ⁇ at the position z of the multi-core optical fiber are calculated from the bending strain ⁇ for each outer core (step S33), and the bending angle ⁇ is differentiated by the arc length.
  • the Frenet-Serret formula calculate the position coordinates in the 3D space of the linear structure whose 3D shape is unknown from the curvature ⁇ and the torsion at the position z, and time. Identifying changes (step S34), I do.
  • FIG. 5 is a diagram illustrating a specific procedure when measuring a linear structure whose three-dimensional shape is unknown by this shape measuring method.
  • the multi-core optical fiber 10 is laid along a linear structure whose three-dimensional shape is known (for example, a water pipe that has already been laid), and the rear Brillouin scattered light of each core in a steady state (reference state). Get the distribution characteristics of. This is used as the reference data.
  • the distribution characteristics of the rear Brillouin scattered light of each core are acquired again. This is used as comparison data.
  • Step C Next, the difference between the comparison data and the reference data is calculated for each core and each position z to derive the difference strain.
  • Step D Next, the differential strain of the central core is subtracted from the differential strain of each outer peripheral core, and the bending strain ⁇ at the position z is derived for each outer peripheral core.
  • Step E Substitute the ⁇ derived by Step D into the relational expression (1) between the bending strain ⁇ of the outer peripheral core, the curvature ⁇ of the multi-core optical fiber 10 and the bending angle ⁇ .
  • ⁇ r ⁇ ⁇ ⁇ cos ( ⁇ - ⁇ )
  • r is the distance between the centers of the central core 11 and the outer peripheral core 12
  • ⁇ is an angle representing the position of the outer peripheral core 12 in the cross section of the multi-core optical fiber 10.
  • r is 125 ⁇ m and ⁇ is 0 °, 120 °, and 240 ° for each outer core.
  • a ternary simultaneous equation can be obtained. From this equation, the curvature ⁇ and the bending angle ⁇ at the position z are calculated.
  • the bending strain ⁇ differs for each core, but since the curvature ⁇ and the bending angle ⁇ are the same in all cores, the curvature ⁇ and the bending angle ⁇ at the position z are determined by the least squares method.
  • Step F Using the Frenet-Serret formula, the position vector (three-dimensional shape) of the multi-core optical fiber 10 is determined from the curvature ⁇ and the bending angle ⁇ for each distance z determined by Step E.
  • FIG. 6 is a diagram illustrating a first example in which the three-dimensional shape of the multi-core optical fiber 10 is measured by the shape measuring method.
  • the longitudinal distance (position z described above) of the multi-core optical fiber 10 is represented by s [m].
  • the three-dimensional space in which the multi-core optical fiber 10 is arranged is represented by the x-axis, the y-axis, and the z-axis. Note that the z-axis here is different from the position z described above.
  • the reference data of Step A in FIG. 5 was acquired in a state where the multi-core optical fiber 10 for sensing was linearly stretched.
  • the vicinity of the center of the multi-core optical fiber 10 was rotated clockwise with a constant curvature, and the comparison data of Step B in FIG. 5 was acquired.
  • 6 (a), (b), and (c) are the results of the reference data of Step A, the comparison data of Step B, and the difference distortion of Step C, respectively.
  • the dotted line, broken line, and alternate long and short dash line in the figure indicate the difference between the outer peripheral cores.
  • 6 (d) and 6 (e) show the evaluation results of the curvature ⁇ and the angle ⁇ of the multi-core optical fiber 10, respectively.
  • FIG. 6F shows the spatial coordinates of the multi-core optical fiber 10 in consideration of the position vector, and the curvature given to the vicinity of the center of the multi-core optical fiber 10 can be detected with high accuracy.
  • FIG. 7 is a diagram illustrating a second example in which the three-dimensional shape of the multi-core optical fiber 10 is measured by the shape measuring method. Also in FIGS. 7A to 7E, the distance in the longitudinal direction of the multi-core optical fiber 10 (position z described above) is represented by s [m]. Also in FIG. 7 (f), the three-dimensional space in which the multi-core optical fiber 10 is arranged is represented by the x-axis, the y-axis, and the z-axis. Note that the z-axis here is different from the position z described above.
  • the reference data of Step A in FIG. 5 was acquired in a state where the multi-core optical fiber 10 for sensing was linearly stretched.
  • the vicinity of the center of the multi-core optical fiber 10 was rotated counterclockwise with a constant curvature, and the comparison data of Step B in FIG. 5 was acquired.
  • 7 (a), (b), and (c) are the results of the reference data of Step A, the comparison data of Step B, and the difference distortion of Step C, respectively.
  • the dotted line, broken line, and alternate long and short dash line in the figure indicate the difference between the outer peripheral cores.
  • 7 (d) and 7 (e) show the evaluation results of the curvature ⁇ and the angle ⁇ of the multi-core optical fiber 10, respectively.
  • FIG. 7 (f) shows the spatial coordinates of the multi-core optical fiber 10 in consideration of the position vector, and the shape change in the direction opposite to that of FIG. 6 (f) (counterclockwise) can be detected with high accuracy.
  • the system configuration of the shape measurement system of the present embodiment is different from the system configuration of the shape measurement system described in FIG. 1 in that the multi-core optical fiber 10 is given a known twist and the analysis procedure of the analysis device 30 is different. ..
  • FIG. 8 is a diagram illustrating a multi-core optical fiber 10.
  • FIG. 8A describes bending strain measured by the multi-core optical fiber 10 in a twisted state
  • FIG. 8B describes bending strain measured by the multi-core optical fiber 10 in a twisted state.
  • the core 0 is the central core
  • the cores 1 to 3 are the outer peripheral cores. If there is no bending strain, the bending strain becomes zero and overlaps in all cores regardless of the presence or absence of twisting. In addition, the bending strain of the central core is zero regardless of the presence or absence of twisting.
  • FIG. 8 (B) bending with a constant curvature is applied as in FIG. 8 (A), but in the section where bending is applied, the positions of the outer peripheral cores are exchanged by twisting, so that the amount of strain is positive or negative. Can be confirmed to be inverted.
  • twisting is required in all the sections arranged in the linear structure.
  • the analyzer 30 has a linear structure whose three-dimensional shape is unknown based on the torsional strain generated in the multi-core optical fiber 10 and the known torsional strain when calculating the position coordinates. Estimate the unintended twist that occurs when arranging along an object, and exclude the effects of the unintended twist.
  • step S33 described in the shape measuring method of FIG. 4 is different from the description of the first embodiment. That is, in step S03, The amount of strain at the position z of each core of the multi-core optical fiber obtained from the rear Brillouin scattered light distribution when the multi-core optical fiber is arranged along a linear structure whose three-dimensional shape is unknown, and the above-mentioned The difference from the strain amount of the position z of each core of the multi-core optical fiber obtained from the rear Brillouin scattered light distribution when the multi-core optical fiber is arranged along a linear structure whose three-dimensional shape is known.
  • Step S31 The bending strain ⁇ bending, i of each of the outer peripheral cores is calculated by subtracting the differential strain of the central core from the differential strain of each of the outer peripheral cores (step S32). Using the relational expression (2) excluding the influence of the unintended twist, the curvature ⁇ and the bending angle ⁇ at the position z of the multi-core optical fiber are calculated from the bending strain ⁇ bending, i for each outer core. (Step S33), and using the Frenet-Serret formula, calculate the position coordinates in the three-dimensional space of the linear structure whose three-dimensional shape is unknown from the curvature ⁇ and the bending angle ⁇ at the position z (step S34). ), I do.
  • FIG. 9 is a diagram illustrating a specific procedure when measuring a linear structure whose three-dimensional shape is unknown by this shape measuring method.
  • "i" is a core number.
  • i of 1 or more is the outer peripheral core.
  • the multi-core optical fiber 10 is laid along a linear structure whose three-dimensional shape is known (for example, a water pipe that has already been laid), and the rear Brillouin scattered light of each core in a steady state (reference state). Get the distribution characteristics of. This is used as the reference data.
  • Step B Next, in a state where the three-dimensional shape of the linear structure has changed (for example, a water pipe that is expected to be deformed due to an earthquake or the like), the distribution characteristics of the rear Brillouin scattered light of each core are acquired again. This is used as comparison data.
  • Step C Next, the difference between the comparison data and the reference data is calculated for each core and each position z to derive the difference strain.
  • Step D Next, the differential strain of the central core is subtracted from the differential strain of each outer peripheral core, and the strains ⁇ mix and i at the position z are derived for each outer peripheral core (i here is 1 or more).
  • the strains ⁇ mix and i include bending strains ⁇ bending and i and torsional strains ⁇ twisting (where i is 1 or more). Since the relative positional relationship of the outer peripheral cores does not change, the torsional strain ⁇ twisting at the position z is used as shown in the following equation by utilizing the fact that the bending strain ⁇ bending and i of all the outer peripheral cores become zero when added. Is calculated. max ⁇ i ⁇ means the maximum core number. The torsional strain ⁇ twisting calculated by the equation (7) is substituted into the equation (5) to calculate the specific torsion angle ⁇ (z) at the position z.
  • FIG. 11 is an image diagram illustrating the specific twist angle ⁇ .
  • the specific twist angle ⁇ is a twist angle per unit length.
  • r is the distance from the center to the center of the outer peripheral core in the cross section of the multi-core optical fiber 10.
  • k2 is a torsion correction coefficient and is expressed by the following equation.
  • p is the spin rate (the amount of twist per unit length) of the outer peripheral core, which is a design value or a value obtained by a measurement method described later.
  • Step E The bending strain ⁇ bending, i of the outer peripheral core i is calculated by subtracting the ⁇ twisting derived by Step D- ⁇ from the strain ⁇ mix, i of each outer peripheral core.
  • the specific twist angle ⁇ (z) calculated in the equation (5) is substituted into the equation (4) to calculate the angle ⁇ i representing the position of the outer peripheral core i at the position z.
  • a i is the initial angle of the outer peripheral core i.
  • the bending strain ⁇ bending, i of the outer peripheral core i and the angle ⁇ i derived by the equation (4) are substituted into the relational expression (2) of the curvature ⁇ and the bending angle ⁇ of the multi-core optical fiber 10.
  • k 1 is a twist correction coefficient represented by the equation (3).
  • ⁇ is the Poisson's ratio of the multi-core optical fiber.
  • a ternary simultaneous equation can be obtained.
  • FIG. 10 is an image diagram illustrating a radius r, a bending direction (curvature ⁇ ), and a bending angle ⁇ of a multi-core optical fiber having an outer peripheral core of 3.
  • 10 (A) is a view seen from the side surface of the multi-core optical fiber
  • FIG. 10 (B) is a cross-sectional view of the multi-core optical fiber.
  • Step F Using the Frenet-Serret formula, the position vector (three-dimensional shape) of the multi-core optical fiber 10 is determined from the curvature ⁇ and the bending angle ⁇ for each distance z determined by Step E. In addition, in order to improve the position accuracy, it is preferable to correct the three-dimensional shape by using a known end point.
  • FIG. 12 is a diagram showing the results of measuring the strain of each core of the twisted multi-core optical fiber in the z direction.
  • the number of twists may differ from the design value due to manufacturing errors and the like. Therefore, the strain period fluctuates in each core. Therefore, the number of twists (spin rate) is calculated from the periodic fluctuation of the strain by the RTM (short-time Fourier transform).
  • the value of the spin rate p calculated from the measured strain period can be used.
  • the shape measurement system according to the present invention can realize a measurement dynamic range of several km to several tens of km by using rear Brillouin scattered light for shape identification of a linear structure. Further, by setting the distance between the center core of the multi-core optical fiber used for detecting the shape change and the center of the outer peripheral core to 120 ⁇ m or more, it is possible to detect a minute change having a curvature of 0.5 [1 / m] or less.
  • Multi-core optical fiber 11 Central core 12: Outer peripheral core 20: Measuring device 30: Analytical device

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Nonlinear Science (AREA)
  • Optics & Photonics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Length Measuring Devices By Optical Means (AREA)
PCT/JP2021/021283 2020-06-05 2021-06-03 形状測定システム及び形状測定方法 Ceased WO2021246497A1 (ja)

Priority Applications (3)

Application Number Priority Date Filing Date Title
JP2022528897A JP7376052B2 (ja) 2020-06-05 2021-06-03 形状測定システム及び形状測定方法
EP21818720.1A EP4163586B1 (en) 2020-06-05 2021-06-03 Shape measurement system and shape measurement method
US17/915,924 US12203847B2 (en) 2020-06-05 2021-06-03 Shape measurement system and shape measurement method

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2020098349 2020-06-05
JP2020-098349 2020-06-05

Publications (1)

Publication Number Publication Date
WO2021246497A1 true WO2021246497A1 (ja) 2021-12-09

Family

ID=78831274

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2021/021283 Ceased WO2021246497A1 (ja) 2020-06-05 2021-06-03 形状測定システム及び形状測定方法

Country Status (4)

Country Link
US (1) US12203847B2 (https=)
EP (1) EP4163586B1 (https=)
JP (1) JP7376052B2 (https=)
WO (1) WO2021246497A1 (https=)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2023120055A1 (ja) * 2021-12-24 2023-06-29 日本電信電話株式会社 形状測定システム及び形状測定方法
CN118149723A (zh) * 2024-05-09 2024-06-07 武汉理工大学 一种长距离水下线缆姿态检测方法及系统
WO2025047077A1 (ja) * 2023-08-28 2025-03-06 日本電信電話株式会社 測定システム、測定方法、及び、解析装置

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115900579B (zh) * 2023-01-06 2023-05-26 山东大学 一种自校正可拼接式光纤位移场传感系统及其校正方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2013505441A (ja) * 2009-09-18 2013-02-14 ルナ イノベーションズ インコーポレイテッド 光学的位置および／または形状センシング
JP2016102691A (ja) * 2014-11-27 2016-06-02 日本電信電話株式会社 光ファイバの曲げ形状測定装置及びその曲げ形状測定方法
JP2018527041A (ja) * 2015-06-15 2018-09-20 コーニンクレッカ フィリップス エヌ ヴェＫｏｎｉｎｋｌｉｊｋｅ Ｐｈｉｌｉｐｓ Ｎ．Ｖ． 後方散乱反射光測定を用いて医療デバイスの位置及び／又は形状を感知するための光学的形状感知システム並びに方法
JP2019522776A (ja) * 2016-05-11 2019-08-15 インテュイティブ サージカル オペレーションズ， インコーポレイテッド 安全性のためのマルチコア光ファイバ内の冗長コア
JP2019531487A (ja) * 2016-08-10 2019-10-31 フラウンホファー ゲセルシャフト ツール フェールデルンク ダー アンゲヴァンテン フォルシュンク エー．ファオ． 光学導波路の屈曲及び／またはひずみを決定するための方法

Family Cites Families (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4985405B2 (ja) * 2005-09-29 2012-07-25 住友電気工業株式会社 センサ及びそれを用いた外乱測定方法
JP2007101508A (ja) * 2005-10-07 2007-04-19 Sumitomo Electric Ind Ltd 温度測定方法及び温度測定装置
JP5769676B2 (ja) * 2012-08-17 2015-08-26 公益財団法人地球環境産業技術研究機構 物質の圧力、温度、ひずみ分布測定システム、これを用いた二酸化炭素地中貯留の監視方法、二酸化炭素注入による地層安定性への影響評価方法、および結氷監視方法
US9366526B2 (en) * 2012-11-30 2016-06-14 Neubrex Co., Ltd. Three-dimensional position measurement system
US10228556B2 (en) * 2014-04-04 2019-03-12 The General Hospital Corporation Apparatus and method for controlling propagation and/or transmission of electromagnetic radiation in flexible waveguide(s)
CN105371781B (zh) * 2015-11-13 2018-09-07 华中科技大学 一种三维形状测量方法
US10145681B2 (en) * 2016-07-19 2018-12-04 Corning Incorporated Brillouin-based distributed bend fiber sensor and method for using same
US10530114B2 (en) * 2017-08-31 2020-01-07 United States Of America As Represented By The Administrator Of Nasa Polarization maintaining, large mode area (PMVLMA) erbium-doped optical fiber and amplifier
CN110243302B (zh) * 2018-03-08 2020-11-06 桂林电子科技大学 反射式多芯循环串接光纤形状传感器
EP3832254B1 (en) * 2018-07-31 2023-08-30 Furukawa Electric Co., Ltd. Cable, cable shape sensing system, sensing system, and method for sensing cable shape
EP3969866B1 (en) * 2019-05-12 2025-12-24 Hampidjan HF. A process for ascertaining the elongation of a load-bearing cable
WO2021086129A1 (ko) * 2019-10-31 2021-05-06 (주)옵토닉스 3차원 곡선 형상을 측정하기 위한 특수 광섬유 및 그의 제조 방법, 그리고 특수 광섬유를 이용하여 3차원 곡선 형상을 측정하는 시스템

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2013505441A (ja) * 2009-09-18 2013-02-14 ルナ イノベーションズ インコーポレイテッド 光学的位置および／または形状センシング
JP2016102691A (ja) * 2014-11-27 2016-06-02 日本電信電話株式会社 光ファイバの曲げ形状測定装置及びその曲げ形状測定方法
JP2018527041A (ja) * 2015-06-15 2018-09-20 コーニンクレッカ フィリップス エヌ ヴェＫｏｎｉｎｋｌｉｊｋｅ Ｐｈｉｌｉｐｓ Ｎ．Ｖ． 後方散乱反射光測定を用いて医療デバイスの位置及び／又は形状を感知するための光学的形状感知システム並びに方法
JP2019522776A (ja) * 2016-05-11 2019-08-15 インテュイティブ サージカル オペレーションズ， インコーポレイテッド 安全性のためのマルチコア光ファイバ内の冗長コア
JP2019531487A (ja) * 2016-08-10 2019-10-31 フラウンホファー ゲセルシャフト ツール フェールデルンク ダー アンゲヴァンテン フォルシュンク エー．ファオ． 光学導波路の屈曲及び／またはひずみを決定するための方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
ROGER G. DUNCAN, HIGH-ACCURACY FIBER-OPTIC SHAPE SENSING, 30 April 2020 (2020-04-30), Retrieved from the Internet <URL:https://doi.org/10.1117/12.720914>
See also references of EP4163586A4
YONG-LAE PARK: "Real-Time Estimation of 3-D Needle Shape and Deflection for MRI-Guided Interventions", IEEE/ASME, vol. 15, no. 6, 2010, XP011339830, DOI: 10.1109/TMECH.2010.2080360

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2023120055A1 (ja) * 2021-12-24 2023-06-29 日本電信電話株式会社 形状測定システム及び形状測定方法
JP7637395B2 (ja) 2021-12-24 2025-02-28 日本電信電話株式会社 形状測定システム及び形状測定方法
WO2025047077A1 (ja) * 2023-08-28 2025-03-06 日本電信電話株式会社 測定システム、測定方法、及び、解析装置
CN118149723A (zh) * 2024-05-09 2024-06-07 武汉理工大学 一种长距离水下线缆姿态检测方法及系统

Also Published As

Publication number Publication date
EP4163586B1 (en) 2025-04-09
EP4163586A1 (en) 2023-04-12
JP7376052B2 (ja) 2023-11-08
US20230147800A1 (en) 2023-05-11
JPWO2021246497A1 (https=) 2021-12-09
EP4163586A4 (en) 2024-02-21
US12203847B2 (en) 2025-01-21

Similar Documents

Publication Publication Date Title
JP7376052B2 (ja) 形状測定システム及び形状測定方法
US8970845B1 (en) In-situ three-dimensional shape rendering from strain values obtained through optical fiber sensors
JP2023112062A (ja) 安全性のためのマルチコア光ファイバ内の冗長コア
CN105865752B (zh) 采用分布式偏振串扰分析仪全面评判保偏光纤特性的方法和装置
JP6346852B2 (ja) 光ファイバの曲げ形状測定装置及びその曲げ形状測定方法
Duncan et al. Characterization of a fiber-optic shape and position sensor
CN111982000A (zh) 一种基于Beta标架的光纤形状重构方法及装置
JP7315009B2 (ja) 光ファイバケーブルセンシング装置、光ファイバケーブルセンシング方法、及びプログラム
JP7306565B2 (ja) テープ心線を用いた測定装置及び測定方法
Liu et al. 3D‐Shape Sensing via Monte Carlo‐Optimized Spatial Resolution in Multi‐Channel Parallel OFDR
Askins et al. Bend and twist sensing in a multi-core optical fiber
JP7637395B2 (ja) 形状測定システム及び形状測定方法
Zhou et al. Scheme and algorithm designs for optical fiber sensors based configuration inversion of umbilical cables in deep-water
JP2022110605A (ja) ケーブル、センシングケーブル、およびセンシングシステム
JP7406767B2 (ja) 光ファイバケーブルセンシングシステム、光ファイバケーブルセンシング方法、及び光ファイバケーブル
Maeda et al. Shape sensing for large-scale line structures with multi-core fibers by using brillouin OTDR
JP7406768B2 (ja) 光ファイバケーブルセンシング装置、光ファイバケーブルセンシング方法、及びプログラム
CN117606377A (zh) 一种基于弱反射光纤光栅的海洋柔性管道曲率监测方法
Nakamoto et al. An approach to low-curvature 2D shape sensing using an indoor optical cable
WO2025047077A1 (ja) 測定システム、測定方法、及び、解析装置
JP7523071B2 (ja) 光ファイバ撚り周期算出システム、光ファイバ撚り周期算出方法、光ファイバ撚り周期算出装置及び光ファイバ撚り周期算出プログラム
JP2026065458A (ja) 解析システム、解析装置、及び、解析方法
OGURA et al. Detection of Buckling in Columnar Structures by Shape Sensing with Multi-Core Optical Fibers
Ogura et al. Shape Sensing of Large-diameter Multi-core Fibers by DSS with Different Spatial Resolutions
WO2024228359A1 (ja) 検出システム、検出方法、及び、解析装置

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 21818720

Country of ref document: EP

Kind code of ref document: A1

ENP Entry into the national phase

Ref document number: 2022528897

Country of ref document: JP

Kind code of ref document: A

WWE Wipo information: entry into national phase

Ref document number: 2021818720

Country of ref document: EP

ENP Entry into the national phase

Ref document number: 2021818720

Country of ref document: EP

Effective date: 20230105

NENP Non-entry into the national phase

Ref country code: DE

WWG Wipo information: grant in national office

Ref document number: 2021818720

Country of ref document: EP