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
Application number
PCT/JP2021/021283
Other languages
English (en)
French (fr)
Inventor
信智 半澤
和秀 中島
隆 松井
英晶 村山
良太 和田
真輝人 小林
Original Assignee
日本電信電話株式会社
国立大学法人東京大学
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 日本電信電話株式会社, 国立大学法人東京大学 filed Critical 日本電信電話株式会社
Priority to JP2022528897A priority Critical patent/JP7376052B2/ja
Priority to EP21818720.1A priority patent/EP4163586A4/en
Priority to US17/915,924 priority patent/US20230147800A1/en
Publication of WO2021246497A1 publication Critical patent/WO2021246497A1/ja

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)

Abstract

線状の被測定物の３次元形状を長距離且つ高分解能で導出することができる形状測定システム及び形状測定方法を提供することを目的とする。　本発明に係る形状測定システムは、断面の中心に配置される中心コア１１、及び中心コア１１の外側かつ同心円状に等間隔で配置される３個以上の外周コア１２を有するマルチコア光ファイバ１０と、マルチコア光ファイバ１０の各コアの伝搬方向における後方ブリルアン散乱光分布を計測する計測装置２０と、３次元形状が未知である線状構造物に沿って配置されたマルチコア光ファイバ１０と３次元形状が既知である線状構造物に沿って配置されたマルチコア光ファイバ１０の前記後方ブリルアン散乱光分布から前記３次元形状が未知である線状構造物の３次元空間における位置座標を算出する解析装置３０と、を備える。

Description

