Location evaluation device, location evaluation method, and program

Abstract

A position evaluation device (50) according to the present disclosure comprises: a vibration characteristics calculation unit (51) in which a signal indicating the characteristic vibration generated at each position of an optical fiber is input from a sensor and sensing data indicating the vibration characteristics of each position of the optical fiber is calculated on the basis of the input signal; a dissimilarity calculation unit (52) that calculates dissimilarity in the sensing data between two adjacent points in the optical fiber; and an environmental change position estimation unit (53) that uses the dissimilarity as a basis to estimate an environmental change position at which the environment of the optical fiber changes.

Description

[Technical Field]

[0001] The present disclosure relates to a location assessment device, a location assessment method, and a computer-readable medium. [Background technology]

[0002] A technology called optical fiber sensing makes it possible to detect vibrations / sounds occurring in any section of an optical fiber. Specifically, in optical fiber sensing, an optical fiber sensor inputs coherent pulsed light into an optical fiber and receives the backscattered light of the pulsed light from the optical fiber. At this time, the optical fiber sensor detects the phase difference between the backscattered light generated at two points on the optical fiber, thereby detecting vibrations / sounds acting on the optical fiber in the phase difference evaluation section (gauge length section) between the two points. Such optical fiber sensors are realized by phase-sensitive OTDRs (Phase-Sensitive Optical Time Domain Reflectometers) or DASs (Distributed Acoustic Sensors), but the following description will be given assuming that the optical fiber sensor is a DAS.

[0003] However, existing optical fiber cables for communication, including existing optical fibers, are sometimes laid in areas away from the ground surface. Examples of such optical fiber cables for communication include optical fiber cables suspended on poles such as utility poles and steel towers, and OPGW (Optical Ground Wire).

[0004] Figure 1 shows an example of the configuration of a sensing system using existing optical fibers suspended on pillars aligned in one dimension. In the sensing system shown in Fig. 1, existing optical fibers are suspended from poles 1 to 3, which are utility poles, steel towers, etc. A DAS is connected to one end of the optical fibers.

[0005] DAS can detect the environment around the optical fiber based on vibration information that indicates vibrations acting on the optical fiber, such as wind and rain hitting the optical fiber, the presence or absence of lightning, the vibration mode of the pole, and the presence or absence of living things.

[0006] DAS can also detect abnormalities occurring around optical fibers based on sound information indicating sounds applied to the optical fibers, such as gunshots, explosions, and abnormal sounds caused by accidents.

[0007] However, if the installation conditions of the optical fiber are not taken into consideration, there is a possibility that the location of the vibration / sound generation may be estimated incorrectly. Examples of installation conditions of optical fiber include excess optical fiber sections generated by optical fiber fusion splicing work, and overhead sections of optical fiber suspended above a pole.

[0008] Figure 2 shows an example of incorrect estimation of the location of vibration due to excess length of optical fiber. As shown in Figure 2, the position information of the vibration occurrence point measured by the DAS is the "length of the optical fiber from the DAS to the vibration point" (hereinafter defined as "DAS coordinates"). For example, the DAS can measure the vibration occurrence point on the DAS coordinates based on the time difference between the time when pulsed light is input and the time when the backscattered light of that pulsed light is received.

[0009] In the example of Figure 2, the excess length of the optical fiber is included in pillar 2. Therefore, if vibration occurs at a position farther away from the DAS than pillar 2, the vibration occurrence point on the distance from the DAS toward the pillar (hereinafter defined as "real world coordinates") will no longer match the vibration occurrence point on the DAS coordinates.

[0010] Therefore, when building a sensing system using optical fibers suspended on pillars, it is necessary to accurately know the environmental change positions where the optical fiber environment changes as DAS coordinates. The environmental change positions of the optical fiber include, for example, the positions of the pillars on which the optical fiber is suspended, the overhead sections of the optical fiber suspended between the pillars, and the excess length sections of the optical fiber.

[0011] An example of a method for matching real-world coordinates and DAS coordinates is to generate an event at a position whose real-world coordinates are known, and then match that position on the real-world coordinates with the vibration occurrence location or vibration occurrence section measured by the DAS for the vibration caused by that event.

[0012] Figure 3 shows an example of a method for aligning the real-world coordinates with the DAS coordinates by artificially vibrating the column that suspends the optical fiber. As shown in Figure 3, when a pillar is artificially vibrated, the DAS measures the vibration occurrence section on the DAS coordinate system. Then, the position of the pillar on the real-world coordinate system is matched with the vibration occurrence section on the DAS coordinate system. In the example of Figure 3, there is an excess length section of optical fiber on the vibrated pillar. Therefore, the vibration occurrence section on the DAS coordinate system corresponds to the excess length section of optical fiber.

[0013] However, the method shown in FIG. 3 has the following problems. - Since each pillar needs to be inspected individually, it takes a lot of man-hours. When a pillar is vibrated, the vibration is transmitted to the overhead optical fiber on both sides of the pillar, causing the vibrating section to expand instantly. This method cannot be applied to poles where it is difficult to artificially generate vibrations (for example, large poles such as steel towers).

[0014] Therefore, recently, a method other than the method shown in FIG. 3 has been proposed for determining the position on the DAS coordinate system of the pillar that suspends the optical fiber. For example, Patent Document 1 discloses a technology in which a section where the strength of vibration detected by optical fiber sensing is above a threshold is determined to be a section where a unique pattern of a utility pole occurs, and the point in that section where the strongest vibration occurs is estimated to be the position of the utility pole. [Prior art documents] [Patent documents]

[0015] [Patent Document 1] International Publication No. 2020 / 044648 Summary of the Invention [Problem to be solved by the invention]

[0016] As described above, the technology described in Patent Document 1 focuses on the threshold of vibration intensity (i.e., the magnitude of vibration), and estimates the point where the vibration with the greatest intensity occurs as the location of the utility pole. Therefore, the technology described in Patent Document 1 cannot estimate the location of a utility pole for points other than the point where the vibration with the greatest intensity occurs. Therefore, it is believed that there is still room for improvement in the accuracy of estimating the position of an optical fiber where an environmental change occurs, such as the position of a pillar.

[0017] In view of the above-mentioned problems, an object of the present disclosure is to provide a position evaluation device, a position evaluation method, and a computer-readable medium that can improve the accuracy of estimating the position of an environmental change in an optical fiber. [Means for solving the problem]

[0018] According to one aspect, a position evaluation device includes: a vibration characteristic calculation unit that receives a signal indicating a natural vibration generated at each position of the optical fiber from the sensor and calculates sensing data indicating the vibration characteristic at each position of the optical fiber based on the received signal; a difference calculation unit that calculates a difference in the sensing data between two adjacent points on the optical fiber; and an environment change position estimation unit that estimates an environment change position where the environment of the optical fiber changes based on the degree of difference.

[0019] A location evaluation method according to one aspect includes: A location estimation method executed by a location estimation device, comprising: a vibration characteristics calculation step of receiving a signal indicating a natural vibration occurring at each position of the optical fiber from a sensor and calculating sensing data indicating the vibration characteristics at each position of the optical fiber based on the received signal; a difference calculation step of calculating a difference in the sensing data between two adjacent points of the optical fiber; and an environment change position estimating step of estimating an environment change position where the environment of the optical fiber changes based on the degree of difference.

[0020] According to one aspect, a computer-readable medium includes: A non-transitory computer-readable medium storing a program to be executed by a computer, The program a vibration characteristics calculation step of receiving a signal indicating a natural vibration occurring at each position of the optical fiber from a sensor and calculating sensing data indicating the vibration characteristics at each position of the optical fiber based on the received signal; a difference calculation step of calculating a difference in the sensing data between two adjacent points of the optical fiber; and an environment change position estimating step of estimating an environment change position where the environment of the optical fiber changes based on the degree of difference. [Effects of the Invention]

[0021] According to the above-described aspects, it is possible to provide a position evaluation device, a position evaluation method, and a computer-readable medium that can improve the accuracy of estimating the position of an environmental change in an optical fiber. [Brief explanation of the drawings]

