孫大鵬，王永亮，王樹樂

(海軍91439部隊96分隊，遼寧 大連 116041)

基于實船抗爆炸沖擊試驗的爆源定位方法

孫大鵬，王永亮，王樹樂

(海軍91439部隊96分隊，遼寧 大連 116041)

對實船抗爆炸沖擊試驗中魚雷爆炸點的幾種定位方法進行了闡述，并結合實船抗某型魚雷爆炸沖擊試驗的實際測量情況，對各種定位方法的應用效果、適用條件進行了綜合比較和分析，提出了實船抗爆炸沖擊試驗爆炸點定位測量的一種有效方法，適用于其他實船抗爆炸沖擊試驗和其他型號魚雷實航打靶試驗，具有通用性和推廣價值。

爆炸試驗;爆源定位;機械拉線定位法;沖擊波零時定位法;沖擊波時延定位法

0 引言

隨著造船工業(yè)的發(fā)展，世界各國對艦船在極端情況下的強度要求越來越高，我國也對艦船抗爆炸沖擊的能力提出了更高的要求。在實船抗爆炸沖擊試驗過程中，由于海上試驗條件的復雜性，靶船位置會受到艦船操控、水流、風力、涌浪等多方面因素的影響，使戰(zhàn)雷攻擊靶船的相對位置產(chǎn)生極大不確定性。戰(zhàn)雷攻擊靶船的位置對于試驗的準確測量、數(shù)據(jù)結果分析、魚雷爆炸威力和艦艇毀傷評定將起到至關重要的作用，也是試驗組織單位、設計生產(chǎn)單位、作戰(zhàn)使用部隊都密切關注的問題。因此，為準確測定實船抗爆炸沖擊試驗中戰(zhàn)雷爆炸的實時位置，確保試驗實施效果，必須對戰(zhàn)雷爆炸點的定位方法進行深入研究。

1 爆源定位原理

在實船抗沖擊試驗中測量的一個關鍵點是確定戰(zhàn)雷的爆炸方位和距靶船的距離，傳統(tǒng)的水下爆炸沖擊波壓力測量中，爆距的確定是通過鋼性固定爆源與測量點的位置來實現(xiàn)，這種機械拉線定位方法對實船抗爆炸沖擊試驗并不適用。對于魚雷實航爆炸時爆源到測點距離的確定，理論上可通過聲測距技術、零時法測距技術及時延定位技術3種技術途徑來實現(xiàn)。

1.1 沖擊波零時法測距技術

這種方法是通過測量沖擊波從爆源傳播到測點的時間來獲取爆距的。測量中的零點信號可利用貼在爆源上的陶瓷片(或利用插入起爆藥內(nèi)的導線)在起爆瞬間給出。計算公式如下［1］:

式中:R為待測距離;R0為爆源的等效藥包半徑;Co為爆源所在水域的聲速;t為傳播時間;K為修正系數(shù)。修正系數(shù)K不僅與測點的沖擊波峰值壓力Pm有關，還與爆炸當量、海域有關。修正系數(shù)必須通過試驗標定來確定。

1.2 沖擊波時延定位技術

沖擊波時延定位在算法上是對沖擊波零時法的增強。此技術是在已知各測點坐標位置的情況下，利用沖擊波到達各個測點的時延差，通過反演技術來實現(xiàn)爆源動態(tài)瞬時定位，它要求傳播到各測點的沖擊波平均速度相同。這一定位方法的精度主要取決于沖擊波到達測點時間的測量精度和反演計算誤差。目前反演算法已經(jīng)很成熟，其計算誤差可通過程序控制，水下爆炸沖擊波到達各測點的時間測量精度可達0.1 ms，其帶來的定位誤差小于0.2 m。因此，這一定位方法的主要誤差是來自沖擊波傳播到各個測點的沖擊波平均速度相等的假定條件，考慮到這一點，沖擊波時延定位技術中的測點布局必須深入研究，在保障現(xiàn)場試驗可操作的情況下，確定測點與靶船的相對距離及各測點間距，并通過現(xiàn)場的小當量水中爆炸試驗來校準驗證布局設計和定位精度［2］。

通過上述基本分析，實船抗沖擊試驗中的爆源定位測量擬采用沖擊波時延定位方法。

2 沖擊波時延定位方法理論分析

沖擊波時延定位技術的前提是假定傳播到各測點的沖擊波傳播速度相等。根據(jù)目前的科技狀況，應用計算機的定位算法及程序已得到廣泛的應用，隨著數(shù)學方法的發(fā)展，建立了爆源參數(shù)的最優(yōu)化修定方法，其中主要包括Geiger法、最小二乘法、阻尼最小二乘法及Powell法。定位計算一般分初定和修定2步進行。初定爆源位置時，僅須根據(jù)測點提供的沖擊波時延，將時延方程線性化，從而解出爆源的位置，即x，y，z和0時 T。

