The invention discloses an automatic calibration and test method of optical fiber length of optical fiber temperature distribution tester. The existing optical fiber length calibration of optical fiber temperature distribution tester mainly depends on manual calibration. Due to many subjective factors in manual analysis, the calibration effect will be affected; and the manual calibration efficiency is very low, which is difficult to meet the production and use of optical fiber temperature distribution tester Requirements in the process. The invention is based on the proposed multiple downsampling composite data length analysis method based on the Stokes and anti Stokes data, through in-depth analysis of the Stokes and anti Stokes data in the optical fiber temperature distribution tester system, realizes the identification calculation of the end of the optical fiber of the Stokes curve and anti Stokes curve, and on this basis realizes the test of the optical fiber temperature distribution The analysis of the length of the optical fiber connected to the instrument shows that the automatic test method saves labor cost, has high calibration accuracy, high calibration efficiency, does not need hardware modification, and does not increase hardware cost.
【技术实现步骤摘要】
光纤温度分布测试仪光纤长度自动校准及自动测试方法
本专利技术涉及光纤温度分布测试仪领域，具体涉及光纤温度分布测试仪光纤长度自动校准及自动测试方法。
技术介绍
光纤温度分布测试仪(又称拉曼光时域反射计，简称ROTDR)由于分布式可连续、测试距离长、耐腐蚀、本质安全、抗电磁干扰、抗雷击、寿命长等特点已广泛用于电力电缆测温度监测、智能管廊温度监测、输油储油设施温度监测及火灾预警等领域，但是在ROTDR的研制、生产及使用过程中，对温度检测及火灾预警的定位是否准确严重依赖于对光纤长度的准确计算与校准。目前对于ROTDR光纤长度校准主要依赖于人工手动校准，通过人为手动对光纤温度曲线、斯托克斯曲线或反斯托克斯曲线进行分析，实现对光纤长度的校准，由于人工分析中主观因素较多，会影响校准效果；而且人工手动校准效率很低，难以满足ROTDR生产、使用过程中的需求；同时人工手动校准也导致人工成本增加。
技术实现思路
针对现有的光纤温度分布测试仪光纤长度校准主要依赖于人工手动校准，效率很低，难以满足光纤温度分布测试仪生产、使用过程中的需求的问题，本专利技术提供了光纤温度分布测试仪光纤长度自动校准及自动测试方法。本专利技术采用以下的技术方案：一种光纤温度分布测试仪光纤长度自动校准及自动测试方法，包括以下步骤：步骤101：启动光纤温度分布测试仪，断开被测光纤；步骤102：启动光纤长度自动校准功能，输入光纤实际长度FL，光纤1550nm折射率IR_IN，1550nm折射率IR_IN可由光纤厂...
【技术保护点】
1.一种光纤温度分布测试仪光纤长度自动校准及自动测试方法，其特征在于，包括以下步骤：/n步骤101：启动光纤温度分布测试仪，断开被测光纤；/n步骤102：启动光纤长度自动校准功能，输入光纤实际长度FL，光纤1550nm折射率IR_IN，1550nm折射率IR_IN可由光纤厂家给出，设置反斯托克斯光折射率IR_AS为IR_IN，设置斯托克斯光折射率IR_S为IR_IN，设置空间分辨率参数SFR为1m，采样分辨率SAMF为0.2m，设置测试量程FR＝光纤长度FL+1km，设置内部光纤斯托克斯补偿数据长度DC_S为0，内部光纤反斯托克斯补偿数据长度DC_AS为0，计算测试数据数目N，N的典型值大于等于(光纤长度FL+1km)÷采样分辨率SAMF；/n步骤103：启动光纤末端长度计算功能，计算内部光纤斯托克斯曲线末端PFS_S为FS_S，内部光纤反斯托克斯曲线末端PFS_AS为FS_AS；/n步骤104：将被测光纤接入光纤温度分布测试仪；/n步骤105：再次启动光纤末端长度计算功能，计算被测光纤斯托克斯曲线末端PFE_S为FS_S，被测光纤反斯托克斯曲线末端PFE_AS为FS_AS；/n步骤1...
【技术特征摘要】
1.一种光纤温度分布测试仪光纤长度自动校准及自动测试方法，其特征在于，包括以下步骤：
步骤101：启动光纤温度分布测试仪，断开被测光纤；
步骤102：启动光纤长度自动校准功能，输入光纤实际长度FL，光纤1550nm折射率IR_IN，1550nm折射率IR_IN可由光纤厂家给出，设置反斯托克斯光折射率IR_AS为IR_IN，设置斯托克斯光折射率IR_S为IR_IN，设置空间分辨率参数SFR为1m，采样分辨率SAMF为0.2m，设置测试量程FR＝光纤长度FL+1km，设置内部光纤斯托克斯补偿数据长度DC_S为0，内部光纤反斯托克斯补偿数据长度DC_AS为0，计算测试数据数目N，N的典型值大于等于(光纤长度FL+1km)÷采样分辨率SAMF；
步骤103：启动光纤末端长度计算功能，计算内部光纤斯托克斯曲线末端PFS_S为FS_S，内部光纤反斯托克斯曲线末端PFS_AS为FS_AS；
步骤104：将被测光纤接入光纤温度分布测试仪；
步骤105：再次启动光纤末端长度计算功能，计算被测光纤斯托克斯曲线末端PFE_S为FS_S，被测光纤反斯托克斯曲线末端PFE_AS为FS_AS；
步骤106：读取斯托克斯距离数据DSX[0～N]，该数据通过斯托克斯光折射率IR_S以及采样分辨率SAMF计算得到；，反斯托克斯距离数据DASX[0～N]，该数据通过反斯托克斯光折射率IR_AS以及采样分辨率SAMF计算得到；计算斯托克斯数据光纤测试长度S_DL为DSX[PFE_S]和DSX[PFS_S]的差值，反斯托克斯数据光纤测试长度AS_DL为DASX[PFE_AS]和DASX[PFS_AS]的差值；
步骤107：计算IR_S为(2*S_DL/FL-1)*IR_IN，IR_AS为(2*AS_DL/FL-1)*IR_IN；
步骤108：根据IR_S与IR_AS重新计算DSX[0]到DSX[N]、DSAX[0]到DSAX[N]；
步骤109：设置DC_S为DSX[PFS_S]，DC_AS为DSAX[PFS_AS]；
步骤110：启动光纤末端长度计算功能，计算被测光纤斯托克斯曲线末端PFE_S为FS_S，被测光纤反斯托克斯曲线末端PFE_AS为FS_AS；
步骤111：判断DSX[PFE_S]与FL差值的绝对值是否小于SFR，如果是，转步骤112，否则转步骤115；
步骤112：判断DSAX[PFE_AS]与FL差值的绝对值是否小于SFR，如果是，转步骤113，否则转步骤115；
步骤113：判断DSX[PFE_S]与DSAX[PFE_AS]差值的绝对值是否小于SFR，如果是，转步骤114，否则转步骤115；
步骤114：校准成功，结束；
步骤115：校准失败，检查光路，结束。
2.根据权利要求1所述的一种光纤温度分布测试仪光纤长度自动校准及自动测试方法，其特征在于，步骤103，步骤105和步骤110的光纤长度末端计算功能具体步骤如下所示，通过对斯托克斯测试数据与反斯托克斯数据分析，实现光纤长度末端的准确计算与识别：
步骤:10301：启动温度分布测试功能，等待测试完毕；
步骤10302：读取斯托克斯测试数据DS[0～N]，反斯托克斯测试数据DAS[0～N]，读取测试数据数目N；
步骤10303：计算斯托克斯测试数据噪声DNS为DS[N-500]到DS[N]的平均值，反斯托克斯测试数据噪声DNAS为DAS[N-500]到DAS[N]的平均值；
步骤10304：计算去除噪声斯托克斯数据DS_SN[0]到DS_SN[N]分别为DS[0]-DNS到DS[N]-DNS的值，去除噪声反斯托克斯数据DAS_SN[0]到DAS_SN[N]分别为DAS[0]-DNAS到DAS[N]-DNAS的值；
步骤10305：根据DS_SN[0]到DS_SN[N]的去除噪声斯托克斯数据，DAS_SN[0]到DAS_SN[N]的去除噪声反斯托克斯数据，计算内部光路斯托克斯曲线末端FS_S，反斯托克斯曲线末端FS_AS；
步骤10306：输出斯托克斯光纤长度FS_S，反斯托克斯光纤长度FS_AS，结束。
3.根据权利要求2所述的一种光纤温度分布测试仪光纤长度自动校准及自动测试方法，其特征在于，步骤10305具体包括以下步骤：
步骤201：读取DS_SN[0]到DS_SN[N]的去除噪声斯托克斯数据，DAS_SN[0]到DAS_SN[N]的去除噪声反斯托克斯数据，读取测试空间分辨率参数SFR，采样分辨率SAMF；
步骤202：计算分解层数SSRN为SFR÷SAMF的值向上取整，设置I为0，M为0，SFR典型值为1m，SAMF典型值为0.2m，且SAMF小于SFR，MM为N/SSRN；
步骤203：计算对数分解斯托克斯数据DS_SN_SSR[I][M]为DS_SN[M*SSRN+I]的对数，计算对数分解反斯托克斯数据DAS_SN_SSR[I][M]为DAS_SN[M*SSRN+I]的对数；
步骤204：判断M是否大于MM，如果是，转步骤206，否则转步骤205；
步骤205：设置M的值加1，并转步骤203再次进行计算；
步骤206：判断I是否大于SSRN，如果是，转步骤208，否则转步骤207；
步骤207：设置l的值加1，并转步骤203再次进行计算；
步骤208：设置J为0，I为0，设置末端判定阈值ETH，典型值为2dB/m*SFR；
步骤209：计算对数分解斯托克斯差分数据DS_SN_SSR_DIFF[I][J]为DS_SN_SSR[I][J]和DS_SN_SSR[I][J+1]的差值，计算对数分解反斯托克斯差分数据DAS_SN_SSR_DIFF[I][J]为DAS_SN_SSR[I][J]和DAS_SN_SSR[I][J+1]的差值；
步骤210：判断J是否大于MM-1，如果是，转步骤212，否则转步骤211；
步骤211：设置J的值加1，并转步骤209再次进行计算；
步骤212：记录DS_SN_SSR_DIFF[I][0]到DS_SN_SSR_DIFF[I][MM-1]中所有大于ETH的值，存入波包数组DS_SN_SSR_FEWAVE[I][0]到DS_SN_SSR_FEWAVE[I][MM-1]中，记录DAS_SN_SSR_DIFF[I][0]到DAS_SN_SSR_DIFF[I][MM-1]中所有大于ETH的值，存入波包数组DAS_SN_SSR_FEWAVE[I][0]到DAS_SN_SSR_FEWAVE[I][MM-1]中；
步骤213：判断I是否大于SSRN，如果是，转步骤215，否则转步骤214；
步骤214：设置l的值加1，并转步骤209再次进行计算；
步骤215：根据DS_SN_SSR_FEWAVE[0][0]到DS_SN_SSR_FEWAVE[SSRN][MM-1]、DS_SN_SSR_DIFF[0][0]到DS_SN_SSR_DIFFSSRN][MM-1]与DS_SN_SSR[0][0]到DS_SN_SSR[SSRN][MM]进行计算，得到光纤长度FS_S；
根据判断条件对DAS_SN_SSR_FEWAVE[0][0]到DAS_SN_SSR_FEWAVE[SSRN][MM-1]、DAS_SN_SSR_DIFF[0][0]到DAS_SN_SSR_DIFF[SSRN][MM-1]与DAS_SN_SSR[0][0]到DAS_SN_SSR[SSRN][MM]进行计算，得到光纤长度FS_AS；
步骤216：输出斯托克斯光纤长度FS_S，反斯托克斯光纤长度FS_AS，结束。
4.根据权利要求3所述的一种光纤温度分布测试仪光纤长度自动校准及自动测试方法，其特征在于，步骤212多重差分波包判定数组计算步骤包括：
步骤21201：设置JJ为0，读取I的值，DS_SN_SSR_FEWAVE[I][0]到DS_SN_SSR_FEWAVE[I][MM]全部设置为0；
步骤21202：判断DS_SN_SSR_DIFF[I][JJ]是否大于等于ETH，如果是，转步骤21203，否则转步骤21209；
步骤21203：设置临时计数变量TEMPJ为JJ+1，并转步骤21204；
步骤21204：判断DS_SN_SSR_DIFF[I][TEMPJ]是否大于等于ETH，如果是，转步骤21207，否则转步骤21205；
步骤21205：判断TEMPJ是否大于等于MM-1，如果是，转步骤21207，否则转步骤21206；
步骤21206：设置TEMPJ为TEMPJ+1，并转步骤21204再次进行判断；
步骤21207：设置DS_SN_SSR_FEWAVE[I][JJ]到DS_SN_SSR_FEWAVE[I][TEMPJ-1]全部为1，并转步骤21208；
步骤21208：设置JJ为TEMPJ，并转步骤21202；
步骤21209：判断DS_SN_SSR_DIFF[I][JJ]是否小于等于-1*ETH，如果是，转步骤21210，否则转步骤21215；
步骤21210：设置TEMPJ为JJ+1，并转步骤21211；
步骤21211：判断DS_SN_SSR_DIFF[I][JJ]是否小于等于-1*ETH，如果是，转步骤21214，否则转步骤21212；
步骤21212：判断TEMPJ是否大于等于MM-1，如果是，转步骤21214，否则转步骤21213；
步骤21213：设置TEMPJ的值加1，并转步骤21211再次进行判断；
步骤21214：设置DS_SN_SSR_FEWAVE[I][JJ]到DS_SN_SSR_FEWAVE[I][TEMPJ-1]全部为-1，并转步骤21208；
步骤21215：判断JJ是否大于等于MM-1，如果是，转步骤21217，否则转步骤21216；
步骤21216：设置JJ的值加1，并转步骤21202；
步骤21217：设置JJ为0，读取I的值，DAS_SN_SSR_FEWAVE[I][0]DS_SN_SSR_FEWAVE[I][MM]全部设置为0；
步骤21218：判断DAS_SN_SSR_DIFF[I][JJ]是否大于等于ETH，如果是，转步骤21219，否则转步骤21225；
步骤21219：设置TEMPJ为JJ+1，并转步骤21220；
步骤21220：判断DAS_SN_SSR_DIFF[I][TEMPJ]是否大于等于ETH，如果是，转步骤21223，否则转步骤21221；
步骤21221：判断TEMPJ是否大于等于MM-1，如果是，转步骤21223，否则转步骤21222；
步骤21222：设置TEMPJ为TEMPJ+1，并转步骤21220再次进行判断；
步骤21223：设置DAS_SN_SSR_FEWAVE[I][JJ]到DAS_SN_SSR_FEWAVE[I][TEMPJ-1]全部为1，并转步骤21224；
步骤21224：设置JJ为TEMPJ，并转步骤21218；
步骤21225：判断DAS_SN_SSR_DIFF[I][JJ]是否小于等于-1*ETH，如果是，转步骤21226，否则转步骤21231；
步骤21226：设置TEMPJ为JJ+1，并转步骤21227；
步骤21227：判断DAS_SN_SSR_DIFF[I][JJ]是否小于等于-1*ETH，如果是，转步骤21230，否则转步骤21228；
步骤21228：判断TEMPJ是否大于等于MM-1，如果是，转步骤21230，否则转步骤21229；
步骤21229：设置TEMPJ的值加1，并转步骤21217再次进行判断；
步骤21230：设置DAS_SN_SSR_FEWAVE[I][JJ]到DAS_SN_SS...
【专利技术属性】
技术研发人员：袁明，张洋，闫继送，徐瑞，龚侃，郭洪龙，李鹏，徐玉华，张志辉，毕宗义，
申请(专利权)人：中国电子科技集团公司第四十一研究所，
类型：发明
国别省市：山东;37
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。