[0022] [Figure 1] FIG. 1 is a diagram showing an example of the configuration of a sensing system using existing optical fibers suspended over pillars aligned in one dimension. [Figure 2] 10A and 10B are diagrams illustrating an example in which the location of vibration occurrence is erroneously estimated due to an excess length of an optical fiber. [Figure 3] FIG. 10 is a diagram showing an example of a method for aligning real-world coordinates with DAS coordinates by artificially vibrating a column that suspends an optical fiber. [Figure 4] 1 is a diagram illustrating an application example of a position evaluation device according to a first embodiment. [Figure 5]1 is a diagram illustrating an example of the configuration of a position evaluation device according to a first embodiment. [Figure 6] 4 is a diagram showing an example of a phase difference signal input as an input signal to a vibration characteristic extraction unit according to the first embodiment. FIG. [Figure 7] 5 is a flowchart showing an example of the flow of operations of the vibration characteristics extraction unit according to the first embodiment. [Figure 8] 10 is a flowchart showing an example of the flow of operations of a dissimilarity calculation unit according to the first embodiment. [Figure 9] 10 is a flowchart showing an example of the flow of operations of a pillar position calculation unit according to the first embodiment. [Figure 10] FIG. 10 is a diagram illustrating an example of a method for analyzing the degree of difference by the pillar position calculation unit according to the first embodiment. [Figure 11] 2 is a diagram showing an example of the configuration of a sensing system assumed in a specific example of the operation of the position evaluation device according to the first embodiment. FIG. [Figure 12] 5 is a diagram showing an example of the degree of difference obtained in a specific example of the operation of the position evaluation device according to the first embodiment. FIG. [Figure 13] FIG. 10 is a diagram illustrating an application example of the position evaluation device according to the second embodiment. [Figure 14] FIG. 10 is a diagram illustrating an example of the configuration of a position evaluation device according to a second embodiment. [Figure 15] 10 is a flowchart showing an example of the flow of operations of a weighted dissimilarity calculation unit according to the second embodiment. [Figure 16] 10 is a flowchart showing an example of the flow of operations of a marginal length section calculation unit according to the second embodiment. [Figure 17] FIG. 10 is a diagram illustrating an example of the configuration of a sensing system assumed in a specific example of the operation of the position evaluation device according to the second embodiment. [Figure 18] 10A and 10B are diagrams illustrating examples of frequency average values, dissimilarity, and weighted dissimilarity of power spectra obtained in a specific example of the operation of the position evaluation device according to the second embodiment. [Figure 19] FIG. 19 is an enlarged view of the X region shown in FIG. [Figure 20] FIG. 10 is a diagram illustrating an application example of a position evaluation device according to a third embodiment. [Figure 21] FIG. 10 is a diagram illustrating an example of the configuration of a position evaluation device according to a third embodiment. [Figure 22] 11 is a flowchart showing an example of the flow of operations of a pillar position calculation unit according to the third embodiment. [Figure 23] FIG. 11 is a diagram showing an example of an analysis target section determined by a pillar position calculation unit according to the third embodiment. [Figure 24] FIG. 11 is a diagram illustrating an example of a window function configured by a pillar position calculation unit according to the third embodiment. [Figure 25] FIG. 11 is a diagram illustrating an example of the relationship between the window function and the dissimilarity when an appropriate offset point is determined by the pillar position calculation unit according to the third embodiment. [Figure 26] FIG. 10 is a diagram illustrating an application example of a position evaluation device according to a fourth embodiment. [Figure 27] FIG. 10 is a diagram illustrating an example of the configuration of a position evaluation device according to a fourth embodiment. [Figure 28] 13 is a flowchart showing an example of the flow of operations of a marginal length section calculation unit according to the fourth embodiment. [Figure 29] 10 is a flowchart showing an example of the flow of operations of a pillar position calculation unit according to the fourth embodiment. [Figure 30] FIG. 13 is a diagram showing an example of an analysis target section determined by a pillar position calculation unit according to the fourth embodiment. [Figure 31] FIG. 13 is a diagram illustrating an example of cross-correlation between a weighted dissimilarity and a window function according to the fourth embodiment. [Figure 32] FIG. 13 is a diagram showing an example of a method for calculating the left extra length by a column position calculation unit according to the fourth embodiment. [Figure 33] FIG. 13 is a diagram illustrating an example of the configuration of a position evaluation device according to a fifth embodiment. [Figure 34] FIG. 2 is a block diagram illustrating an example of a hardware configuration of a computer that realizes the position evaluation device according to each embodiment. DETAILED DESCRIPTION OF THE INVENTION

[0023] Hereinafter, embodiments of the present disclosure will be described with reference to the drawings. Note that the following description and drawings have been omitted and simplified as appropriate for clarity of explanation. In addition, in each of the following drawings, the same elements are given the same reference numerals, and duplicate explanations are omitted as necessary. Furthermore, specific numerical values ​​etc. shown below are merely examples to facilitate understanding of the present disclosure, and are not limited thereto.

[0024] In each embodiment of the present disclosure described below, the following (1) and / or (2) are performed. (1) Estimation of the position of the pole that suspends the optical fiber / Estimation of the aerial section of the optical fiber Specifically, we estimate the position of each pillar that suspends the optical fiber on the DAS coordinate system, and also estimate the distance on the DAS coordinate system from the start to the end of the overhead section of the optical fiber that runs between two pillars.

[0025] (2) Estimation of excess length of optical fiber Specifically, the extent of the excess optical fiber section (optical fiber localized at a certain location in real space) is estimated on the DAS coordinate system.

[0026] The basic concept for carrying out (1) and / or (2) above in each embodiment will be described below. (1) Estimation of the position of the pole that suspends the optical fiber / Estimation of the aerial section of the optical fiber The motion of the optical fiber in the aerial section can be regarded as the motion of a string, and in each embodiment, this is utilized to estimate the aerial section of the optical fiber on the DAS coordinate system. For example, if the total weight of the optical fiber is large, the force suspending the optical fiber also becomes large. As a result, the optical fiber is constantly subjected to a large tension, making it prone to dynamic distortion. Furthermore, optical fibers are constantly swaying due to wind. In addition, the transmission speed of a signal transmitted through an optical fiber is determined by the linear density (ρ), tension (T), and span length (L) of the optical fiber. TIFF0007758178000001.tif721 is characterized. In addition, the optical fiber has a steady-state nth natural vibration mode (frequency f n ) vibration is observed. TIFF0007758178000002.tif734

[0027] On the other hand, the position of the pillar that suspends the optical fiber is the boundary part of each aerial section of the optical fiber. In each embodiment, this fact is utilized to estimate the position of the pillar on the DAS coordinate system. For example, vibrations traveling through an optical fiber are reflected at the boundary of a pillar. In addition, the vibration characteristics change for each overhead section, with the pillars as the boundary.

[0028] (2) Estimation of excess length of optical fiber In most cases, the excess length of optical fiber is glued or fixed to the pole, making it less susceptible to dynamic distortion. Therefore, the excess length of optical fiber is insensitive to background noise such as wind and vibration of the pole itself (i.e., the noise level is small). In each embodiment, this fact is utilized to estimate the excess length of optical fiber on the DAS coordinates. Each embodiment of the present disclosure will be described in detail below.

[0029] <First Embodiment> In the first embodiment, as shown in FIG. 4, the position of a pillar is estimated and output as a DAS coordinate value. To redefine the DAS coordinate, the DAS coordinate is the length of the optical fiber measured by the DAS relative to a certain point.

[0030] An example of the configuration of the position evaluation device 10 according to the first embodiment will be described with reference to FIG. As shown in FIG. 5, the position evaluation device 10 according to the first embodiment includes a vibration characteristics extraction unit 11, a dissimilarity calculation unit 12, and a pillar position calculation unit 13.

[0031] The vibration characteristic extraction unit 11 is connected to a DAS (not shown) and receives a phase difference signal of backscattered light obtained from the DAS as an input signal. The vibration characteristic extraction unit 11 then extracts the input signal in a certain time interval, calculates a power spectrum for the extracted input signal, and extracts and outputs a power spectrum in a frequency band including the fundamental vibration mode in the imaginary interval from the calculated power spectrum.

[0032] The dissimilarity calculation unit 12 determines an evaluation interval for evaluating the dissimilarity. Furthermore, the dissimilarity calculation unit 12 calculates and outputs the dissimilarity between the two power spectra between the evaluation intervals for all DAS coordinate points, based on the power spectra obtained by the vibration characteristic extraction unit 11.

[0033] Based on the dissimilarity at all DAS coordinate points obtained by the dissimilarity calculation unit 12, the pillar position calculation unit 13 estimates the position where the dissimilarity takes a maximum peak value as the position of the pillar, and outputs an output signal representing the DAS coordinate value corresponding to the estimated position.

[0034] The position evaluation device 10 according to the first embodiment will be described in more detail below. First, the input signal will be described. The DAS inputs pulsed light into an optical fiber and receives backscattered light (Rayleigh scattered light) from the optical fiber in response to the input pulsed light. The DAS also detects the phase difference between the backscattered light generated at two points on the optical fiber, and generates a phase difference signal indicating the detected phase difference. TIFF0007758178000003.tif517 can be obtained. This phase difference signal is proportional to the dynamic strain of the optical fiber in the phase difference evaluation section (gauge length section), which is the section between the two points. The vibration characteristic extraction unit 11 receives this phase difference signal as an input signal.

[0035] The phase difference signal will be described with reference to FIG. In the phase difference signal, d indicates the distance from the DAS to the measurement point in the longitudinal direction of the optical fiber, and is expressed as follows: TIFF0007758178000004.tif569 where p is the DAS coordinate label (integer). ADC is the frequency of the ADC (analog-digital converter) installed in the DAS. c is the speed of light in the optical fiber, expressed as c = c0 / n. c0 is the speed of light in a vacuum, and n is the refractive index of the optical fiber core (approximately 1.46 for a silica glass core). d unit is the spacing between discrete points in the spatial direction. For example, f ADC If is 125MHz, d unit is approximately 0.82 m.

[0036] Furthermore, t in the phase difference signal indicates the measurement time, and is expressed as follows: TIFF0007758178000005.tif1228 where q is the time interval label (an integer). Pulse is the frequency at which the DAS launches pulsed light into the optical fiber.

[0037] The gauge length is given as follows depending on the DAS settings: TIFF0007758178000006.tif523 where g is an integer value. It should be noted that the smaller the gauge length, the higher the spatial resolution that can be measured, and therefore the smaller the gauge length, the more preferable it is.

[0038] From the above, the phase difference signal is expressed as follows: TIFF0007758178000007.tif5125 where, TIFF0007758178000008.tif513 is a vector representing time series data at a certain DAS coordinate.

