高曉文 張自麗 葉博 葛輝良

（第七一五研究所，杭州，310023）

光纖水聽器具有傳感與傳輸一體、濕端無電子器件，不易受電磁干擾、適用于大規(guī)模陣列復(fù)用等優(yōu)點(diǎn)，已應(yīng)用于拖曳線列陣聲吶、寬孔徑陣舷側(cè)陣聲吶、海底布放岸基聲吶等[1-3]。光纖水聽器一般采用干涉式結(jié)構(gòu)，利用待測(cè)聲壓對(duì)干涉儀兩臂的調(diào)制將聲壓信號(hào)轉(zhuǎn)換為光學(xué)相位信息，并通過光學(xué)相干檢測(cè)的方法獲取聲壓信息，因此，信號(hào)檢測(cè)技術(shù)是光纖水聽器的關(guān)鍵技術(shù)之一[4,5]。目前，國內(nèi)的光纖水聽器研制單位普遍采用相位調(diào)制載波（Phase Generated Carrier，PGC）調(diào)制解調(diào)方案作為光纖水聽器的信號(hào)檢測(cè)手段，通過引入載波對(duì)待測(cè)信號(hào)進(jìn)行調(diào)制，使得待測(cè)信號(hào)位于調(diào)制信號(hào)的邊帶上，再利用微分交叉相乘（Differentiate Cross Multiply，DCM）算法得到待測(cè)信號(hào)[6-8]。為了避免邊帶間的信號(hào)頻率混疊，待測(cè)信號(hào)的頻率以及動(dòng)態(tài)范圍受制于載波的頻率。與之相比，3×3耦合器解調(diào)方案不需要引入載波調(diào)制，利用3路輸出信號(hào)之間的相位關(guān)系結(jié)合數(shù)字信號(hào)處理技術(shù)即可獲取待測(cè)信號(hào)，實(shí)現(xiàn)簡(jiǎn)單，運(yùn)算量小，待測(cè)信號(hào)的頻率以及動(dòng)態(tài)范圍受制于采樣頻率，動(dòng)態(tài)范圍較PGC解調(diào)方案大[9,10]。

1 3×3耦合器解調(diào)算法優(yōu)化

1.1 解調(diào)原理

3×3耦合器解調(diào)方案一般采用圖1中的Michelson干涉儀結(jié)構(gòu)實(shí)現(xiàn)。激光器輸出光經(jīng)環(huán)形器進(jìn)入3×3耦合器，耦合器的兩臂輸出分別作為參考臂及信號(hào)臂，分別經(jīng)法拉第鏡反射，返回信號(hào)再次通過3×3耦合器分3路輸出，由3路光電探測(cè)器進(jìn)行光電轉(zhuǎn)換。

圖1 3×3耦合器解調(diào)方案的Michelson實(shí)現(xiàn)

理論上，3×3耦合器輸出的3路信號(hào)僅存在120°的相位差，可表示為：

式中，A為輸出信號(hào)的直流量；B為輸出信號(hào)的交流量系數(shù)；為待測(cè)聲信號(hào)，φn為噪聲，φ0初始相位。

將3路輸出信號(hào)相加得到干涉信號(hào)的直流量：

分別將3路信號(hào)減去直流量，得到：

求上式中各項(xiàng)平方和，得到交流量系數(shù)B：

利用交流量系數(shù)對(duì)式（3）進(jìn)行歸一化，得到3路信號(hào)：

將式（5）的兩個(gè)信號(hào)微分相減，并與另一個(gè)信號(hào)相乘，得到：

同理分別得到另兩路信號(hào)：

式（6）～（8）相加得到：

積分后經(jīng)高通，即得到所需信號(hào)：

3×3耦合器解調(diào)方案的實(shí)現(xiàn)可分為兩部分，一是信號(hào)調(diào)理，即對(duì)3路輸出信號(hào)進(jìn)行去直流及交流量歸一化處理；二是利用3路輸出信號(hào)的相位關(guān)系進(jìn)行待測(cè)信號(hào)的解調(diào)。

1.2 解調(diào)算法的缺陷分析及優(yōu)化

