CN102095538A - Data demodulation method for polarization maintaining fiber stress sensing - Google Patents

Data demodulation method for polarization maintaining fiber stress sensing Download PDF

Info

Publication number
CN102095538A
CN102095538A CN 201110045099 CN201110045099A CN102095538A CN 102095538 A CN102095538 A CN 102095538A CN 201110045099 CN201110045099 CN 201110045099 CN 201110045099 A CN201110045099 A CN 201110045099A CN 102095538 A CN102095538 A CN 102095538A
Authority
CN
China
Prior art keywords
coupling
signal
data
polarization
component
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN 201110045099
Other languages
Chinese (zh)
Other versions
CN102095538B (en
Inventor
张红霞
叶雯婷
陈信伟
贾大功
温国强
刘铁根
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Tianjin University
Original Assignee
Tianjin University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Tianjin University filed Critical Tianjin University
Priority to CN2011100450993A priority Critical patent/CN102095538B/en
Publication of CN102095538A publication Critical patent/CN102095538A/en
Application granted granted Critical
Publication of CN102095538B publication Critical patent/CN102095538B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Testing Of Optical Devices Or Fibers (AREA)

Abstract

The invention discloses a data demodulation method for polarization maintaining fiber stress sensing. The method comprises the following steps of: acquiring an interference signal by using a polarization maintaining fiber stress sensing test system to acquire photovoltage data; pre-processing the photovoltage data, averaging every n points to acquire new photovoltage data, and data volume; decomposing a photovoltage data signal into an intrinsic mode function (IMF) component and an allowance by an experimental mode method; finding out a basic component and identifying a small coupling point; reconstructing coupling intensity of a signal, and calculating the stress magnitude through intensity of each coupling point so as to realize the distributed detection of a system; and calculating indexes which can show the system test accuracy, such as a signal-to-noise ratio and the like. The method can effectively improve the signal-to-noise ratio of the polarization maintaining fiber stress sensing system, and improves the identification capacity of weak coupling points and the sensitivity of a coupling test. Under the environment of relatively low signal-to-noise ratio, microstress can be detected well.

Description

The data demodulation method of Stress Profile for Polarization-Maintaining sensing
Technical field
The present invention relates to polarizer stress parameters measuring method, belong to field of optical measuring technologies, relate in particular to a kind of data demodulation method of protecting inclined to one side stresses of parts sensing based on interferometric method optics.
Background technology
During the disturbance of factor such as polarization maintaining optical fibre is subjected to external force, magnetic field, bending, reverse the polarization coupled phenomenon can take place, by can realize the sensing of disturbance physical quantity to the measurement of polarization maintaining optical fibre polarization coupled.The distributed polarization maintaining optical fibre sensor of white light interference polarization coupled is widely used in stress, strain, temperature, gas and position sensing owing to its high spatial resolution and wide dynamic range.The characteristics of its maximum are to realize distributed sensing on simple optical fiber, and any position on the optical fiber can be as sensing unit.Because polarization maintaining optical fibre has the polarization hold facility, making does not need to consider the Polarization Control problem in measuring process.
This system can be widely used in the structural health detection.In addition, an important use is that it can intrinsic Coupling point of detection fiber gyro ring and polarization interference.
Utilize polarization maintaining optical fibre polarization coupled principle to realize that stress sensing has polarization time domain reflectometry (POTDR), Michelson white light interference method, carrier being frequency modulated (FMCW) method, modelocked fiber cyclic laser method, detects interferometric method and synthesizing optical coherence function method based on Kerr effect detection method, wavelength.Wherein the white light michelson interferometry is simple in structure with it, and measuring accuracy is high and be widely used.
The factor that influences measuring accuracy has: the AD converter noise of the stability of mechanical scanning step-length, power supply, the noise of detector and resolution characteristic, capture card etc.In order to improve system signal noise ratio, need adopt a series of measures, cooperate corresponding data processing to improve experimental result, such as, data pre-service, Fourier smothing filtering, small echo denoising etc.When the frequency band of the useful information of signal and noise did not overlap, the Fourier smothing filtering just can be obtained effect preferably, yet, a specific cutoff frequency is set and is not suitable for nonlinear noise and has the signal structure of sudden change.In the wavelet transformation, in a single day selected when wavelet basis, in whole decomposition, restructuring procedure, just can't change, and the signal part is not had adaptivity.The essence of wavelet transformation is the Fourier transform with flexible time frequency window.The empirical modal method is suitable for handling the time dependent non-linear non-stationary signal of a large amount of frequencies based on signal local characteristics and adaptivity.
Patent 201010183969.9 sends the in-situ stresses strain data by the timer poll, shows in real time after arithmetic mean filtering.Patent 200920108831.5 usefulness digital signal processing chips and extensive FPGA (Field Programmable Gate Array) array chip demodulate the dependent variable in the fiber grating sensing system as signal processing unit with improved binomial/trinomial fitting algorithm.Most patents all are to do innovation on system architecture, and this patent is in conjunction with the dependent variable of empirical modal method demodulation polarization maintaining optical fibre, and is simple in structure, are easy to realize.
Summary of the invention
The present invention seeks to be to improve the recognition capability of weak coupling point in the Stress Profile for Polarization-Maintaining sensing and the signal to noise ratio (S/N ratio) of raising stiffness of coupling test macro, propose a kind of data demodulation method of the Stress Profile for Polarization-Maintaining sensing based on the empirical modal method.
The present invention has taken into full account the coupling effect that many different power on the optical polarization device cause Coupling point, by data pre-service and improved empirical modal algorithm, make measurement result accurately, reliable.
The data demodulation method of Stress Profile for Polarization-Maintaining sensing provided by the invention comprises:
1st, use the Stress Profile for Polarization-Maintaining sensing testing system to gather interference signal, obtain the photovoltage data
Figure 2011100450993100002DEST_PATH_IMAGE001
2nd, with the photovoltage data
Figure 39129DEST_PATH_IMAGE001
Carry out the data pre-service, every n point obtains new photovoltage data after doing on average
Figure 391613DEST_PATH_IMAGE002
,
Figure 2011100450993100002DEST_PATH_IMAGE003
,
Figure 209527DEST_PATH_IMAGE004
For
Figure 647462DEST_PATH_IMAGE001
Number (value of n need be selected for use by the flexible in size of the data volume that collects, n is big more, then average noise reduction is good more, but if original data volume is little, then n should correspondingly get smaller, for subsequent algorithm provides enough data), purpose is the random noise of effectively eliminating in the acquired signal, and data volume can be dwindled n doubly;
3rd, use the empirical modal method with the photovoltage data-signal
Figure 379663DEST_PATH_IMAGE002
Be decomposed into
Figure 2011100450993100002DEST_PATH_IMAGE005
Individual IMF component and a surplus, promptly
Figure 914550DEST_PATH_IMAGE006
, wherein,
Figure 841049DEST_PATH_IMAGE002
Be the photovoltage data after average, Be the IMF component, reflected different frequencies and amplitude composition,
Figure 828596DEST_PATH_IMAGE005
Be the number of IMF component,
Figure 52904DEST_PATH_IMAGE008
,
Figure 425986DEST_PATH_IMAGE004
For
Figure 772653DEST_PATH_IMAGE002
The data volume size, it be empirical modal decompose satisfy cut-off condition after, the signal number that satisfies the IMF condition of being stripped from, it is with the difference of original signal, the result can be different.
Figure 2011100450993100002DEST_PATH_IMAGE009
Be a surplus, the IMF component that the empirical modal algorithm can be high with frequency at first extracts;
The basic step of empirical modal is as follows:
3.1st, find signal
Figure 998229DEST_PATH_IMAGE002
All Local Extremum, connect maximum value and minimal value respectively, and obtain the coenvelope of original signal with the cubic spline curve match
Figure 26228DEST_PATH_IMAGE010
And lower envelope
3.2nd, calculate the average of coenvelope and lower envelope
Figure 581712DEST_PATH_IMAGE012
, and ask signal
Figure 833702DEST_PATH_IMAGE002
Difference with the gained average
Figure 2011100450993100002DEST_PATH_IMAGE013
3.3rd,
Figure 546574DEST_PATH_IMAGE014
Regard new raw data as, repeat the 3.1st step and the 3.2nd step (soon Regard new raw data as, find All Local Extremum, connect maximum value and minimal value respectively, and obtain its coenvelope with the cubic spline curve match And lower envelope
Figure 148829DEST_PATH_IMAGE016
,Calculate the average of coenvelope and lower envelope , and ask signal
Figure 411314DEST_PATH_IMAGE014
Difference with the gained average
Figure 781115DEST_PATH_IMAGE018
, by that analogy), satisfy two conditions of IMF up to the signal that generates, with the signal definition that satisfies condition be
Figure 2011100450993100002DEST_PATH_IMAGE019
, the definition surplus is
Figure 45612DEST_PATH_IMAGE020
,
Figure 2011100450993100002DEST_PATH_IMAGE021
(each IMF needs two conditions satisfying as follows:
1, local maximum and local minimizing number sum will equate with the number of striding zero point or differ one at most.
2, on the time point in office, be zero by the mean value of defined coenvelope of local maximum and the defined lower envelope of local minimum.)
3.4th, surplus
Figure 904984DEST_PATH_IMAGE020
Regard new raw data as, repeat the 3.1st step to the 3.3rd step up to extracting all IMF components; When satisfying stop condition, iteration stopping (stop condition utilizes the component standard deviation (SD) of double The selection result as criterion,
Figure 326869DEST_PATH_IMAGE022
, typical SD size is at 0.2 to 0.3);
4th, find out substrate component and identify little Coupling point from the IMF component that decomposites, concrete grammar is as follows:
With each IMF component of decompositing respectively with the surplus addition, promptly
Figure 2011100450993100002DEST_PATH_IMAGE023
, the stiffness of coupling of obtaining separately distributes
Figure 562678DEST_PATH_IMAGE024
Be used to judge substrate and little Coupling point; Generally the stiffness of coupling computing method of Cai Yonging are:
Figure 2011100450993100002DEST_PATH_IMAGE025
, wherein
Figure 947261DEST_PATH_IMAGE026
With
Figure 2011100450993100002DEST_PATH_IMAGE027
Being respectively Coupling point interferes envelope and principal maximum to interfere the amplitude of envelope.Because surplus is the trend term of signal, DC component that can representation signal; In conjunction with Stress Profile for Polarization-Maintaining sensor-based system characteristic, the substrate component is the useful signal that identifies at first, is defined as
Figure 728266DEST_PATH_IMAGE028
, can think to identify first of little Coupling point
Figure 2011100450993100002DEST_PATH_IMAGE029
Component; Physical characteristics by Coupling point, stiffness of coupling is higher than average level 10dB and promptly is judged to be little Coupling point, therefore can from figure, identify little Coupling point and reach thereof easily, the scope that contains little Coupling point in this IMF component is superimposed to the substrate component with reconstruction signal, promptly
Figure 198300DEST_PATH_IMAGE030
, wherein
Figure 2011100450993100002DEST_PATH_IMAGE031
It is the medium and small Coupling point part of each IMF;
5th, reconstruction signal
Figure 785270DEST_PATH_IMAGE032
Stiffness of coupling, stress intensity can calculate by the each point stiffness of coupling, thereby has realized the Distributed Detection of system.
Can set up stiffness of coupling and the optic fibre force size and Orientation has following relation according to the theoretical model of three sections polarization maintaining optical fibres of polarization coupled serial connection strain gauge:
?(1)
Wherein
Figure 712775DEST_PATH_IMAGE034
Be beat length of polarization maintaining optical fiber,
Figure 2011100450993100002DEST_PATH_IMAGE035
Be stress intensity,
Figure 225533DEST_PATH_IMAGE036
Be stress direction,
Figure 622011DEST_PATH_IMAGE038
Be the action length of power, can obtain the optic fibre force size by the stool and urine of polarization coupled point stiffness of coupling.
?
Advantage of the present invention and beneficial effect:
The present invention can effectively improve the signal to noise ratio (S/N ratio) of Stress Profile for Polarization-Maintaining sensor-based system, increases the recognition capability of weak coupling point, improves the coupling measurement sensitivity.Under the lower environment of signal to noise ratio (S/N ratio), also can quite good detecting go out microstress.
Description of drawings
Fig. 1 is the process flow diagram of the data demodulates algorithm of Stress Profile for Polarization-Maintaining sensing;
Fig. 2 is a system construction drawing of making the Stress Profile for Polarization-Maintaining sensing of the white light interference method;
Output interferogram when Fig. 3 is a Coupling point;
Fig. 4 is the original signal that computer acquisition is arrived;
Fig. 5 is the stiffness of coupling figure of Fig. 4 correspondence;
Fig. 6 carries out empirical modal to Fig. 4 to decompose the first seven IMF component and surplus that obtains;
Fig. 7 be with the first six IMF respectively with the surplus addition after the stiffness of coupling figure that tries to achieve;
Fig. 8 is with the later signal of demodulating algorithm reconstruct;
Fig. 9 (promptly uses a diameter 125 under embodiment 1 situation
Figure 2011100450993100002DEST_PATH_IMAGE039
, the optical fiber of clapping long 3.1mm experimentizes, and the action direction of power is 15 °, action length is the polarization maintaining optical fibre of 12.5mm), the size of stiffness of coupling and the relation between the stress intensity;
Figure 10 is the stress envelope that this stress sensing system demodulation is come out under embodiment 1 situation.
Among the figure, the 1st, SLD light source, the 2nd, the polarizer, the 3rd, polarization maintaining optical fibre, the 4th, collimation lens, the 5th, turntable, the 6th, Glan prism, the 7th, beam splitter, the 8th, quiet arm, the 9th, swing arm, the 10th, convergent lens, the 11st, stepper motor, the 12nd, capture card, the 13rd, computing machine;
Figure 245628DEST_PATH_IMAGE026
For exciting mould and coupled modes to interfere the amplitude of envelope,
Figure 230901DEST_PATH_IMAGE027
For exciting mould to interfere the amplitude of envelope,
Figure 602977DEST_PATH_IMAGE040
Be Michelson interferometer two arm optical path differences.
Embodiment
Embodiment 1:
Fig. 1 is the process flow diagram of the data demodulates algorithm of Stress Profile for Polarization-Maintaining sensing; Fig. 2 is a hardware configuration of the present invention, wherein, the 1st, SLD light source, the 2nd, the polarizer, the 3rd, polarization maintaining optical fibre, the 4th, collimation lens, the 5th, turntable, the 6th, Glan prism, the 7th, beam splitter, the 8th, quiet arm, the 9th, swing arm, the 10th, convergent lens, the 11st, stepper motor, the 12nd, capture card, the 13rd, computing machine.
The principle and the course of work of the inventive method are as follows:
The distributed Stress Profile for Polarization-Maintaining sensing testing system of white light interference polarization coupled is based on polarization maintaining optical fibre polarization coupled principle.It utilizes external disturbance to cause polarization maintaining optical fibre generation polarization coupled, obtains the position and the size of disturbance point by the Michelson interferometer demodulation, and its structure as shown in Figure 2.
The light that SLD light source 1 sends enters polarization maintaining optical fibre 3 behind the optical fiber polarizer 2, the optical fiber polarizer docks with polarization maintaining optical fibre by adapter, causes that one excites mould to transmit in optical fiber.Do the time spent when certain position of polarization maintaining optical fibre is subjected to external force, can make this point that polarization coupled takes place and inspire one with the coupled modes of importing light wave polarization direction quadrature.The light wave that carries disturbance information incides on the turntable 4, by entering Michelson interferometer 7,8,9,10 behind the prism 6.Because the existence of optic fiber polarization mould dispersion makes two patterns have delay inequality at the optical fiber exit end
Figure 2011100450993100002DEST_PATH_IMAGE041
, obtaining interference signal by 9 these delay inequalities of compensation of the index glass in the mobile Michelson interferometer, it can be expressed as form:
Figure 486750DEST_PATH_IMAGE042
(2)
In the formula
Figure 2011100450993100002DEST_PATH_IMAGE043
Be the interference DC component,
Figure 914058DEST_PATH_IMAGE044
For two feature axis propagation constants of polarization maintaining optical fibre poor,
Figure 2011100450993100002DEST_PATH_IMAGE045
Be the coherence length of laser, Be Michelson interferometer two arm optical path differences, Be the light velocity,
Figure 2011100450993100002DEST_PATH_IMAGE047
Be the Coupling point stiffness of coupling,
Figure 776207DEST_PATH_IMAGE038
Be the length of Coupling point apart from the optical fiber exit end.Fig. 3 represents to have in the optical fiber interferogram of a Coupling point.
The amplitude of stiffness of coupling and interference signal envelope has following relation:
(3)
In the formula For exciting mould and coupled modes to interfere the amplitude of envelope,
Figure 41469DEST_PATH_IMAGE027
Interfere the amplitude of envelope for exciting mould.Distributed sensing system not only is embodied in the position that can locate each disturbance point, goes back its size of energy measurement.Native system is according to interferometer two arm optical path differences
Figure 899835DEST_PATH_IMAGE040
Realize the location positioning of stress point, by measuring the size that the Coupling point stiffness of coupling comes deterministic force.Stress point stiffness of coupling and optic fibre force size and Orientation have following relation:
Figure 295044DEST_PATH_IMAGE033
(1)
Wherein
Figure 906154DEST_PATH_IMAGE034
Be beat length of polarization maintaining optical fiber,
Figure 994196DEST_PATH_IMAGE035
Be stress intensity,
Figure 572814DEST_PATH_IMAGE036
Be stress direction,
Figure 771714DEST_PATH_IMAGE038
Action length for power.
The spatial resolution of sensor-based system
Figure 237330DEST_PATH_IMAGE048
For:
Figure 2011100450993100002DEST_PATH_IMAGE049
(4)
In the formula
Figure 309322DEST_PATH_IMAGE050
Be the light source spectrum width, Be the light source center wavelength.Under the situation that the optical fiber bat is grown and the light source center wavelength is certain, the spatial resolution of system depends primarily on the light source spectrum width.
 
The concrete steps of Stress Profile for Polarization-Maintaining sensing data demodulation method provided by the invention are as follows:
Step1 uses the Stress Profile for Polarization-Maintaining sensing testing system to gather interference signal, obtains the photovoltage data.
In the embodiment of the invention, light source is selected the SLD-561 light source of SUPERLUM for use, and can export with 1310 nm is the light signal of centre wavelength.Follow the tracks of the steady road of voltage stabilizing power supply with the DH1718D-4 two-way and provide 6 V DC voltage, carry out data acquisition with the USB6251 of NI for modulator.Regulate SLD output photocurrent at 136.0 nA, the DC terms of the photovoltage that obtains is 0.74 V, photodetector output connects the analog differential input port of capture card, and sweep velocity is set to 0.2 mm/s, and the capture card sampling rate is set to 200 points of each step motor drive pulse collection.Use a diameter 125
Figure 460991DEST_PATH_IMAGE039
, the optical fiber of clapping long 3.1mm experimentizes, and the action direction of power is 15 °, and action length is 12.5mm, and the signal to noise ratio (S/N ratio) of trying to achieve the signal noise section that collects is 56.28 dB.
Step2 carries out the data pre-service, asks and obtains new photovoltage data after average for per 100, and as shown in Figure 4, Fig. 5 is its corresponding stiffness of coupling.
Step3 has become 16 IMF components and a surplus with the empirical modal method with signal decomposition, and Fig. 6 has shown preceding 7 natural mode of vibration component IMF1-IMF7 and surplus res.
Step4 finds out the substrate component and discerns little Coupling point.
With each IMF component respectively with the surplus addition, promptly
Figure 198003DEST_PATH_IMAGE023
, the stiffness of coupling of obtaining separately distributes
Figure 783705DEST_PATH_IMAGE024
, Fig. 7 has shown the first six component
Figure 213549DEST_PATH_IMAGE052
From figure, can obviously find out,
Figure 2645DEST_PATH_IMAGE014
Comprising three little Coupling points, is the useful signal that identifies at first, is basal signal so IMF1 adds res, Compare
Figure 917249DEST_PATH_IMAGE014
Comprised a new Coupling point, the scope that can be got new Coupling point by algorithm routine is the 14254-14744 point.Be left
Figure 560720DEST_PATH_IMAGE024
All no longer comprise new Coupling point.So reconstruction signal
Figure 895886DEST_PATH_IMAGE030
, promptly the IMF1 of signal length, res and index are that the IMF2 summation of 14254-14744 gets final product.
Step5 asks reconstruction signal
Figure 687125DEST_PATH_IMAGE032
Stiffness of coupling, as shown in Figure 8.Stress intensity can calculate by the each point stiffness of coupling, and Fig. 9 causes the situation of change of Coupling point stiffness of coupling with the size of power for polarization maintaining optical fibre power under this performance.Stress distribution and size that Figure 10 comes out for sensor-based system demodulation under this performance.Calculating signal to noise ratio (S/N ratio) is 72.04dB.Demodulation method of the present invention can effectively improve the signal to noise ratio (S/N ratio) of Stress Profile for Polarization-Maintaining sensor-based system.Under the high-noise environment, stress position and size that this Stress Profile for Polarization-Maintaining sensor-based system still can be exerted oneself demodulation and cause Coupling point have increased the recognition capability of weak coupling point, have improved the coupling measurement sensitivity.
After using this stress sensing demodulating algorithm, the weak coupling point that is difficult to detect can be detected, by Fig. 5 and Fig. 9 as can be seen, experiment causes Coupling point with containing four power in the polarization maintaining optical fibre, because ambient noise interference, the signal that hardware system collects are as shown in Figure 5, promptly detection is exerted oneself and is caused the existence of Coupling point, by obtaining Fig. 9 behind the demodulation method of the present invention, the size and the position probing of four Coupling points can be come out.Be submerged in the stronger noise even signal is described, this algorithm also to a certain extent the power that obtains in the polarization maintaining optical fibre of demodulation cause Coupling point.

Claims (2)

1. the data demodulation method of a Stress Profile for Polarization-Maintaining sensing is characterized in that, this method comprises:
1st, use the Stress Profile for Polarization-Maintaining sensing testing system to gather interference signal, obtain the photovoltage data
Figure 2011100450993100001DEST_PATH_IMAGE001
2nd, with the photovoltage data
Figure 713934DEST_PATH_IMAGE001
Carry out the data pre-service, every n point obtains new photovoltage data after doing on average
Figure 292552DEST_PATH_IMAGE002
,
Figure 2011100450993100001DEST_PATH_IMAGE003
,
Figure 553769DEST_PATH_IMAGE004
For
Figure 35697DEST_PATH_IMAGE001
The data volume size, purpose is the random noise of effectively eliminating in the acquired signal, and data volume can be dwindled n doubly;
3rd, use the empirical modal method with the photovoltage data-signal
Figure 91377DEST_PATH_IMAGE002
Be decomposed into
Figure 2011100450993100001DEST_PATH_IMAGE005
Individual IMF component and a surplus, promptly
Figure 157291DEST_PATH_IMAGE006
, wherein,
Figure 956620DEST_PATH_IMAGE002
Be the photovoltage data after average,
Figure 2011100450993100001DEST_PATH_IMAGE007
Be the IMF component, reflected different frequencies and amplitude composition,
Figure 293055DEST_PATH_IMAGE005
Be the number of IMF component,
Figure 785216DEST_PATH_IMAGE008
, For
Figure 613549DEST_PATH_IMAGE002
The data volume size,
Figure 2011100450993100001DEST_PATH_IMAGE009
Be a surplus, the IMF component that the empirical modal algorithm can be high with frequency at first extracts;
4th, find out substrate component and identify little Coupling point from the IMF component that decomposites, concrete grammar is as follows:
With each IMF component of decompositing respectively with the surplus addition, promptly
Figure 319337DEST_PATH_IMAGE010
, the stiffness of coupling of obtaining separately distributes Be used to judge substrate and little Coupling point; Generally the stiffness of coupling computing method of Cai Yonging are:
Figure 529870DEST_PATH_IMAGE012
, wherein
Figure 2011100450993100001DEST_PATH_IMAGE013
With
Figure 570376DEST_PATH_IMAGE014
Being respectively Coupling point interferes envelope and principal maximum to interfere the amplitude of envelope; Because surplus is the trend term of signal, DC component that can representation signal; In conjunction with Stress Profile for Polarization-Maintaining sensor-based system characteristic, the substrate component is the useful signal that identifies at first, is defined as
Figure 2011100450993100001DEST_PATH_IMAGE015
, can think to identify first of little Coupling point
Figure 524556DEST_PATH_IMAGE016
Component; Physical characteristics by Coupling point, stiffness of coupling is higher than average level 10dB and promptly is judged to be little Coupling point, therefore can from figure, identify little Coupling point and reach thereof easily, the scope that contains little Coupling point in this IMF component is superimposed to the substrate component with reconstruction signal, promptly
Figure 2011100450993100001DEST_PATH_IMAGE017
, wherein
Figure 84851DEST_PATH_IMAGE018
It is the medium and small Coupling point part of each IMF;
5th, reconstruction signal
Figure 2011100450993100001DEST_PATH_IMAGE019
Stiffness of coupling, stress intensity can calculate by the each point stiffness of coupling, thereby has realized the Distributed Detection of system.
2. method according to claim 1 is characterized in that the computing method of described stress intensity of the 5th step are as follows:
Can set up stiffness of coupling and the optic fibre force size and Orientation has following relation according to the theoretical model of three sections polarization maintaining optical fibres of polarization coupled serial connection strain gauge:
Figure 168082DEST_PATH_IMAGE020
(1)
Wherein
Figure 2011100450993100001DEST_PATH_IMAGE021
Be beat length of polarization maintaining optical fiber,
Figure 994087DEST_PATH_IMAGE022
Be stress intensity,
Figure 2011100450993100001DEST_PATH_IMAGE023
Be stress direction,
Figure 2011100450993100001DEST_PATH_IMAGE025
Be the action length of power, can obtain the optic fibre force size by the stool and urine of polarization coupled point stiffness of coupling.
CN2011100450993A 2011-02-25 2011-02-25 Data demodulation method for polarization maintaining fiber stress sensing Expired - Fee Related CN102095538B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2011100450993A CN102095538B (en) 2011-02-25 2011-02-25 Data demodulation method for polarization maintaining fiber stress sensing

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2011100450993A CN102095538B (en) 2011-02-25 2011-02-25 Data demodulation method for polarization maintaining fiber stress sensing

Publications (2)

Publication Number Publication Date
CN102095538A true CN102095538A (en) 2011-06-15
CN102095538B CN102095538B (en) 2012-02-22

Family

ID=44128723

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2011100450993A Expired - Fee Related CN102095538B (en) 2011-02-25 2011-02-25 Data demodulation method for polarization maintaining fiber stress sensing

Country Status (1)

Country Link
CN (1) CN102095538B (en)

Cited By (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102384799A (en) * 2011-09-29 2012-03-21 国电南京自动化股份有限公司 Frequency sweeping and data processing method based on Brillouin distributed fiber sensing system correlation detection scheme
CN102706538A (en) * 2012-05-29 2012-10-03 天津大学 Polarization maintaining optical fiber extinction-ratio self-adaptive demodulation method
CN103808337A (en) * 2012-11-12 2014-05-21 中国航空工业第六一八研究所 White light interference optical sensor noise signal processing method
CN104573248A (en) * 2015-01-16 2015-04-29 东南大学 EMD based fiber-optic gyroscope temperature drift multi-scale extreme learning machine training method
CN104951798A (en) * 2015-06-10 2015-09-30 上海大学 Method for predicting non-stationary fluctuating wind speeds by aid of LSSVM (least square support vector machine) on basis of EMD (empirical mode decomposition)
CN105205495A (en) * 2015-09-02 2015-12-30 上海大学 Non-stationary fluctuating wind speed forecasting method based on EMD-ELM
CN107246952A (en) * 2017-05-19 2017-10-13 北京邮电大学 A kind of polarization mode dispersion measurement precision improvement method and system
CN108489647A (en) * 2018-01-16 2018-09-04 天津大学 A kind of method of dynamic stress frequency in demodulation polarization maintaining optical fibre
CN108489640A (en) * 2018-01-17 2018-09-04 天津大学 A kind of distributed dynamic stress frequency measurement method based on white light interference
CN108955984A (en) * 2018-05-31 2018-12-07 沈阳工业大学 A kind of principal direction of stress determination method that can directly indicate direction
CN110231815A (en) * 2019-05-31 2019-09-13 上海波汇科技有限公司 A kind of DAS system signal detection and the method for assessment
CN110851897A (en) * 2019-10-15 2020-02-28 天津大学 Aqueduct stress-strain prediction method under multi-factor correlation
CN111102994A (en) * 2020-01-19 2020-05-05 枣庄学院 Method for simultaneously measuring polarization coupling distribution of optical fiber ring and beat length of ring-wound optical fiber
CN111103122A (en) * 2019-12-03 2020-05-05 天津大学 Method for extracting weak coupling point for polarization maintaining fiber distributed polarization coupling detection
CN111238772A (en) * 2019-12-31 2020-06-05 中国航空工业集团公司西安飞机设计研究所 Optical fiber ring detection device and detection method based on polarization crosstalk principle
CN111555803A (en) * 2020-05-22 2020-08-18 中天宽带技术有限公司 Bidirectional multi-core optical fiber crosstalk calculation method and device and computer readable storage medium

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20060066839A1 (en) * 2004-08-06 2006-03-30 Payton Robert M Natural fiber span reflectometer providing a virtual signal sensing array capability
JP2008082921A (en) * 2006-09-28 2008-04-10 Oki Electric Ind Co Ltd Optical fiber sensor system
CN101319921A (en) * 2008-07-15 2008-12-10 浙江大学 Method and equipment for optical fibre optical grating sensing network demodulation
CN101358886A (en) * 2008-09-25 2009-02-04 西安科技大学 Grating anchor rod force-measuring device and method for monitoring anchor rod stress variation

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20060066839A1 (en) * 2004-08-06 2006-03-30 Payton Robert M Natural fiber span reflectometer providing a virtual signal sensing array capability
JP2008082921A (en) * 2006-09-28 2008-04-10 Oki Electric Ind Co Ltd Optical fiber sensor system
CN101319921A (en) * 2008-07-15 2008-12-10 浙江大学 Method and equipment for optical fibre optical grating sensing network demodulation
CN101358886A (en) * 2008-09-25 2009-02-04 西安科技大学 Grating anchor rod force-measuring device and method for monitoring anchor rod stress variation

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
《9th International Conference on Optical Communications and Networks》 20101027 Xinwei Chen, etc Spectral domain demodulation of fibre optics position and displacement sensors by Fourier-transform spectral interferogram 第65-69页 1-2 , 2 *
《中国博士学位论文全文数据库 工程科技II辑》 20101015 党淑雯 光纤陀螺的信号分析及滤波技术研究 第99-105页 1-2 , 2 *
《中国激光》 20090331 张锦龙,等 基于保偏光纤的高精度光纤光栅传感解调方法 第728-731页 1-2 第36卷, 第3期 2 *

Cited By (25)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102384799A (en) * 2011-09-29 2012-03-21 国电南京自动化股份有限公司 Frequency sweeping and data processing method based on Brillouin distributed fiber sensing system correlation detection scheme
CN102384799B (en) * 2011-09-29 2014-01-08 国电南京自动化股份有限公司 Frequency sweeping and data processing method based on Brillouin distributed fiber sensing system correlation detection scheme
CN102706538A (en) * 2012-05-29 2012-10-03 天津大学 Polarization maintaining optical fiber extinction-ratio self-adaptive demodulation method
CN102706538B (en) * 2012-05-29 2014-10-15 天津大学 Polarization maintaining optical fiber extinction-ratio self-adaptive demodulation method
CN103808337A (en) * 2012-11-12 2014-05-21 中国航空工业第六一八研究所 White light interference optical sensor noise signal processing method
CN104573248A (en) * 2015-01-16 2015-04-29 东南大学 EMD based fiber-optic gyroscope temperature drift multi-scale extreme learning machine training method
CN104951798A (en) * 2015-06-10 2015-09-30 上海大学 Method for predicting non-stationary fluctuating wind speeds by aid of LSSVM (least square support vector machine) on basis of EMD (empirical mode decomposition)
CN105205495A (en) * 2015-09-02 2015-12-30 上海大学 Non-stationary fluctuating wind speed forecasting method based on EMD-ELM
CN107246952A (en) * 2017-05-19 2017-10-13 北京邮电大学 A kind of polarization mode dispersion measurement precision improvement method and system
CN108489647A (en) * 2018-01-16 2018-09-04 天津大学 A kind of method of dynamic stress frequency in demodulation polarization maintaining optical fibre
CN108489647B (en) * 2018-01-16 2020-10-02 天津大学 Method for demodulating dynamic stress frequency in polarization maintaining optical fiber
CN108489640A (en) * 2018-01-17 2018-09-04 天津大学 A kind of distributed dynamic stress frequency measurement method based on white light interference
CN108489640B (en) * 2018-01-17 2020-10-02 天津大学 Distributed dynamic stress frequency measurement method based on white light interference
CN108955984A (en) * 2018-05-31 2018-12-07 沈阳工业大学 A kind of principal direction of stress determination method that can directly indicate direction
CN108955984B (en) * 2018-05-31 2020-07-14 沈阳工业大学 Principal stress direction judgment method capable of directly indicating direction
CN110231815A (en) * 2019-05-31 2019-09-13 上海波汇科技有限公司 A kind of DAS system signal detection and the method for assessment
CN110851897A (en) * 2019-10-15 2020-02-28 天津大学 Aqueduct stress-strain prediction method under multi-factor correlation
CN110851897B (en) * 2019-10-15 2024-03-22 天津大学 Aqueduct stress-strain prediction method under multi-factor association
CN111103122A (en) * 2019-12-03 2020-05-05 天津大学 Method for extracting weak coupling point for polarization maintaining fiber distributed polarization coupling detection
CN111103122B (en) * 2019-12-03 2024-04-30 天津大学 Polarization maintaining fiber distributed polarization coupling detection weak coupling point extraction method
CN111238772A (en) * 2019-12-31 2020-06-05 中国航空工业集团公司西安飞机设计研究所 Optical fiber ring detection device and detection method based on polarization crosstalk principle
CN111102994A (en) * 2020-01-19 2020-05-05 枣庄学院 Method for simultaneously measuring polarization coupling distribution of optical fiber ring and beat length of ring-wound optical fiber
CN111102994B (en) * 2020-01-19 2021-08-20 枣庄学院 Method for simultaneously measuring polarization coupling distribution of optical fiber ring and beat length of ring-wound optical fiber
CN111555803A (en) * 2020-05-22 2020-08-18 中天宽带技术有限公司 Bidirectional multi-core optical fiber crosstalk calculation method and device and computer readable storage medium
CN111555803B (en) * 2020-05-22 2021-07-27 中天宽带技术有限公司 Bidirectional multi-core optical fiber crosstalk calculation method and device and computer readable storage medium

Also Published As

Publication number Publication date
CN102095538B (en) 2012-02-22

Similar Documents

Publication Publication Date Title
CN102095538B (en) Data demodulation method for polarization maintaining fiber stress sensing
CN112097813B (en) 3X 3 coupler photoelectric detection method and device based on optical frequency modulation
CN102384799B (en) Frequency sweeping and data processing method based on Brillouin distributed fiber sensing system correlation detection scheme
CN102650526B (en) Open-loop detecting circuit for frequency modulated continuous wave optical fiber gyroscope based on phase comparison
CN101270991A (en) System for measuring interfere type optic fiber gyroscope eigenfrequency and half-wave voltage adopting square wave modulation
CN107179097A (en) The optical fiber interference type sensing modulation-demo-demodulation method and device modulated based on phase shift signal
CN104567959A (en) Large-dynamic interference type optical fiber sensor based on two-channel unbalanced interferometer
CN107389978A (en) A kind of weak reflective Bragg gratings accelerometer and its method for sensing
CN102607620A (en) Reflection-type interference optical fiber sensor system
CN102419247B (en) High-precision detection device and method of reflection type optical fiber turbulence
CN104006948A (en) Method for demodulating polarization maintaining optical fiber polarization coupling point position based on multimodal division cycle
CN103335600B (en) Ratio Fringe-counting method and displacement transducer demodulating system thereof based on dual F-P interferometer system
CN105277270A (en) Dual-mode vibration detection system based on fiber grating sensing
CN103196541A (en) Object vibration spectrum measurement system based on laser
Yu et al. Distributed optical fiber vibration sensing using phase-generated carrier demodulation algorithm
CN105092015A (en) Non-contact fiber vibration sensing system and method
Bai et al. Noise reduction method of Φ-OTDR system based on EMD-TFPF algorithm
CN104076165B (en) Contactless transient speed method for sensing
CN108204827A (en) A kind of phase-shifted fiber grating demodulating system
CN104048617B (en) Sensing modulating method for high-precision fiber bragg grating insensitive to polarization state changes
CN203443662U (en) Device for improving signal-to-noise ratio of laser micro-vibration sensing system
CN202648677U (en) Reflecting interference type optical fiber sensor system
CN103454447B (en) Michelson optical fiber accelerometer based on frequency modulation and modulation method thereof
CN102221356B (en) Device and method for measuring laser incident angle by sinusoidally modulating multi-beam laser heterodyne secondary harmonics with Doppler galvanometer
CN108489640B (en) Distributed dynamic stress frequency measurement method based on white light interference

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20120222

Termination date: 20210225

CF01 Termination of patent right due to non-payment of annual fee