[0039] Next, an example of the flow of operations of the vibration characteristic extraction unit 11 will be described with reference to FIG. 7, first, the vibration characteristic extraction unit 11 extracts an input signal in a certain time interval (step S11). The input signal in a certain time interval is the above-mentioned phase difference signal obtained from the DAS. This refers to data extracted for a desired time interval from TIFF0007758178000009.tif517. If the extracted data has N components in the time direction, the phase difference signal can be expressed as an N-dimensional vector as follows: TIFF0007758178000010.tif513

[0040] Next, the vibration characteristic extraction unit 11 calculates a power spectrum for the input signal extracted in step S11 (step S12). Specifically, the vibration characteristic extraction unit 11 calculates a power spectrum for the input signal extracted in step S11 (step S12). Perform a Fourier transform on TIFF0007758178000011.tif513 and obtain the absolute value of the resulting Fourier component (power spectrum) Calculate TIFF0007758178000012.tif613.

[0041] Thereafter, the vibration characteristic extraction unit 11 extracts a power spectrum of a certain frequency band from the power spectrum calculated in step S12, and outputs the extracted power spectrum to the dissimilarity calculation unit 12 (step S13). The certain frequency band is a frequency band that includes the fundamental vibration mode in the aerial section of the optical fiber. For example, if a maximum peak value of the power spectrum is observed around 40 Hz, the vibration characteristic extraction unit 11 extracts a power spectrum of a frequency band of 30-50 Hz.

[0042] Next, an example of the flow of operations of the dissimilarity calculation unit 12 will be described with reference to FIG. 8, first, the dissimilarity calculation unit 12 determines an evaluation interval for evaluating the dissimilarity (step S21). The evaluation interval for evaluating the dissimilarity is an interval for evaluating the dissimilarity of the power spectrum between two pieces of data.

[0043] Thereafter, the dissimilarity calculation unit 12 calculates the dissimilarity between the two power spectra between the evaluation intervals determined in step S21 for all DAS coordinate points based on the power spectra obtained by the vibration characteristic extraction unit 11, and outputs the calculated dissimilarity to the pillar position calculation unit 13 (step S22).

[0044] Next, an example of the flow of operations of the pillar position calculation unit 13 will be described with reference to FIG. 9, first, the pillar position calculation unit 13 detects the maximum peak value of the dissimilarity based on the dissimilarity at all DAS coordinate points obtained by the dissimilarity calculation unit 12 (step S31). Note that any method may be used to detect the maximum peak value. Furthermore, the number of maximum peak values ​​to be detected may be determined, for example, depending on the number of pillars whose positions are to be estimated.

[0045] Thereafter, the pillar position calculation unit 13 estimates the position of the maximum peak value detected in step S31 as the position of the pillar, extracts the DAS coordinate value corresponding to the estimated position, and outputs an output signal representing the extracted DAS coordinate value (step S32).

[0046] Next, an example of a method for analyzing the degree of difference by the pillar position calculation unit 13 will be described with reference to FIG. Here, the dissimilarity calculation unit 12 calculates the power spectrum TIFF0007758178000013.tif613 and power spectrum The dissimilarity D(p) between TIFF0007758178000014.tif622 and the image is calculated as follows: TIFF0007758178000015.tif1480Here, a is the evaluation interval for the dissimilarity.

[0047] for example, TIFF0007758178000016.tif613 and If both TIFF0007758178000017.tif622 and are power spectra of the fictitious part and do not cross a pillar, the power spectra of both are similar, and D(p) takes a value close to 0.

[0048] on the other hand, TIFF0007758178000018.tif613 and When TIFF0007758178000019.tif622 and TIFF0007758178000019.tif622 cross the imaginary part and the pillar, the similarity of the power spectra of the two becomes small, and D(p) takes a value close to 1.

[0049] Therefore, the pillar position calculation unit 13 analyzes the degree of difference from the above viewpoints. FIG. 10 shows an example in which the evaluation interval a is 1. In the example of FIG. 10, the pillar position calculation unit 13 estimates the position of p=N where D(p) is closest to 1, that is, the position of p=N where the dissimilarity takes the maximum peak value, as the position of the pillar, and outputs the DAS coordinate value corresponding to the estimated position.

[0050] Next, a specific example of the operation of the position evaluation device 10 according to the first embodiment will be described. In this specific example, as shown in FIG. 11, a sensing system is assumed in which an optical fiber of approximately 30 m is suspended from a pillar and a DAS is connected to the left end of the optical fiber. Under this assumption, in this specific example, the left part of the overhead section of the optical cable is set to 0 m in DAS coordinate value (length of optical fiber), and the DAS coordinate value of the pole position is estimated.

[0051] First, the vibration characteristic extraction unit 11 extracts the power spectrum of the frequency band from 0 to 50 Hz for the input signal for 10 seconds. unit Set =0.82m and g=2. Next, the dissimilarity calculation unit 12 calculates the dissimilarity using these power spectra. Here, a=1 is set.

[0052] As a result, the dissimilarity shown in FIG. 12 is obtained. In this case, the pillar position calculation unit 13 estimates the position where the dissimilarity takes the maximum peak as the position of the pillar, and outputs a DAS coordinate value of 15.2 m corresponding to the estimated position.

[0053] As described above, according to the first embodiment, the vibration characteristic extraction unit 11 receives a phase difference signal of backscattered light as an input signal, extracts the input signal in a certain time interval, calculates a power spectrum for the extracted input signal, and extracts a power spectrum in a certain frequency band from the calculated power spectrum. The dissimilarity calculation unit 12 calculates the dissimilarity between two power spectra between evaluation intervals for all DAS coordinate points. The pillar position calculation unit 13 estimates the position where the dissimilarity takes a maximum peak value as the position of the pillar, and extracts and outputs the DAS coordinate value corresponding to the estimated position.

[0054] Therefore, even if the point is not the point where the strongest vibration occurs, as long as the difference in the power spectrum between evaluation intervals is maximized, the point can be estimated as the position of the pillar. This makes it possible to improve the accuracy of estimating the position of the pillar, which is the position where the optical fiber environment changes.

[0055] <Embodiment 2> In the second embodiment, the position of an extra length section of optical fiber is estimated and output as a DAS coordinate value, as shown in Fig. 13. Specifically, the left and right ends of the extra length of optical fiber are expressed as the extra length left and the extra length right, respectively, and the DAS coordinate values ​​of the positions of the extra length left and the extra length right are estimated and output.

[0056] An example of the configuration of the position evaluation device 20 according to the second embodiment will be described with reference to FIG. As shown in FIG. 14, a position evaluation device 20 according to the second embodiment includes a vibration characteristic extraction unit 21, a weighted dissimilarity calculation unit 22, and a marginal length calculation unit .

[0057] The input signal input to the vibration characteristic extraction unit 21 is the same as the input signal according to the first embodiment described above. The vibration characteristic extraction unit 21 is similar to the vibration characteristic extraction unit 11 according to the first embodiment described above.

[0058] The weighted dissimilarity calculation unit 22 calculates the frequency average value of the power spectrum at all DAS coordinate points based on the power spectrum obtained by the vibration characteristic extraction unit 21. The weighted dissimilarity calculation unit 22 also determines an evaluation interval for evaluating the dissimilarity. The weighted dissimilarity calculation unit 22 also calculates the dissimilarity between two power spectra between evaluation intervals at all DAS coordinate points based on the power spectrum obtained by the vibration characteristic extraction unit 21. The weighted dissimilarity calculation unit 22 also calculates and outputs weighted dissimilarity at all DAS coordinate points, in which the dissimilarity is weighted by the frequency average value of the power spectrum, based on the obtained dissimilarity at all DAS coordinate points and the frequency average value of the power spectrum. The weighted dissimilarity calculation unit 22 also outputs the frequency average value of the power spectrum at all DAS coordinate points.

[0059] The excess section calculation unit 23 sets the excess section range on the DAS coordinates based on the average frequency values ​​of the power spectrum at all DAS coordinate points obtained by the weighted dissimilarity calculation unit 22. Furthermore, based on the weighted dissimilarity at all DAS coordinate points obtained by the weighted dissimilarity calculation unit 22, the excess section calculation unit 23 estimates positions within the set excess section range where the weighted dissimilarity takes a maximum peak value as the left and right excess section positions, extracts DAS coordinate values ​​corresponding to the estimated positions, and outputs output signals representing the extracted DAS coordinate values.

[0060] The position evaluation device 20 according to the second embodiment will be described in more detail below. First, an example of the flow of operations of the weighted dissimilarity calculation section 22 will be described with reference to FIG. As shown in FIG. 15, first, the weighted dissimilarity calculation unit 22 calculates the average frequency value of the power spectrum at all DAS coordinate points based on the power spectrum obtained by the vibration characteristic extraction unit 21 (step S41).

[0061] Next, the weighted dissimilarity calculation unit 22 determines an evaluation interval for evaluating the dissimilarity (step S42). Next, the weighted dissimilarity calculation unit 22 calculates the dissimilarity between the two power spectra between the evaluation intervals determined in step S42 for all DAS coordinate points based on the power spectra obtained by the vibration characteristic extraction unit 21 (step S43).

[0062] Thereafter, the weighted dissimilarity calculation unit 22 calculates weighted dissimilarity for all DAS coordinate points based on the average frequency value of the power spectrum and the dissimilarity at all DAS coordinate points obtained in steps S41 and S43, and outputs the calculated weighted dissimilarity to the extra length interval calculation unit 23. The weighted dissimilarity is obtained by weighting the dissimilarity by the average frequency value of the power spectrum and is obtained by multiplying the average frequency value of the power spectrum by the dissimilarity. The weighted dissimilarity can characterize the state of points where the dissimilarity is high and the vibration intensity is small. Furthermore, the weighted dissimilarity calculation unit 22 outputs the average frequency value of the power spectrum for all DAS coordinate points obtained in step S41 to the extra length interval calculation unit 23 (step S44).