形状測定システム及び形状測定方法
　本開示は、光ファイバを用いたリフレクトメトリ技術により、パイプラインや海底ケーブル等の被測定対象の３次元形状を導出する装置および方法に関する。
　マルチコア光ファイバの各コアの周波数領域の反射スペクトルをＯＦＤＲ（Ｏｐｔｉｃａｌ　Ｆｒｅｑｕｅｎｃｙ　Ｄｏｍａｉｎ　Ｒｅｆｌｅｃｔｏｍｅｔｒｙ）で測定及び解析し、被測定対象の３次元形状を導出する技術が提案されている（例えば、非特許文献１を参照。）。また、センシング媒体（光ファイバ）の全長にＦＢＧ（Ｆｉｂｅｒ　Ｂｒａｇｇ　Ｇｒａｔｉｎｇ）を付与し、測定分解能を改善する手法も提案されている（例えば、非特許文献２を参照。）。
Ｒｏｇｅｒ　Ｇ．　Ｄｕｎｃａｎ，　ｅｔｃ．，　"Ｈｉｇｈ－ａｃｃｕｒａｃｙ　ｆｉｂｅｒ－ｏｐｔｉｃ　ｓｈａｐｅ　ｓｅｎｓｉｎｇ"，　ｈｔｔｐｓ：／／ｄｏｉ．ｏｒｇ／１０．１１１７／１２．７２０９１４，　２０２０年４月３０日検索 Ｙｏｎｇ－Ｌａｅ　Ｐａｒｋ，　ｅｔｃ．，　"Ｒｅａｌ－Ｔｉｍｅ　Ｅｓｔｉｍａｔｉｏｎ　ｏｆ　３－Ｄ　Ｎｅｅｄｌｅ　Ｓｈａｐｅ　ａｎｄ　Ｄｅｆｌｅｃｔｉｏｎ　ｆｏｒ　ＭＲＩ－Ｇｕｉｄｅｄ　Ｉｎｔｅｒｖｅｎｔｉｏｎｓ"，　ＩＥＥＥ／ＡＳＭＥ　ＴＲＡＮＳＡＣＴＩＯＮＳ　ＯＮ　ＭＥＣＨＡＴＲＯＮＩＣＳ，　ＶＯＬ．　１５，　ＮＯ．　６，　ＤＥＣＥＭＢＥＲ　２０１０
　ＯＦＤＲで被測定対象の３次元形状を導出しようとした場合、次のような困難性が発生する。
（１）ＯＦＤＲは、数１０ｍｍオーダーの高分解能を実現できるが、測定距離が数１０ｍ程度に限定され、長距離の３次元形状を導出することが困難である。
（２）光ファイバにＦＢＧを付与すれば、その光ファイバの製造容易性や経済性を向上させることが困難となる。
（３）ＦＢＧの付与など準分布測定により形状同定を行う技術は、測定点数に制約があり、長距離にわたる形状同定が困難である。
（４）ＦＢＧを付与した光ファイバを非線形材料の被測定対象に埋め込んだ形態の場合、温度や応力でＦＢＧによる反射の中心周波数が変化し、計測する光ファイバの歪と曲率にも非線形性が生じるため、被測定物の形状変化を特定することが困難である。
　そこで、本発明は、前記課題を解決するために、線状の被測定物の３次元形状を長距離且つ高分解能で導出することができる形状測定システム及び形状測定方法を提供することを目的とする。
　上記目的を達成するために、本発明に係る形状測定システムは、所定のコア配置を持つマルチコア光ファイバをセンシング媒体として用い、ＢＯＴＤＲ（Ｂｒｉｌｌｏｕｉｎ　Ｏｐｔｉｃａｌ　Ｔｉｍｅ　Ｄｏｍａｉｎ　Ｒｅｆｌｅｃｔｏｍｅｔｒｙ）で取得したデータを解析することとした。
　具体的には、本発明に係る形状測定システムは、
　断面の中心に配置される中心コア、及び前記中心コアの外側かつ同心円状に等間隔で配置される３個以上の外周コアを有するマルチコア光ファイバと、
　前記マルチコア光ファイバの各コアの伝搬方向における後方ブリルアン散乱光分布を計測する計測装置と、
　３次元形状が未知である線状構造物に沿って配置された前記マルチコア光ファイバと３次元形状が既知である線状構造物に沿って配置された前記マルチコア光ファイバの前記後方ブリルアン散乱光分布から前記３次元形状が未知である線状構造物の３次元空間における位置座標を算出し、その時間変化を同定する解析装置と、
を備える。
　また、本発明に係る形状測定方法は、
　断面の中心に配置される中心コア、及び前記中心コアの外側かつ同心円状に等間隔で配置される３個以上６個以下の外周コアを有するマルチコア光ファイバを線状構造物に沿って配置すること、
　前記マルチコア光ファイバの各コアの伝搬方向における後方ブリルアン散乱光分布を計測すること、及び
　３次元形状が未知である線状構造物に沿って配置された前記マルチコア光ファイバと３次元形状が既知である線状構造物に沿って配置された前記マルチコア光ファイバの前記後方ブリルアン散乱光分布から前記３次元形状が未知である線状構造物の３次元空間における位置座標を算出し、その時間変化を同定すること、
を行う。
　マルチコア光ファイバのコア毎に後方ブリルアン散乱光を測定することで長距離測定を可能とする。また、マルチコア光ファイバの中心コアと外周コアとの中心間距離、並びに隣接する外周コア間の距離を保つことで被測定物の微小変動を検知することができる。従って、本発明は、線状の被測定物の３次元形状を長距離且つ高分解能で導出することができる形状測定システム及び形状測定方法を提供することができる。
　具体的な解析手法は次の通りである。
　前記マルチコア光ファイバの長手方向の位置をｚとすると、
　前記解析装置は、
　前記３次元形状が未知である線状構造物に沿って前記マルチコア光ファイバを配置したときの前記後方ブリルアン散乱光分布から得られた前記マルチコア光ファイバの各コアの位置ｚの歪量と、前記３次元形状が既知である線状構造物に沿って前記マルチコア光ファイバを配置したときの前記後方ブリルアン散乱光分布から得られた前記マルチコア光ファイバの各コアの位置ｚの歪量との差分である位置ｚの差分歪を計算すること、
　前記外周コアそれぞれの前記差分歪から前記中心コアの前記差分歪を減算して前記外周コアそれぞれの曲げ歪εを計算すること、
　関係式（１）を用い、前記外周コア毎の前記曲げ歪εから前記マルチコア光ファイバの位置ｚにおける曲率κと曲げ角度βを計算し、曲げ角度βを弧長で微分して捩率を計算すること、及び
　フレネ・セレの公式を用い、位置ｚにおける曲率κと捩率から前記３次元形状が未知である線状構造物の３次元空間における位置座標を算出し、その時間変化を同定すること、
を行うことを特徴とする。
尚、曲げ角度βはＭＣＦ断面において、ＭＣＦ中心と曲率中心を結ぶベクトルと基準方向ベクトル、例えば、ＭＣＦ中心とある特定の一つのコアを結ぶベクトルのなす角度を表す。
［関係式（１）］
　ε＝r・κ・ｃｏｓ（α－β）
ただし、rは前記中心コアと前記外周コアとの中心間距離、αは前記マルチコア光ファイバの断面における前記外周コアの位置を表す角度である。
　本手法において、構造物にマルチコア光ファイバを沿わせる際に、意図しない捩じりが生じることにより、計測結果に誤差が生じる可能性がある。このような場合、
　前記マルチコア光ファイバには、既知の捩じりが与えられていること、及び
　前記解析装置は、前記位置座標を算出するときに、前記マルチコア光ファイバに生じている捩じり歪と前記既知の捩じりによる歪とに基づいて前記３次元形状が未知である線状構造物に沿って配置するときに発生した意図しない捩じりを推定し、前記意図しない捩じりによる影響を除外すること、
が好ましい。
　この場合の具体的な解析手法は次の通りである。
　前記マルチコア光ファイバの長手方向の位置をｚとすると、
　前記解析装置は、
　前記３次元形状が未知である線状構造物に沿って前記マルチコア光ファイバを配置したときの前記後方ブリルアン散乱光分布から得られた前記マルチコア光ファイバの各コアの位置ｚの歪量と、前記３次元形状が既知である線状構造物に沿って前記マルチコア光ファイバを配置したときの前記後方ブリルアン散乱光分布から得られた前記マルチコア光ファイバの各コアの位置ｚの歪量との差分である位置ｚの差分歪を計算すること、
　前記外周コアそれぞれの前記差分歪から前記中心コアの前記差分歪を減算して前記外周コアそれぞれの曲げ歪εｂｅｎｄｉｎｇ，ｉを計算すること、
　前記意図しない捩じりによる影響を除外する関係式（２）を用い、前記外周コア毎の前記曲げ歪εｂｅｎｄｉｎｇ，ｉから前記マルチコア光ファイバの位置ｚにおける曲率κと曲げ角度βを計算すること、及び
　フレネ・セレの公式を用い、位置ｚにおける曲率κと曲げ角度βから前記３次元形状が未知である線状構造物の３次元空間における位置座標を算出すること、
を行うことを特徴とする。
Figure JPOXMLDOC01-appb-M000011
 ただし、rは前記中心コアと前記外周コアとの中心間距離、ｋ１は式（３）で表される捩じりの補正係数、νは前記マルチコア光ファイバのポアソン比、ａｉはコアｉの初期角度、ωｉは式（４）で表される前記マルチコア光ファイバの位置ｚでの断面における前記外周コアの位置を表す角度、ｐは前記外周コアのスピンレート、εｔｗｉｓｔｉｎｇは前記マルチコア光ファイバの位置ｚに発生している捩じり歪、φ（ｚ）は式（５）で表される前記マルチコア光ファイバの位置ｚでの比捩じり角、及びｋ２は式（６）で表される捩じりの補正係数である。
