Field of the Invention
The embodiments of the present invention relate to a data processor and a data processing method for processing strain distribution data which is measured by a distributed optical fiber sensor.
Priority is claimed on Japanese Patent Application No. 2015-256344, filed on Dec. 28, 2015, the contents of which are incorporated herein by reference.
Description of the Related Art
A Fiber Bragg Grating (FBG) sensor is generally used in the field of an optical fiber sensing for measuring a strain by using an optical fiber sensor. The FBG sensor is disclosed in, for example, Non-Patent Document “The Optical Fiber Sensors Guide—Fundamentals & Applications—, Micron Optics” which is available at http://www.micronoptics.com/wp-content/uploads/2015/06/Updated_Optical_Fiber_Sensors_Guide_130529.pdf. The FBG sensor is the sensor in which plural FBG elements (a diffraction grating) are formed inside a core of the optical fiber. The strain of the position where the FBG elements are formed in the optical fiber can be measured by analyzing an optical spectrum of a reflected beam of light of a Bragg wavelength inside the optical fiber by using the FBG sensor. The FBG sensor is referred to as a quasi-distributed optical fiber sensor.
The FBG elements are discretized in the FBG sensor, and the position of the FBG element cannot be changed by workers. Therefore, the strain distribution data at an arbitrary position cannot be acquired with the FBG sensor.
In contrast, a set of strain distribution data can be acquired with the distributed optical fiber sensor using such as the Brillouin scattering phenomenon by measuring the strain at each position while switching (with random access) the position to be measured. Such the optical fiber sensor is referred to as a random access optical fiber sensor.
The random access optical fiber sensor can measure the strain at only one point of position at a certain point of time. It means the strain distribution data of each point acquired by the random access optical fiber sensor has mutually different measured times. Thus, the random access optical fiber sensor is not used to observe the strain distribution at the same time. Furthermore, the random access optical fiber sensor is not used to observe the vibration mode of a measurement target even if a displacement of the measurement target is acquired from the measured strain.
The present invention relates to a data processor and a data processing method which is capable of acquiring strain distribution data at the same time.
One aspect of the present invention, a data processor having an acquirer for acquiring strain distribution data in which a distributed optical fiber sensor measures strains at a set of positions of an optical fiber at different times; and an interpolator for calculating the strains at the set of positions of the optical fiber at the same time by interpolating the positions and the times of the strain distribution data acquired by the acquirer.
The data processor according to the aspect of the present invention, the interpolator is for calculating the strain at the set of positions of the optical fiber at the same time by interpolating the times and the positions of the strain distribution data acquired by the acquirer.
The data processor according to the aspect of the present invention, the interpolator is for calculating the strain at the set of times of the optical fiber at the same position by interpolating the strains and the positions of the strain distribution data acquired by the acquirer.
The data processor according to the aspect of the present invention further includes a calculator for calculating a temporal change of the strain distribution of a measurement target having the optical fiber attached thereto from the strain distribution data interpolated by the interpolator.
The data processor according to the aspect of the present invention includes an output unit for outputting the strain distribution data interpolated by the interpolator.
One aspect of the present invention, a data processing method having acquiring strain distribution data in which a distributed optical fiber sensor measures a strain at a set of positions of an optical fiber at different times; and calculating the strains at the set of positions of the optical fiber at the same time by interpolating positions and times of the strain distribution data.
The data processing method according to the aspect of the present invention acquiring strain distribution data in which a distributed optical fiber sensor measures a strain at a set of positions of an optical fiber at different times; and calculating the strains at the set of positions of the optical fiber at the same time by interpolating positions and times of the strain distribution data.
The data processing method according to the aspect of the present invention acquiring strain distribution data in which a distributed optical fiber sensor measures a strain at a set of positions of an optical fiber at different times; and calculating the strains at the set of times of the optical fiber at the same position by interpolating positions and times of the strain distribution data.
According to the present invention, it is capable of acquiring the strain distribution data in the same or identical time.
Hereinafter, a detailed description of a preferred embodiment of a data processor, a data processing method according to the present invention will be described in detail below with reference to the drawings. In all figures, the same reference numerals are attached to such as the same or corresponding configuration and the same or corresponding part.
The optical fiber sensor 2 measures a strain at plural positions of a measurement target 4 at mutually different times. The optical fiber sensor 2 is a distributed optical fiber sensor. The distributed optical fiber sensor 2 includes such as the BOCDA (Brillouin Optical Correlation Domain Analysis), the BOCDR (Brillouin Optical Correlation Domain Reflectometry), and the BOTDR (Brillouin Optical Time Domain Reflectometer) using the Brillouin scattering phenomenon. The BOCDA and the BOCDR are referred to as a random access optical fiber sensor, and are particularly suitable for the present embodiment. The optical fiber sensor 2 includes an optical fiber 21 and a measurer 22.
The optical fiber 21 is attached to the measurement target 4. The strain of the measurement target 4 occurs in the optical fiber 21 in a pseudo manner. The optical fiber 21 is preferably brought into close contact with the measurement target 4 and installed in a linear shape or a loop shape in order to improve a measurement sensitivity of the strain.
The measurer 22 measures a strain of the optical fiber 21. The measurer 22 measures a strain at the plural positions of the optical fiber 21 by using such as the stimulated Brillouin scattering phenomenon. It is known that the frequency difference between an incident beam of light and a Brillouin backscattering beam of light depends on a strain amount at where the stimulated Brillouin scattering beam of light generated. For example, using BOCDA to detect the frequency difference at a specific position, the measurer 22 generates the stimulated Brillouin scattering beam of light by the pump beam of light incident on the optical fiber 21, and modulates a frequency of the light source (not shown). The measurer 22 selects the stimulated Brillouin scattering beam of light at a specific position and observes the stimulated Brillouin scattering beam of light. Then, the measurer 22 detects the stimulated Brillouin scattering beam of light which is locally generated, and calculates the strain of the optical fiber 21 at the position where the stimulated Brillouin scattering beam of light is generated.
However, the measurer 22 can measure only one point of the strain at one time. Therefore, the measurer 22 continuously changes the position at which the stimulated Brillouin scattering is locally generated by changing the modulation frequency of the light source. According to the above operation, the measurer 22 can measure the strain distribution over the entire length of the optical fiber 21. That is, the optical fiber sensor 2 measures the strain of the plural positions of the optical fiber 21 at plural different times. As a result, it can be the acquired strain distribution data which is acquired by measuring the strains at the plural positions of the optical fiber 21 at the mutually plural different times. The strain distribution data is composed of the plural strain measurement values measured by the measurer 22. Specifically, the strain distribution data is composed of at least mutually different four or more strain measurement values acquired by measuring the strains of the two or more positions of the optical fiber 21 at the mutually different two or more times. In this case, if it is acquired the mutually different six or more strain measurement values which are acquired by measuring the strains at the two or more positions of the optical fiber 21 at the mutually different three or more times, it is suitable for observing the vibration mode of the measurement target 4 described below.
The strain distribution data processor 3 is constructed by a computer device such as a personal computer. The strain distribution data processor 3 is capable of connecting with the optical fiber sensor 2, and calculates the strain at the plural positions of the measurement target 4 at the same (identical) time based on the strain distribution data acquired from the optical fiber sensor 2.
The acquirer 31 acquires the strain distribution data of the optical fiber 21 acquired by the optical fiber sensor 2. Then, the acquirer 31 transmits the acquired strain distribution data to the interpolator 32.
The interpolator 32 calculates the strain of the plural positions of the optical fiber 21 at the same time by the interpolation (a surface interpolation) of the time and the position of the strain distribution data acquired by the acquirer 31. This interpolation can be carried out by various known methods, for example, it can be carried out by following formula 1 or formula 2. The formula 1 is an Exact Method that interpolates so that the observation value is always passed, and the formula 2 is an Approximate Method that the observation value is not passed (Lam, Nina Siu-Ngan. “Spatial interpolation methods: are view.” The American Cartographer 10.2 (1983): 129-150.). In the formula 1 and the formula 2, the “p” represents the strain, the “x” represents the position, the “y” represents the time, the “a” represents the coefficient, and the “z” represents the observed strain. In the formula 1 and the formula 2, the “N” is an integer of 2 or more, and preferably an integer of 3 or more.
Hereinafter, an interpolation processing by the interpolator 32 will be described with reference to
As shown in
Therefore, the interpolator 32 calculates or acquires the strain (the distribution of strain) of the optical fiber 21 at every predetermined time interval and at every predetermined position by interpolating the strain distribution data acquired by the acquirer 31 in terms of the time and the position (a time axis direction and/or a position axis direction). The interpolator 32 calculates the strain of the mutually different plural positions of the optical fiber 21 at the same time (hereinafter, including the meaning of the specific time).
Furthermore, the interpolator 32 calculates the strain at the mutually different plural times of the optical fiber 21 at the same position (hereinafter, including the meaning of the specific position) by interpolation (the surface interpolation) of the strain and the position of the strain distribution data of which the acquirer 31 acquired. This interpolation can be carried out by various known methods and can be carried out, for example, by the above formula 1 or formula 2.
The calculator 33 calculates the temporal change of the strain distribution of the measurement target 4 having the optical fiber 21 attached thereto from the strain distribution data (interpolated strain distribution data) generated or interpolated by the interpolator 32. As described above, the interpolator 32 interpolates the strain distribution data at the plural positions and the mutually different plural times of the measurement target 4. The calculator 33 calculates the temporal change of the strain distribution of the measurement target 4 from the strain distribution data. The temporal change of the strain distribution of the measurement target 4 can be represented by, for example, the natural frequency of the measurement target 4 and the vibration mode of the measurement target 4. The natural frequency and the vibration mode can be calculated by a well-known method.
The output unit 34 outputs the strain at the mutually different plural positions at the same time of the optical fiber 21 and the strain at the mutually different plural times at the same position of the optical fiber 21 as the strain distribution data. That is, the output unit 34 outputs the strain distribution data generated or interpolated by the interpolator 32 (see
Next, a strain distribution data processing method will be described.
First, in the measuring step S1, the measurer 22 measures the strain at the plural positions of the optical fiber 21 at the mutually different plural times. The optical fiber sensor 2 acquires the strain distribution data acquired by measuring the strain at the plural positions of the optical fiber 21 at the mutually different plural times.
Next, in the acquiring step S2, the acquirer 31 acquires the strain distribution data acquired in the measuring step S1. In other words, the acquirer 31 acquires the strain distribution data in which the optical fiber sensor 2 measures the strain at the plural positions of the optical fiber 21 at the mutually different plural times respectively.
Next, in the interpolating step S3, the interpolator 32 interpolates the strain distribution data acquired in the acquiring step S2 in the time direction and in the position direction, and calculates the strain at the mutually different plural positions of the optical fiber 21 at the same time. Furthermore, the interpolator 32 interpolates the strain distribution data acquired in the acquisition step S2 in the time direction and in the position direction, and calculates the strain at the mutually different plural times of the optical fiber 21 at the same position.
Next, in the calculating step S4, the calculator 33 calculates the temporal change or the positional change of the strain distribution of the measurement target 4 based on the strain distribution data generated or interpolated in the interpolating step S3.
Next, in the outputting step S5, the output unit 34 outputs the strain of the mutually different plural positions at the same time of the optical fiber 21 calculated in the interpolating step S3 and the strain at the mutually different plural times at the same position of the optical fiber 21 calculated in the interpolating step S3. In the outputting step S5, it may output the temporal change or the positional change of the strain distribution of the measurement target 4 which is calculated in the calculating step S4.
As described above, in the present embodiment, the distributed optical fiber sensor 2 measures the strain at the plural positions of the optical fiber 21 at the mutually different plural times, and the acquirer 31 acquires the strain distribution data based on the data of the strain. The interpolator 32 calculates the strain at the mutually different plural positions of the optical fiber 21 at the same time by interpolating the strain distribution data in the time direction and in the position direction. Thus, it is possible to acquire the strain distribution data at the plural positions of the optical fiber 21 at the same time. Further, it is possible to acquire the strain distribution data at the same time of the measurement target 4 by attaching the one or more optical fibers 21 to one or more positions of the measurement target 4. As a result, it is possible to observe the strain distribution of the measurement target 4 at the same time. In this case, it is possible to observe the vibration mode of the measurement target 4 by calculating the displacement from the strain.
The interpolator 32 interpolates the strain distribution data in the position direction and the time direction, and calculates the strain at the mutually different plural times at the same position of the optical fiber 21. The interpolator 32 can acquire the temporal change of the strain at the same position of the optical fiber 21. The interpolator 32 can acquire the temporal change of the strain at the same position of the measurement target 4 by attaching the one or more optical fibers 21 to plural points of the measurement target 4.
The calculator 33 calculates the temporal change of the strain distribution of the measurement target 4 from the strain distribution data generated or interpolated by the interpolator 32. The worker can easily observe such as the vibration mode of the measurement target 4.
The output unit 34 outputs the strain distribution data calculated by the interpolator 32. The worker can easily observe the strain distribution data at the same time of the measurement target 4 and the temporal change of the strain at the same position. The output unit 34 outputs the displacement calculated from the strain. The worker can easily observe the vibration mode of the measurement target 4.
Although preferred embodiments of the present invention have been described above, the present invention is not limited to the above-described embodiments.
Next, examples of a data processor according to the present invention will be described. However, the present invention is not limited to the following examples.
The strain is measured by each strain gauge 16 at a rate of 200 Hz synchronously. Distribution data of the measurement value is referred to as measured strain distribution data of the strain gauge 16.
The strain at the position corresponding to each strain gauge 16 was measured while traveling the positions at a rate of 70 Hz by using the optical fiber sensor 17 simultaneously. The distribution data of the measurement value is referred to as a measured strain distribution data of the optical fiber sensor 17.
The strain at the position corresponding to each strain gauge 16 of the optical fiber 18 at the same time is calculated by interpolating the strain distribution data acquired by the optical fiber sensor 17 in the time direction and the position direction. The calculated strain distribution data is referred to as interpolated strain distribution data of the optical fiber sensor 17.
As shown in
On the other hand, as shown in
As shown in
As shown in
On the other hand, as shown in
As shown in
Next, the table 2 shows the compared result acquired by comparing the natural frequency of the beam member 14 calculated based on the measured strain distribution data of the strain gauge 16 and the natural vibration of the beam member 14 calculated based on the interpolated strain distribution data of the optical fiber sensor 17. The comparison result shows the natural vibration of the beam member 14 in 10 seconds after the beam member 14 was vibrated. The natural vibration of the beam member 14 corresponds to the temporal change of the strain of the entire beam member 14.
As shown in the table 2, the value of the natural vibration of the beam member 14 calculated based on the measured strain distribution data of the strain gauge 16 is same in the conditions 2 and in the conditions 3, and the value of the natural vibration of the beam member 14 calculated based on the interpolated strain distribution of the optical fiber sensor 17 is same in the conditions 2 and in the conditions 3. For this reason, it is difficult to distinguish the conditions 2 from the conditions 3 by the natural vibration. However, the value of the natural vibration of the beam member 14 calculated based on the measured strain distribution data of the strain gauge 16 and the value of the natural vibration of the beam member 14 calculated based on the interpolated strain distribution data of the optical fiber sensor 17 are very close together under all conditions. Therefore, it is understood that the interpolated strain distribution data of the optical fiber sensor 17 has sufficient accuracy as the parameter representing the change in the natural vibration.
A table 3 shows the compared result acquired by comparing a damping ratio of the beam member 14 calculated based on the measured strain distribution data of the strain gauge 16 and the damping ratio of the beam member 14 calculated based on the interpolated strain distribution data of the optical fiber sensor 17. The comparison result shows the damping ratio of the beam member 14 in 10 seconds after the beam member 14 was vibrated.
As shown in the table 3, the damping ratio of the beam member 14 calculated based on the measured strain distribution data of the strain gauge 16 and the damping ratio of the beam member 14 calculated based on the interpolated strain distribution data of the optical fiber sensor 17 goes away as the conditions progresses from the conditions 1 to the conditions 4. However, as described above, since this condition reproduces the situation where the bridge model deteriorates according to the conditions 1 go to the conditions 4, it is possible to estimate the chronological change tendency of the damping ratio in the bridge model from this comparison result.
Next,
As shown in
On the other hand, as shown in
Then, as shown in
Next, the measured strain distribution data of the strain gauge 16 and the interpolated strain distribution data of the optical fiber sensor 17 are compared.
As shown in
Next, in order to ascertain the interpolation accuracy of the interpolated strain distribution data of the optical fiber sensor 17, it is calculated that the correlation coefficient R2 between the measured strain distribution data of the strain gauge 16 and the measured strain distribution data of the optical fiber sensor 17, and the correlation coefficient R2 between the measured strain distribution data of the strain gauge 16 and the interpolated strain distribution data of the optical fiber sensor 17. The calculated results are shown in a table 4.
Further, the two correlation coefficients R2, the discrepancy in natural frequency between the natural frequency based on the measured strain distribution data of the strain gauge 16 and the natural frequency based on the interpolated strain distribution data of the optical fiber sensor 17, and the discrepancy in the damping ratio between the damping ratio based on the measured strain distribution data of the strain gauge 16 and the damping ratio based on the measured strain distribution data of the optical fiber sensor 17 are plotted on
As shown in the table 4 and
Number | Date | Country | Kind |
---|---|---|---|
2015-256344 | Dec 2015 | JP | national |