[0063] Next, an example of the flow of operations of the extra length section calculation unit 23 will be described with reference to FIG. 16, first, the extra length section calculation unit 23 sets an extra length section range on the DAS coordinates based on the frequency average value of the power spectrum at all DAS coordinate points obtained by the weighted dissimilarity calculation unit 22 (step S51). Specifically, the extra length section calculation unit 23 sets the range on the DAS coordinates where the frequency average value of the power spectrum is below a threshold as the extra length section range. This sets a rough range that can be considered as an extra length section of the optical fiber. The threshold is a percentile value with respect to the frequency average value of the power spectrum.

[0064] Thereafter, the excess length section calculation unit 23 estimates, within the excess length section range set in step S51, positions where the weighted dissimilarity takes a maximum peak value as the positions of the left and right excess length of the optical fiber, based on the weighted dissimilarity at all DAS coordinate points obtained by the weighted dissimilarity calculation unit 22. Then, the excess length section calculation unit 23 extracts DAS coordinate values ​​corresponding to the positions estimated to be the left and right excess length positions, and outputs output signals representing the extracted DAS coordinate values ​​(step S52).

[0065] Next, a specific example of the operation of the position evaluation device 20 according to the second embodiment will be described. In this specific example, we assume a sensing system in which three pillars contain approximately 50 m of extra optical fiber, approximately 30 m of optical fiber is suspended between each of the three pillars, and a DAS is connected to the left end of the optical fiber, as shown in Figure 17. Therefore, the total length of the optical fiber is approximately 210 m.

[0066] Based on this assumption, in this specific example, the left portion of the excess length of the optical fiber included in the left-hand pillar of the three pillars mentioned above is set to 0 m in DAS coordinate value (length of optical fiber), and the DAS coordinate values ​​of the left and right positions of the excess length of the optical fiber included in the central pillar (approximately 80 m to 130 m in DAS coordinate value) are estimated.

[0067] First, the vibration characteristic extraction unit 21 extracts the power spectrum of the frequency band from 0 to 50 Hz for the input signal for 10 seconds. unit Set =0.82m and g=2. Next, the weighted dissimilarity calculation unit 22 uses these power spectra to calculate the average frequency values ​​of the power spectra, dissimilarity, and weighted dissimilarity at all DAS coordinate points, where a=2 is set.

[0068] As a result, it is assumed that the frequency average values, dissimilarity, and weighted dissimilarity of the power spectrum are obtained as shown in Fig. 18 and Fig. 19. Note that Fig. 19 is an enlarged view of the X region shown in Fig. 18.

[0069] Here, the extra length section of the optical fiber is hardly affected by the wind, and therefore the frequency average value of the power spectrum of the extra length section is small. Therefore, first, the extra length section calculation unit 23 sets the range on the DAS coordinates where the frequency average value of the power spectrum is below a threshold as the extra length section range, which is a rough range that can be considered as an extra length section of the optical fiber. Here, the threshold is set as the 70th percentile.

[0070] On the other hand, at the position where the optical fiber changes from the slack section to the overhead section, the vibration applied to the optical fiber changes, and the characteristics of the backscattered light change, so the dissimilarity increases. At this time, the weighted dissimilarity, which is obtained by weighting the dissimilarity with the frequency average value of the power spectrum, shows maximum peaks at the left and right parts of the slack section of the optical fiber.

[0071] Therefore, the excess length section calculation unit 23 next estimates the positions within the excess length section range where the weighted dissimilarity takes the maximum peak value as the left and right positions of the excess length of the optical fiber, and outputs DAS coordinate values ​​of 80.36 m and 129.2 m corresponding to the estimated positions.

[0072] As described above, according to the second embodiment, the vibration characteristic extraction unit 21 receives a phase difference signal of backscattered light as an input signal, extracts the input signal in a certain time interval, calculates a power spectrum for the extracted input signal, and extracts a power spectrum in a certain frequency band from the calculated power spectrum. The weighted dissimilarity calculation unit 22 calculates the frequency average value of the power spectrum at all DAS coordinate points, calculates the dissimilarity between two power spectra between evaluation intervals, and calculates a weighted dissimilarity by weighting the dissimilarity by the frequency average value of the power spectrum. The extra length section calculation unit 23 sets an extra length section range on the DAS coordinates, estimates positions within the set extra length section range where the weighted dissimilarity takes a maximum peak value as the left and right positions of the extra length of the optical fiber, and extracts and outputs DAS coordinate values ​​corresponding to the estimated positions.

[0073] Therefore, even if the points are not the points where the strongest vibration occurs, as long as the weighted difference in the power spectrum between evaluation intervals is maximized, the points can be estimated as the left and right positions of the excess length of the optical fiber. This makes it possible to improve the accuracy of estimating the positions of the excess length sections, which are the positions of the optical fiber where the environment changes.

[0074] <Third Embodiment> In the third embodiment, as shown in FIG. 20, the positions of the pillars are estimated and output as DAS coordinate values ​​by referring to the distance between the pillars given in advance. The inter-pillar distance indicates the distance between adjacent pillars, and is assumed to be given in advance from the following information, for example: · Positioning of the distance between pillars using GPS (Global Positioning System) information and map information The optical fiber length recorded when the optical fiber was suspended between two posts

[0075] An example of the configuration of a position evaluation device 30 according to the third embodiment will be described with reference to FIG. As shown in FIG. 21, a position evaluation device 30 according to the third embodiment includes a vibration characteristics extraction unit 31, a dissimilarity calculation unit 32, and a pillar position calculation unit 33.

[0076] The input signal input to the vibration characteristic extraction unit 31 is the same as the input signal according to the first and second embodiments described above. The vibration characteristic extraction unit 31 is similar to the vibration characteristic extraction units 11 and 21 according to the first and second embodiments described above. The dissimilarity calculation section 32 is similar to the dissimilarity calculation section 12 according to the first embodiment described above.

[0077] The pillar position calculation unit 33 determines, as the analysis target section, a section that includes all pillars whose positions are to be estimated. The pillar position calculation unit 33 also configures a window function for the analysis target section based on the pre-given inter-pillar distance. The pillar position calculation unit 33 also searches for a DAS coordinate value (offset point) at which the value of the cross-correlation function between the window function and the dissimilarity takes the maximum value, based on the dissimilarity at all DAS coordinate points obtained by the dissimilarity calculation unit 32 and the window function configured above. The pillar position calculation unit 33 also outputs the position of each pillar as a DAS coordinate value, based on the pre-given inter-pillar distance and the offset point searched above. Specifically, the pillar position calculation unit 33 estimates the position of each pillar to be the position of the offset point and the position obtained by adding the offset point to the inter-pillar distance, and outputs a DAS coordinate value corresponding to the estimated position.

[0078] The position evaluation device 30 according to the third embodiment will be described in more detail below. First, an example of the flow of operations of the pillar position calculation unit 33 will be described with reference to FIG. 22, first, the pillar position calculation unit 33 determines an analysis target section (step S61). The analysis target section is a section of DAS coordinates that includes all pillars whose positions are to be estimated.

[0079] Next, the pillar position calculation unit 33 constructs a window function for the analysis target section based on the distance between pillars given in advance (step S62). The window function will be described in detail later. Next, based on the dissimilarity at all DAS coordinate points obtained by the dissimilarity calculation unit 32 and the window function configured in step S62, the pillar position calculation unit 33 searches for the DAS coordinate value (offset point) in the analysis section at which the value of the cross-correlation function between the window function and the dissimilarity is maximum (step S63).

[0080] Thereafter, the pillar position calculation unit 33 adds the inter-pillar distance to the offset point based on the pre-given inter-pillar distance and the offset point searched for in step S63, thereby outputting the position of each pillar as a DAS coordinate value (step S64). Specifically, the pillar position calculation unit 33 estimates the position of each pillar as the position of the offset point and the position obtained by adding the inter-pillar distance to the offset point, and outputs DAS coordinate values ​​corresponding to the estimated positions.

[0081] Next, an example of an analysis target section determined by the pillar position calculation unit 33 will be described with reference to FIG. 23, the pillar position calculation unit 33 sets a label for each pillar and determines the analysis target section. As described above, the analysis target section is a section of DAS coordinates that includes all pillars whose positions are to be estimated (here, pillars 0 to N).

[0082] Here, the distances between the columns that correspond to the labels of each column are defined as d1, d2, ... In addition, the DAS coordinate value from the DAS to the left end of the analysis section (point p=0) is set to d0. This DAS coordinate value is set to an arbitrary value so as to include column 0.

[0083] Next, an example of the window function configured by the pillar position calculation unit 33 will be described with reference to FIG. As shown in FIG. 24, the window function is constructed from inter-column distances p1, p2, . . . as follows: TIFF0007758178000020.tif16137 Also, the window width p l can be set appropriately as a value of about a. In addition, in order to improve the accuracy of the analysis results, the window width p l The analysis may be performed while changing the value of