根據(jù)1.1節(jié)分析， 3×3耦合器的理想解調(diào)基于以下3個(gè)前提條件：（1）3路信號(hào)相位差120°；（2）3路信號(hào)直流量系數(shù)相等；（3）3路信號(hào)交流量系數(shù)相等。在實(shí)際應(yīng)用中，由于工藝限制，無法保證3×3耦合器結(jié)構(gòu)的理想對(duì)稱，即無法實(shí)現(xiàn)通道間精確的120°相位差以及1：1：1的分光比。此外，由于整個(gè)系統(tǒng)采用了非互易結(jié)構(gòu)光路，且不采用全保偏光路，使得3路光信號(hào)在系統(tǒng)傳輸過程中的偏振態(tài)及傳輸損耗的變化互不相同，即無法保證干涉條紋對(duì)稱性及穩(wěn)定性。因此，實(shí)現(xiàn)3×3耦合器理想解調(diào)的3個(gè)條件均不成立[11]。

在實(shí)際應(yīng)用中，需要對(duì)3×3耦合器的解調(diào)方案作優(yōu)化，以近似滿足3個(gè)前提條件。本文采用一種改進(jìn)的3×3耦合器解調(diào)方案，利用3×3耦合器輸出信號(hào)間的相位差，任意兩路信號(hào)可以形成利薩如圖，根據(jù)利薩如圖特點(diǎn)構(gòu)建數(shù)學(xué)模型，求取3路輸出通道的各項(xiàng)系數(shù)，完成解調(diào)前信號(hào)調(diào)節(jié)步驟。圖2為改進(jìn)后的3×3耦合器解調(diào)方案原理框圖。

圖2 改進(jìn)的3×3耦合器解調(diào)方案原理框圖

實(shí)際光路輸出信號(hào)為：

式中，A1、A2、A3為3×3耦合器3路輸出信號(hào)的直流量；B1、B2、B3為3×3耦合器3路輸出信號(hào)的交流量系數(shù)；β1為通道1和通道2兩路信號(hào)之間的相位差，β2為通道1和通道3兩路信號(hào)之間的相位差。聯(lián)立通道1的輸出I1與通道2的輸出I2，消除φ得到：

上式是一個(gè)橢圓方程表達(dá)式，對(duì)比標(biāo)準(zhǔn)橢圓方程：

可以得到通道1和通道2的輸出干涉信號(hào)中各參數(shù)的表達(dá)式與橢圓擬合式中各系數(shù)的對(duì)應(yīng)關(guān)系為：

通過最小二乘法擬合的方式獲取I1與I2聯(lián)立的橢圓曲線，根據(jù)擬合曲線得到橢圓系數(shù)代入式（14）得到I1與I2的各項(xiàng)系數(shù)[12]。分別對(duì)3×3耦合器的3路干涉信號(hào)兩兩組合進(jìn)行上述計(jì)算，即可得到各項(xiàng)干涉系數(shù)、直信號(hào)，利用交流量系數(shù)進(jìn)行歸一化：

至此，得到了如式（1）中的3路信號(hào)，依次通過式（2）～（10）的解算，得到聲信號(hào)φS。

2 仿真

為了驗(yàn)證解調(diào)算法的有效性，利用Matlab進(jìn)行了3×3耦合器的解調(diào)仿真，仿真信號(hào)設(shè)置如下：

式中，Ni為輸出信號(hào)上噪聲擾動(dòng)，Pi為干涉系統(tǒng)的相位噪聲，為模擬聲信號(hào)，C為模擬聲信號(hào)幅值，f為聲信號(hào)頻率。仿真信號(hào)的輸出如圖3所示。

圖3 3×3耦合器輸出信號(hào)仿真圖

根據(jù)算法要求，3路輸出兩兩組合進(jìn)行橢圓擬合，其中通道1和通道2輸出信號(hào)的擬合結(jié)果如圖4所示。

表1是經(jīng)橢圓擬合方法得到的信號(hào)系數(shù)與仿真信號(hào)系數(shù)的關(guān)系?？梢钥吹剑?jīng)過橢圓擬合得到的通道參數(shù)，其誤差值均小于0.1%。

表1 通道參數(shù)擬合結(jié)果

經(jīng)解調(diào)得到的信號(hào)與仿真信號(hào)的對(duì)比如圖5所示，幅值差異小于0.1%。

圖5 仿真信號(hào)與解調(diào)信號(hào)對(duì)比