在實船抗沖擊試驗中，由于定位陣測點的坐標會因海流的作用而隨機變化，不能用絕對坐標位置來描述魚雷爆炸點的定位問題。因此，在定位陣上建立動態(tài)坐標系，通過定位來確定相對動態(tài)坐標系的瞬時爆炸位置、沖擊波傳播速度V和起爆時刻T。

2.1 實船抗爆炸沖擊試驗戰(zhàn)雷爆炸點的初定

坐標系中魚雷的爆炸深度z以設定值為標準。在初定爆炸點參數(shù)的基礎上，用廣義逆優(yōu)化方法對參數(shù)進行修定。

設某次爆炸試驗接收到的測點有n個(n≥5)。其中任一測點的坐標為(xi，yi)，相應的爆炸沖擊波時延為Ti。待求的爆源位置為x，y，爆炸深度為z，起爆時刻為T。按照爆炸點到各測點的沖擊波傳播速度相等的假設，可列出下列方程［3－4］:

式中:V為由爆炸點傳播到測點的沖擊波速度。

式(4)為非線性方程，為便于求解，首先應將其線性化，變?yōu)榫€性方程。為此展開式(4)得到:

對于第1個測點(i=1)可得:

將式(5)中 i分別代以 2，3，…，n，與式(6)相減得到:

式(7)構成了初定方程，將初定方程用矩陣的形式表達，即

展開，可得

2.2 戰(zhàn)雷爆炸點參數(shù)的修定

為了取得更精確的定位值，必須對初定結果進行修定。我們采用Powell方法，該方法是一種直接尋找目標函數(shù)極小值的有效方法。

設第i個測點沖擊波時延的觀測值為Toi，由給定的初定爆源參數(shù)求出其理論時延為τi(X)，由此定義的目標函數(shù)為二次型函數(shù):

解決二次型的極值問題，共扼梯度法是一個有效的方法。作為對共扼梯度法改進的Powell方法，求式(10)目標函數(shù)的極小值是很有用的［5－6］。

在一維尋找時，其步長STEP的選取與具體問題中的自變量有關。對于水中爆炸，其單位為米。結束收斂的準則可以采用迭代前后2次的目標函數(shù)的差值達到預定的限差

3 魚雷爆源定位分析與計算

3.1 試驗概況

實船抗某型魚雷爆炸沖擊試驗的測量布陣如圖1所示，共布設3個測量陣，每個測量陣相互獨立，其中包括6個測點組成的定位陣和10個測點組成的自由場沖擊波超壓測量陣。定位陣中心點離爆炸目標區(qū)中心點的垂直距離不小于20 m，初定為30 m，定位陣的長度為12 m，沖擊波超壓測量陣的最近測點離爆炸目標區(qū)中心點的距離不小于10 m，初定為10 m，測量陣的長度為20 m。

圖1 艦船抗某型魚雷爆炸沖擊試驗的測量布陣Fig.1 Ship resistance measurements lineup of a torpedo impulse test explosion

3.2 模擬計算的定位結果

應用前述定位算法，設計計算程序。首先對定位方法和定位程序進行了檢驗，假定動態(tài)坐標系建立在定位陣上，給定定位陣各測點的坐標分別為(0，0，－4)，(0，3，－4)，(0，6，－4)，(0，8，－4)(0，10，－4)，(0，12，－4)。爆源位置初定為(30，5，－4)，沖擊波傳播速度為c=1 500 m/s，起爆時刻為0.0 s。用理論計算得到的各測點沖擊波時延作為觀測時延，對爆源位置進行反演，定位結果如圖2所示。表明此定位算法合理，計算程序可行。

圖2 沖擊波時延定位計算結果Fig.2 Calculation results of shock delay positioning

其次對定位誤差進行分析，對沖擊波時延引入誤差，對各測點測得的沖擊波到時疊加時延的千分之一，得到圖3的定位結果，其定位誤差小于2%。數(shù)值實驗結果表明定位精度受時延誤差的影響較大。

圖3 沖擊波時延定位疊加隨機誤差后的定位結果Fig.3 Shock delay positioning result of composition random errors

從模擬計算的結果可知，本文給出的定位算法是合理可行的，但是對于其定位的真正精度必須通過試驗數(shù)據(jù)來驗證。

4 對魚雷爆炸點的定位測量

以前，國內(nèi)也曾組織過魚雷實航打靶試驗，但試驗中沒有進行爆源定位、壓力等各種參數(shù)的測量工作，只是宏觀上考察了魚雷爆炸效果。因此，沒有可供參考的經(jīng)驗。而近年進行的某型艦船抗爆炸沖擊試驗中針對爆源的定位測量是1次重要的嘗試，科學分析試驗數(shù)據(jù)，對于實船抗沖擊試驗中的戰(zhàn)雷爆源定位方法研究有重要的借鑒作用。

該項試驗各次試驗的定位結果見表1。本文利用前面討論的沖擊波時延定位算法，對爆距重新進行了計算，一并列入表1。

?

機械拉線法用于試驗時對爆源的布放，爆源與被試艦舷側的水平距離通過拉線粗略地估算后進行布放，就基本上滿足要求，但爆源與舷側中心點的橫向偏差較難把握，每次布放時不是偏向艦首就是偏向艦尾，如圖4所示。

由于實船抗爆炸沖擊試驗參試設備多，試驗環(huán)節(jié)多，實施復雜，從爆源布放到起爆時間間隔較長。在此期間，被試艦受海水流速流向影響，其位置可能發(fā)生很大變化，導致爆距產(chǎn)生較大漂移。因此，機械拉線法僅用于爆源布放，并采取有效措施盡量使被試艦艦位穩(wěn)定，才能提高定位精度。該方法更適于湖上靜水試驗條件，并不適用于海上實船抗沖擊試驗。

沖擊波零時法對爆源位置的準確定位，適用于艦船(靜態(tài))抗沖擊試驗的近距離測量。該方法誤差源為爆距測量誤差、陣元坐標測量誤差和解算誤差等，其中以爆距測量誤差的影響最大，并隨著爆距的增加而增大。例如第1次試驗時，沖擊波零時法得到的爆距與布放距離相差較大(布放要求100 m，計算結果為143 m)，其原因是試驗時海浪較大，布放時海流又大，被試艦發(fā)生脫錨，因此船位漂移了40多m。該方法只適用于艦船處于拋錨狀態(tài)，即相對靜止狀態(tài)。當艦船處于航行狀態(tài)時，無法拖帶電纜航行，自由場壓力傳感器也不能布放。

沖擊波時延定位法在零時法的基礎上加以改進，利用Powell方法來求解目標函數(shù)極小值，定位精度有所提高。對于魚雷實航打靶試驗，也同樣可以應用，這一點從表1的計算結果得到驗證。

5 結語

通過對幾種定位方法的比較分析，我們認為對于實船抗爆炸沖擊試驗中魚雷爆炸時刻的位置測量，可以采用沖擊波零時法和沖擊波時延定位法。

在對沖擊波時延定位法的理論計算和實際測量結果的對比中，我們認為該方法在實船抗爆炸沖擊試驗中應用是切實可行的，定位精度滿足工程要求。同時利用實船抗沖擊試驗中的試驗數(shù)據(jù)證明了作者的研究是科學合理的，對于其他實船抗爆炸沖擊試驗和其他型號魚雷實航打靶試驗具有通用性和推廣價值。

［1］庫爾.水下爆炸［M］.羅耀杰，等譯.北京:國防工業(yè)出版社，1960.

［2］程素秋，董俊哲，張姝紅.實船爆炸試驗測量規(guī)模需求綜述［J］.水雷戰(zhàn)與艦船防護，2007，(2):12－15.

CHENG Su-qiu，DONG Jun-ze，ZHANG Shu-hong.Measure scope demand of ship shock trials summarize［J］.Mine Warfare and Ship Self-defence，2007，(2):12 －15.

［3］楊德森，等.基于倒譜分析的水聲信號被動定位時延估計算法研究［J］.系統(tǒng)仿真學報，2009，(2):298 －300.

YANG De-sen，et al.Research of time-delay estimation algorithm for passive localization based on cepstrum analysi［J］.Joumal of system simulation，2009，(2):298 －300.

［4］程素秋，邵宗戰(zhàn)，譚捷.遭受沉底爆炸載荷作用下的圓柱結構沖擊損傷計算［J］.艦船科學技術，2008，30(6):131－135.

CHENG Su-qiu，SHAO Zhong-zhan，TAN Jie.Shock damnification and simulation ofcylindersubjected to seabed underwater explosion loads［J］.Ship Science and Technology，2008，30(6):131 －135.

［5］數(shù)學手冊編寫組.數(shù)學手冊［M］.北京:高等教育出版社，1979.

［6］彭永鎮(zhèn)，曹岳海，陳剛.基于多基元的水下目標高精度測速方法［J］.聲學技術，2008，(6):46－48.

PENG Yong-zhen，CAO Yue-hai，CHEN Gang.A method of high precision underwater target velocity measurement［J］.Thchnical Acoustics，2008，(6):46 －48.

The positioning method of explosion source based on ship resistance impulse test

SUN Da-peng，WANG Yong-liang，WANG Shu-le
(No 91439 Unit of PLA，Dalian 116041，China)

Several positioning methods of torpedo stopping point in impulse test of ship resistance of underwater explosion are briefly described in this paper，and combined with some type of ship resistance to impulse test of a torpedo explosion of practical measurement situation，application effect of various positioning method，and that the applicable conditions comprehensive are compared and analyzed，then we put forward the ship impulse test and defused positioning measurement anti-explosion the general method.It is suitable for other ship impulse test and other models anti-explosion real navigation targeting test，which has versatile and popularize value.

explosion test;explosion source positioning;mechanical cable positioning method;shockwave flinging positioning method;shockwave delay positioning method

O382+.1

A

1672－7649(2012)04－0036－04

10.3404/j.issn.1672－7649.2012.04.008

2011－06－14;

2011－07－18

孫大鵬(1962－)，男，高級工程師，從事水中兵器試驗測量工作。