[0084] Next, an example of the cross-correlation function calculated by the pillar position calculation unit 33 will be described. In the analysis interval, the cross-correlation function between the window function W(p) and the dissimilarity D(p) is given by: TIFF0007758178000021.tif1781Here, the dissimilarity D(p) becomes larger near the position of the pillar. Therefore, the pillar position calculation unit 33 calculates the DAS coordinate value (offset point) p where the value of the cross-correlation function is the maximum value in the analysis section. offset Explore. The pillar position calculation unit 33 calculates the offset point p offset After determining the coordinates, the following DAS coordinate values ​​representing the position of each pillar N' are output. TIFF0007758178000022.tif1654

[0085] FIG. 25 shows the calculation of the appropriate offset point p offset 25 shows an example of the relationship between the window function and the dissimilarity when the offset point p offset is the DAS coordinate value of column 0, and the distance between columns given in advance is the offset point p offset The values ​​added together are the DAS coordinate values ​​of pillar 1 to pillar N, respectively.

[0086] As described above, according to the third embodiment, the vibration characteristic extraction unit 31 receives a phase difference signal of backscattered light as an input signal, extracts the input signal in a certain time interval, calculates a power spectrum for the extracted input signal, and extracts a power spectrum in a certain frequency band from the calculated power spectrum. The dissimilarity calculation unit 32 calculates the dissimilarity between two power spectra between evaluation intervals at all DAS coordinate points. The pillar position calculation unit 33 determines the analysis interval and configures a window function for the analysis interval based on a predetermined pillar-to-pillar distance. The pillar position calculation unit 33 also searches for a DAS coordinate value (offset point) at which the value of the cross-correlation function between the window function and the dissimilarity is maximized, and outputs a DAS coordinate value corresponding to the position of each pillar by adding the predetermined pillar-to-pillar distance to the offset point. That is, the pillar position calculation unit 33 estimates the position of each pillar as the position of the offset point and the position obtained by adding the offset point to the pillar-to-pillar distance, and outputs a DAS coordinate value corresponding to the estimated position.

[0087] Therefore, it is possible to estimate the position of each pillar even at points other than the point where the strongest vibration occurs. This makes it possible to improve the accuracy of estimating the position of each pillar, which is the position where the optical fiber's environment changes. It also makes it possible to estimate the position of each pillar simultaneously.

[0088] <Fourth Embodiment> In the fourth embodiment, as shown in FIG. 26, when the pillars containing the extra length section of the optical fiber are known, the position of each pillar and the position of the extra length section are estimated and output as DAS coordinate values ​​by referring to the distance between pillars given in advance.

[0089] An example of the configuration of a position evaluation device 40 according to the fourth embodiment will be described with reference to FIG. As shown in FIG. 27, the position evaluation device 40 according to the fourth embodiment includes a vibration characteristic extraction unit 41, a weighted dissimilarity calculation unit 42, an excess length section calculation unit 43, and a pillar position calculation unit 44.

[0090] The input signal input to the vibration characteristic extraction unit 41 is the same as the input signal according to the first, second and third embodiments described above. The vibration characteristic extraction unit 41 is similar to the vibration characteristic extraction units 11, 21, and 31 according to the first, second, and third embodiments described above. The weighted dissimilarity calculation unit 42 is similar to the weighted dissimilarity calculation unit 22 according to the second embodiment described above.

[0091] The extra length section calculation unit 43 sets an extra length section range on the DAS coordinates based on the average frequency values ​​of the power spectrum at all the DAS coordinate points obtained by the weighted dissimilarity calculation unit 42.

[0092] The pillar position calculation unit 44 determines the section from a pillar including an extra length section to the next pillar including an extra length section as the analysis target section. The pillar position calculation unit 44 also configures a window function for the analysis target section based on a pre-given inter-pillar distance. The pillar position calculation unit 44 also searches for a DAS coordinate value (offset point) at which the value of the cross-correlation function between the window function and the weighted dissimilarity is maximized, based on the weighted dissimilarity at all DAS coordinate points obtained by the weighted dissimilarity calculation unit 42 and the window function configured above. The pillar position calculation unit 44 also outputs DAS coordinate values ​​corresponding to the position of each pillar based on the pre-given inter-pillar distance and the offset point searched above. That is, the pillar position calculation unit 44 estimates the position of each pillar as the position of the offset point and the position obtained by adding the offset point to the inter-pillar distance, and outputs DAS coordinate values ​​corresponding to the estimated positions.

[0093] Furthermore, the pole position calculation unit 44 outputs a DAS coordinate value corresponding to the position to the right of the excess length of optical fiber, based on the offset point searched for as described above. That is, the pole position calculation unit 44 estimates the position of the offset point as the position to the right of the excess length, and outputs a DAS coordinate value corresponding to the estimated position. Also, the pole position calculation unit 44 outputs a DAS coordinate value corresponding to the position to the left of the excess length of optical fiber, based on the weighted dissimilarity of each position on the DAS side of the position of the pole including the excess length section. That is, the pole position calculation unit 44 estimates the position where the weighted dissimilarity takes a maximum peak value on the DAS side of the position of the pole including the excess length section as the position to the left of the excess length, and outputs a DAS coordinate value corresponding to the estimated position.

[0094] The position evaluation device 40 according to the fourth embodiment will be described in more detail below. First, an example of the flow of operations of the extra length section calculation unit 43 will be described with reference to FIG. 28, the extra length section calculation unit 43 sets an extra length section range on the DAS coordinates based on the frequency average value of the power spectrum at all DAS coordinate points obtained by the weighted dissimilarity calculation unit 42 (step S71). Specifically, the extra length section calculation unit 43 sets the range on the DAS coordinates where the frequency average value of the power spectrum is below the threshold as the extra length section range. This sets an extra length section range, which is a rough range that can be considered as an extra length section of optical fiber. Note that the threshold is a percentile value with respect to the frequency average value of the power spectrum. Furthermore, the number of extra length section ranges is equal to the number of columns that include extra length sections.

[0095] Next, an example of the flow of operations of the pillar position calculation unit 44 will be described with reference to FIG. 29, first, the column position calculation unit 44 determines an analysis target section (step S81). The analysis target section is a section of DAS coordinates from a column including an extra length section to the next column including an extra length section.

[0096] Next, the pillar position calculation unit 44 configures a window function for the analysis target section based on the distance between pillars given in advance (step S82). The window function is the same as the window function according to the third embodiment described above.

[0097] Next, the pillar position calculation unit 44 searches for the DAS coordinate value (offset point) in the analysis section where the value of the cross-correlation function between the window function and the weighted dissimilarity is maximum, based on the weighted dissimilarity at all DAS coordinate points obtained by the weighted dissimilarity calculation unit 42 and the window function configured in step S82 (step S83).

[0098] Next, the pillar position calculation unit 44 adds the offset points to the inter-pillar distance based on the pre-given inter-pillar distance and the offset points searched for in step S83, and outputs DAS coordinate values ​​corresponding to the position of each pillar (step S84). That is, the pillar position calculation unit 44 estimates the positions of the offset points and the positions obtained by adding the inter-pillar distance to the offset points as the positions of each pillar, and outputs DAS coordinate values ​​corresponding to the estimated positions.

[0099] Next, the pole position calculation unit 44 outputs the DAS coordinate value of the offset point found in step S83 as the DAS coordinate value corresponding to the position to the right of the excess length of the optical fiber (step S85). That is, the pole position calculation unit 44 estimates the position of the offset point as the position to the right of the excess length, and outputs the DAS coordinate value corresponding to the estimated position.

[0100] Thereafter, the pole position calculation unit 44 outputs the DAS coordinate value corresponding to the position where the weighted dissimilarity takes the maximum peak value based on the weighted dissimilarity of each position on the DAS side of the pole position including the extra length section, as the DAS coordinate value corresponding to the position to the left of the extra length of the optical fiber (step S86). That is, the pole position calculation unit 44 estimates the position where the weighted dissimilarity takes the maximum peak value on the DAS side of the pole position including the extra length section as the position to the left of the extra length, and outputs the DAS coordinate value corresponding to the estimated position.

[0101] Next, an example of an analysis target section determined by the pillar position calculation unit 44 will be described with reference to FIG. 30, the pillar position calculation unit 44 sets a label for each pillar and determines the analysis target section. As described above, the analysis target section is the section of DAS coordinates from a pillar including an extra length section to the next pillar including an extra length section.

[0102] 30, column 0 includes an excess length section, and the next column that includes an excess length section is column N+1. Therefore, the column position calculation unit 44 determines the section from column 0 to column N, which is immediately before column N+1, as analysis section 1, and determines the section from column N+1 onwards as analysis section 2.

[0103] Here, the distances between the columns that correspond to the labels of each column are defined as d1, d2, ... In addition, the DAS coordinate value from the DAS to the left end of the analysis section (point p=0) is set to d0. This DAS coordinate value is set to an arbitrary value so as to include column 0.

[0104] Next, an example of the cross-correlation between the weighted dissimilarity and the window function will be described with reference to FIG. As described above, the weighted dissimilarity is obtained by weighting the dissimilarity by the average frequency value of the power spectrum, and is obtained by multiplying the average frequency value of the power spectrum by the dissimilarity.

[0105] As shown in Fig. 31, when the window function is configured based on the right excess length, the value of the cross-correlation function between the weighted dissimilarity and the window function is larger than when the window function is configured based on the left excess length. Therefore, the column position calculation unit 44 configures the window function based on the right excess length. Then, the column position calculation unit 44 calculates the DAS coordinate value (offset point) p where the value of the cross-correlation function is maximum. offset The DAS coordinate value of the offset point corresponds to the position to the right of the excess length of the optical fiber. Thereafter, the pole position calculation unit 44 estimates the DAS coordinate value of the position of each pole, in the same way as in the third embodiment described above.

[0106] Next, an example of a method for calculating the left extra length by the column position calculation unit 44 will be described with reference to FIG. As shown in Figure 32, the pole position calculation unit 44 estimates the DAS coordinate value at which the weighted dissimilarity takes the maximum peak value among the DAS coordinate values ​​of each position on the DAS side of the pole position including the excess length section as the DAS coordinate value corresponding to the position to the left of the excess length of the optical fiber.

[0107] As mentioned above, the offset value p offset The DAS coordinate value corresponds to the position to the right of the excess length of the optical fiber. Therefore, in the example of FIG. 32, the pillar position calculation unit 44 calculates 0 <p<p offset Among the DAS coordinate values, the DAS coordinate value with the maximum weighted dissimilarity is estimated as the DAS coordinate value corresponding to the left position of the excess length of the optical fiber.

[0108] As described above, according to the fourth embodiment, the vibration characteristic extraction unit 41 receives a phase difference signal of backscattered light as an input signal, extracts the input signal in a certain time interval, calculates a power spectrum for the extracted input signal, and extracts a power spectrum in a certain frequency band from the calculated power spectrum. The weighted dissimilarity calculation unit 42 calculates the frequency average value of the power spectrum at all DAS coordinate points, calculates the dissimilarity between two power spectra between evaluation intervals, and calculates a weighted dissimilarity by weighting the dissimilarity by the frequency average value of the power spectrum. The marginal interval calculation unit 43 sets the marginal interval range on the DAS coordinate. The pillar position calculation unit 44 determines the analysis target interval and configures a window function for the analysis target interval based on a predetermined pillar-to-pillar distance. The pillar position calculation unit 44 also searches for a DAS coordinate value (offset point) at which the value of the cross-correlation function between the window function and the weighted dissimilarity takes the maximum value, and outputs a DAS coordinate value corresponding to the position of each pillar by adding the predetermined pillar-to-pillar distance to the offset point. The pole position calculation unit 44 also outputs the offset point as the DAS coordinate value of the right side of the excess length of optical fiber. The pole position calculation unit 44 also outputs, among the DAS coordinate values ​​of each position on the DAS side of the right side of the excess length, the DAS coordinate value at which the weighted dissimilarity takes a maximum peak value as the DAS coordinate value of the position on the left side of the excess length of optical fiber. That is, the pole position calculation unit 44 estimates the position of the offset point and the position obtained by adding the offset point to the distance between poles as the position of each pole, estimates the position of the offset point as the position of the right side of the excess length, and estimates the position at which the weighted dissimilarity takes a maximum peak value on the DAS side of the position of the pole including the excess length section as the position of the left side of the excess length. The pole position calculation unit 44 then outputs DAS coordinate values ​​corresponding to the estimated positions of each pole, the right side of the excess length, and the left side of the excess length.

[0109] Therefore, points other than the point where the strongest vibration occurs can be estimated as the position of each pillar or the position of the extra length of optical fiber. This makes it possible to improve the accuracy of estimating the position of each pillar or the position of the extra length of optical fiber, which are the positions where the optical fiber's environment changes. It also makes it possible to estimate the position of each pillar simultaneously.

[0110] <Fifth Embodiment> The fifth embodiment corresponds to an embodiment that is a superordinate concept of the first to fourth embodiments described above. An example of the configuration of a position evaluation device 50 according to the fifth embodiment will be described with reference to FIG. As shown in FIG. 33, a position evaluation device 50 according to the fifth embodiment includes a vibration characteristics calculation unit 51, a dissimilarity calculation unit 52, and an environment change position estimation unit 53.

[0111] The vibration characteristic calculation unit 51 receives a signal indicating the natural vibration occurring at each position of the optical fiber from the sensor, and calculates sensing data indicating the vibration characteristic at each position of the optical fiber based on the received signal. The vibration characteristic calculation unit 51 corresponds to the vibration characteristic extraction units 11, 21, 31, and 41 according to the above-described first, second, third, and fourth embodiments. The sensor corresponds to a phase-sensing OTDR or a DAS.

[0112] The dissimilarity calculation unit 52 calculates the dissimilarity of sensing data between two adjacent points on the optical fiber. The dissimilarity calculation unit 52 corresponds to the dissimilarity calculation units 12 and 32 according to the first and third embodiments and the weighted dissimilarity calculation units 22 and 42 according to the second and fourth embodiments.

[0113] The environment change position estimation unit 53 estimates the environment change position where the environment of the optical fiber changes based on the degree of difference. The environment change position estimation unit 53 corresponds to the pole position calculation units 13, 33, 44 according to the above-mentioned first, third, and fourth embodiments and the extra length section calculation units 23, 43 according to the above-mentioned second and fourth embodiments.

[0114] Since the fifth embodiment is configured as described above, it is possible to estimate a point other than the point where the maximum vibration occurs as the environmental change position where the environment of the optical fiber changes, thereby improving the accuracy of estimating the environmental change position of the optical fiber.

[0115] The optical fiber may be an optical fiber suspended above a pole. In this case, the vibration characteristic calculation unit 51 may calculate a power spectrum in a predetermined frequency band as sensing data. Furthermore, the dissimilarity calculation unit 52 may calculate the dissimilarity of the power spectrum between two adjacent points on the optical fiber. Furthermore, the environment change position estimation unit 53 may estimate the position where the dissimilarity takes a maximum peak value as the position of the pillar suspending the optical fiber.

[0116] Alternatively, the vibration characteristic calculation unit 51 may calculate a power spectrum in a predetermined frequency band as sensing data. Furthermore, the dissimilarity calculation unit 52 may calculate the dissimilarity of the power spectrum between two adjacent points on the optical fiber, and may also calculate a weighted dissimilarity by weighting the calculated dissimilarity by the frequency average value of the power spectrum. Furthermore, the environmental change position estimation unit 53 may set the range in which the frequency average value of the power spectrum is below a threshold as the range of the extra length section of the optical fiber, and estimate the positions within the range of the extra length section where the weighted dissimilarity takes a maximum peak value as the left and right ends of the extra length section.

[0117] Alternatively, the vibration characteristic calculation unit 51 may calculate a power spectrum in a predetermined frequency band as sensing data. Furthermore, the dissimilarity calculation unit 52 may calculate the dissimilarity of power spectra between two adjacent points on the optical fiber. Furthermore, the environmental change position estimation unit 53 may determine an analysis target section that includes all pillars whose positions are to be estimated, configure a window function based on a pre-given inter-pillar distance indicating the distance between adjacent pillars in the analysis target section, search for an offset point where the correlation function between the window function and the dissimilarity takes a maximum value, and estimate the position of the offset point and the position obtained by adding the offset point to the inter-pillar distance as the position of the pillar suspending the optical fiber.

[0118] Alternatively, the vibration characteristic calculation unit 51 may calculate a power spectrum in a predetermined frequency band as the sensing data. The dissimilarity calculation unit 52 may calculate the dissimilarity of the power spectrum between two adjacent points on the optical fiber and calculate a weighted dissimilarity by weighting the calculated dissimilarity by the frequency average value of the power spectrum. The environmental change position estimation unit 53 may determine an analysis target section from a pillar including an extra length section of the optical fiber to the next pillar including the extra length section, configure a window function based on a predetermined pillar-to-pillar distance indicating the distance between adjacent pillars in the analysis target section, search for an offset point where the correlation function between the window function and the weighted dissimilarity takes a maximum value, estimate the position of the offset point and a position obtained by adding the offset point to the pillar-to-pillar distance as the positions of the pillars suspending the optical fiber, and estimate the position of the offset point and positions on the sensor side of the offset point where the weighted dissimilarity takes a maximum peak value as the left and right ends of the extra length section.

[0119] In this way, in the case of an optical fiber suspended from a pole, the environmental change position of the optical fiber estimated in the fifth embodiment is, for example, the position of the pole suspending the optical fiber and the extra length section of the optical fiber.

[0120] However, the environmental change position of the optical fiber estimated in this fifth embodiment is not limited to this. This fifth embodiment may also be used to estimate the position of a boundary point where physical properties change significantly (for example, a fixed point of the optical fiber) as the environmental change position. For example, an optical fiber included in an optical submarine cable has a boundary point between a section buried in the seabed and a section exposed in the sea and swaying due to waves. This fifth embodiment may also be used to estimate the position of such a boundary point as the environmental change position.

[0121] <Hardware configuration of the position evaluation device according to the embodiment> An example of the hardware configuration of the computer 90 that realizes the position evaluation devices 10, 20, 30, 40, and 50 according to the first, second, third, fourth, and fifth embodiments will be described with reference to FIG.

[0122] 34, a computer 90 includes a processor 91, a memory 92, a storage 93, an input / output interface (input / output I / F) 94, and a communication interface (communication I / F) 95. The processor 91, the memory 92, the storage 93, the input / output interface 94, and the communication interface 95 are connected by a data transmission path for transmitting and receiving data to and from each other.

[0123] The processor 91 is, for example, an arithmetic processing device such as a CPU (Central Processing Unit) or a GPU (Graphics Processing Unit). The memory 92 is, for example, a memory such as a RAM (Random Access Memory) or a ROM (Read Only Memory). The storage 93 is, for example, a storage device such as an HDD (Hard Disk Drive), an SSD (Solid State Drive), or a memory card. The storage 93 may also be a memory such as a RAM or a ROM.

