CN110455438A - Fiber optic temperature distribution detector fiber lengths are calibrated automatically and automatic test approach - Google Patents
Fiber optic temperature distribution detector fiber lengths are calibrated automatically and automatic test approach Download PDFInfo
- Publication number
- CN110455438A CN110455438A CN201910723408.4A CN201910723408A CN110455438A CN 110455438 A CN110455438 A CN 110455438A CN 201910723408 A CN201910723408 A CN 201910723408A CN 110455438 A CN110455438 A CN 110455438A
- Authority
- CN
- China
- Prior art keywords
- ssr
- das
- fewave
- value
- goes
- 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
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01K—MEASURING TEMPERATURE; MEASURING QUANTITY OF HEAT; THERMALLY-SENSITIVE ELEMENTS NOT OTHERWISE PROVIDED FOR
- G01K15/00—Testing or calibrating of thermometers
- G01K15/005—Calibration
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Optical Communication System (AREA)
- Data Exchanges In Wide-Area Networks (AREA)
Abstract
It is calibrated automatically the invention discloses a kind of fiber optic temperature distribution detector fiber lengths and automatic test approach, existing fiber optic temperature distribution detector fiber lengths calibration, which depends on, manually calibrates, since subjective factor is more in manual analysis, calibration effect will affect；And manually calibration efficiency is very low, it is difficult to meet the production of fiber optic temperature distribution detector, in use process the needs of.Multiple down-sampled complex data length analysis method based on Stokes and Anti-Stokes data of the present invention to propose, by being analysed in depth to Stokes in fiber optic temperature distribution detector system with Anti-Stokes data, realize that the identification to Stokes curve and anti-Stokes curve optical fiber connector calculates, and realize analysis to fiber optic temperature distribution detector incoming fiber optic length based on this, which saves labour turnover, calibration accuracy is high, calibration efficiency is high, without hardware modifications, do not increase hardware cost.
Description
Technical field
The present invention relates to fiber optic temperature distribution detector fields, and in particular to fiber optic temperature distribution detector fiber lengths are certainly
Dynamic calibration and automatic test approach.
Background technique
Fiber optic temperature distribution detector (also known as Raman optical time domain reflectometer, abbreviation ROTDR) can continuous, survey due to distribution
The features such as trying distance, corrosion-resistant, essential safety, electromagnetism interference, anti-lightning, long service life is widely used in power cable thermometric
The fields such as degree monitoring, intelligent piping lane temperature monitoring, the monitoring of oil transportation oil storage temperature and fire alarm, but grinding in ROTDR
In system, production and use process, whether the positioning of temperature detection and fire alarm is accurately depended critically upon to fiber lengths
It is accurate to calculate and calibrate.
The calibration of ROTDR fiber lengths is depended at present and is manually calibrated, by artificially manually to optical fiber temperature
Write music line, Stokes curve or anti-Stokes curve is analyzed, and the calibration to fiber lengths is realized, due to manual analysis
Middle subjective factor is more, will affect calibration effect；And manually calibration efficiency is very low, it is difficult to meet ROTDR production, use
Demand in the process；Manually calibrating simultaneously also causes cost of labor to increase.
Summary of the invention
It depends on for the calibration of existing fiber optic temperature distribution detector fiber lengths and manually calibrates, efficiency is very
It is low, it is difficult to meet the production of fiber optic temperature distribution detector, in use process the needs of the problem of, the present invention provides fiber optic temperatures
Distribution detector fiber lengths are calibrated automatically and automatic test approach.
The following technical solution is employed by the present invention:
A kind of fiber optic temperature distribution detector fiber lengths are calibrated automatically and automatic test approach, comprising the following steps:
Step 101: starting fiber optic temperature distribution detector disconnects tested optical fiber；
Step 102: starting fiber lengths automatic calibration function, input optical fibre physical length FL, optical fiber 1550nm refractive index
IR_IN, 1550nm refractive index IR_IN can be provided by optical fiber producer, and setting anti-Stokes optical index IR_AS is IR_IN, if
Setting Stokes optical index IR_S is IR_IN, and installation space resolution parameter SFR is 1m, and sampling resolution SAMF is 0.2m,
Testing range FR=fiber lengths FL+1km is set, and setting internal optical fiber Stokes offset data length DC_S is 0, interior lights
Fine anti-Stokes offset data length DC_AS is 0, calculates test data number N, and the representative value of N is more than or equal to (fiber lengths
FL+1km) ÷ sampling resolution SAMF；
Step 103: starting optical fiber connector length computation function, calculating internal optical fiber Stokes curve end PFS_S is
FS_S, internal optical fiber anti-Stokes curve end PFS_AS are FS_AS；
Step 104: by tested optical fiber incoming fiber optic Temperature Distribution tester；
Step 105: being again started up optical fiber connector length computation function, calculate tested optical fiber Stokes curve end PFE_
S is FS_S, and tested optical fiber anti-Stokes curve end PFE_AS is FS_AS；
Step 106: reading Stokes range data DSX [0~N], which passes through Stokes optical index IR_S
And sampling resolution SAMF is calculated；, anti-Stokes range data DASX [0~N], which passes through anti-Stokes
Optical index IR_AS and sampling resolution SAMF are calculated；Calculating Stokes data test optical fiber length S_DL is DSX
The difference of [PFE_S] and DSX [PFS_S], Anti-Stokes data test optical fiber length AS_DL are DASX [PFE_AS] and DASX
The difference of [PFS_AS]；
Step 107: calculating IR_S is (2*S_DL/FL-1) * IR_IN, and IR_AS is (2*AS_DL/FL-1) * IR_IN；
Step 108: DSX [0] to DSX [N], DSAX [0] to DSAX [N] are recalculated according to IR_S and IR_AS；
Step 109: setting DC_S is DSX [PFS_S], and DC_AS is DSAX [PFS_AS]；
Step 110: starting optical fiber connector length computation function, calculating tested optical fiber Stokes curve end PFE_S is
FS_S, tested optical fiber anti-Stokes curve end PFE_AS are FS_AS；
Step 111: judge whether DSX [PFE_S] and the absolute value of FL difference are less than SFR, if so, 112 are gone to step, it is no
Then go to step 115；
Step 112: judge whether DSAX [PFE_AS] and the absolute value of FL difference are less than SFR, if so, 113 are gone to step,
Otherwise 115 are gone to step；
Step 113: judging whether DSX [PFE_S] and the absolute value of DSAX [PFE_AS] difference are less than SFR, if so, turning
Step 114,115 are otherwise gone to step；
Step 114: calibrating successfully, terminate；
Step 115: calibration failure checks optical path, terminates.
Preferably, step 103, the fiber lengths end computing function of step 105 and step 110 specific step is as follows institute
Show, by analyzing Stokes test data and Anti-Stokes data, realizing the accurate calculating of fiber lengths end and knowing
Other:
Step: 10301: start-up temperature distribution tests function, waiting are completed；
Step 10302: reading Stokes test data DS [0~N], anti-Stokes test data DAS [0~N], read
Take test data number N；
Step 10303: calculating Stokes test data noise DNS is the average value that DS [N-500] arrives DS [N], it is anti-this
Lentor test data noise DNAS is the average value that DAS [N-500] arrives DAS [N]；
Step 10304: calculating removal noise Stokes data DS_SN [0] to DS_SN [N] is respectively that DS [0]-DNS is arrived
The value of DS [N]-DNS, removing noise Anti-Stokes data DAS_SN [0] to DAS_SN [N] is respectively that DAS [0]-DNAS is arrived
The value of DAS [N]-DNAS；
Step 10305: the removal noise Stokes data of DS_SN [N] are arrived according to DS_SN [0], DAS_SN [0] is arrived
The removal noise Anti-Stokes data of DAS_SN [N] calculates internal optical path Stokes curve end FS_S, anti-Stokes
Curve end FS_AS；
Step 10306: output Stokes optical fiber length FS_S, anti-Stokes optical fiber length FS_AS terminate.
Preferably, step 10305 specifically includes the following steps:
Step 201: reading the removal noise Stokes data that DS_SN [0] arrives DS_SN [N], DAS_SN [0] arrives DAS_
The removal noise Anti-Stokes data of SN [N], read test spatial resolution parameter SFR, sampling resolution SAMF；
Step 202: calculating the value that Decomposition order SSRN is SFR ÷ SAMF and round up, setting I is 0, M 0, and SFR is typical
Value is 1m, and SAMF representative value is 0.2m, and SAMF is less than SFR, MM N/SSRN；
Step 203: calculating pair that logarithmic resolution Stokes data DS_SN_SSR [I] [M] is DS_SN [M*SSRN+I]
Number calculates the logarithm that logarithmic resolution Anti-Stokes data DAS_SN_SSR [I] [M] is DAS_SN [M*SSRN+I]；
Step 204: judging whether M is greater than MM, if so, going to step 206, otherwise go to step 205；
Step 205: the value that M is arranged adds 1, and goes to step 203 and calculated again；
Step 206: judging whether I is greater than SSRN, if so, going to step 208, otherwise go to step 207；
Step 207: the value that l is arranged adds 1, and goes to step 203 and calculated again；
Step 208: setting J is 0, I 0, and end decision threshold ETH, representative value 2dB/m*SFR is arranged；
Step 209: calculating logarithmic resolution Stokes differential data DS_SN_SSR_DIFF [I] [J] is DS_SN_SSR
The difference of [I] [J] and DS_SN_SSR [I] [J+1] calculate logarithmic resolution anti-Stokes differential data DAS_SN_SSR_DIFF
[I] [J] is the difference of DAS_SN_SSR [I] [J] and DAS_SN_SSR [I] [J+1]；
Step 210: judging whether J is greater than MM-1, if so, going to step 212, otherwise go to step 211；
Step 211: the value that J is arranged adds 1, and goes to step 209 and calculated again；
Step 212: DS_SN_SSR_DIFF [I] [0] is all in DS_SN_SSR_DIFF [I] [MM-1] is greater than for record
The value of ETH, deposit wave packet array DS_SN_SSR_FEWAVE [I] [0] are arrived in DS_SN_SSR_FEWAVE [I] [MM-1], record
DAS_SN_SSR_DIFF [I] [0] is stored in wave packet array to all values greater than ETH in DAS_SN_SSR_DIFF [I] [MM-1]
DAS_SN_SSR_FEWAVE [I] [0] is in DAS_SN_SSR_FEWAVE [I] [MM-1]；
Step 213: judging whether I is greater than SSRN, if so, going to step 215, otherwise go to step 214；
Step 214: the value that l is arranged adds 1, and goes to step 209 and calculated again；
Step 215: DS_SN_SSR_FEWAVE [SSRN] [MM-1], DS_ are arrived according to DS_SN_SSR_FEWAVE [0] [0]
SN_SSR_DIFF [0] [0] arrives DS_SN_SSR_DIFFSSRN] [MM-1] and DS_SN_SSR [0] [0] to DS_SN_SSR
[SSRN] [MM] is calculated, and fiber lengths FS_S is obtained；
According to Rule of judgment to DAS_SN_SSR_FEWAVE [0] [0] to DAS_SN_SSR_FEWAVE [SSRN] [MM-1],
DAS_SN_SSR_DIFF [0] [0] arrives DAS_SN_ to DAS_SN_SSR_DIFF [SSRN] [MM-1] and DAS_SN_SSR [0] [0]
SSR [SSRN] [MM] is calculated, and fiber lengths FS_AS is obtained；
Step 216: output Stokes optical fiber length FS_S, anti-Stokes optical fiber length FS_AS terminate.
Preferably, the multiple difference wave packet of step 212 determines that array calculates step and includes:
Step 21201: setting JJ is 0, reads the value of I, and DS_SN_SSR_FEWAVE [I] [0] arrives DS_SN_SSR_
FEWAVE [I] [MM] is all set to 0；
Step 21202: judging whether DS_SN_SSR_DIFF [I] [JJ] is more than or equal to ETH, if so, going to step
21203, otherwise go to step 21209；
Step 21203: it is JJ+1 that interim counting variable TEMPJ, which is arranged, and goes to step 21204；
Step 21204: judging whether DS_SN_SSR_DIFF [I] [TEMPJ] is more than or equal to ETH, if so, going to step
21207, otherwise go to step 21205；
Step 21205: judging whether TEMPJ is more than or equal to MM-1, if so, going to step 21207, otherwise go to step
21206；
Step 21206: setting TEMPJ is TEMPJ+1, and goes to step 21204 and judged again；
Step 21207: setting DS_SN_SSR_FEWAVE [I] [JJ] to DS_SN_SSR_FEWAVE [I] [TEMPJ-1] is complete
Portion is 1, and goes to step 21208；
Step 21208: setting JJ is TEMPJ, and goes to step 21202；
Step 21209: judging whether DS_SN_SSR_DIFF [I] [JJ] is less than or equal to -1*ETH, if so, going to step
21210, otherwise go to step 21215；
Step 21210: setting TEMPJ is JJ+1, and goes to step 21211；
Step 21211: judging whether DS_SN_SSR_DIFF [I] [JJ] is less than or equal to -1*ETH, if so, going to step
21214, otherwise go to step 21212；
Step 21212: judging whether TEMPJ is more than or equal to MM-1, if so, going to step 21214, otherwise go to step
21213；
Step 21213: the value that TEMPJ is arranged adds 1, and goes to step 21211 and judged again；
Step 21214: setting DS_SN_SSR_FEWAVE [I] [JJ] to DS_SN_SSR_FEWAVE [I] [TEMPJ-1] is complete
Portion is -1, and goes to step 21208；
Step 21215: judging whether JJ is more than or equal to MM-1, if so, going to step 21217, otherwise go to step 21216；
Step 21216: the value that JJ is arranged adds 1, and goes to step 21202；
Step 21217: setting JJ is 0, reads the value of I, DAS_SN_SSR_FEWAVE [I] [0] DS_SN_SSR_FEWAVE
[I] [MM] is all set to 0；
Step 21218: judging whether DAS_SN_SSR_DIFF [I] [JJ] is more than or equal to ETH, if so, going to step
21219, otherwise go to step 21225；
Step 21219: setting TEMPJ is JJ+1, and goes to step 21220；
Step 21220: judging whether DAS_SN_SSR_DIFF [I] [TEMPJ] is more than or equal to ETH, if so, going to step
21223, otherwise go to step 21221；
Step 21221: judging whether TEMPJ is more than or equal to MM-1, if so, going to step 21223, otherwise go to step
21222；
Step 21222: setting TEMPJ is TEMPJ+1, and goes to step 21220 and judged again；
Step 21223: DAS_SN_SSR_FEWAVE [I] [TEMPJ-1] is arrived in setting DAS_SN_SSR_FEWAVE [I] [JJ]
All 1, and go to step 21224；
Step 21224: setting JJ is TEMPJ, and goes to step 21218；
Step 21225: judging whether DAS_SN_SSR_DIFF [I] [JJ] is less than or equal to -1*ETH, if so, going to step
21226, otherwise go to step 21231；
Step 21226: setting TEMPJ is JJ+1, and goes to step 21227；
Step 21227: judging whether DAS_SN_SSR_DIFF [I] [JJ] is less than or equal to -1*ETH, if so, going to step
21230, otherwise go to step 21228；
Step 21228: judging whether TEMPJ is more than or equal to MM-1, if so, going to step 21230, otherwise go to step
21229；
Step 21229: the value that TEMPJ is arranged adds 1, and goes to step 21217 and judged again；
Step 21230: DAS_SN_SSR_FEWAVE [I] [TEMPJ-1] is arrived in setting DAS_SN_SSR_FEWAVE [I] [JJ]
All -1, and go to step 21224；
Step 21231: judging whether JJ is more than or equal to MM-1, if so, going to step 21233, otherwise go to step 21232；
Step 21232: the value that JJ is arranged adds 1, and goes to step 21218；
Step 21233: the value of DS_SN_SSR_FEWAVE [I] [MM] is arrived in output DS_SN_SSR_FEWAVE [I] [0],
The value that DAS_SN_SSR_FEWAVE [I] [0] arrives DAS_SN_SSR_FEWAVE [I] [MM] is exported, is terminated.
Preferably, the step of step 215 calculates comprehensive weight accumulation wave packet array according to many times decomposition difference wave packet array
Include:
Step 21501: reading DS_SN_SSR_FEWAVE [0] [0] to DS_SN_SSR_FEWAVE [SSRN] [MM-1],
DS_SN_SSR_DIFF [0] [0] to DS_SN_SSR_DIFF [SSRN] [MM-1], DS_SN_SSR [0] [0] DS_SN_SSR are arrived
The value of [SSRN] [MM],
It reads DAS_SN_SSR_FEWAVE [0] [0] and arrives DAS_SN_SSR_FEWAVE [SSRN] [MM-1], DAS_SN_
SSR_DIFF [0] [0] arrives DAS_SN_SSR to DAS_SN_SSR_DIFF [SSRN] [MM-1], DAS_SN_SSR [0] [0]
The value of [SSRN] [MM]；
Step 21502: setting JJ is 0, II 0, the comprehensive interim array DS_SN_ of weight accumulation wave packet of Stokes data
SSR_FEWAVE_T [0~MM] is all set to 0；
Step 21503: the value for calculating DS_SN_SSR_FEWAVE_T [JJ] is DS_SN_SSR_FEWAVE_T [JJ] plus DS_
SN_SSR_FEWAVE[II][JJ]；
Step 21504: judging whether II is more than or equal to SSRN, if so, going to step 21506, otherwise go to step 21505；
Step 21505: the value that II is arranged adds 1, and goes to step 21503；
Step 21506: judging whether DS_SN_SSR_FEWAVE_T [JJ-1] meets: 0 < DS_SN_SSR_FEWAVE_T
[JJ-1] < SSRN-1 or 0 > DS_SN_SSR_FEWAVE_T [JJ-1] > 1-SSRN, and DS_SN_SSR_FEWAVE_T
[JJ] is 0, and DS_SN_SSR_FEWAVE_T [JJ-2] is 0, if so, going to step 21507, otherwise goes to step 21508；
Step 21507: setting DS_SN_SSR_FEWAVE_T [JJ-1] is 0, and goes to step 21508；
Step 21508: judging whether JJ is more than or equal to MM-1, if so, going to step 21511, otherwise go to step 21509；
Step 21509: the value that JJ is arranged adds 1, and goes to step 21510；
Step 21510: setting II is 0, and goes to step 21503 and carry out assignment again；
Step 21511: setting JJ is 0, II 0, the comprehensive interim array DAS_ of weight accumulation wave packet of Anti-Stokes data
SN_SSR_FEWAVE_T [0] to DAS_SN_SSR_FEWAVE_T [MM] is all set to 0, goes to step 21512；
Step 21512: setting DAS_SN_SSR_FEWAVE_T [JJ] value be DAS_SN_SSR_FEWAVE_T [JJ]+
DS_SN_SSR_FEWAVE [II] [JJ], goes to step 21513；
Step 21513: judging whether II is more than or equal to SSRN, if so, going to step 21515, otherwise go to step 21514；
Step 21514: the value that II is arranged adds 1, and goes to step 21512 and carry out assignment again；
Step 21515:: judge whether DAS_SN_SSR_FEWAVE_T [JJ-1] meets: 0 < DAS_SN_SSR_
FEWAVE_T [JJ-1] < SSRN-1 or 0 > DAS_SN_SSR_FEWAVE_T [JJ-1] > 1-SSRN, and DAS_SN_
SSR_FEWAVE_T [JJ] is 0, and DAS_SN_SSR_FEWAVE_T [JJ-2] is 0, if so, going to step 21516, otherwise
Go to step 21517；
Step 21516: setting DAS_SN_SSR_FEWAVE_T [JJ-1] is 0, goes to step 21517；
Step 21517: judging whether JJ is more than or equal to MM-1, if so, going to step 21520, otherwise go to step 21518；
Step 21518: the value that JJ is arranged adds 1, goes to step 21519；
Step 21519: setting II is 0, and goes to step 21512 and carry out assignment again；
Step 21520: calculating and the comprehensive weight accumulation wave packet of Stokes is set and veritify data DS_SN_SSR_FEWAVE_
T_HY [JJ+II] is DS_SN_SSR_FEWAVE_T [II] [JJ], and the comprehensive weight accumulation wave packet of anti-Stokes veritifies data
DAS_SN_SSR_FEWAVE_T_HY [JJ+II] be DAS_SN_SSR_FEWAVE_T [II] [JJ], wherein the range of II be 0~
The range of MM-1, JJ are 0~SSRN；
Step 21521: DS_SN_SSR_FEWAVE_T_HY [N], DS_ are arrived according to DS_SN_SSR_FEWAVE_T_HY [0]
The value of SN [0] to DS_SN [N] calculate Stokes optical fiber length FS_S, are arrived according to DAS_SN_SSR_FEWAVE_T_HY [0]
The value of DAS_SN_SSR_FEWAVE_T_HY [N], DAS_SN [0] to DAS_SN [N] calculate anti-Stokes optical fiber length FS_
AS terminates.
Preferably, step 21521 includes: the step of calculating optical fiber connector according to comprehensive weight accumulation wave packet array
Step 301: reading DS_SN_SSR_FEWAVE_T_HY [0] and arrive DS_SN_SSR_FEWAVE_T_HY [N], DS_SN
[0] value of DS_SN [N] is arrived, and calculates FS_S, DAS_SN_SSR_FEWAVE_T_HY [0] is read and arrives DAS_SN_SSR_
FEWAVE_T_HY [N], DAS_SN [0] to DAS_SN [N] value and calculate FS_AS, setting signal decision threshold STH, representative value
For 20dB；
Step 302: setting JJ is N, JAE N, JAL 0, JAR 0；
Step 303: judging whether DS_SN_SSR_FEWAVE_T_HY [JJ] is greater than 0, if so, going to step 308, otherwise
Go to step 304；
Step 304: judging whether JJ is greater than 0, if so, going to step 305, otherwise go to step 306；
Step 305: the value that JJ is arranged subtracts 1, and goes to step 303 and judged again；
Step 306: judging whether JAE is less than N, if so, going to step 314, otherwise go to step 307；
Step 307: feedback error information can not analyze end, terminate；
Step 308: setting II is JJ-1, goes to step 309；
Step 309: judging whether DS_SN_SSR_FEWAVE_T_HY [II] is 0, if so, going to step 315, is otherwise turned
Step 310；
Step 310: judging whether II is greater than 0, if so, going to step 311, otherwise go to step 312；
Step 311: the value that II is arranged subtracts 1, goes to step 309 and is judged again；
Step 312: it is corresponding for DS_SN_SSR_FEWAVE_T_HY [0~N] that the value of interim wave bag data TW [0~N] is set
Value, the interim value for denoising data TD [0~N] is DS_SN [0~N] corresponding value, and TS 0, TE JJ calculate JAR；
Step 313: judging whether JAR is less than STH, if so, going to step 314, otherwise go to step 307；
Step 314: setting FS_S is 0, goes to step 319；
Step 315: the value of setting TW [0] to TW [N] is that DS_SN_SSR_FEWAVE_T_HY [0] arrives DS_SN_SSR_
The corresponding value of FEWAVE_T_HY [N], TD [0] to TD [N] are that DS_SN [0] arrives DS_SN [N] corresponding value, and TS II+1, TE are
JJ calculates JAL and JAR；
Step 316: if JAL is more than or equal to STH and JAR is less than STH, going to step 318, otherwise go to step 317；
Step 317: setting JJ is II-1, and JAE II-1 goes to step 303, judged again；
Step 318: setting FS_S is II-1, goes to step 319；
Step 319: setting JJ is N, and JAE N, JAL 0, JAR 0 go to step 320；
Step 320: judging whether DAS_SN_SSR_FEWAVE_T_HY [JJ] is greater than 0, if so, going to step 324, otherwise
Go to step 321；
Step 321: judging whether JJ is greater than 0, if so, going to step 322, otherwise go to step 323；
Step 322: the value that JJ is arranged subtracts 1, goes to step 320 and is judged again；
Step 323: judging whether JAE is less than N, if so, going to step 330, otherwise go to step 307；
Step 324: setting II is JJ-1, goes to step 325；
Step 325: judging whether DAS_SN_SSR_FEWAVE_T_HY [II] is 0, if so, going to step 331, is otherwise turned
Step 326；
Step 326: judging whether II is greater than 0, if so, going to step 327, otherwise go to step 328；
Step 327: the value that II is arranged subtracts 1, goes to step 325 and is judged again；
Step 328: the value of setting TW [0] to TW [N] is that DAS_SN_SSR_FEWAVE_T_HY [0] arrives DAS_SN_SSR_
The corresponding value of FEWAVE_T_HY [N], TD [0] to TD [N] are that DAS_SN [0] arrives DAS_SN [N] corresponding value, and TS 0, TE are
JJ calculates JAR, is detailed in step 31501；
Step 329: judging whether JAR is less than STH, if so, going to step 330, otherwise go to step 307；
Step 330: setting FS_AS is 0, goes to step 335；
Step 331: the value of setting TW [0] to TW [N] is that DAS_SN_SSR_FEWAVE_T_HY [0] arrives DAS_SN_SSR_
The corresponding value of FEWAVE_T_HY [N], TD [0] to TD [N] are that DAS_SN [0] arrives DAS_SN [N] corresponding value, TS II+1, TE
For JJ, JAL and JAR is calculated；
Step 332: if JAL is more than or equal to STH and JAR is less than STH, going to step 334, otherwise go to step 333；
Step 333: setting JJ is II-1, and JAE II-1 goes to step 320 and judged again；
Step 334: setting FS_AS is II-1, goes to step 335；
Step 335: the value of output FS_S and FS_AS terminates.
Preferably, the calculating step of JAR and JAL includes:
Step 31501: reading TW [0] and arrive TW [N], TD [0] arrives the value of TD [N], reads the value of TS and TE；
Step 31502: setting II is TE+1, goes to step 31503；
Step 31503: judging whether TW [II] is greater than 0, if so, going to step 31506, otherwise go to step 31504；
Step 31504: judging whether II is less than N, if so, going to step 31505, otherwise go to step 31507；
Step 31505: the value that II is arranged adds 1, goes to step 31502 and carries out assignment again；
Step 31506: setting JAR is the average value that TD [TE] arrives TD [II-1], goes to step 31508；
Step 31507: setting JAR is the average value that TD [TE+1] arrives TD [N], goes to step 31508；
Step 31508: judging whether TS is 0, if so, going to step 31509, otherwise goes to step 31510；
Step 31509: exporting the value of JAR, go to step 31517；
Step 31510: setting II is TS-1, goes to step 31511；
Step 31511: judging whether TW [II] is greater than 0, if so, going to step 31514, otherwise go to step 31512；
Step 31512: judging whether II is greater than 0, if so, going to step 31513, otherwise go to step 31515；
Step 31513: the value that II is arranged subtracts 1, goes to step 31510 and carries out assignment again；
Step 31514: setting JAL is the average value that TD [II+1] arrives TD [TS-1], goes to step 31516；
Step 31515: setting JAL is the average value that TD [0] arrives TD [TS-1], goes to step 31516；
Step 31516: the value of output JAL and JAR；
Step 31517: analysis terminates.
The invention has the advantages that:
Fiber optic temperature distribution detector fiber lengths provided by the invention are calibrated automatically and automatic test approach, with proposition
Multiple down-sampled complex data length analysis method based on Stokes and Anti-Stokes data, by fiber optic temperature point
Stokes is analysed in depth with Anti-Stokes data in cloth tester system, is realized to Stokes curve and this anti-support
The identification of gram this curve optical fiber connector calculates, and is realized based on this to fiber optic temperature distribution detector incoming fiber optic length
Analysis, which saves labour turnover, calibration accuracy is high, calibration efficiency is high, without hardware modifications, do not increase hardware
Cost can effectively improve the lengthy calibration speed in the production of fiber optic temperature distribution detector and use process.
Detailed description of the invention
Fig. 1 is calibrated automatically for fiber optic temperature distribution detector fiber lengths and the flow chart of automatic test approach.
Fig. 2 is fiber lengths end computing function idiographic flow schematic diagram.
Fig. 3 is Stokes and anti-Stokes optical fiber size calculating step schematic diagram.
Fig. 4 is that multiple difference wave packet determines array schematic diagram of calculation flow.
Fig. 5 is that comprehensive weight accumulation wave packet array flow diagram is calculated according to many times decomposition difference wave packet array.
Fig. 6 is to calculate optical fiber connector flow diagram according to comprehensive weight accumulation wave packet array.
Fig. 7 is JAR and JAL schematic diagram of calculation flow.
Specific embodiment
A specific embodiment of the invention is described further in the following with reference to the drawings and specific embodiments:
In conjunction with Fig. 1 to Fig. 7, a kind of fiber optic temperature distribution detector fiber lengths are calibrated automatically and automatic test approach, packet
Include following steps:
Step 101: starting fiber optic temperature distribution detector disconnects tested optical fiber；
Step 102: starting fiber lengths automatic calibration function, input optical fibre physical length FL, optical fiber 1550nm refractive index
IR_IN, 1550nm refractive index IR_IN can be provided by optical fiber producer, and setting anti-Stokes optical index IR_AS is IR_IN, if
Setting Stokes optical index IR_S is IR_IN, and installation space resolution parameter SFR is 1m, and sampling resolution SAMF is 0.2m,
Testing range FR=fiber lengths FL+1km is set, and setting internal optical fiber Stokes offset data length DC_S is 0, interior lights
Fine anti-Stokes offset data length DC_AS is 0, calculates test data number N, and the representative value of N is more than or equal to (fiber lengths
FL+1km) ÷ sampling resolution SAMF；
Step 103: starting optical fiber connector length computation function, calculating internal optical fiber Stokes curve end PFS_S is
FS_S, internal optical fiber anti-Stokes curve end PFS_AS are FS_AS；
Step 104: by tested optical fiber incoming fiber optic Temperature Distribution tester；
Step 105: being again started up optical fiber connector length computation function, calculate tested optical fiber Stokes curve end PFE_
S is FS_S, and tested optical fiber anti-Stokes curve end PFE_AS is FS_AS；
Step 106: reading Stokes range data DSX [0~N], which passes through Stokes optical index IR_S
And sampling resolution SAMF is calculated；, anti-Stokes range data DASX [0~N], which passes through anti-Stokes
Optical index IR_AS and sampling resolution SAMF are calculated；Calculating Stokes data test optical fiber length S_DL is DSX
The difference of [PFE_S] and DSX [PFS_S], Anti-Stokes data test optical fiber length AS_DL are DASX [PFE_AS] and DASX
The difference of [PFS_AS]；
Step 107: calculating IR_S is (2*S_DL/FL-1) * IR_IN, and IR_AS is (2*AS_DL/FL-1) * IR_IN；
Step 108: DSX [0] to DSX [N], DSAX [0] to DSAX [N] are recalculated according to IR_S and IR_AS；
Step 109: setting DC_S is DSX [PFS_S], and DC_AS is DSAX [PFS_AS]；
Step 110: starting optical fiber connector length computation function, calculating tested optical fiber Stokes curve end PFE_S is
FS_S, tested optical fiber anti-Stokes curve end PFE_AS are FS_AS；
Step 111: judge whether DSX [PFE_S] and the absolute value of FL difference are less than SFR, if so, 112 are gone to step, it is no
Then go to step 115；
Step 112: judge whether DSAX [PFE_AS] and the absolute value of FL difference are less than SFR, if so, 113 are gone to step,
Otherwise 115 are gone to step；
Step 113: judging whether DSX [PFE_S] and the absolute value of DSAX [PFE_AS] difference are less than SFR, if so, turning
Step 114,115 are otherwise gone to step；
Step 114: calibrating successfully, terminate；
Step 115: calibration failure checks optical path, terminates.
It is as shown in Figure 2:
Step 103, the fiber lengths end computing function of step 105 and step 110 is shown in specific step is as follows, by right
Stokes test data and Anti-Stokes data are analyzed, and realize the accurate calculating and identification of fiber lengths end:
Step: 10301: start-up temperature distribution tests function, waiting are completed；
Step 10302: reading Stokes test data DS [0~N], anti-Stokes test data DAS [0~N], read
Take test data number N；
Step 10303: calculating Stokes test data noise DNS is the average value that DS [N-500] arrives DS [N], it is anti-this
Lentor test data noise DNAS is the average value that DAS [N-500] arrives DAS [N]；
Step 10304: calculating removal noise Stokes data DS_SN [0] to DS_SN [N] is respectively that DS [0]-DNS is arrived
The value of DS [N]-DNS, removing noise Anti-Stokes data DAS_SN [0] to DAS_SN [N] is respectively that DAS [0]-DNAS is arrived
The value of DAS [N]-DNAS；
Step 10305: the removal noise Stokes data of DS_SN [N] are arrived according to DS_SN [0], DAS_SN [0] is arrived
The removal noise Anti-Stokes data of DAS_SN [N] calculates internal optical path Stokes curve end FS_S, anti-Stokes
Curve end FS_AS；
Step 10306: output Stokes optical fiber length FS_S, anti-Stokes optical fiber length FS_AS terminate.
As shown in figure 3, Stokes and anti-Stokes optical fiber length computation process are as follows:
Step 10305 specifically includes the following steps:
Step 201: reading the removal noise Stokes data that DS_SN [0] arrives DS_SN [N], DAS_SN [0] arrives DAS_
The removal noise Anti-Stokes data of SN [N], read test spatial resolution parameter SFR, sampling resolution SAMF；
Step 202: calculating the value that Decomposition order SSRN is SFR ÷ SAMF and round up, setting I is 0, M 0, and SFR is typical
Value is 1m, and SAMF representative value is 0.2m, and SAMF is less than SFR, MM N/SSRN；
Step 203: calculating pair that logarithmic resolution Stokes data DS_SN_SSR [I] [M] is DS_SN [M*SSRN+I]
Number calculates the logarithm that logarithmic resolution Anti-Stokes data DAS_SN_SSR [I] [M] is DAS_SN [M*SSRN+I]；
Step 204: judging whether M is greater than MM, if so, going to step 206, otherwise go to step 205；
Step 205: the value that M is arranged adds 1, and goes to step 203 and calculated again；
Step 206: judging whether I is greater than SSRN, if so, going to step 208, otherwise go to step 207；
Step 207: the value that l is arranged adds 1, and goes to step 203 and calculated again；
Step 208: setting J is 0, I 0, and end decision threshold ETH, representative value 2dB/m*SFR is arranged；
Step 209: calculating logarithmic resolution Stokes differential data DS_SN_SSR_DIFF [I] [J] is DS_SN_SSR
The difference of [I] [J] and DS_SN_SSR [I] [J+1] calculate logarithmic resolution anti-Stokes differential data DAS_SN_SSR_DIFF
[I] [J] is the difference of DAS_SN_SSR [I] [J] and DAS_SN_SSR [I] [J+1]；
Step 210: judging whether J is greater than MM-1, if so, going to step 212, otherwise go to step 211；
Step 211: the value that J is arranged adds 1, and goes to step 209 and calculated again；
Step 212: DS_SN_SSR_DIFF [I] [0] is all in DS_SN_SSR_DIFF [I] [MM-1] is greater than for record
The value of ETH, deposit wave packet array DS_SN_SSR_FEWAVE [I] [0] are arrived in DS_SN_SSR_FEWAVE [I] [MM-1], record
DAS_SN_SSR_DIFF [I] [0] is stored in wave packet array to all values greater than ETH in DAS_SN_SSR_DIFF [I] [MM-1]
DAS_SN_SSR_FEWAVE [I] [0] is in DAS_SN_SSR_FEWAVE [I] [MM-1]；
Step 213: judging whether I is greater than SSRN, if so, going to step 215, otherwise go to step 214；
Step 214: the value that l is arranged adds 1, and goes to step 209 and calculated again；
Step 215: DS_SN_SSR_FEWAVE [SSRN] [MM-1], DS_ are arrived according to DS_SN_SSR_FEWAVE [0] [0]
SN_SSR_DIFF [0] [0] arrives DS_SN_SSR_DIFFSSRN] [MM-1] and DS_SN_SSR [0] [0] to DS_SN_SSR
[SSRN] [MM] is calculated, and fiber lengths FS_S is obtained；
According to Rule of judgment to DAS_SN_SSR_FEWAVE [0] [0] to DAS_SN_SSR_FEWAVE [SSRN] [MM-1],
DAS_SN_SSR_DIFF [0] [0] arrives DAS_SN_ to DAS_SN_SSR_DIFF [SSRN] [MM-1] and DAS_SN_SSR [0] [0]
SSR [SSRN] [MM] is calculated, and fiber lengths FS_AS is obtained；
Step 216: output Stokes optical fiber length FS_S, anti-Stokes optical fiber length FS_AS terminate.
It is as shown in Figure 4:
The multiple difference wave packet of step 212 determines that array calculates step and includes:
Step 21201: setting JJ is 0, reads the value of I, and DS_SN_SSR_FEWAVE [I] [0] arrives DS_SN_SSR_
FEWAVE [I] [MM] is all set to 0；
Step 21202: judging whether DS_SN_SSR_DIFF [I] [JJ] is more than or equal to ETH, if so, going to step
21203, otherwise go to step 21209；
Step 21203: it is JJ+1 that interim counting variable TEMPJ, which is arranged, and goes to step 21204；
Step 21204: judging whether DS_SN_SSR_DIFF [I] [TEMPJ] is more than or equal to ETH, if so, going to step
21207, otherwise go to step 21205；
Step 21205: judging whether TEMPJ is more than or equal to MM-1, if so, going to step 21207, otherwise go to step
21206；
Step 21206: setting TEMPJ is TEMPJ+1, and goes to step 21204 and judged again；
Step 21207: setting DS_SN_SSR_FEWAVE [I] [JJ] to DS_SN_SSR_FEWAVE [I] [TEMPJ-1] is complete
Portion is 1, and goes to step 21208；
Step 21208: setting JJ is TEMPJ, and goes to step 21202；
Step 21209: judging whether DS_SN_SSR_DIFF [I] [JJ] is less than or equal to -1*ETH, if so, going to step
21210, otherwise go to step 21215；
Step 21210: setting TEMPJ is JJ+1, and goes to step 21211；
Step 21211: judging whether DS_SN_SSR_DIFF [I] [JJ] is less than or equal to -1*ETH, if so, going to step
21214, otherwise go to step 21212；
Step 21212: judging whether TEMPJ is more than or equal to MM-1, if so, going to step 21214, otherwise go to step
21213；
Step 21213: the value that TEMPJ is arranged adds 1, and goes to step 21211 and judged again；
Step 21214: setting DS_SN_SSR_FEWAVE [I] [JJ] to DS_SN_SSR_FEWAVE [I] [TEMPJ-1] is complete
Portion is -1, and goes to step 21208；
Step 21215: judging whether JJ is more than or equal to MM-1, if so, going to step 21217, otherwise go to step 21216；
Step 21216: the value that JJ is arranged adds 1, and goes to step 21202；
Step 21217: setting JJ is 0, reads the value of I, DAS_SN_SSR_FEWAVE [I] [0] DS_SN_SSR_FEWAVE
[I] [MM] is all set to 0；
Step 21218: judging whether DAS_SN_SSR_DIFF [I] [JJ] is more than or equal to ETH, if so, going to step
21219, otherwise go to step 21225；
Step 21219: setting TEMPJ is JJ+1, and goes to step 21220；
Step 21220: judging whether DAS_SN_SSR_DIFF [I] [TEMPJ] is more than or equal to ETH, if so, going to step
21223, otherwise go to step 21221；
Step 21221: judging whether TEMPJ is more than or equal to MM-1, if so, going to step 21223, otherwise go to step
21222；
Step 21222: setting TEMPJ is TEMPJ+1, and goes to step 21220 and judged again；
Step 21223: DAS_SN_SSR_FEWAVE [I] [TEMPJ-1] is arrived in setting DAS_SN_SSR_FEWAVE [I] [JJ]
All 1, and go to step 21224；
Step 21224: setting JJ is TEMPJ, and goes to step 21218；
Step 21225: judging whether DAS_SN_SSR_DIFF [I] [JJ] is less than or equal to -1*ETH, if so, going to step
21226, otherwise go to step 21231；
Step 21226: setting TEMPJ is JJ+1, and goes to step 21227；
Step 21227: judging whether DAS_SN_SSR_DIFF [I] [JJ] is less than or equal to -1*ETH, if so, going to step
21230, otherwise go to step 21228；
Step 21228: judging whether TEMPJ is more than or equal to MM-1, if so, going to step 21230, otherwise go to step
21229；
Step 21229: the value that TEMPJ is arranged adds 1, and goes to step 21217 and judged again；
Step 21230: DAS_SN_SSR_FEWAVE [I] [TEMPJ-1] is arrived in setting DAS_SN_SSR_FEWAVE [I] [JJ]
All -1, and go to step 21224；
Step 21231: judging whether JJ is more than or equal to MM-1, if so, going to step 21233, otherwise go to step 21232；
Step 21232: the value that JJ is arranged adds 1, and goes to step 21218；
Step 21233: the value of DS_SN_SSR_FEWAVE [I] [MM] is arrived in output DS_SN_SSR_FEWAVE [I] [0],
The value that DAS_SN_SSR_FEWAVE [I] [0] arrives DAS_SN_SSR_FEWAVE [I] [MM] is exported, is terminated.
It is as shown in Figure 5:
Step 215 calculates the step of integrating weight accumulation wave packet array according to many times decomposition difference wave packet array
Step 21501: reading DS_SN_SSR_FEWAVE [0] [0] to DS_SN_SSR_FEWAVE [SSRN] [MM-1],
DS_SN_SSR_DIFF [0] [0] to DS_SN_SSR_DIFF [SSRN] [MM-1], DS_SN_SSR [0] [0] DS_SN_SSR are arrived
The value of [SSRN] [MM],
It reads DAS_SN_SSR_FEWAVE [0] [0] and arrives DAS_SN_SSR_FEWAVE [SSRN] [MM-1], DAS_SN_
SSR_DIFF [0] [0] arrives DAS_SN_SSR to DAS_SN_SSR_DIFF [SSRN] [MM-1], DAS_SN_SSR [0] [0]
The value of [SSRN] [MM]；
Step 21502: setting JJ is 0, II 0, the comprehensive interim array DS_SN_ of weight accumulation wave packet of Stokes data
SSR_FEWAVE_T [0~MM] is all set to 0；
Step 21503: the value for calculating DS_SN_SSR_FEWAVE_T [JJ] is DS_SN_SSR_FEWAVE_T [JJ] plus DS_
SN_SSR_FEWAVE[II][JJ]；
Step 21504: judging whether II is more than or equal to SSRN, if so, going to step 21506, otherwise go to step 21505；
Step 21505: the value that II is arranged adds 1, and goes to step 21503；
Step 21506: judging whether DS_SN_SSR_FEWAVE_T [JJ-1] meets: 0 < DS_SN_SSR_FEWAVE_T
[JJ-1] < SSRN-1 or 0 > DS_SN_SSR_FEWAVE_T [JJ-1] > 1-SSRN, and DS_SN_SSR_FEWAVE_T
[JJ] is 0, and DS_SN_SSR_FEWAVE_T [JJ-2] is 0, if so, going to step 21507, otherwise goes to step 21508；
Step 21507: setting DS_SN_SSR_FEWAVE_T [JJ-1] is 0, and goes to step 21508；
Step 21508: judging whether JJ is more than or equal to MM-1, if so, going to step 21511, otherwise go to step 21509；
Step 21509: the value that JJ is arranged adds 1, and goes to step 21510；
Step 21510: setting II is 0, and goes to step 21503 and carry out assignment again；
Step 21511: setting JJ is 0, II 0, the comprehensive interim array DAS_ of weight accumulation wave packet of Anti-Stokes data
SN_SSR_FEWAVE_T [0] to DAS_SN_SSR_FEWAVE_T [MM] is all set to 0, goes to step 21512；
Step 21512: setting DAS_SN_SSR_FEWAVE_T [JJ] value be DAS_SN_SSR_FEWAVE_T [JJ]+
DS_SN_SSR_FEWAVE [II] [JJ], goes to step 21513；
Step 21513: judging whether II is more than or equal to SSRN, if so, going to step 21515, otherwise go to step 21514；
Step 21514: the value that II is arranged adds 1, and goes to step 21512 and carry out assignment again；
Step 21515:: judge whether DAS_SN_SSR_FEWAVE_T [JJ-1] meets: 0 < DAS_SN_SSR_
FEWAVE_T [JJ-1] < SSRN-1 or 0 > DAS_SN_SSR_FEWAVE_T [JJ-1] > 1-SSRN, and DAS_SN_
SSR_FEWAVE_T [JJ] is 0, and DAS_SN_SSR_FEWAVE_T [JJ-2] is 0, if so, going to step 21516, otherwise
Go to step 21517；
Step 21516: setting DAS_SN_SSR_FEWAVE_T [JJ-1] is 0, goes to step 21517；
Step 21517: judging whether JJ is more than or equal to MM-1, if so, going to step 21520, otherwise go to step 21518；
Step 21518: the value that JJ is arranged adds 1, goes to step 21519；
Step 21519: setting II is 0, and goes to step 21512 and carry out assignment again；
Step 21520: calculating and the comprehensive weight accumulation wave packet of Stokes is set and veritify data DS_SN_SSR_FEWAVE_
T_HY [JJ+II] is DS_SN_SSR_FEWAVE_T [II] [JJ], and the comprehensive weight accumulation wave packet of anti-Stokes veritifies data
DAS_SN_SSR_FEWAVE_T_HY [JJ+II] be DAS_SN_SSR_FEWAVE_T [II] [JJ], wherein the range of II be 0~
The range of MM-1, JJ are 0~SSRN；
Step 21521: DS_SN_SSR_FEWAVE_T_HY [N], DS_ are arrived according to DS_SN_SSR_FEWAVE_T_HY [0]
The value of SN [0] to DS_SN [N] calculate Stokes optical fiber length FS_S, are arrived according to DAS_SN_SSR_FEWAVE_T_HY [0]
The value of DAS_SN_SSR_FEWAVE_T_HY [N], DAS_SN [0] to DAS_SN [N] calculate anti-Stokes optical fiber length FS_
AS terminates.
It is as shown in Figure 6:
Step 21521 includes: the step of calculating optical fiber connector according to comprehensive weight accumulation wave packet array
Step 301: reading DS_SN_SSR_FEWAVE_T_HY [0] and arrive DS_SN_SSR_FEWAVE_T_HY [N], DS_SN
[0] value of DS_SN [N] is arrived, and calculates FS_S, DAS_SN_SSR_FEWAVE_T_HY [0] is read and arrives DAS_SN_SSR_
FEWAVE_T_HY [N], DAS_SN [0] to DAS_SN [N] value and calculate FS_AS, setting signal decision threshold STH, representative value
For 20dB；
Step 302: setting JJ is N, JAE N, JAL 0, JAR 0；
Step 303: judging whether DS_SN_SSR_FEWAVE_T_HY [JJ] is greater than 0, if so, going to step 308, otherwise
Go to step 304；
Step 304: judging whether JJ is greater than 0, if so, going to step 305, otherwise go to step 306；
Step 305: the value that JJ is arranged subtracts 1, and goes to step 303 and judged again；
Step 306: judging whether JAE is less than N, if so, going to step 314, otherwise go to step 307；
Step 307: feedback error information can not analyze end, terminate；
Step 308: setting II is JJ-1, goes to step 309；
Step 309: judging whether DS_SN_SSR_FEWAVE_T_HY [II] is 0, if so, going to step 315, is otherwise turned
Step 310；
Step 310: judging whether II is greater than 0, if so, going to step 311, otherwise go to step 312；
Step 311: the value that II is arranged subtracts 1, goes to step 309 and is judged again；
Step 312: it is corresponding for DS_SN_SSR_FEWAVE_T_HY [0~N] that the value of interim wave bag data TW [0~N] is set
Value, the interim value for denoising data TD [0~N] is DS_SN [0~N] corresponding value, and TS 0, TE JJ calculate JAR；
Step 313: judging whether JAR is less than STH, if so, going to step 314, otherwise go to step 307；
Step 314: setting FS_S is 0, goes to step 319；
Step 315: the value of setting TW [0] to TW [N] is that DS_SN_SSR_FEWAVE_T_HY [0] arrives DS_SN_SSR_
The corresponding value of FEWAVE_T_HY [N], TD [0] to TD [N] are that DS_SN [0] arrives DS_SN [N] corresponding value, and TS II+1, TE are
JJ calculates JAL and JAR；
Step 316: if JAL is more than or equal to STH and JAR is less than STH, going to step 318, otherwise go to step 317；
Step 317: setting JJ is II-1, and JAE II-1 goes to step 303, judged again；
Step 318: setting FS_S is II-1, goes to step 319；
Step 319: setting JJ is N, and JAE N, JAL 0, JAR 0 go to step 320；
Step 320: judging whether DAS_SN_SSR_FEWAVE_T_HY [JJ] is greater than 0, if so, going to step 324, otherwise
Go to step 321；
Step 321: judging whether JJ is greater than 0, if so, going to step 322, otherwise go to step 323；
Step 322: the value that JJ is arranged subtracts 1, goes to step 320 and is judged again；
Step 323: judging whether JAE is less than N, if so, going to step 330, otherwise go to step 307；
Step 324: setting II is JJ-1, goes to step 325；
Step 325: judging whether DAS_SN_SSR_FEWAVE_T_HY [II] is 0, if so, going to step 331, is otherwise turned
Step 326；
Step 326: judging whether II is greater than 0, if so, going to step 327, otherwise go to step 328；
Step 327: the value that II is arranged subtracts 1, goes to step 325 and is judged again；
Step 328: the value of setting TW [0] to TW [N] is that DAS_SN_SSR_FEWAVE_T_HY [0] arrives DAS_SN_SSR_
The corresponding value of FEWAVE_T_HY [N], TD [0] to TD [N] are that DAS_SN [0] arrives DAS_SN [N] corresponding value, and TS 0, TE are
JJ calculates JAR, is detailed in step 31501；
Step 329: judging whether JAR is less than STH, if so, going to step 330, otherwise go to step 307；
Step 330: setting FS_AS is 0, goes to step 335；
Step 331: the value of setting TW [0] to TW [N] is that DAS_SN_SSR_FEWAVE_T_HY [0] arrives DAS_SN_SSR_
The corresponding value of FEWAVE_T_HY [N], TD [0] to TD [N] are that DAS_SN [0] arrives DAS_SN [N] corresponding value, TS II+1, TE
For JJ, JAL and JAR is calculated；
Step 332: if JAL is more than or equal to STH and JAR is less than STH, going to step 334, otherwise go to step 333；
Step 333: setting JJ is II-1, and JAE II-1 goes to step 320 and judged again；
Step 334: setting FS_AS is II-1, goes to step 335；
Step 335: the value of output FS_S and FS_AS terminates.
It is as shown in Figure 7:
The calculating step of JAR and JAL includes:
Step 31501: reading TW [0] and arrive TW [N], TD [0] arrives the value of TD [N], reads the value of TS and TE；
Step 31502: setting II is TE+1, goes to step 31503；
Step 31503: judging whether TW [II] is greater than 0, if so, going to step 31506, otherwise go to step 31504；
Step 31504: judging whether II is less than N, if so, going to step 31505, otherwise go to step 31507；
Step 31505: the value that II is arranged adds 1, goes to step 31502 and carries out assignment again；
Step 31506: setting JAR is the average value that TD [TE] arrives TD [II-1], goes to step 31508；
Step 31507: setting JAR is the average value that TD [TE+1] arrives TD [N], goes to step 31508；
Step 31508: judging whether TS is 0, if so, going to step 31509, otherwise goes to step 31510；
Step 31509: exporting the value of JAR, go to step 31517；
Step 31510: setting II is TS-1, goes to step 31511；
Step 31511: judging whether TW [II] is greater than 0, if so, going to step 31514, otherwise go to step 31512；
Step 31512: judging whether II is greater than 0, if so, going to step 31513, otherwise go to step 31515；
Step 31513: the value that II is arranged subtracts 1, goes to step 31510 and carries out assignment again；
Step 31514: setting JAL is the average value that TD [II+1] arrives TD [TS-1], goes to step 31516；
Step 31515: setting JAL is the average value that TD [0] arrives TD [TS-1], goes to step 31516；
Step 31516: the value of output JAL and JAR；
Step 31517: analysis terminates.
Certainly, the above description is not a limitation of the present invention, and the present invention is also not limited to the example above, this technology neck
The variations, modifications, additions or substitutions that the technical staff in domain is made within the essential scope of the present invention also should belong to of the invention
Protection scope.
Claims (7)
1. a kind of fiber optic temperature distribution detector fiber lengths are calibrated automatically and automatic test approach, which is characterized in that including with
Lower step:
Step 101: starting fiber optic temperature distribution detector disconnects tested optical fiber；
Step 102: starting fiber lengths automatic calibration function, input optical fibre physical length FL, optical fiber 1550nm refractive index IR_
IN, 1550nm refractive index IR_IN can be provided by optical fiber producer, and setting anti-Stokes optical index IR_AS is IR_IN, setting
Stokes optical index IR_S is IR_IN, and installation space resolution parameter SFR is 1m, and sampling resolution SAMF is 0.2m, if
Testing range FR=fiber lengths FL+1km is set, setting internal optical fiber Stokes offset data length DC_S is 0, internal optical fiber
Anti-Stokes offset data length DC_AS is 0, calculates test data number N, and the representative value of N is more than or equal to (fiber lengths FL
+ 1km) ÷ sampling resolution SAMF；
Step 103: starting optical fiber connector length computation function, calculating internal optical fiber Stokes curve end PFS_S are FS_S,
Internal optical fiber anti-Stokes curve end PFS_AS is FS_AS；
Step 104: by tested optical fiber incoming fiber optic Temperature Distribution tester；
Step 105: being again started up optical fiber connector length computation function, calculating tested optical fiber Stokes curve end PFE_S is
FS_S, tested optical fiber anti-Stokes curve end PFE_AS are FS_AS；
Step 106: read Stokes range data DSX [0~N], the data by Stokes optical index IR_S and
Sampling resolution SAMF is calculated；, anti-Stokes range data DASX [0~N], the data pass through anti-Stokes light roll over
It penetrates rate IR_AS and sampling resolution SAMF is calculated；Calculating Stokes data test optical fiber length S_DL is DSX
The difference of [PFE_S] and DSX [PFS_S], Anti-Stokes data test optical fiber length AS_DL are DASX [PFE_AS] and DASX
The difference of [PFS_AS]；
Step 107: calculating IR_S is (2*S_DL/FL-1) * IR_IN, and IR_AS is (2*AS_DL/FL-1) * IR_IN；
Step 108: DSX [0] to DSX [N], DSAX [0] to DSAX [N] are recalculated according to IR_S and IR_AS；
Step 109: setting DC_S is DSX [PFS_S], and DC_AS is DSAX [PFS_AS]；
Step 110: starting optical fiber connector length computation function, calculating tested optical fiber Stokes curve end PFE_S are FS_S,
Tested optical fiber anti-Stokes curve end PFE_AS is FS_AS；
Step 111: judging whether DSX [PFE_S] and the absolute value of FL difference are less than SFR, if so, going to step 112, otherwise turn
Step 115；
Step 112: judging whether DSAX [PFE_AS] and the absolute value of FL difference are less than SFR, if so, going to step 113, otherwise
Go to step 115；
Step 113: judging whether DSX [PFE_S] and the absolute value of DSAX [PFE_AS] difference are less than SFR, if so, going to step
114, otherwise go to step 115；
Step 114: calibrating successfully, terminate；
Step 115: calibration failure checks optical path, terminates.
2. a kind of fiber optic temperature distribution detector fiber lengths according to claim 1 are calibrated automatically and automatic test side
Method, which is characterized in that step 103, the fiber lengths end computing function of step 105 and step 110 is shown in specific step is as follows,
By analyzing Stokes test data and Anti-Stokes data, the accurate calculating and identification of realization fiber lengths end:
Step: 10301: start-up temperature distribution tests function, waiting are completed；
Step 10302: reading Stokes test data DS [0~N], anti-Stokes test data DAS [0~N], read and survey
Try data number N；
Step 10303: calculating Stokes test data noise DNS is the average value that DS [N-500] arrives DS [N], anti-stoke
This test data noise DNAS is the average value that DAS [N-500] arrives DAS [N]；
Step 10304: calculating removal noise Stokes data DS_SN [0] to DS_SN [N] is respectively that DS [0]-DNS arrives DS
The value of [N]-DNS, removal noise Anti-Stokes data DAS_SN [0] to DAS_SN [N] are respectively DAS [0]-DNAS to DAS
The value of [N]-DNAS；
Step 10305: the removal noise Stokes data of DS_SN [N] are arrived according to DS_SN [0], DAS_SN [0] arrives DAS_SN
The removal noise Anti-Stokes data of [N] calculates internal optical path Stokes curve end FS_S, anti-Stokes curve end
Hold FS_AS；
Step 10306: output Stokes optical fiber length FS_S, anti-Stokes optical fiber length FS_AS terminate.
3. a kind of fiber optic temperature distribution detector fiber lengths according to claim 2 are calibrated automatically and automatic test side
Method, which is characterized in that step 10305 specifically includes the following steps:
Step 201: reading the removal noise Stokes data that DS_SN [0] arrives DS_SN [N], DAS_SN [0] arrives DAS_SN [N]
Removal noise Anti-Stokes data, read test spatial resolution parameter SFR, sampling resolution SAMF；
Step 202: calculating the value that Decomposition order SSRN is SFR ÷ SAMF and round up, setting I is 0, M 0, and SFR representative value is
1m, SAMF representative value are 0.2m, and SAMF is less than SFR, MM N/SSRN；
Step 203: calculating the logarithm that logarithmic resolution Stokes data DS_SN_SSR [I] [M] is DS_SN [M*SSRN+I], meter
Calculate the logarithm that logarithmic resolution Anti-Stokes data DAS_SN_SSR [I] [M] is DAS_SN [M*SSRN+I]；
Step 204: judging whether M is greater than MM, if so, going to step 206, otherwise go to step 205；
Step 205: the value that M is arranged adds 1, and goes to step 203 and calculated again；
Step 206: judging whether I is greater than SSRN, if so, going to step 208, otherwise go to step 207；
Step 207: the value that l is arranged adds 1, and goes to step 203 and calculated again；
Step 208: setting J is 0, I 0, and end decision threshold ETH, representative value 2dB/m*SFR is arranged；
Step 209: calculating logarithmic resolution Stokes differential data DS_SN_SSR_DIFF [I] [J] is DS_SN_SSR [I] [J]
With the difference of DS_SN_SSR [I] [J+1], calculate logarithmic resolution anti-Stokes differential data DAS_SN_SSR_DIFF [I] [J]
For the difference of DAS_SN_SSR [I] [J] and DAS_SN_SSR [I] [J+1]；
Step 210: judging whether J is greater than MM-1, if so, going to step 212, otherwise go to step 211；
Step 211: the value that J is arranged adds 1, and goes to step 209 and calculated again；
Step 212: all greater than ETH's in record DS_SN_SSR_DIFF [I] [0] to DS_SN_SSR_DIFF [I] [MM-1]
Value, deposit wave packet array DS_SN_SSR_FEWAVE [I] [0] are arrived in DS_SN_SSR_FEWAVE [I] [MM-1], and DAS_ is recorded
SN_SSR_DIFF [I] [0] is stored in wave packet array DAS_ to all values greater than ETH in DAS_SN_SSR_DIFF [I] [MM-1]
SN_SSR_FEWAVE [I] [0] is in DAS_SN_SSR_FEWAVE [I] [MM-1]；
Step 213: judging whether I is greater than SSRN, if so, going to step 215, otherwise go to step 214；
Step 214: the value that l is arranged adds 1, and goes to step 209 and calculated again；
Step 215: DS_SN_SSR_FEWAVE [SSRN] [MM-1], DS_SN_ are arrived according to DS_SN_SSR_FEWAVE [0] [0]
SSR_DIFF [0] [0] arrives DS_SN_SSR_DIFFSSRN] [MM-1] and DS_SN_SSR [0] [0] to DS_SN_SSR [SSRN]
[MM] is calculated, and fiber lengths FS_S is obtained；
DAS_SN_SSR_FEWAVE [SSRN] [MM-1], DAS_ are arrived to DAS_SN_SSR_FEWAVE [0] [0] according to Rule of judgment
SN_SSR_DIFF [0] [0] arrives DAS_SN_SSR to DAS_SN_SSR_DIFF [SSRN] [MM-1] and DAS_SN_SSR [0] [0]
[SSRN] [MM] is calculated, and fiber lengths FS_AS is obtained；
Step 216: output Stokes optical fiber length FS_S, anti-Stokes optical fiber length FS_AS terminate.
4. a kind of fiber optic temperature distribution detector fiber lengths according to claim 3 are calibrated automatically and automatic test side
Method, which is characterized in that the multiple difference wave packet of step 212 determines that array calculates step and includes:
Step 21201: setting JJ is 0, reads the value of I, and DS_SN_SSR_FEWAVE [I] [0] arrives DS_SN_SSR_FEWAVE [I]
[MM] is all set to 0；
Step 21202: judge whether DS_SN_SSR_DIFF [I] [JJ] is more than or equal to ETH, if so, 21203 are gone to step, it is no
Then go to step 21209；
Step 21203: it is JJ+1 that interim counting variable TEMPJ, which is arranged, and goes to step 21204；
Step 21204: judge whether DS_SN_SSR_DIFF [I] [TEMPJ] is more than or equal to ETH, if so, 21207 are gone to step,
Otherwise 21205 are gone to step；
Step 21205: judging whether TEMPJ is more than or equal to MM-1, if so, going to step 21207, otherwise go to step 21206；
Step 21206: setting TEMPJ is TEMPJ+1, and goes to step 21204 and judged again；
Step 21207: setting DS_SN_SSR_FEWAVE [I] [JJ] to DS_SN_SSR_FEWAVE [I] [TEMPJ-1] is all
1, and go to step 21208；
Step 21208: setting JJ is TEMPJ, and goes to step 21202；
Step 21209: judge whether DS_SN_SSR_DIFF [I] [JJ] is less than or equal to -1*ETH, if so, 21210 are gone to step,
Otherwise 21215 are gone to step；
Step 21210: setting TEMPJ is JJ+1, and goes to step 21211；
Step 21211: judge whether DS_SN_SSR_DIFF [I] [JJ] is less than or equal to -1*ETH, if so, 21214 are gone to step,
Otherwise 21212 are gone to step；
Step 21212: judging whether TEMPJ is more than or equal to MM-1, if so, going to step 21214, otherwise go to step 21213；
Step 21213: the value that TEMPJ is arranged adds 1, and goes to step 21211 and judged again；
Step 21214: setting DS_SN_SSR_FEWAVE [I] [JJ] to DS_SN_SSR_FEWAVE [I] [TEMPJ-1] is all
It is -1, and goes to step 21208；
Step 21215: judging whether JJ is more than or equal to MM-1, if so, going to step 21217, otherwise go to step 21216；
Step 21216: the value that JJ is arranged adds 1, and goes to step 21202；
Step 21217: setting JJ is 0, reads the value of I, DAS_SN_SSR_FEWAVE [I] [0] DS_SN_SSR_FEWAVE [I]
[MM] is all set to 0；
Step 21218: judge whether DAS_SN_SSR_DIFF [I] [JJ] is more than or equal to ETH, if so, 21219 are gone to step, it is no
Then go to step 21225；
Step 21219: setting TEMPJ is JJ+1, and goes to step 21220；
Step 21220: judging whether DAS_SN_SSR_DIFF [I] [TEMPJ] is more than or equal to ETH, if so, going to step
21223, otherwise go to step 21221；
Step 21221: judging whether TEMPJ is more than or equal to MM-1, if so, going to step 21223, otherwise go to step 21222；
Step 21222: setting TEMPJ is TEMPJ+1, and goes to step 21220 and judged again；
Step 21223: setting DAS_SN_SSR_FEWAVE [I] [JJ] to DAS_SN_SSR_FEWAVE [I] [TEMPJ-1] is all
It is 1, and goes to step 21224；
Step 21224: setting JJ is TEMPJ, and goes to step 21218；
Step 21225: judging whether DAS_SN_SSR_DIFF [I] [JJ] is less than or equal to -1*ETH, if so, going to step
21226, otherwise go to step 21231；
Step 21226: setting TEMPJ is JJ+1, and goes to step 21227；
Step 21227: judging whether DAS_SN_SSR_DIFF [I] [JJ] is less than or equal to -1*ETH, if so, going to step
21230, otherwise go to step 21228；
Step 21228: judging whether TEMPJ is more than or equal to MM-1, if so, going to step 21230, otherwise go to step 21229；
Step 21229: the value that TEMPJ is arranged adds 1, and goes to step 21217 and judged again；
Step 21230: setting DAS_SN_SSR_FEWAVE [I] [JJ] to DAS_SN_SSR_FEWAVE [I] [TEMPJ-1] is all
It is -1, and goes to step 21224；
Step 21231: judging whether JJ is more than or equal to MM-1, if so, going to step 21233, otherwise go to step 21232；
Step 21232: the value that JJ is arranged adds 1, and goes to step 21218；
Step 21233: the value of DS_SN_SSR_FEWAVE [I] [MM] is arrived in output DS_SN_SSR_FEWAVE [I] [0],
The value that DAS_SN_SSR_FEWAVE [I] [0] arrives DAS_SN_SSR_FEWAVE [I] [MM] is exported, is terminated.
5. a kind of fiber optic temperature distribution detector fiber lengths according to claim 3 are calibrated automatically and automatic test side
Method, which is characterized in that the step of step 215 calculates comprehensive weight accumulation wave packet array according to many times decomposition difference wave packet array is wrapped
It includes:
Step 21501: reading DS_SN_SSR_FEWAVE [0] [0] and arrive DS_SN_SSR_FEWAVE [SSRN] [MM-1], DS_SN_
SSR_DIFF [0] [0] arrives DS_SN_SSR_DIFF [SSRN] [MM-1], DS_SN_SSR [0] [0] DS_SN_SSR to [SSRN]
The value of [MM],
It reads DAS_SN_SSR_FEWAVE [0] [0] and arrives DAS_SN_SSR_FEWAVE [SSRN] [MM-1], DAS_SN_SSR_
DIFF [0] [0] arrives DAS_SN_SSR [SSRN] to DAS_SN_SSR_DIFF [SSRN] [MM-1], DAS_SN_SSR [0] [0]
The value of [MM]；
Step 21502: setting JJ is 0, II 0, the comprehensive interim array DS_SN_SSR_ of weight accumulation wave packet of Stokes data
FEWAVE_T [0~MM] is all set to 0；
Step 21503: the value for calculating DS_SN_SSR_FEWAVE_T [JJ] is DS_SN_SSR_FEWAVE_T [JJ] plus DS_SN_
SSR_FEWAVE[II][JJ]；
Step 21504: judging whether II is more than or equal to SSRN, if so, going to step 21506, otherwise go to step 21505；
Step 21505: the value that II is arranged adds 1, and goes to step 21503；
Step 21506: judging whether DS_SN_SSR_FEWAVE_T [JJ-1] meets: 0 < DS_SN_SSR_FEWAVE_T [JJ-
1] < SSRN-1 or 0 > DS_SN_SSR_FEWAVE_T [JJ-1] > 1-SSRN, and DS_SN_SSR_FEWAVE_T [JJ]
It is 0, and DS_SN_SSR_FEWAVE_T [JJ-2] is 0, if so, going to step 21507, otherwise goes to step 21508；
Step 21507: setting DS_SN_SSR_FEWAVE_T [JJ-1] is 0, and goes to step 21508；
Step 21508: judging whether JJ is more than or equal to MM-1, if so, going to step 21511, otherwise go to step 21509；
Step 21509: the value that JJ is arranged adds 1, and goes to step 21510；
Step 21510: setting II is 0, and goes to step 21503 and carry out assignment again；
Step 21511: setting JJ is 0, II 0, the comprehensive interim array DAS_SN_ of weight accumulation wave packet of Anti-Stokes data
SSR_FEWAVE_T [0] to DAS_SN_SSR_FEWAVE_T [MM] is all set to 0, goes to step 21512；
Step 21512: the value of setting DAS_SN_SSR_FEWAVE_T [JJ] is DAS_SN_SSR_FEWAVE_T [JJ]+DS_SN_
SSR_FEWAVE [II] [JJ], goes to step 21513；
Step 21513: judging whether II is more than or equal to SSRN, if so, going to step 21515, otherwise go to step 21514；
Step 21514: the value that II is arranged adds 1, and goes to step 21512 and carry out assignment again；
Step 21515:: judge whether DAS_SN_SSR_FEWAVE_T [JJ-1] meets: 0 < DAS_SN_SSR_FEWAVE_T
[JJ-1] < SSRN-1 or 0 > DAS_SN_SSR_FEWAVE_T [JJ-1] > 1-SSRN, and DAS_SN_SSR_FEWAVE_
T [JJ] is 0, and DAS_SN_SSR_FEWAVE_T [JJ-2] is 0, if so, going to step 21516, otherwise goes to step 21517；
Step 21516: setting DAS_SN_SSR_FEWAVE_T [JJ-1] is 0, goes to step 21517；
Step 21517: judging whether JJ is more than or equal to MM-1, if so, going to step 21520, otherwise go to step 21518；
Step 21518: the value that JJ is arranged adds 1, goes to step 21519；
Step 21519: setting II is 0, and goes to step 21512 and carry out assignment again；
Step 21520: calculating and the comprehensive weight accumulation wave packet of Stokes is set and veritify data DS_SN_SSR_FEWAVE_T_HY
[JJ+II] is DS_SN_SSR_FEWAVE_T [II] [JJ], and the comprehensive weight accumulation wave packet of anti-Stokes veritifies data DAS_SN_
SSR_FEWAVE_T_HY [JJ+II] is DAS_SN_SSR_FEWAVE_T [II] [JJ], and wherein the range of II is 0~MM-1, JJ
Range be 0~SSRN；
Step 21521: DS_SN_SSR_FEWAVE_T_HY [N], DS_SN [0] are arrived according to DS_SN_SSR_FEWAVE_T_HY [0]
Value to DS_SN [N] calculates Stokes optical fiber length FS_S, arrives DAS_SN_ according to DAS_SN_SSR_FEWAVE_T_HY [0]
The value of SSR_FEWAVE_T_HY [N], DAS_SN [0] to DAS_SN [N] calculate anti-Stokes optical fiber length FS_AS, terminate.
6. a kind of fiber optic temperature distribution detector fiber lengths according to claim 5 are calibrated automatically and automatic test side
Method, which is characterized in that step 21521 includes: the step of calculating optical fiber connector according to comprehensive weight accumulation wave packet array
Step 301: reading DS_SN_SSR_FEWAVE_T_HY [0] and arrived to DS_SN_SSR_FEWAVE_T_HY [N], DS_SN [0]
The value of DS_SN [N], and FS_S is calculated, it reads DAS_SN_SSR_FEWAVE_T_HY [0] and arrives DAS_SN_SSR_FEWAVE_T_HY
The value of [N], DAS_SN [0] to DAS_SN [N] simultaneously calculate FS_AS, setting signal decision threshold STH, representative value 20dB；
Step 302: setting JJ is N, JAE N, JAL 0, JAR 0；
Step 303: judging whether DS_SN_SSR_FEWAVE_T_HY [JJ] is greater than 0, if so, going to step 308, otherwise turn to walk
Rapid 304；
Step 304: judging whether JJ is greater than 0, if so, going to step 305, otherwise go to step 306；
Step 305: the value that JJ is arranged subtracts 1, and goes to step 303 and judged again；
Step 306: judging whether JAE is less than N, if so, going to step 314, otherwise go to step 307；
Step 307: feedback error information can not analyze end, terminate；
Step 308: setting II is JJ-1, goes to step 309；
Step 309: judging whether DS_SN_SSR_FEWAVE_T_HY [II] is 0, if so, going to step 315, is otherwise gone to step
310；
Step 310: judging whether II is greater than 0, if so, going to step 311, otherwise go to step 312；
Step 311: the value that II is arranged subtracts 1, goes to step 309 and is judged again；
Step 312: the value that interim wave bag data TW [0~N] is arranged is that DS_SN_SSR_FEWAVE_T_HY [0~N] is corresponding
Value, the interim value for denoising data TD [0~N] are DS_SN [0~N] corresponding value, and TS 0, TE JJ calculate JAR；
Step 313: judging whether JAR is less than STH, if so, going to step 314, otherwise go to step 307；
Step 314: setting FS_S is 0, goes to step 319；
Step 315: the value of setting TW [0] to TW [N] is that DS_SN_SSR_FEWAVE_T_HY [0] arrives DS_SN_SSR_FEWAVE_
The corresponding value of T_HY [N], TD [0] to TD [N] are that DS_SN [0] arrives DS_SN [N] corresponding value, and TS II+1, TE JJ is calculated
JAL and JAR；
Step 316: if JAL is more than or equal to STH and JAR is less than STH, going to step 318, otherwise go to step 317；
Step 317: setting JJ is II-1, and JAE II-1 goes to step 303, judged again；
Step 318: setting FS_S is II-1, goes to step 319；
Step 319: setting JJ is N, and JAE N, JAL 0, JAR 0 go to step 320；
Step 320: judging whether DAS_SN_SSR_FEWAVE_T_HY [JJ] is greater than 0, if so, going to step 324, otherwise turn to walk
Rapid 321；
Step 321: judging whether JJ is greater than 0, if so, going to step 322, otherwise go to step 323；
Step 322: the value that JJ is arranged subtracts 1, goes to step 320 and is judged again；
Step 323: judging whether JAE is less than N, if so, going to step 330, otherwise go to step 307；
Step 324: setting II is JJ-1, goes to step 325；
Step 325: judging whether DAS_SN_SSR_FEWAVE_T_HY [II] is 0, if so, going to step 331, is otherwise gone to step
326；
Step 326: judging whether II is greater than 0, if so, going to step 327, otherwise go to step 328；
Step 327: the value that II is arranged subtracts 1, goes to step 325 and is judged again；
Step 328: the value of setting TW [0] to TW [N] is that DAS_SN_SSR_FEWAVE_T_HY [0] arrives DAS_SN_SSR_
The corresponding value of FEWAVE_T_HY [N], TD [0] to TD [N] are that DAS_SN [0] arrives DAS_SN [N] corresponding value, and TS 0, TE are
JJ calculates JAR, is detailed in step 31501；
Step 329: judging whether JAR is less than STH, if so, going to step 330, otherwise go to step 307；
Step 330: setting FS_AS is 0, goes to step 335；
Step 331: the value of setting TW [0] to TW [N] is that DAS_SN_SSR_FEWAVE_T_HY [0] arrives DAS_SN_SSR_
The corresponding value of FEWAVE_T_HY [N], TD [0] to TD [N] are that DAS_SN [0] arrives DAS_SN [N] corresponding value, TS II+1, TE
For JJ, JAL and JAR is calculated；
Step 332: if JAL is more than or equal to STH and JAR is less than STH, going to step 334, otherwise go to step 333；
Step 333: setting JJ is II-1, and JAE II-1 goes to step 320 and judged again；
Step 334: setting FS_AS is II-1, goes to step 335；
Step 335: the value of output FS_S and FS_AS terminates.
7. a kind of fiber optic temperature distribution detector fiber lengths according to claim 6 are calibrated automatically and automatic test side
Method, which is characterized in that the calculating step of JAR and JAL includes:
Step 31501: reading TW [0] and arrive TW [N], TD [0] arrives the value of TD [N], reads the value of TS and TE；
Step 31502: setting II is TE+1, goes to step 31503；
Step 31503: judging whether TW [II] is greater than 0, if so, going to step 31506, otherwise go to step 31504；
Step 31504: judging whether II is less than N, if so, going to step 31505, otherwise go to step 31507；
Step 31505: the value that II is arranged adds 1, goes to step 31502 and carries out assignment again；
Step 31506: setting JAR is the average value that TD [TE] arrives TD [II-1], goes to step 31508；
Step 31507: setting JAR is the average value that TD [TE+1] arrives TD [N], goes to step 31508；
Step 31508: judging whether TS is 0, if so, going to step 31509, otherwise goes to step 31510；
Step 31509: exporting the value of JAR, go to step 31517；
Step 31510: setting II is TS-1, goes to step 31511；
Step 31511: judging whether TW [II] is greater than 0, if so, going to step 31514, otherwise go to step 31512；
Step 31512: judging whether II is greater than 0, if so, going to step 31513, otherwise go to step 31515；
Step 31513: the value that II is arranged subtracts 1, goes to step 31510 and carries out assignment again；
Step 31514: setting JAL is the average value that TD [II+1] arrives TD [TS-1], goes to step 31516；
Step 31515: setting JAL is the average value that TD [0] arrives TD [TS-1], goes to step 31516；
Step 31516: the value of output JAL and JAR；
Step 31517: analysis terminates.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910723408.4A CN110455438B (en) | 2019-08-07 | 2019-08-07 | Automatic calibration and automatic test method for optical fiber length of optical fiber temperature distribution tester |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910723408.4A CN110455438B (en) | 2019-08-07 | 2019-08-07 | Automatic calibration and automatic test method for optical fiber length of optical fiber temperature distribution tester |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110455438A true CN110455438A (en) | 2019-11-15 |
| CN110455438B CN110455438B (en) | 2020-12-01 |
Family
ID=68485186
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910723408.4A Active CN110455438B (en) | 2019-08-07 | 2019-08-07 | Automatic calibration and automatic test method for optical fiber length of optical fiber temperature distribution tester |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110455438B (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111043967A (en) * | 2019-12-13 | 2020-04-21 | 北京航天控制仪器研究所 | Method and device for monitoring optical fiber length of distributed optical fiber temperature sensing system in real time |
| CN112050967A (en) * | 2020-08-05 | 2020-12-08 | 中电科仪器仪表有限公司 | Optical fiber temperature automatic calibration and compensation method of optical fiber temperature distribution tester |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2007106827A1 (en) * | 2006-03-13 | 2007-09-20 | Sensortran, Inc. | Dual source calibration for distributed temperature systems |
| CN101794506A (en) * | 2010-02-09 | 2010-08-04 | 上海波汇通信科技有限公司 | Method and device used for data calibration in distributed type optical fiber temperature sensor |
| CN103364112A (en) * | 2013-07-31 | 2013-10-23 | 北京航空航天大学 | Parameter calibration and automatic correction method for distributed optical fiber temperature measurement system |
| CN103424137A (en) * | 2013-08-06 | 2013-12-04 | 神华集团有限责任公司 | Position calibration method and device for distributed optical fiber sensing system |
-
2019
- 2019-08-07 CN CN201910723408.4A patent/CN110455438B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2007106827A1 (en) * | 2006-03-13 | 2007-09-20 | Sensortran, Inc. | Dual source calibration for distributed temperature systems |
| CN101794506A (en) * | 2010-02-09 | 2010-08-04 | 上海波汇通信科技有限公司 | Method and device used for data calibration in distributed type optical fiber temperature sensor |
| CN103364112A (en) * | 2013-07-31 | 2013-10-23 | 北京航空航天大学 | Parameter calibration and automatic correction method for distributed optical fiber temperature measurement system |
| CN103424137A (en) * | 2013-08-06 | 2013-12-04 | 神华集团有限责任公司 | Position calibration method and device for distributed optical fiber sensing system |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111043967A (en) * | 2019-12-13 | 2020-04-21 | 北京航天控制仪器研究所 | Method and device for monitoring optical fiber length of distributed optical fiber temperature sensing system in real time |
| CN111043967B (en) * | 2019-12-13 | 2021-10-01 | 北京航天控制仪器研究所 | Method and device for monitoring optical fiber length of distributed optical fiber temperature sensing system in real time |
| CN112050967A (en) * | 2020-08-05 | 2020-12-08 | 中电科仪器仪表有限公司 | Optical fiber temperature automatic calibration and compensation method of optical fiber temperature distribution tester |
| CN112050967B (en) * | 2020-08-05 | 2022-04-05 | 中电科思仪科技股份有限公司 | Optical fiber temperature automatic calibration and compensation method of optical fiber temperature distribution tester |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110455438B (en) | 2020-12-01 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109474333B (en) | Electric power optical cable real-time monitoring protection system based on deep learning | |
| CN108414113B (en) | Fire alarm system and method for predicting optical fiber temperature by using multipoint temperature discrete coefficients | |
| US20170023504A1 (en) | Method For Lightning Stroke Identification And Location On Optical Fiber Composite Overhead Ground Wire | |
| CN109282837B (en) | Demodulation method of Bragg fiber grating staggered spectrum based on LSTM network | |
| CN110455438A (en) | Fiber optic temperature distribution detector fiber lengths are calibrated automatically and automatic test approach | |
| CN111006786B (en) | Double-path high-precision temperature demodulation method based on distributed optical fiber Raman sensing system | |
| CN105509979A (en) | Fiber optic negative pressure wave-based oil and gas pipeline leakage monitoring positioning system and method | |
| WO2017076189A1 (en) | Otdr event analysis algorithm based on difference window and template matching | |
| CN113358239A (en) | FBG-based wavelength feature identification method | |
| CN113281609B (en) | Active traveling wave positioning method and system for power distribution network fault and storage medium | |
| CN112683418B (en) | Raman scattered light double-path demodulation method for optical fiber distributed temperature measurement | |
| CN111141413B (en) | Distributed optical fiber temperature measurement system and method based on double-end temperature calibration | |
| CN102539113B (en) | A kind of tracing analysis method based on Brillouin optical time domain analysis instrument BOTDA | |
| CN116221625A (en) | Natural gas pipeline monitoring and positioning system and method based on optical fiber acoustic wave technology | |
| CN202631153U (en) | Single-port distributed optic fiber temperature sensor with automatic compensation function | |
| CN102208190A (en) | Measurement method and device for inhibiting unstable noise equipment convergence time | |
| CN115615332A (en) | State on-line monitoring system of optical fiber composite overhead ground wire | |
| CN110763259B (en) | Time domain distributed optical fiber sensing method with variable spatial resolution | |
| CN205280311U (en) | Oil gas pipeline leak detection positioning system based on optic fibre rarefaction wave | |
| CN114742103A (en) | City monitoring data processing method and device based on Internet of things and storage medium | |
| CN110672777B (en) | Catalytic combustion type methane sensor abnormal data identification and analysis method and system | |
| CN108303173B (en) | Distributed optical fiber sensing pipeline disturbance event detection method | |
| CN112327216B (en) | Cable health state assessment device and method based on edge calculation | |
| CN114965325B (en) | Chloride ion detection method and device and computer readable storage medium | |
| CN112050967B (en) | Optical fiber temperature automatic calibration and compensation method of optical fiber temperature distribution tester |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |