发明名称 一种脉冲光的光纤布里渊谱拟合方法
摘要 一种脉冲光的光纤布里渊谱拟合方法,采用布里渊光时域反射计传感测量系统测得传感光纤的布里渊谱频域采样值;对所测布里渊谱采样值{<img file="dest_path_image001.GIF" wi="61" he="27" />,<i>k</i>=1,2…N}采用如下(1)式函数进行拟合,<img file="427151dest_path_image002.GIF" wi="587" he="139" />。本发明可提高布里渊频移的检测精度,可同时获得光纤中声子衰减系数、光纤布里渊频移及光纤的散射脉冲光的相干时间三个参数。
申请公布号 CN104535094B 申请公布日期 2017.02.22
申请号 CN201510020266.7 申请日期 2015.01.14
申请人 华北电力大学(保定) 发明人 杨志;李永倩;尚秋峰;赵丽娟
分类号 G01D5/353(2006.01)I 主分类号 G01D5/353(2006.01)I
代理机构 石家庄冀科专利商标事务所有限公司 13108 代理人 李羡民;高锡明
主权项 一种脉冲光的光纤布里渊谱拟合方法,其特征是,所述方法包括以下步骤:步骤A:采用布里渊光时域反射计(BOTDR)传感测量系统测得传感光纤的布里渊谱频域采样值{(f<sub>k</sub>,S<sub>k</sub>),k=1,2…N},N为采样点数,f<sub>k</sub>为布里渊谱的扫描检测频点,S<sub>k</sub>为对应频点f<sub>k</sub>处的布里渊谱样值,采样频率范围(Δf=f<sub>N</sub>‑f<sub>1</sub>)至少为3倍的布里渊谱半值全宽(即3dB带宽),并使布里渊谱峰值频率位于采样频率范围中心区域,确保采样值包含有一个完整的3dB带宽布里渊谱样值,其中△f为采样频率范围,f<sub>N</sub>为布里渊谱的第N个扫描检测频点,f<sub>1</sub>为布里渊谱的第1个扫描检测频点;步骤B:对所测布里渊谱采样值{(f<sub>k</sub>,S<sub>k</sub>),k=1,2…N}采用下式函数进行拟合,<maths num="0001"><math><![CDATA[<mfenced open = "" close = ""><mtable><mtr><mtd><mrow><msub><mi>S</mi><mrow><mi>S</mi><mi>P</mi></mrow></msub><mrow><mo>(</mo><mi>f</mi><mo>)</mo></mrow><mo>=</mo><mi>A</mi><mfrac><mi>&Gamma;</mi><mrow><msup><mrow><mo>(</mo><mi>&Gamma;</mi><mo>/</mo><mn>2</mn><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><mn>4</mn><msup><mi>&pi;</mi><mn>2</mn></msup><msup><mrow><mo>(</mo><mi>f</mi><mo>-</mo><msub><mi>f</mi><mi>B</mi></msub><mo>)</mo></mrow><mn>2</mn></msup></mrow></mfrac><mo>+</mo><mi>A</mi><mfrac><mrow><mn>2</mn><msup><mi>e</mi><mrow><mo>-</mo><mi>&Gamma;</mi><mo>&CenterDot;</mo><mi>T</mi><mo>/</mo><mn>2</mn></mrow></msup><mo>&lsqb;</mo><msup><mrow><mo>(</mo><mi>&Gamma;</mi><mo>/</mo><mn>2</mn><mo>)</mo></mrow><mn>2</mn></msup><mo>-</mo><mn>4</mn><msup><mi>&pi;</mi><mn>2</mn></msup><msup><mrow><mo>(</mo><mi>f</mi><mo>-</mo><msub><mi>f</mi><mi>B</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>&rsqb;</mo><mi>c</mi><mi>o</mi><mi>s</mi><mo>&lsqb;</mo><mn>2</mn><mi>&pi;</mi><mrow><mo>(</mo><mi>f</mi><mo>-</mo><msub><mi>f</mi><mi>B</mi></msub><mo>)</mo></mrow><mi>T</mi><mo>&rsqb;</mo></mrow><mrow><mi>T</mi><msup><mrow><mo>&lsqb;</mo><msup><mrow><mo>(</mo><mi>&Gamma;</mi><mo>/</mo><mn>2</mn><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><mn>4</mn><msup><mi>&pi;</mi><mn>2</mn></msup><msup><mrow><mo>(</mo><mi>f</mi><mo>-</mo><msub><mi>f</mi><mi>B</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>&rsqb;</mo></mrow><mn>2</mn></msup></mrow></mfrac></mrow></mtd></mtr><mtr><mtd><mrow><mo>-</mo><mi>A</mi><mfrac><mrow><mn>4</mn><msup><mi>&pi;e</mi><mrow><mo>-</mo><mi>&Gamma;</mi><mo>&CenterDot;</mo><mi>T</mi><mo>/</mo><mn>2</mn></mrow></msup><mi>&Gamma;</mi><mrow><mo>(</mo><mi>f</mi><mo>-</mo><msub><mi>f</mi><mi>B</mi></msub><mo>)</mo></mrow><mi>sin</mi><mo>&lsqb;</mo><mn>2</mn><mi>&pi;</mi><mrow><mo>(</mo><mi>f</mi><mo>-</mo><msub><mi>f</mi><mi>B</mi></msub><mo>)</mo></mrow><mi>T</mi><mo>&rsqb;</mo><mo>+</mo><mn>2</mn><mo>&lsqb;</mo><msup><mrow><mo>(</mo><mi>&Gamma;</mi><mo>/</mo><mn>2</mn><mo>)</mo></mrow><mn>2</mn></msup><mo>-</mo><mn>4</mn><msup><mi>&pi;</mi><mn>2</mn></msup><msup><mrow><mo>(</mo><mi>f</mi><mo>-</mo><msub><mi>f</mi><mi>B</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>&rsqb;</mo></mrow><mrow><mi>T</mi><msup><mrow><mo>&lsqb;</mo><msup><mrow><mo>(</mo><mi>&Gamma;</mi><mo>/</mo><mn>2</mn><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><mn>4</mn><msup><mi>&pi;</mi><mn>2</mn></msup><msup><mrow><mo>(</mo><mi>f</mi><mo>-</mo><msub><mi>f</mi><mi>B</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>&rsqb;</mo></mrow><mn>2</mn></msup></mrow></mfrac></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0001135366370000011.GIF" wi="1826" he="428" /></maths>其中,A、Γ、f<sub>B</sub>和T为待拟合参量,A为幅度系数、Γ为声子衰减系数、f<sub>B</sub>为布里渊频移,T为散射脉冲光的相干时间,f为布里渊散射的斯托克斯光或反斯托克斯光与光源频率之差;上式右边第一项为洛仑兹函数,上式右边第二、三项为脉冲光带来的光谱展宽部分;可应用Levenberg‑Marquardt算法进行拟合;步骤C:将拟合完成后所得到的A、Γ、f<sub>B</sub>和T的值代入上式,即为脉冲光的布里渊谱,Γ、f<sub>B</sub>和T的值分别为所测得的光纤中的声子衰减系数、布里渊频移和散射脉冲光的相干时间。
地址 071003 河北省保定市永华北大街619号