[0124] A program is stored in the storage 93. When the program is loaded into a computer, it includes a set of instructions (or software code) that causes the computer 90 to perform one or more functions of the above-described position evaluation devices 10, 20, 30, 40, and 50. The components of the above-described position evaluation devices 10, 20, 30, 40, and 50 may be realized by the processor 91 reading and executing the program stored in the storage 93. Furthermore, the storage function of the above-described position evaluation devices 10, 20, 30, 40, and 50 may be realized by the memory 92 or the storage 93.

[0125] The above-described programs may also be stored on non-transitory computer-readable media or tangible storage media. By way of example and not limitation, computer-readable media or tangible storage media include RAM, ROM, flash memory, SSD or other memory technology, CD (Compact Disc)-ROM, DVD (Digital Versatile Disc), Blu-ray® disc or other optical disk storage, magnetic cassette, magnetic tape, magnetic disk storage or other magnetic storage devices. The programs may also be transmitted on a transitory computer-readable medium or communication medium. By way of example and not limitation, transitory computer-readable media or communication media include electrical, optical, acoustic, or other forms of propagated signals.

[0126] The input / output interface 94 is connected to a display device 941, an input device 942, a sound output device 943, etc. The display device 941 is a device that displays a screen corresponding to drawing data processed by the processor 91, such as an LCD (Liquid Crystal Display), a CRT (Cathode Ray Tube) display, or a monitor. The input device 942 is a device that accepts operational inputs from an operator, such as a keyboard, a mouse, or a touch sensor. The display device 941 and the input device 942 may be integrated and realized as a touch panel. The sound output device 943 is a device that outputs sound corresponding to audio data processed by the processor 91, such as a speaker.

[0127] The communication interface 95 transmits and receives data to and from an external device. For example, the communication interface 95 communicates with the external device via a wired communication path or a wireless communication path.

[0128] Although the present disclosure has been described above with reference to the embodiments, the present disclosure is not limited to the above-described embodiments. Various modifications that can be understood by those skilled in the art can be made to the configuration and details of the present disclosure within the scope of the present disclosure. For example, some or all of the above-described embodiments may be used in combination with each other.

[0129] Furthermore, some or all of the above-described embodiments can be described as, but are not limited to, the following supplementary notes. (Appendix 1) a vibration characteristic calculation unit that receives a signal indicating a natural vibration generated at each position of the optical fiber from the sensor and calculates sensing data indicating the vibration characteristic at each position of the optical fiber based on the received signal; a difference calculation unit that calculates a difference in the sensing data between two adjacent points on the optical fiber; and an environment change position estimation unit that estimates an environment change position where the environment of the optical fiber changes based on the degree of difference. (Appendix 2) The optical fiber is a pillar-mounted optical fiber, the vibration characteristic calculation unit calculates a power spectrum in a predetermined frequency band as the sensing data, the dissimilarity calculation unit calculates a dissimilarity of the power spectrum between two adjacent points on the optical fiber; the environment change position estimation unit estimates a position where the degree of difference takes a maximum peak value as the position of a pillar suspending the optical fiber; 2. The position evaluation device of claim 1. (Appendix 3) The optical fiber is a pillar-mounted optical fiber, the vibration characteristic calculation unit calculates a power spectrum in a predetermined frequency band as the sensing data, the dissimilarity calculation unit calculates a dissimilarity of the power spectrum between two adjacent points of the optical fiber, and calculates a weighted dissimilarity by weighting the calculated dissimilarity with a frequency average value of the power spectrum; the environmental change position estimating unit sets a range in which the frequency average value of the power spectrum is below a threshold as a range of the extra length section of the optical fiber, and estimates positions within the range of the extra length section where the weighted dissimilarity takes a maximum peak value as the left and right ends of the extra length section. 2. The position evaluation device of claim 1. (Appendix 4) The optical fiber is a pillar-mounted optical fiber, the vibration characteristic calculation unit calculates a power spectrum in a predetermined frequency band as the sensing data, the dissimilarity calculation unit calculates a dissimilarity of the power spectrum between two adjacent points on the optical fiber; The environment change position estimation unit The section including all the pillars whose positions are to be estimated is determined as the analysis section. A window function is configured based on a predetermined column distance indicating a distance between adjacent columns in the analysis section; searching for an offset point in the analysis interval where a correlation function between the window function and the dissimilarity takes a maximum value; In the analysis section, the position of the offset point and the position obtained by adding the distance between the pillars to the offset point are estimated as the positions of the pillars suspending the optical fiber. 2. The position evaluation device of claim 1. (Appendix 5) The optical fiber is a pillar-mounted optical fiber, the vibration characteristic calculation unit calculates a power spectrum in a predetermined frequency band as the sensing data, the dissimilarity calculation unit calculates a dissimilarity of the power spectrum between two adjacent points of the optical fiber, and calculates a weighted dissimilarity by weighting the calculated dissimilarity with a frequency average value of the power spectrum; The environment change position estimation unit determining a section from a pole including the extra length section of the optical fiber to a next pole including the extra length section as an analysis section; A window function is configured based on a predetermined column distance indicating a distance between adjacent columns in the analysis section; searching for an offset point in the analysis interval where a correlation function between the window function and the weighted dissimilarity takes a maximum value; In the analysis section, the position of the offset point and the position obtained by adding the distance between the pillars to the offset point are estimated as the positions of the pillars suspending the optical fiber; In the analysis target section, the position of the offset point and a position on the sensor side of the offset point where the weighted difference takes a maximum peak value are estimated to be the left and right ends of the extra length section. 2. The position evaluation device of claim 1. (Appendix 6) the predetermined frequency band is a frequency band including a fundamental vibration mode in the aerial section of the optical fiber; 6. A position evaluation device according to any one of appendices 2 to 5. (Appendix 7) A location estimation method executed by a location estimation device, comprising: a vibration characteristics calculation step of receiving a signal indicating a natural vibration occurring at each position of the optical fiber from a sensor and calculating sensing data indicating the vibration characteristics at each position of the optical fiber based on the received signal; a difference calculation step of calculating a difference in the sensing data between two adjacent points of the optical fiber; and an environment change position estimating step of estimating an environment change position where the environment of the optical fiber changes based on the degree of difference. (Appendix 8) The optical fiber is a pillar-mounted optical fiber, In the vibration characteristic calculation step, a power spectrum in a predetermined frequency band is calculated as the sensing data; In the dissimilarity calculation step, a dissimilarity of the power spectrum between two adjacent points on the optical fiber is calculated; In the environmental change position estimating step, a position where the difference degree takes a maximum peak value is estimated to be the position of a pillar suspending the optical fiber. 7. The location evaluation method described in Appendix 7. (Appendix 9) The optical fiber is a pillar-mounted optical fiber, In the vibration characteristic calculation step, a power spectrum in a predetermined frequency band is calculated as the sensing data; In the dissimilarity calculation step, a dissimilarity of the power spectrum between two adjacent points of the optical fiber is calculated, and a weighted dissimilarity is calculated by weighting the calculated dissimilarity by a frequency average value of the power spectrum; In the environmental change position estimating step, a range in which the frequency average value of the power spectrum is below a threshold is set as a range of the extra length section of the optical fiber, and positions within the range of the extra length section where the weighted dissimilarity takes a maximum peak value are estimated to be the left and right ends of the extra length section. 7. The location evaluation method described in Appendix 7. (Appendix 10) The optical fiber is a pillar-mounted optical fiber, In the vibration characteristic calculation step, a power spectrum in a predetermined frequency band is calculated as the sensing data; In the dissimilarity calculation step, a dissimilarity of the power spectrum between two adjacent points on the optical fiber is calculated; In the environment change position estimation step, The section including all the pillars whose positions are to be estimated is determined as the analysis section. A window function is configured based on a predetermined column distance indicating a distance between adjacent columns in the analysis section; searching for an offset point in the analysis interval where a correlation function between the window function and the dissimilarity takes a maximum value; In the analysis section, the position of the offset point and the position obtained by adding the distance between the pillars to the offset point are estimated as the positions of the pillars suspending the optical fiber. 7. The location evaluation method described in Appendix 7. (Appendix 11) The optical fiber is a pillar-mounted optical fiber, In the vibration characteristic calculation step, a power spectrum in a predetermined frequency band is calculated as the sensing data; In the dissimilarity calculation step, a dissimilarity of the power spectrum between two adjacent points of the optical fiber is calculated, and a weighted dissimilarity is calculated by weighting the calculated dissimilarity by a frequency average value of the power spectrum; In the environment change position estimation step, determining a section from a pole including the extra length section of the optical fiber to a next pole including the extra length section as an analysis section; A window function is configured based on a predetermined column distance indicating a distance between adjacent columns in the analysis section; searching for an offset point in the analysis interval where a correlation function between the window function and the weighted dissimilarity takes a maximum value; In the analysis section, the position of the offset point and the position obtained by adding the distance between the pillars to the offset point are estimated as the positions of the pillars suspending the optical fiber; In the analysis target section, the position of the offset point and a position on the sensor side of the offset point where the weighted difference takes a maximum peak value are estimated to be the left and right ends of the extra length section. 7. The location evaluation method described in Appendix 7. (Appendix 12) the predetermined frequency band is a frequency band including a fundamental vibration mode in the aerial section of the optical fiber; 12. The position evaluation method according to any one of appendices 8 to 11. (Appendix 13) A non-transitory computer-readable medium storing a program to be executed by a computer, The program a vibration characteristics calculation step of receiving a signal indicating a natural vibration occurring at each position of the optical fiber from a sensor and calculating sensing data indicating the vibration characteristics at each position of the optical fiber based on the received signal; a difference calculation step of calculating a difference in the sensing data between two adjacent points of the optical fiber; and an environmental change position estimating step of estimating an environmental change position where an environment of the optical fiber changes based on the degree of difference. (Appendix 14) The optical fiber is a pillar-mounted optical fiber, In the vibration characteristic calculation step, a power spectrum in a predetermined frequency band is calculated as the sensing data; In the dissimilarity calculation step, a dissimilarity of the power spectrum between two adjacent points on the optical fiber is calculated; In the environmental change position estimating step, a position where the difference degree takes a maximum peak value is estimated to be the position of a pillar suspending the optical fiber. 14. The computer-readable medium of claim 13. (Appendix 15) The optical fiber is a pillar-mounted optical fiber, In the vibration characteristic calculation step, a power spectrum in a predetermined frequency band is calculated as the sensing data; In the dissimilarity calculation step, a dissimilarity of the power spectrum between two adjacent points of the optical fiber is calculated, and a weighted dissimilarity is calculated by weighting the calculated dissimilarity by a frequency average value of the power spectrum; In the environmental change position estimating step, a range in which the frequency average value of the power spectrum is below a threshold is set as a range of the extra length section of the optical fiber, and positions within the range of the extra length section where the weighted dissimilarity takes a maximum peak value are estimated to be the left and right ends of the extra length section. 14. The computer-readable medium of claim 13. (Appendix 16) The optical fiber is a pillar-mounted optical fiber, In the vibration characteristic calculation step, a power spectrum in a predetermined frequency band is calculated as the sensing data; In the dissimilarity calculation step, a dissimilarity of the power spectrum between two adjacent points on the optical fiber is calculated; In the environment change position estimation step, The section including all the pillars whose positions are to be estimated is determined as the analysis section. A window function is configured based on a predetermined column distance indicating a distance between adjacent columns in the analysis section; searching for an offset point in the analysis interval where a correlation function between the window function and the dissimilarity takes a maximum value; In the analysis section, the position of the offset point and the position obtained by adding the distance between the pillars to the offset point are estimated as the positions of the pillars suspending the optical fiber. 14. The computer-readable medium of claim 13. (Appendix 17) The optical fiber is a pillar-mounted optical fiber, In the vibration characteristic calculation step, a power spectrum in a predetermined frequency band is calculated as the sensing data; In the dissimilarity calculation step, a dissimilarity of the power spectrum between two adjacent points of the optical fiber is calculated, and a weighted dissimilarity is calculated by weighting the calculated dissimilarity by a frequency average value of the power spectrum; In the environment change position estimation step, determining a section from a pole including the extra length section of the optical fiber to a next pole including the extra length section as an analysis section; A window function is configured based on a predetermined column distance indicating a distance between adjacent columns in the analysis section; searching for an offset point in the analysis interval where a correlation function between the window function and the weighted dissimilarity takes a maximum value; In the analysis section, the position of the offset point and the position obtained by adding the distance between the pillars to the offset point are estimated as the positions of the pillars suspending the optical fiber; In the analysis target section, the position of the offset point and a position on the sensor side of the offset point where the weighted difference takes a maximum peak value are estimated to be the left and right ends of the extra length section. 14. The computer-readable medium of claim 13. (Appendix 18) the predetermined frequency band is a frequency band including a fundamental vibration mode in the aerial section of the optical fiber; 18. The computer-readable medium of any one of claims 14 to 17. [Explanation of symbols]

