綦 磊，樊 帆，孫立臣，李明利，于兆吉

(1.北京衛(wèi)星環(huán)境工程研究所,北京100094；2.中國空間技術(shù)研究院,北京100094；3.天津大學(xué)精密測試技術(shù)及儀器國家重點實驗室,天津300072)

1 引言

隨著在軌航天器數(shù)量的增加,空間碎片的數(shù)量也與日俱增,一旦與航天器發(fā)生碰撞,必將造成航天器損傷,甚至造成航天器艙體泄漏,嚴(yán)重威脅航天器的正常運(yùn)行和航天員的生命安全［1-2］。因此,盡早判斷碰撞及泄漏是否發(fā)生,并對碰撞及泄漏源進(jìn)行定位具有重要意義［3］。

近年來,國內(nèi)外許多學(xué)者對航天器在軌碰撞及泄漏定位展開研究,并取得了一系列成果。龐寶君等［4-5］研究了空間碎片碰撞產(chǎn)生的聲發(fā)射信號的模態(tài)特征,并利用模態(tài)聲達(dá)時差法實現(xiàn)碎片碰撞的定位以及碰撞損傷等級的識別,但該方法無法解決泄漏定位問題。劉武鋼等［6-7］研究了多孔鋁板高速撞擊聲發(fā)射定位方法,并分析了高速撞擊聲發(fā)射信號在鋁板中的衰減規(guī)律,實現(xiàn)了碰撞的有效定位。孫偉等［8］研制了一種基于可視化超聲的便攜式在軌檢漏儀,可以對泄漏源進(jìn)行定位,但該儀器需要航天員操作,無法實現(xiàn)實時監(jiān)測。Holland等［9-10］提出了一種基于相關(guān)法的泄漏定位方法,該方法可以有效降低定位算法對硬件的要求,但檢測時間較長,無法對碰撞這種突發(fā)信號源進(jìn)行定位。針對上述問題,本文提出一種基于兩級陣列傳感器的碰撞及泄漏定位方法,以實現(xiàn)對碰撞及泄漏的實時監(jiān)測,并具有較高的定位精度。

2 碰撞及泄漏定位原理

本方法應(yīng)用兩級傳感器陣列,一級傳感器陣列是由多個陣元組成,將一級傳感器陣列按照一定的拓?fù)浣Y(jié)構(gòu)粘貼于航天器艙壁上組成二級傳感器陣列,如圖1所示。每一個一級傳感器陣列只有1號陣元一直處于工作狀態(tài),其余陣元處于待機(jī)狀態(tài),只有1號陣元采集到過閾值的信號才工作,這樣可以有效降低檢測系統(tǒng)的功耗。

碰撞產(chǎn)生的聲發(fā)射信號為突發(fā)型信號,如圖2所示,每一個碰撞信號都包含縱波和橫波,縱波波速較快,幅值較小,橫波波速較慢,幅值較大,經(jīng)過一定傳播距離后,縱波和橫波分離開來。因此,合理設(shè)置信號門限值,可以得到縱波或者橫波的到達(dá)時刻,從而根據(jù)相應(yīng)的波速利用聲達(dá)時差法(TDOA)對碰撞點進(jìn)行定位。

TDOA定位原理如圖3所示［11］：假設(shè)三個一級傳感器陣列的1號陣元分別位于A(xA,yA)、B(xB,yB)和C(xC,yC)三個位置,同步接收碰撞點T(xT,yT)發(fā)出的聲發(fā)射信號,接收時刻分別記為tAT、tBT、tCT,計算出時間差 tAB＝ ｜tAT?tBT｜和 tAC＝ ｜tAT?tCT｜,由時間差乘以聲波的傳播速度v,就可以得到距離差dAB和dAC如式(1)～(2)。由于已知tAT、tBT、tCT三者的大小關(guān)系,因此式(1)和(2)各代表雙曲線的一支MN和RS,聯(lián)立式(1)和(2),可以求解碰撞點的坐標(biāo)(xT,yT)。

泄漏信號屬于連續(xù)聲信號,無法通過設(shè)置信號門限值來直接獲得到達(dá)時間,傳統(tǒng)TDOA不能對此種信號進(jìn)行定位［12］。但連續(xù)信號聲源指向性特征明顯,可以基于一級傳感器陣列應(yīng)用波束形成的方法對泄漏源進(jìn)行定向,并利用二級傳感器陣列進(jìn)行泄漏定位,其原理如圖4所示。

波束形成器的原理是對傳感器陣列中各陣元采集的信號進(jìn)行時域平移再求和,得到增強(qiáng)信號輸出。同一泄漏源產(chǎn)生的泄漏信號具有一致性,而噪聲信號具有隨機(jī)性,因此波束形成器可通過信號疊加技術(shù)來提高信噪比［10］,流程如圖5所示。

假設(shè)一級傳感器陣列由n個陣元組成,1號陣元為參考陣元,其接收到的泄漏聲信號為ψ(t),那么陣列中其他陣元i接收到的信號如式(3)：

式中,α(R)是與距離R有關(guān)的衰減因子；ΣΨm代表噪聲及反射回波；Δti代表第i陣元相對于參考陣元的時間時延。將各陣元信號進(jìn)行合成并求聲功率可以得到式(4)：

式中,時延Δti取決于傳感器陣列的空間分布、聲波入射到陣列的角度θ及聲波傳播速度v。對于特定的傳感器陣列及航天器艙體結(jié)構(gòu),傳感器陣列各陣元位置參數(shù)及波速v均為固定值,因此,聲功率僅為角度 θ的函數(shù)［13］。 θ從0°到360°掃描,可以獲得角度與聲功率關(guān)系曲線,如圖6所示。聲功率最大處對應(yīng)的角度即為泄漏聲波向傳感器陣列入射的角度。

利用兩個一級傳感器陣列,分別計算出泄漏聲源入射一級陣列傳感器的角度θ1和θ2,結(jié)合圖4中的幾何關(guān)系,泄漏點坐標(biāo)(x,y)滿足式(5),求解可得到泄漏點坐標(biāo)。

3 試驗與數(shù)據(jù)分析

利用激光碎片發(fā)生器產(chǎn)生5 km/s的微小碎片,碰撞到100 cm×100 cm金屬鋁平板上模擬空間碎片碰撞航天器的情況,如圖7所示。利用三個一級陣列傳感器中的組成二級傳感器陣列進(jìn)行定位,以金屬鋁板中心為坐標(biāo)原點,建立坐標(biāo)系,三個一級陣列傳感器所在坐標(biāo)分別為(－40 cm,－40 cm),(40 cm,－40 cm),(0 cm,40 cm),取門限電平為20 mV(本底噪聲為5 mV),定位結(jié)果如表1所示。

表1 碰撞定位結(jié)果Table 1 Results of impact locating

從上表可以看出,在1 m3檢測范圍內(nèi),碰撞定位的相對誤差小于2%。

在100 cm×100 cm金屬鋁平板上事先制作若干直徑1 mm、1.5 mm、2.0 mm的漏孔,真空泵連接波紋管對漏孔抽氣模擬航天器在軌泄漏的情況,如圖8所示。利用二個一級陣列傳感器中的組成二級傳感器陣列進(jìn)行定位,以金屬鋁板中心為坐標(biāo)原點,建立坐標(biāo)系,二個一級陣列傳感器所在坐標(biāo)分別為(－40 cm,－40 cm),(40 cm,－40 cm),定位結(jié)果如表2所示。

表2 泄漏定位結(jié)果Table 2 Results of leak locating

從上表可以看出,在1 m2檢測范圍內(nèi),泄漏定位的相對誤差小于10%。

4 結(jié)論

1)這種方法基于聲發(fā)射原理,可以實時檢測碰撞及泄漏事件的發(fā)生；

2)這種方法能夠?qū)崿F(xiàn)在1 m2的范圍內(nèi),碰撞定位相對誤差小于2%,直徑大于1 mm漏孔泄漏定位相對誤差小于10%,滿足在軌航天器碰撞及泄漏檢測的需求。