Figure JPOXMLDOC01-appb-M000012
Figure JPOXMLDOC01-appb-M000013
Figure JPOXMLDOC01-appb-M000014
Figure JPOXMLDOC01-appb-M000015
　この場合、前記解析装置は、前記外周コアそれぞれの曲げ歪εｂｅｎｄｉｎｇ，ｉの周期的な変動から前記外周コアのスピンレートｐを算出することが好ましい。マルチコア光ファイバの製造誤差などにより捩じり回数の誤差をも考慮でき、被測定物の形状変化を高い精度で特定することができる。
　前記マルチコア光ファイバは、前記外周コアの数は３以上であり、クラッド直径が３７５μｍ以上であることが好ましい。
　また、前記マルチコア光ファイバは、前記中心コアと前記外周コアとの中心間距離が１２０μｍ以上であることが好ましい。
　既存のシングルモードファイバ（クラッド直径が１２５μｍ）をバンドルすることでファンイン／ファンアウト（ＦＩＦＯ）デバイスを容易に実現できる。また、中心コアと外周コアの中心間の間隔が広いほど外周コアに生じる歪が大きくなり、高ダイナミックレンジ化することができる。
　なお、上記各発明は、可能な限り組み合わせることができる。
　本発明は、線状の被測定物の３次元形状を長距離且つ高分解能で導出することができる形状測定システム及び形状測定方法を提供することができる。
本発明に係る形状測定システムの構成を説明する図である。 本発明に係る形状測定システムのマルチコア光ファイバの断面を説明する図である。 本発明に係る形状測定システムの効果を説明する図である。 本発明に係る形状測定方法を説明する図である。 本発明に係る形状測定方法を説明する図である。 本発明に係る形状測定システムでの計測例を説明する図である。 本発明に係る形状測定システムでの計測例を説明する図である。 本発明に係る形状測定システムに備わるマルチコア光ファイバを説明する図である。 本発明に係る形状測定方法を説明する図である。 本発明に係る形状測定システムに備わるマルチコア光ファイバを説明する図である。 マルチコア光ファイバの比捩じり角を説明する図である。 マルチコア光ファイバのスピンレート実測方法を説明する図である。
　添付の図面を参照して本発明の実施形態を説明する。以下に説明する実施形態は本発明の実施例であり、本発明は、以下の実施形態に制限されるものではない。なお、本明細書及び図面において符号が同じ構成要素は、相互に同一のものを示すものとする。
［実施形態１］
　図１は、本実施形態の形状測定システムを説明する図である。本形状測定システムは、
　断面の中心に配置される中心コア１１、及び中心コア１１の外側かつ同心円状に等間隔で配置される３個以上の外周コア１２を有するマルチコア光ファイバ１０と、
　マルチコア光ファイバ１０の各コアの伝搬方向における後方ブリルアン散乱光分布を計測する計測装置２０と、
　３次元形状が未知である線状構造物に沿って配置されたマルチコア光ファイバ１０と３次元形状が既知である線状構造物に沿って配置されたマルチコア光ファイバ１０の前記後方ブリルアン散乱光分布から前記３次元形状が未知である線状構造物の３次元空間における位置座標を算出し、その時間変化を同定する解析装置３０と、
を備える。
　尚、図１では、マルチコア光ファイバ１０の特定の１コアと、後方ブリルアン散乱光分布を計測する計測装置２０とが接続される例を示している。他の構成として、マルチコア光ファイバの各コアを単一コア光ファイバに分離するファンアウト機構、および光スイッチを介してマルチコア光ファイバ１０と計測装置２０とを接続する形態でも構わない。
　本形状測定システムは、センシング媒体であるマルチコア光ファイバ１０と、マルチコア光ファイバ１０の各コアの伝搬方向における後方ブリルアン散乱光を検出する計測装置（ＢＯＴＤＲ）２０と、ＢＯＴＤＲ１２で取得された測定データを解析する解析装置３０とを備える。
　マルチコア光ファイバ１０は、被測定対象である線状構造物の長手方向に沿って設置される。線状構造物とは、例えば、パイプライン、ライザー（海洋掘削や海洋生産において海底から海面上の設備まで流体が通れるパイプ）、海底ケーブル等である。図２は、マルチコア光ファイバ１０の断面を説明する図である。マルチコア光ファイバ１０は、断面の中心にある中心コア、および断面中心から概ね同一円周上に概ね等間隔で配置された外周コアの計４個のコアを有する。本実施形態では、４個のコアは概ね同等の屈折率分布および光学特性を有するものを使用したが、各コアの屈折率分布や光学特性が意図的に異なるように配置される構造でも構わない。
　また、コアの数の制限については次の通りである。コア数が３以下（外周コア２個以下）では、形状の同定を行うことができない。一方、コア数が５以上（外周コア４個以上）の場合、測定精度の面では向上するが、測定時間はコア数分増える。
　また、ある半径の円は、それに接する同じ半径の６個の円で取り囲むことができる。この場合、円の数はトータルで７個である。つまり、同じ直径のシングルコアの光ファイバを７本を最密で配置したときのコア位置のマルチコア光ファイバとすることができ、ＦＩＦＯを容易に作ることができる。
　一方、外周の円が７個以上（円の数はトータルで８個以上）になると、１つの円の周りに７つ以上の円を配置すると中心の円と外周の円との間に隙間が生まれる。つまり、同じ直径のシングルコアの光ファイバを８本以上では最密で配置することができず、マルチコア光ファイバのクラッド径が大きくなり、ファンイン／ファンアウト（ＦＩＦＯ）の作成が難しくなる。
　従って、マルチコア光ファイバの外周コア数は３個以上６個以下が望ましい。
　マルチコア光ファイバ１０は、一般的な通信用光ファイバより大きなクラッド外径（＞１２５μｍ）を有する。光ファイバの中心に配置された中心コア１１と、中心コア１１の外周に配置された外周コア１２との中心間距離は３０μｍより大きく設定する。その理由は、クラッド外径１２５μｍの標準的なマルチコア光ファイバにおいて漏洩損失などの影響を抑えて配置可能な中心コアと外周コアの中心間距離が３０μｍ程度であるからである。
　図３は、関係式（１）において（α－β）が０としたマルチコア光ファイバ１０の中心の曲率に対する前記中心間距離と外周コア１２の歪量との関係を説明する図である。図３より、中心コア１１と外周コア１２の中心間距離を１２０μｍ以上で、曲率κ＝０．５［１／ｍ］においてもクラッド外径１２５μｍの標準的なマルチコア光ファイバで得られる歪量の５倍以上の歪量を得ることができることがわかる。これは、中心コア１１と外周コア１２の中心間距離を長くすることで、マルチコア光ファイバ１０に発生した小さい曲率の形状変化に対して感度を向上できることを意味する。
　以上説明したように、既存のシングルモードファイバ（クラッド直径１２５μｍ）をバンドルすることでＦＩＦＯデバイスを容易に実現できること、並びに既存のシングルモードファイバで作成されたＦＩＦＯを用いればマルチコア光ファイバ１０の中心コア１１と外周コア１２の中心間距離が１２５μｍとなり小さな形状変化に対しても十分な感度を得られること、から、本実施形態では、マルチコア光ファイバ１０のクラッド外径を３７５μｍと（＝１２５μｍ×３）として検討を行った。
　次に、図４を用いてマルチコア光ファイバ１０を用いて被測定物の形状変化を測定する形状測定方法を説明する。本形状測定方法は、
　マルチコア光ファイバ１０を線状構造物に沿って配置すること（ステップＳ０１）、
　マルチコア光ファイバ１０の各コアの伝搬方向における後方ブリルアン散乱光分布を計測すること（ステップＳ０２）、及び
　３次元形状が未知である線状構造物に沿って配置されたマルチコア光ファイバ１０と３次元形状が既知である線状構造物に沿って配置されたマルチコア光ファイバ１０の後方ブリルアン散乱光分布から３次元形状が未知である線状構造物の３次元空間における位置座標を算出し、時間変化を同定すること（ステップＳ０３）、
を行う。なお、計測装置２０がステップＳ０２を行い、解析装置３０がステップＳ０３を行う。ステップＳ０２では、ブリルアン周波数シフトの温度依存性を補正することが望ましい。
　ステップＳ０３についてより詳細に説明する、測定装置２０からマルチコア光ファイバ１０の長手方向の位置をｚとする。
　ステップＳ０３では、
　３次元形状が未知である線状構造物に沿って前記マルチコア光ファイバを配置したときの前記後方ブリルアン散乱光分布から得られた前記マルチコア光ファイバの各コアの位置ｚの歪量と、前記３次元形状が既知である線状構造物に沿って前記マルチコア光ファイバを配置したときの前記後方ブリルアン散乱光分布から得られた前記マルチコア光ファイバの各コアの位置ｚの歪量との差分である位置ｚの差分歪を計算すること（ステップＳ３１）、
　前記外周コアそれぞれの前記差分歪から前記中心コアの前記差分歪を減算して前記外周コアそれぞれの曲げ歪εを計算すること（ステップＳ３２）、
　関係式（１）を用い、前記外周コア毎の前記曲げ歪εから前記マルチコア光ファイバの位置ｚにおける曲率κと曲げ角度βを計算すること（ステップＳ３３）、曲げ角度βを弧長で微分して捩率を計算すること、及び
　フレネ・セレの公式を用い、位置ｚにおける曲率κと捩率から前記３次元形状が未知である線状構造物の３次元空間における位置座標を算出し、時間変化を同定すること（ステップＳ３４）、
を行う。
　図５は、本形状測定方法で３次元形状が未知である線状構造物を計測するときの具体的手順を説明する図である。
［Ｓｔｅｐ　Ａ］
　まず、３次元形状が既知である線状構造物（例えば、すでに敷設されている水道管等）に沿ってマルチコア光ファイバ１０を敷設し、定常状態（基準状態）における各コアの後方ブリルアン散乱光の分布特性を取得する。これを基準データとする。
［Ｓｔｅｐ　Ｂ］
　次に、線状構造物の３次元形状が変化した状態（例えば、地震などで変形が想定される水道管等）で、改めて各コアの後方ブリルアン散乱光の分布特性を取得する。これを比較データとする。
［Ｓｔｅｐ　Ｃ］
　次に、コア毎且つ位置ｚ毎に比較データと基準データの差分を計算して差分歪を導出する。
［Ｓｔｅｐ　Ｄ］
　次に、それぞれの外周コアの差分歪から中心コアの差分歪を減算し、外周コア毎に位置ｚにおける曲げ歪εを導出する。
［Ｓｔｅｐ　Ｅ］
　外周コアの曲げ歪εと、マルチコア光ファイバ１０の曲率κと曲げ角度βの関係式（１）に、Ｓｔｅｐ　Ｄで導出したεを代入する。
［関係式（１）］
　ε＝ｒ・κ・ｃｏｓ（α－β）
ただし、ｒは中心コア１１と外周コア１２との中心間距離、αはマルチコア光ファイバ１０の断面における外周コア１２の位置を表す角度である。例えば、ｒは１２５μｍ、αは外周コア毎に０°、１２０°、２４０°である。
　本実施形態では、外周コアが３つあるので３元連立方程式が得られる。この式により位置ｚにおける曲率κと曲げ角度βを算出する。ここで、曲げ歪εはコアごとに異なるが、曲率κと曲げ角度βはいずれのコアにおいても等しいため、最小二乗法で位置ｚにおける曲率κと曲げ角度βを決定する。
［Ｓｔｅｐ　Ｆ］
　フレネ・セレの公式を用い、Ｓｔｅｐ　Ｅで決定した距離ｚ毎の曲率κと曲げ角度βから、マルチコア光ファイバ１０の位置ベクトル（３次元形状）を決定する。なお、位置精度を向上するために、既知の終点を用いて３次元形状を補正することが好ましい。
（実施例１）
　図６は、形状測定方法でマルチコア光ファイバ１０の３次元形状を測定した第１の例を説明する図である。図６（ａ）～（ｅ）では、マルチコア光ファイバ１０の長手方向の距離（上で説明した位置ｚ）をｓ［ｍ］で表している。図６（ｆ）では、マルチコア光ファイバ１０が配置された３次元空間をｘ軸、ｙ軸、ｚ軸で表している。ここでのｚ軸は上で説明した位置ｚと異なることに注意する。
　本実施例では、センシング用のマルチコア光ファイバ１０を直線状に延伸した状態で、図５のＳｔｅｐ　Ａの基準データを取得した。次に、マルチコア光ファイバ１０の中央付近を時計回りに一定曲率で回転させ、図５のＳｔｅｐ　Ｂの比較データを取得した。図６（ａ）、（ｂ）、（ｃ）は、それぞれ図５のＳｔｅｐ　Ａの基準データ、Ｓｔｅｐ　Ｂの比較データ、Ｓｔｅｐ　Ｃの差分歪の結果である。図中の点線、破線、一点鎖線は、それぞれ外周コアの違いを表す。図６（ｄ）、（ｅ）は、それぞれマルチコア光ファイバ１０の曲率κおよび角度βの評価結果を表す。図６（ｆ）は、位置ベクトルを考慮した、マルチコア光ファイバ１０の空間座標を表しており、マルチコア光ファイバ１０の中央付近に与えた曲率が精度良く検出できている。
（実施例２）
　図７は、形状測定方法でマルチコア光ファイバ１０の３次元形状を測定した第２の例を説明する図である。図７（ａ）～（ｅ）でも、マルチコア光ファイバ１０の長手方向の距離（上で説明した位置ｚ）をｓ［ｍ］で表している。図７（ｆ）でも、マルチコア光ファイバ１０が配置された３次元空間をｘ軸、ｙ軸、ｚ軸で表している。ここでのｚ軸は上で説明した位置ｚと異なることに注意する。
　本実施例では、センシング用のマルチコア光ファイバ１０を直線状に延伸した状態で、図５のＳｔｅｐ　Ａの基準データを取得した。次に、マルチコア光ファイバ１０の中央付近を反時計回りに一定曲率で回転させ、図５のＳｔｅｐ　Ｂの比較データを取得した。図７（ａ）、（ｂ）、（ｃ）は、それぞれ図５のＳｔｅｐ　Ａの基準データ、Ｓｔｅｐ　Ｂの比較データ、Ｓｔｅｐ　Ｃの差分歪の結果である。図中の点線、破線、一点鎖線は、それぞれ外周コアの違いを表す。図７（ｄ）、（ｅ）は、それぞれマルチコア光ファイバ１０の曲率κおよび角度βの評価結果を表す。図７（ｆ）は、位置ベクトルを考慮した、マルチコア光ファイバ１０の空間座標を表しており、図６（ｆ）と逆方向（反時計回り）の形状変化が精度良く検出できている。
［実施形態２］
　センシング媒体として線状構造物にマルチコア光ファイバを沿わせる際に、意図しない捩じりが生じ、計測結果に誤差が生じる可能性がある。そこで、本実施形態では、センシング媒体に使用するマルチコア光ファイバに意図的に捩じりを加えることとした。線状構造物の形状同定の際にその捩じり歪以下の意図しない捩じりによる歪分を計算で排除することができ、形状同定の精度を向上させることができる。
　本実施形態の形状測定システムのシステム構成は図１で説明した形状測定システムのシステム構成に対し、マルチコア光ファイバ１０に既知の捩じれが与えられていること、及び解析装置３０の解析手順が相違する。
　図８は、マルチコア光ファイバ１０について説明する図である。図８（Ａ）は捩じれが無い状態のマルチコア光ファイバ１０で測定される曲げ歪、図８（Ｂ）は捩じれがある状態のマルチコア光ファイバ１０で測定される曲げ歪を説明している。コア０が中心コア、コア１からコア３が外周コアである。
　なお、曲げ歪が無い場合、曲げ歪は捩じれの有無によらず全てのコアでゼロとなり重なる。また、中心コアについては捩じれの有無によらず、曲げ歪はゼロとなる。捩じれが無い状態のマルチコア光ファイバ１０の場合、曲率が一定の曲げが付与されると測定される曲げ歪は図８（Ａ）の軸方向の位置ｚ＝２ｍからｚ＝１２ｍの間のようにコアごとに異なる歪量であるが一定の歪量が観測される。
　また、図８（Ｂ）は、軸方向の位置ｚ＝２ｍからｚ＝１２ｍの間で捩じれを付与したマルチコア光ファイバ１０の曲げ歪を説明している。図８（Ｂ）では、図８（Ａ）と同様に曲率が一定の曲げを付与しているが、曲げが付与された区間において、外周コアの位置が捻じれによって入れ替わるため、歪量の正負が反転していることが確認できる。マルチコア光ファイバ１０をセンシング媒体として線状構造物に配置する場合、線状構造物に配置する区間全てに捩じりが必要である。
　解析装置３０は、前記位置座標を算出するときに、マルチコア光ファイバ１０に生じている捩じり歪と前記既知の捩じりによる歪とに基づいて前記３次元形状が未知である線状構造物に沿って配置するときに発生した意図しない捩じりを推定し、前記意図しない捩じりによる影響を除外する。
　本形状測定方法は、図４の形状測定方法で説明したステップＳ３３が実施形態１の説明と異なる。すなわち、ステップＳ０３では、
　前記３次元形状が未知である線状構造物に沿って前記マルチコア光ファイバを配置したときの前記後方ブリルアン散乱光分布から得られた前記マルチコア光ファイバの各コアの位置ｚの歪量と、前記３次元形状が既知である線状構造物に沿って前記マルチコア光ファイバを配置したときの前記後方ブリルアン散乱光分布から得られた前記マルチコア光ファイバの各コアの位置ｚの歪量との差分である位置ｚの差分歪を計算すること（ステップＳ３１）、
　前記外周コアそれぞれの前記差分歪から前記中心コアの前記差分歪を減算して前記外周コアそれぞれの曲げ歪εｂｅｎｄｉｎｇ，ｉを計算すること（ステップＳ３２）、
　前記意図しない捩じりによる影響を除外する関係式（２）を用い、前記外周コア毎の前記曲げ歪εｂｅｎｄｉｎｇ，ｉから前記マルチコア光ファイバの位置ｚにおける曲率κと曲げ角度βを計算すること（ステップＳ３３）、及び
　フレネ・セレの公式を用い、位置ｚにおける曲率κと曲げ角度βから前記３次元形状が未知である線状構造物の３次元空間における位置座標を算出すること（ステップＳ３４）、
を行う。
　図９は、本形状測定方法で３次元形状が未知である線状構造物を計測するときの具体的手順を説明する図である。ここで、“ｉ”はコア番号とする。なお、ｉ＝０は中心コア、１以上のｉは外周コアである。
［Ｓｔｅｐ　Ａ］
　まず、３次元形状が既知である線状構造物（例えば、すでに敷設されている水道管等）に沿ってマルチコア光ファイバ１０を敷設し、定常状態（基準状態）における各コアの後方ブリルアン散乱光の分布特性を取得する。これを基準データとする。
［Ｓｔｅｐ　Ｂ］
　次に、線状構造物の３次元形状が変化した状態（例えば、地震などで変形が想定される水道管等）で、改めて各コアの後方ブリルアン散乱光の分布特性を取得する。これを比較データとする。
［Ｓｔｅｐ　Ｃ］
　次に、コア毎且つ位置ｚ毎に比較データと基準データの差分を計算して差分歪を導出する。
［Ｓｔｅｐ　Ｄ］
　次に、それぞれの外周コアの差分歪から中心コアの差分歪を減算し、外周コア毎に位置ｚにおける歪εｍｉｘ，ｉを導出する（ここのｉは１以上）。
［Ｓｔｅｐ　Ｄ－α］
　歪εｍｉｘ，ｉには、曲げ歪εｂｅｎｄｉｎｇ，ｉと捩じり歪εｔｗｉｓｔｉｎｇが含まれる（ここのｉは１以上）。外周コアの相対的な位置関係は変化しないので、全ての外周コアｉの曲げ歪εｂｅｎｄｉｎｇ，ｉを加算するとゼロになることを利用し、次式のように位置ｚにおける捩じり歪εｔｗｉｓｔｉｎｇを算出する。ｍａｘ｛ｉ｝は最大のコア番号を意味する。
Figure JPOXMLDOC01-appb-M000016
 　式（７）で算出した捩じり歪εｔｗｉｓｔｉｎｇを式（５）に代入して位置ｚにおける比捩じり角φ（ｚ）を算出する。図１１は比捩じり角φを説明するイメージ図である。比捩じり角φは、単位長当たりの捩じれ角度である。