[0130] 10, 20, 30, 40, 50 position evaluation device 11, 21, 31, 41 Vibration characteristics extraction section 12,32,52 Dissimilarity calculation unit 13,33,44 Column position calculation part 22,42 Weighted Dissimilarity Calculation Unit 23,43 Extra length calculation section 51 Vibration characteristics calculation unit 53 Environmental change position estimation unit 90 Computer 91 processors 92 memory 93 Storage 94 Input / Output Interface 941 Display device 942 Input Device 943 Sound Output Device 95 Communication Interface

Claims

1. a vibration characteristic calculation unit that receives a signal indicating a natural vibration generated at each position of the optical fiber from the sensor and calculates sensing data indicating the vibration characteristic at each position of the optical fiber based on the received signal; a difference calculation unit that calculates a difference in the sensing data between two adjacent points on the optical fiber; and an environment change position estimation unit that estimates an environment change position where the environment of the optical fiber changes based on the degree of difference.

2. The optical fiber is a pillar-mounted optical fiber, the vibration characteristic calculation unit calculates a power spectrum in a predetermined frequency band as the sensing data, the dissimilarity calculation unit calculates a dissimilarity of the power spectrum between two adjacent points on the optical fiber; the environment change position estimation unit estimates a position where the degree of difference takes a maximum peak value as the position of a pillar suspending the optical fiber; The position evaluation device according to claim 1 .

3. The optical fiber is a pillar-mounted optical fiber, the vibration characteristic calculation unit calculates a power spectrum in a predetermined frequency band as the sensing data, the dissimilarity calculation unit calculates a dissimilarity of the power spectrum between two adjacent points of the optical fiber, and calculates a weighted dissimilarity by weighting the calculated dissimilarity with a frequency average value of the power spectrum; the environmental change position estimating unit sets a range in which the frequency average value of the power spectrum is below a threshold as a range of the extra length section of the optical fiber, and estimates positions within the range of the extra length section where the weighted dissimilarity takes a maximum peak value as the left and right ends of the extra length section. The position evaluation device according to claim 1 .

4. The optical fiber is a pillar-mounted optical fiber, the vibration characteristic calculation unit calculates a power spectrum in a predetermined frequency band as the sensing data, the dissimilarity calculation unit calculates a dissimilarity of the power spectrum between two adjacent points on the optical fiber; The environment change position estimation unit The section including all the pillars whose positions are to be estimated is determined as the analysis section. A window function is configured based on a predetermined column distance indicating a distance between adjacent columns in the analysis section; searching for an offset point in the analysis interval where a correlation function between the window function and the dissimilarity takes a maximum value; In the analysis section, the position of the offset point and the position obtained by adding the distance between the pillars to the offset point are estimated as the positions of the pillars suspending the optical fiber. The position evaluation device according to claim 1 .

5. The optical fiber is a pillar-mounted optical fiber, the vibration characteristic calculation unit calculates a power spectrum in a predetermined frequency band as the sensing data, the dissimilarity calculation unit calculates a dissimilarity of the power spectrum between two adjacent points of the optical fiber, and calculates a weighted dissimilarity by weighting the calculated dissimilarity with a frequency average value of the power spectrum; The environment change position estimation unit determining a section from a pole including the extra length section of the optical fiber to a next pole including the extra length section as an analysis section; A window function is configured based on a predetermined column distance indicating a distance between adjacent columns in the analysis section; searching for an offset point in the analysis interval where a correlation function between the window function and the weighted dissimilarity takes a maximum value; In the analysis section, the position of the offset point and the position obtained by adding the distance between the pillars to the offset point are estimated as the positions of the pillars suspending the optical fiber; In the analysis target section, the position of the offset point and a position on the sensor side of the offset point where the weighted difference takes a maximum peak value are estimated to be the left and right ends of the extra length section. The position evaluation device according to claim 1 .

6. the predetermined frequency band is a frequency band including a fundamental vibration mode in the aerial section of the optical fiber; The position evaluation device according to any one of claims 2 to 5.

7. A location estimation method executed by a location estimation device, comprising: a vibration characteristics calculation step of receiving a signal indicating a natural vibration occurring at each position of the optical fiber from a sensor and calculating sensing data indicating the vibration characteristics at each position of the optical fiber based on the received signal; a difference calculation step of calculating a difference in the sensing data between two adjacent points of the optical fiber; and an environment change position estimating step of estimating an environment change position where the environment of the optical fiber changes based on the degree of difference.

8. The optical fiber is a pillar-mounted optical fiber, In the vibration characteristic calculation step, a power spectrum in a predetermined frequency band is calculated as the sensing data; In the dissimilarity calculation step, a dissimilarity of the power spectrum between two adjacent points on the optical fiber is calculated, In the environmental change position estimating step, a position where the difference degree takes a maximum peak value is estimated to be the position of a pillar suspending the optical fiber. The location evaluation method according to claim 7 .

9. The optical fiber is a pillar-mounted optical fiber, In the vibration characteristic calculation step, a power spectrum in a predetermined frequency band is calculated as the sensing data; In the dissimilarity calculation step, a dissimilarity of the power spectrum between two adjacent points of the optical fiber is calculated, and a weighted dissimilarity is calculated by weighting the calculated dissimilarity by a frequency average value of the power spectrum; In the environmental change position estimating step, a range in which the frequency average value of the power spectrum is below a threshold is set as a range of the extra length section of the optical fiber, and positions within the range of the extra length section where the weighted dissimilarity takes a maximum peak value are estimated to be the left and right ends of the extra length section. The location evaluation method according to claim 7 .

10. On the computer, a vibration characteristics calculation step of receiving a signal indicating a natural vibration occurring at each position of the optical fiber from a sensor and calculating sensing data indicating the vibration characteristics at each position of the optical fiber based on the received signal; a difference calculation step of calculating a difference in the sensing data between two adjacent points of the optical fiber; an environment change position estimating step of estimating an environment change position where the environment of the optical fiber changes based on the degree of difference.

Images