主权项 |
一种脉冲光的光纤布里渊谱拟合方法,其特征是,所述方法包括以下步骤:步骤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>Γ</mi><mrow><msup><mrow><mo>(</mo><mi>Γ</mi><mo>/</mo><mn>2</mn><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><mn>4</mn><msup><mi>π</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>Γ</mi><mo>·</mo><mi>T</mi><mo>/</mo><mn>2</mn></mrow></msup><mo>[</mo><msup><mrow><mo>(</mo><mi>Γ</mi><mo>/</mo><mn>2</mn><mo>)</mo></mrow><mn>2</mn></msup><mo>-</mo><mn>4</mn><msup><mi>π</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>]</mo><mi>c</mi><mi>o</mi><mi>s</mi><mo>[</mo><mn>2</mn><mi>π</mi><mrow><mo>(</mo><mi>f</mi><mo>-</mo><msub><mi>f</mi><mi>B</mi></msub><mo>)</mo></mrow><mi>T</mi><mo>]</mo></mrow><mrow><mi>T</mi><msup><mrow><mo>[</mo><msup><mrow><mo>(</mo><mi>Γ</mi><mo>/</mo><mn>2</mn><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><mn>4</mn><msup><mi>π</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>]</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>πe</mi><mrow><mo>-</mo><mi>Γ</mi><mo>·</mo><mi>T</mi><mo>/</mo><mn>2</mn></mrow></msup><mi>Γ</mi><mrow><mo>(</mo><mi>f</mi><mo>-</mo><msub><mi>f</mi><mi>B</mi></msub><mo>)</mo></mrow><mi>sin</mi><mo>[</mo><mn>2</mn><mi>π</mi><mrow><mo>(</mo><mi>f</mi><mo>-</mo><msub><mi>f</mi><mi>B</mi></msub><mo>)</mo></mrow><mi>T</mi><mo>]</mo><mo>+</mo><mn>2</mn><mo>[</mo><msup><mrow><mo>(</mo><mi>Γ</mi><mo>/</mo><mn>2</mn><mo>)</mo></mrow><mn>2</mn></msup><mo>-</mo><mn>4</mn><msup><mi>π</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>]</mo></mrow><mrow><mi>T</mi><msup><mrow><mo>[</mo><msup><mrow><mo>(</mo><mi>Γ</mi><mo>/</mo><mn>2</mn><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><mn>4</mn><msup><mi>π</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>]</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的值分别为所测得的光纤中的声子衰减系数、布里渊频移和散射脉冲光的相干时间。 |