Figure JPOXMLDOC01-appb-M000017
 　ここで、ｒはマルチコア光ファイバ１０の断面において中心から外周コアの中心までの距離である。ｋ２は捩じりの補正係数であり、次式で表される。ｐは外周コアのスピンレート（単位長あたりの捩じれ量）であり、設計上の値もしくは後述する測定手法で得た値である。
Figure JPOXMLDOC01-appb-M000018
 ［Ｓｔｅｐ　Ｅ］
　各外周コアの歪εｍｉｘ，ｉからＳｔｅｐ　Ｄ－αで導出したεｔｗｉｓｔｉｎｇを減算して外周コアｉの曲げ歪εｂｅｎｄｉｎｇ，ｉを算出する。
　そして、式（４）に式（５）で算出した比捩じり角φ（ｚ）を代入し、位置ｚにおける外周コアｉの位置を表わす角度ωｉを算出する。ａｉは外周コアｉの初期角度である。
Figure JPOXMLDOC01-appb-M000019
 　さらに、外周コアｉの曲げ歪εｂｅｎｄｉｎｇ，ｉと、式（４）で導出した角度ωｉをマルチコア光ファイバ１０の曲率κと曲げ角度βの関係式（２）に代入する。
Figure JPOXMLDOC01-appb-M000020
 ただし、ｋ１は式（３）で表される捩じりの補正係数である。