3 試驗(yàn)

3.1 測(cè)試系統(tǒng)設(shè)置

在仿真測(cè)試的基礎(chǔ)上，搭建了一套光纖水聽器3×3耦合器數(shù)字解調(diào)系統(tǒng)，系統(tǒng)設(shè)置如圖6所示。光源采用RIO公司生產(chǎn)的超窄線寬激光器，線寬小于3 kHz，輸出經(jīng)環(huán)形器接入3×3耦合器，耦合器的兩臂分別為水聽器的參考臂及傳感臂，其中傳感臂光纖繞制在PZT環(huán)上，由信號(hào)源加載調(diào)制以產(chǎn)生模擬聲信號(hào)，參考臂作去敏處理，兩臂光程差5 m。使用New_Focus公司的光電探測(cè)器進(jìn)行光電轉(zhuǎn)換，并由上位機(jī)采集數(shù)據(jù)，采樣率500 kHz。

圖6 3×3耦合器數(shù)字解調(diào)系統(tǒng)原理框圖

為降低環(huán)境干擾的影響，光源及干涉儀的兩臂放置于隔振隔聲桶內(nèi)。試驗(yàn)系統(tǒng)如圖7所示。

圖7 3×3耦合器數(shù)字解調(diào)試驗(yàn)系統(tǒng)

3.2 測(cè)試結(jié)果

通過信號(hào)源對(duì)干涉儀加載正弦波模擬聲信號(hào)，頻率為1 300 Hz，信號(hào)幅值29.2 rad，對(duì)應(yīng)功率譜密度為26.2 dB/Hz。圖8給出了經(jīng)3×3耦合器解調(diào)得到的模擬聲信號(hào)的頻譜結(jié)果，在2 600 Hz、3 900 Hz處均出現(xiàn)了諧波，諧波幅值分別為-50.8 dB、-46.4 dB，2 kHz處噪聲本底為-95 dB，噪聲水平符合理論預(yù)期。

圖8 3×3耦合器數(shù)字解調(diào)實(shí)測(cè)結(jié)果

由于3×3耦合器的非完美對(duì)稱性，解調(diào)時(shí)存在信號(hào)諧波，定義10次諧波的合功率與信號(hào)功率的比值為信號(hào)失真度，用以判斷是否正常解調(diào)。當(dāng)信號(hào)失真度小于-50 dB時(shí)，判斷為正常解調(diào)。

根據(jù)失真度定義進(jìn)行了動(dòng)態(tài)范圍測(cè)試，得到各頻率下的動(dòng)態(tài)范圍與理論動(dòng)態(tài)范圍，如圖9所示。以2 kHz頻率點(diǎn)為例，可檢測(cè)信號(hào)幅值29.5 dB，噪聲本底以-95 dB計(jì)算，動(dòng)態(tài)范圍為124.5 dB。

圖9 0～4 kHz內(nèi)動(dòng)態(tài)范圍實(shí)測(cè)結(jié)果

4 結(jié)論

本文開展了光纖水聽器3×3耦合器解調(diào)方案的研究，優(yōu)化了光纖水聽器3×3耦合器解調(diào)方案，增加了信號(hào)預(yù)處理調(diào)節(jié)步驟，該方法利用耦合器3路輸出信號(hào)間的利薩如圖，結(jié)合橢圓擬合的方法獲取各通道的系數(shù)，完成對(duì)信號(hào)的預(yù)處理，最后利用微分交叉相乘法提取相位信號(hào)。仿真結(jié)果驗(yàn)證了該方案的有效性，耦合器3路輸出通道的參數(shù)解算結(jié)果與仿真信號(hào)的數(shù)值誤差優(yōu)于0.1%，信號(hào)解調(diào)誤差優(yōu)于0.1%。搭建了實(shí)驗(yàn)室環(huán)境下光纖水聽器3×3耦合器數(shù)字解調(diào)系統(tǒng)，經(jīng)測(cè)試，在水聽器兩臂光程差為5 m，系統(tǒng)采樣率為500 kHz的條件下，實(shí)現(xiàn)2 kHz處信號(hào)解調(diào)的動(dòng)態(tài)范圍124.5 dB，有效提升了光纖水聽器的信號(hào)獲取能力。