Figure JPOXMLDOC01-appb-M000021
 ここで、νは前記マルチコア光ファイバのポアソン比である。
　本実施形態では、外周コアが３つあるので３元連立方程式が得られる。関係式（２）により位置ｚにおける曲率κと曲げ角度βを算出する。ここで、曲げ歪εはコアごとに異なるが、曲率κと曲げ角度βはいずれのコアにおいても等しいため、最小二乗法で位置ｚにおける曲率κと曲げ角度βを決定する。
　なお、図１０は、外周コアが３であるマルチコア光ファイバの半径ｒ、曲げ方向（曲率κ）、曲げ角度βを説明するイメージ図である。図１０（Ａ）はマルチコア光ファイバの側面から見た図、図１０（Ｂ）はマルチコア光ファイバの断面図である。
［Ｓｔｅｐ　Ｆ］
　フレネ・セレの公式を用い、Ｓｔｅｐ　Ｅで決定した距離ｚ毎の曲率κと曲げ角度βから、マルチコア光ファイバ１０の位置ベクトル（３次元形状）を決定する。なお、位置精度を向上するために、既知の終点を用いて３次元形状を補正することが好ましい。
（スピンレート測定方法）
　ここで、スピンレートを実測する方法について説明する。
　図１２は、捩じりを加えたマルチコア光ファイバの各コアの歪をｚ方向に測定した結果を示す図である。マルチコア光ファイバを捩じった場合、製造誤差などにより捩じり回数（スピンレート）に設計値との誤差が生じる。このため、各コアとも、歪の周期が変動する。そこで、ＳＴＦＴ（ｓｈｏｒｔ－ｔｉｍｅ　Ｆｏｕｒｉｅｒ　ｔｒａｎｓｆｏｒｍ）により周期的な歪の変動から捩じり回数（スピンレート）を算出する。上述したステップＤ－αでは、この実測した歪の周期から算出したスピンレートｐの値を利用することができる。
（効果）
　本発明に係る形状測定システムは、線状構造物の形状同定に後方ブリルアン散乱光を用いることで、数ｋｍから数１０ｋｍの測定ダイナミックレンジを実現できる。また、形状変化の検出に用いるマルチコア光ファイバの中心コアと外周コアの中心間距離を１２０μｍ以上としたことで、曲率０．５［１／ｍ］以下の微小変化を検出できる。
１０：マルチコア光ファイバ
１１：中心コア
１２：外周コア
２０：計測装置
３０：解析装置

Claims (14)

  1. 　断面の中心に配置される中心コア、及び前記中心コアの外側かつ同心円状に等間隔で配置される３個以上の外周コアを有するマルチコア光ファイバと、
    　前記マルチコア光ファイバの各コアの伝搬方向における後方ブリルアン散乱光分布を計測する計測装置と、
    　３次元形状が未知である線状構造物に沿って配置された前記マルチコア光ファイバと３次元形状が既知である線状構造物に沿って配置された前記マルチコア光ファイバの前記後方ブリルアン散乱光分布から前記３次元形状が未知である線状構造物の３次元空間における位置座標を算出する解析装置と、
    を備える形状測定システム。
  2. 　前記マルチコア光ファイバの長手方向の位置をｚとすると、
    　前記解析装置は、
    　前記３次元形状が未知である線状構造物に沿って前記マルチコア光ファイバを配置したときの前記後方ブリルアン散乱光分布から得られた前記マルチコア光ファイバの各コアの位置ｚの歪量と、前記３次元形状が既知である線状構造物に沿って前記マルチコア光ファイバを配置したときの前記後方ブリルアン散乱光分布から得られた前記マルチコア光ファイバの各コアの位置ｚの歪量との差分である位置ｚの差分歪を計算すること、
    　前記外周コアそれぞれの前記差分歪から前記中心コアの前記差分歪を減算して前記外周コアそれぞれの曲げ歪εを計算すること、
    　関係式（１）を用い、前記外周コア毎の前記曲げ歪εから前記マルチコア光ファイバの位置ｚにおける曲率κと曲げ角度βを計算すること、及び
    　フレネ・セレの公式を用い、位置ｚにおける曲率κと曲げ角度βから前記３次元形状が未知である線状構造物の３次元空間における位置座標を算出すること、
    を行うことを特徴とする請求項１に記載の形状測定システム。
    ［関係式（１）］
    　ε＝ｒ・κ・ｃｏｓ（α－β）
    ただし、ｒは前記中心コアと前記外周コアとの中心間距離、αは前記マルチコア光ファイバの断面における前記外周コアの位置を表す角度である。
  3. 　前記マルチコア光ファイバは、既知の捩じりが与えられていること、及び
    　前記解析装置は、前記位置座標を算出するときに、前記マルチコア光ファイバに生じている捩じり歪と前記既知の捩じりによる歪とに基づいて前記３次元形状が未知である線状構造物に沿って配置するときに発生した意図しない捩じりを推定し、前記意図しない捩じりによる影響を除外すること、
    を特徴とする請求項１に記載の形状測定システム。
  4. 　前記マルチコア光ファイバの長手方向の位置をｚとすると、
    　前記解析装置は、
    　前記３次元形状が未知である線状構造物に沿って前記マルチコア光ファイバを配置したときの前記後方ブリルアン散乱光分布から得られた前記マルチコア光ファイバの各コアの位置ｚの歪量と、前記３次元形状が既知である線状構造物に沿って前記マルチコア光ファイバを配置したときの前記後方ブリルアン散乱光分布から得られた前記マルチコア光ファイバの各コアの位置ｚの歪量との差分である位置ｚの差分歪を計算すること、
    　前記外周コアそれぞれの前記差分歪から前記中心コアの前記差分歪を減算して前記外周コアそれぞれの曲げ歪εｂｅｎｄｉｎｇ，ｉを計算すること、
    　前記意図しない捩じりによる影響を除外する関係式（２）を用い、前記外周コア毎の前記曲げ歪εｂｅｎｄｉｎｇ，ｉから前記マルチコア光ファイバの位置ｚにおける曲率κと曲げ角度βを計算すること、及び
    　フレネ・セレの公式を用い、位置ｚにおける曲率κと曲げ角度βから前記３次元形状が未知である線状構造物の３次元空間における位置座標を算出すること、
    を行うことを特徴とする請求項３に記載の形状測定システム。
    Figure JPOXMLDOC01-appb-M000001
     ただし、rは前記中心コアと前記外周コアとの中心間距離、ｋ１は式（３）で表される捩じりの補正係数、νは前記マルチコア光ファイバのポアソン比、ａｉはコアｉの初期角度、ωｉは式（４）で表される前記マルチコア光ファイバの位置ｚでの断面における前記外周コアの位置を表す角度、ｐは前記外周コアのスピンレート、εｔｗｉｓｔｉｎｇは前記マルチコア光ファイバの位置ｚに発生している捩じり歪、φ（ｚ）は式（５）で表される前記マルチコア光ファイバの位置ｚでの比捩じり角、及びｋ２は式（６）で表される捩じりの補正係数である。
    Figure JPOXMLDOC01-appb-M000002
    Figure JPOXMLDOC01-appb-M000003
    Figure JPOXMLDOC01-appb-M000004
    Figure JPOXMLDOC01-appb-M000005
  5. 　前記解析装置は、前記外周コアそれぞれの曲げ歪εｂｅｎｄｉｎｇ，ｉの周期的な変動から前記外周コアのスピンレートｐを算出することを特徴とする請求項４に記載の形状測定システム。
  6. 　前記マルチコア光ファイバは、前記外周コアの数は３であり、クラッド直径が３７５μｍ以上であることを特徴とする請求項１から５のいずれかに記載の測定システム。
  7. 　前記マルチコア光ファイバは、前記中心コアと前記外周コアとの中心間距離が１２０μｍ以上であることを特徴とする請求項１から６のいずれかに記載の測定システム。
  8. 　断面の中心に配置される中心コア、及び前記中心コアの外側かつ同心円状に等間隔で配置される３個以上の外周コアを有するマルチコア光ファイバを線状構造物に沿って配置すること、
    　前記マルチコア光ファイバの各コアの伝搬方向における後方ブリルアン散乱光分布を計測すること、及び
    　３次元形状が未知である線状構造物に沿って配置された前記マルチコア光ファイバと３次元形状が既知である線状構造物に沿って配置された前記マルチコア光ファイバの前記後方ブリルアン散乱光分布から前記３次元形状が未知である線状構造物の３次元空間における位置座標を算出すること、
    を行う形状測定方法。
  9. 　前記３次元形状が未知である線状構造物の３次元空間における位置座標を算出するときに、前記マルチコア光ファイバの長手方向の位置をｚとし、
    　前記３次元形状が未知である線状構造物に沿って前記マルチコア光ファイバを配置したときの前記後方ブリルアン散乱光分布から得られた前記マルチコア光ファイバの各コアの位置ｚの歪量と、前記３次元形状が既知である線状構造物に沿って前記マルチコア光ファイバを配置したときの前記後方ブリルアン散乱光分布から得られた前記マルチコア光ファイバの各コアの位置ｚの歪量との差分である位置ｚの差分歪を計算すること、
    　前記外周コアそれぞれの前記差分歪から前記中心コアの前記差分歪を減算して前記外周コアそれぞれの曲げ歪εを計算すること、
    　関係式（１）を用い、前記外周コア毎の前記曲げ歪εから前記マルチコア光ファイバの位置ｚにおける曲率κと曲げ角度βを計算すること、及び
    　フレネ・セレの公式を用い、位置ｚにおける曲率κと曲げ角度βから前記３次元形状が未知である線状構造物の３次元空間における位置座標を算出すること、
    を行うことを特徴とする請求項８に記載の形状測定方法。
    ［関係式（１）］
    　ε＝ｒ・κ・ｃｏｓ（α－β）
    ただし、ｒは前記中心コアと前記外周コアとの中心間距離、αは前記マルチコア光ファイバの断面における前記外周コアの位置を表す角度である。
  10. 　前記マルチコア光ファイバには、既知の捩じりが与えられていること、及び
    　前記位置座標を算出するときに、前記マルチコア光ファイバに生じている捩じり歪と前記既知の捩じりによる歪とに基づいて前記３次元形状が未知である線状構造物に沿って配置するときに発生した意図しない捩じりを推定し、前記意図しない捩じりによる影響を除外すること、
    を特徴とする請求項８に記載の形状測定方法。
  11. 　前記３次元形状が未知である線状構造物の３次元空間における位置座標を算出するときに、前記マルチコア光ファイバの長手方向の位置をｚとし、
    　前記３次元形状が未知である線状構造物に沿って前記マルチコア光ファイバを配置したときの前記後方ブリルアン散乱光分布から得られた前記マルチコア光ファイバの各コアの位置ｚの歪量と、前記３次元形状が既知である線状構造物に沿って前記マルチコア光ファイバを配置したときの前記後方ブリルアン散乱光分布から得られた前記マルチコア光ファイバの各コアの位置ｚの歪量との差分である位置ｚの差分歪を計算すること、
    　前記外周コアそれぞれの前記差分歪から前記中心コアの前記差分歪を減算して前記外周コアそれぞれの曲げ歪εｂｅｎｄｉｎｇ，ｉを計算すること、
    　前記意図しない捩じりによる影響を除外する関係式（２）を用い、前記外周コア毎の前記曲げ歪εｂｅｎｄｉｎｇ，ｉから前記マルチコア光ファイバの位置ｚにおける曲率κと曲げ角度βを計算すること、及び
    　フレネ・セレの公式を用い、位置ｚにおける曲率κと曲げ角度βから前記３次元形状が未知である線状構造物の３次元空間における位置座標を算出すること、
    を行うことを特徴とする請求項１０に記載の形状測定方法。
    Figure JPOXMLDOC01-appb-M000006
     ただし、ｒは前記中心コアと前記外周コアとの中心間距離、ｋ１は式（３）で表される捩じりの補正係数、νは前記マルチコア光ファイバのポアソン比、ａｉはコアｉの初期角度、ωｉは式（４）で表される前記マルチコア光ファイバの位置ｚでの断面における前記外周コアの位置を表す角度、ｐは前記外周コアのスピンレート、εｔｗｉｓｔｉｎｇは前記マルチコア光ファイバの位置ｚに発生している捩じり歪、φ（ｚ）は式（５）で表される前記マルチコア光ファイバの位置ｚでの比捩じり角、及びｋ２は式（６）で表される捩じりの補正係数である。
    Figure JPOXMLDOC01-appb-M000007
    Figure JPOXMLDOC01-appb-M000008
    Figure JPOXMLDOC01-appb-M000009
    Figure JPOXMLDOC01-appb-M000010
  12. 　前記外周コアそれぞれの曲げ歪εｂｅｎｄｉｎｇ，ｉの周期的な変動から前記外周コアのスピンレートｐを算出することを特徴とする請求項１１に記載の形状測定方法。
  13. 　前記マルチコア光ファイバは、前記外周コアの数は３であり、クラッド直径が３７５μｍ以上であることを特徴とする請求項８から１２のいずれかに記載の形状測定方法。
  14. 　前記マルチコア光ファイバは、前記中心コアと前記外周コアとの中心間距離が１２０μｍ以上であることを特徴とする請求項８から１３のいずれかに記載の形状測定方法。
PCT/JP2021/021283 2020-06-05 2021-06-03 形状測定システム及び形状測定方法 WO2021246497A1 (ja)

Priority Applications (3)

Application Number Priority Date Filing Date Title
JP2022528897A JP7376052B2 (ja) 2020-06-05 2021-06-03 形状測定システム及び形状測定方法
EP21818720.1A EP4163586A4 (en) 2020-06-05 2021-06-03 FORM MEASURING SYSTEM AND FORM MEASURING METHOD
US17/915,924 US20230147800A1 (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
JP2020-098349 2020-06-05
JP2020098349 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 WO2021246497A1 (ja) 2020-06-05 2021-06-03 形状測定システム及び形状測定方法

Country Status (4)

Country Link
US (1) US20230147800A1 (ja)
EP (1) EP4163586A4 (ja)
JP (1) JP7376052B2 (ja)
WO (1) WO2021246497A1 (ja)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2023120055A1 (ja) * 2021-12-24 2023-06-29 日本電信電話株式会社 形状測定システム及び形状測定方法

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 (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105371781B (zh) * 2015-11-13 2018-09-07 华中科技大学 一种三维形状测量方法
CN110243302B (zh) * 2018-03-08 2020-11-06 桂林电子科技大学 反射式多芯循环串接光纤形状传感器

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 (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2023120055A1 (ja) * 2021-12-24 2023-06-29 日本電信電話株式会社 形状測定システム及び形状測定方法

Also Published As

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

Similar Documents

Publication Publication Date Title
CN107003473B (zh) 用于应变和温度分离的多芯光纤中的不同的纤芯
CN108603977B (zh) 具有用于安全性的冗余纤芯的多纤芯光学纤维
CN105865752B (zh) 采用分布式偏振串扰分析仪全面评判保偏光纤特性的方法和装置
WO2014083989A1 (ja) ３次元位置計測装置
US20100215311A1 (en) Method and Apparatus for Shape and End Position Determination Using an Optical Fiber
US20080212082A1 (en) Fiber optic position and/or shape sensing based on rayleigh scatter
Askins et al. Bend and twist sensing in a multiple-core optical fiber
Duncan et al. Characterization of a fiber-optic shape and position sensor
JP6346852B2 (ja) 光ファイバの曲げ形状測定装置及びその曲げ形状測定方法
WO2021246497A1 (ja) 形状測定システム及び形状測定方法
CN111982000A (zh) 一种基于Beta标架的光纤形状重构方法及装置
Marković et al. Application of fiber-optic curvature sensor in deformation measurement process
Kang et al. Strain measurements on a cantilever beam with fiber Bragg grating sensors using a pair of collimators
Askins et al. Bend and twist sensing in a multi-core optical fiber
JP7406767B2 (ja) 光ファイバケーブルセンシングシステム、光ファイバケーブルセンシング方法、及び光ファイバケーブル
JP2016102690A (ja) 光ファイバの曲げ形状測定装置及びその曲げ形状測定方法
JP7315009B2 (ja) 光ファイバケーブルセンシング装置、光ファイバケーブルセンシング方法、及びプログラム
WO2023120055A1 (ja) 形状測定システム及び形状測定方法
Maeda et al. Shape Sensing for Large-scale Line Structures with Multi-core Fibers by Using Brillouin OTDR
CN115667840A (zh) 用于沿着三维空间测量微弯和任意微变形的系统
Lebang et al. Detection of displacement using glass optical fiber sensor with various configuration
JP7406768B2 (ja) 光ファイバケーブルセンシング装置、光ファイバケーブルセンシング方法、及びプログラム
US20230152170A1 (en) Measuring device and measuring method using tape core wire
JP2023014654A (ja) 光ファイバ曲げ方向算出システム、光ファイバ曲げ方向算出方法、光ファイバ曲げ方向算出装置及び光ファイバ曲げ方向算出プログラム
OGURA et al. Detection of Buckling in Columnar Structures by Shape Sensing with Multi-Core Optical Fibers

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