发明名称 基于射线声学的三维倾斜海底参数快速测量方法
摘要 本发明提供的是一种基于射线声学的三维倾斜海底参数快速测量方法。声源发出的声信号S(t)经由海底反射,声信号传播到水听器,4个水听器的接收信号为:A1(t)、A2(t)、A3(t)、A4(t),4个水听器信号的到达时刻为:t1、t2、t3、t4,将每个水听器的到达时刻及深度信息代入各自满足的方程,共计28个方程,求解方程组,获得海底参数、即cb;α;h;n;θ,求解的另外24个变量为声线轨迹的位置信息。本发明方法可以实现三维倾斜海底情况下海底参数测量,在实际测量中仅需测量若干条声线到达时间,较其他方法的反演模型简单,计算速度快,精度高,实验方法简便可行。
申请公布号 CN102508247B 申请公布日期 2013.07.31
申请号 CN201110320196.9 申请日期 2011.10.20
申请人 哈尔滨工程大学 发明人 杨士莪;宋扬;黄益旺;朴胜春;朴健浩
分类号 G01S15/02(2006.01)I;G01S5/18(2006.01)I 主分类号 G01S15/02(2006.01)I
代理机构 代理人
主权项 一种基于射线声学的三维倾斜海底参数快速测量方法,其特征是:声源发出的声信号S(t)经由海底反射,声信号传播到水听器,4个水听器的接收信号为:A1(t)、A2(t)、A3(t)、A4(t),4个水听器信号的到达时刻为:t1、t2、t3、t4,将每个水听器的到达时刻及深度信息代入各自满足的方程,共计28个方程,求解方程组,获得海底参数、即cb;α;h;n;θ,求解的另外24个变量为声线轨迹的位置信息;所述各自满足的方程包括:k1ρ1cosβ1=rcosθ+k2ρ2cosβ2          (1)k1ρ1sinβ1=rsinθ+k2ρ2sinβ2          (2) <mrow> <mfrac> <msup> <msub> <mi>&rho;</mi> <mn>1</mn> </msub> <mn>2</mn> </msup> <mrow> <msup> <msub> <mi>&rho;</mi> <mn>1</mn> </msub> <mn>2</mn> </msup> <mo>+</mo> <msup> <mi>H</mi> <mn>2</mn> </msup> </mrow> </mfrac> <mo>=</mo> <mfrac> <mrow> <msup> <mi>n</mi> <mn>2</mn> </msup> <msup> <mrow> <mo>(</mo> <msub> <mi>k</mi> <mn>1</mn> </msub> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> <mn>2</mn> </msup> <msup> <msub> <mi>&rho;</mi> <mn>1</mn> </msub> <mn>2</mn> </msup> </mrow> <mrow> <msup> <mrow> <mo>(</mo> <msub> <mi>k</mi> <mn>1</mn> </msub> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> <mn>2</mn> </msup> <msup> <msub> <mi>&rho;</mi> <mn>1</mn> </msub> <mn>2</mn> </msup> <mo>+</mo> <msup> <mrow> <mo>(</mo> <msub> <mi>k</mi> <mn>1</mn> </msub> <msub> <mi>&rho;</mi> <mn>1</mn> </msub> <mi>sin</mi> <msub> <mi>&beta;</mi> <mn>1</mn> </msub> <mi>tan</mi> <mi>&alpha;</mi> <mo>+</mo> <mi>h</mi> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>3</mn> <mo>)</mo> </mrow> </mrow> <mrow> <mfrac> <msup> <msub> <mi>&rho;</mi> <mn>2</mn> </msub> <mn>2</mn> </msup> <mrow> <msup> <msub> <mi>&rho;</mi> <mn>2</mn> </msub> <mn>2</mn> </msup> <mo>+</mo> <msup> <mi>Z</mi> <mn>2</mn> </msup> </mrow> </mfrac> <mo>=</mo> <mfrac> <mrow> <msup> <mi>n</mi> <mn>2</mn> </msup> <msup> <mrow> <mo>(</mo> <msub> <mi>k</mi> <mn>2</mn> </msub> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> <mn>2</mn> </msup> <msup> <msub> <mi>&rho;</mi> <mn>2</mn> </msub> <mn>2</mn> </msup> </mrow> <mrow> <msup> <mrow> <mo>(</mo> <msub> <mi>k</mi> <mn>2</mn> </msub> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> <mn>2</mn> </msup> <msup> <msub> <mi>&rho;</mi> <mn>2</mn> </msub> <mn>2</mn> </msup> <mo>+</mo> <msup> <mrow> <mo>(</mo> <msub> <mi>k</mi> <mn>1</mn> </msub> <msub> <mi>&rho;</mi> <mn>1</mn> </msub> <mi>sin</mi> <msub> <mi>&beta;</mi> <mn>1</mn> </msub> <mi>tan</mi> <mi>&alpha;</mi> <mo>+</mo> <mi>h</mi> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>4</mn> <mo>)</mo> </mrow> </mrow> <mrow> <mfrac> <mrow> <mrow> <mo>(</mo> <msub> <mi>k</mi> <mn>1</mn> </msub> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> <msub> <mi>&rho;</mi> <mn>1</mn> </msub> <mi>sin</mi> <msub> <mi>&beta;</mi> <mn>1</mn> </msub> <mi>tan</mi> <mi>&alpha;</mi> <mo>-</mo> <mrow> <mo>(</mo> <msub> <mi>k</mi> <mn>1</mn> </msub> <msub> <mi>&rho;</mi> <mn>1</mn> </msub> <mi>sin</mi> <msub> <mi>&beta;</mi> <mn>1</mn> </msub> <mi>tan</mi> <mi>&alpha;</mi> <mo>+</mo> <mi>h</mi> <mo>)</mo> </mrow> </mrow> <msup> <mrow> <mo>[</mo> <msup> <mrow> <mo>(</mo> <msub> <mi>k</mi> <mn>1</mn> </msub> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> <mn>2</mn> </msup> <msup> <msub> <mi>&rho;</mi> <mn>1</mn> </msub> <mn>2</mn> </msup> <mo>+</mo> <msup> <mrow> <mo>(</mo> <msub> <mi>k</mi> <mn>1</mn> </msub> <msub> <mi>&rho;</mi> <mn>1</mn> </msub> <mi>sin</mi> <msub> <mi>&beta;</mi> <mn>1</mn> </msub> <mi>tan</mi> <mi>&alpha;</mi> <mo>+</mo> <mi>h</mi> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>]</mo> </mrow> <mrow> <mn>1</mn> <mo>/</mo> <mn>2</mn> </mrow> </msup> </mfrac> <mo>=</mo> <mfrac> <mrow> <mrow> <mo>(</mo> <msub> <mi>k</mi> <mn>2</mn> </msub> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> <msub> <mi>&rho;</mi> <mn>2</mn> </msub> <mi>sin</mi> <msub> <mi>&beta;</mi> <mn>2</mn> </msub> <mi>tan</mi> <mi>&alpha;</mi> <mo>-</mo> <mrow> <mo>(</mo> <msub> <mi>k</mi> <mn>1</mn> </msub> <msub> <mi>&rho;</mi> <mn>1</mn> </msub> <mi>sin</mi> <msub> <mi>&beta;</mi> <mn>1</mn> </msub> <mi>tan</mi> <mi>&alpha;</mi> <mo>+</mo> <mi>h</mi> <mo>)</mo> </mrow> </mrow> <msup> <mrow> <mo>[</mo> <msup> <mrow> <mo>(</mo> <msub> <mi>k</mi> <mn>2</mn> </msub> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> <mn>2</mn> </msup> <msup> <msub> <mi>&rho;</mi> <mn>2</mn> </msub> <mn>2</mn> </msup> <mo>+</mo> <msup> <mrow> <mo>(</mo> <msub> <mi>k</mi> <mn>1</mn> </msub> <msub> <mi>&rho;</mi> <mn>1</mn> </msub> <mi>sin</mi> <msub> <mi>&beta;</mi> <mn>1</mn> </msub> <mi>tan</mi> <mi>&alpha;</mi> <mo>+</mo> <mi>h</mi> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>]</mo> </mrow> <mrow> <mn>1</mn> <mo>/</mo> <mn>2</mn> </mrow> </msup> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>5</mn> <mo>)</mo> </mrow> </mrow> <mrow> <mi>tan</mi> <mi>&alpha;</mi> <mrow> <mo>(</mo> <msub> <mi>k</mi> <mn>1</mn> </msub> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> <msqrt> <msup> <mi>n</mi> <mn>2</mn> </msup> <mrow> <mo>(</mo> <msup> <msub> <mi>&rho;</mi> <mn>1</mn> </msub> <mn>2</mn> </msup> <mo>+</mo> <msup> <mi>H</mi> <mn>2</mn> </msup> <mo>)</mo> </mrow> <mo>-</mo> <msup> <msub> <mi>&rho;</mi> <mn>1</mn> </msub> <mn>2</mn> </msup> </msqrt> <mrow> <mo>(</mo> <mi>r</mi> <mi>cos</mi> <mi>&theta;</mi> <mo>+</mo> <msub> <mi>&rho;</mi> <mn>2</mn> </msub> <mi>cos</mi> <msub> <mi>&beta;</mi> <mn>2</mn> </msub> <mo>-</mo> <msub> <mi>&rho;</mi> <mn>1</mn> </msub> <mi>cos</mi> <msub> <mi>&beta;</mi> <mn>1</mn> </msub> <mo>)</mo> </mrow> <mo>=</mo> <mrow> <mo>(</mo> <msub> <mi>k</mi> <mn>1</mn> </msub> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> <mrow> <mo>(</mo> <msub> <mi>k</mi> <mn>2</mn> </msub> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> <msub> <mi>&rho;</mi> <mn>1</mn> </msub> <msub> <mi>&rho;</mi> <mn>2</mn> </msub> <mi>sin</mi> <mrow> <mo>(</mo> <msub> <mi>&beta;</mi> <mn>2</mn> </msub> <mo>-</mo> <msub> <mi>&beta;</mi> <mn>1</mn> </msub> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>6</mn> <mo>)</mo> </mrow> </mrow> <mrow> <mi>t</mi> <mo>=</mo> <mfrac> <mn>1</mn> <msub> <mi>c</mi> <mn>0</mn> </msub> </mfrac> <mo>{</mo> <msqrt> <msubsup> <mi>&rho;</mi> <mn>1</mn> <mn>2</mn> </msubsup> <mo>+</mo> <msup> <mi>H</mi> <mn>2</mn> </msup> </msqrt> <mo>+</mo> <msqrt> <msubsup> <mi>&rho;</mi> <mn>2</mn> <mn>2</mn> </msubsup> <mo>+</mo> <msup> <mi>Z</mi> <mn>2</mn> </msup> </msqrt> <mo>}</mo> <mo>+</mo> <mfrac> <mn>1</mn> <msub> <mi>c</mi> <mi>b</mi> </msub> </mfrac> <mo>{</mo> <msqrt> <msup> <mrow> <mo>(</mo> <msub> <mi>k</mi> <mn>1</mn> </msub> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> <mn>2</mn> </msup> <msup> <msub> <mi>&rho;</mi> <mn>1</mn> </msub> <mn>2</mn> </msup> <mo>+</mo> <msup> <mrow> <mo>(</mo> <msub> <mi>k</mi> <mn>1</mn> </msub> <msub> <mi>&rho;</mi> <mn>1</mn> </msub> <mi>sin</mi> <msub> <mi>&beta;</mi> <mn>1</mn> </msub> <mi>tan</mi> <mi>&alpha;</mi> <mo>+</mo> <mi>h</mi> <mo>)</mo> </mrow> <mn>2</mn> </msup> </msqrt> </mrow>(7) <mrow> <mo>+</mo> <msqrt> <msup> <mrow> <mo>(</mo> <msub> <mi>k</mi> <mn>2</mn> </msub> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> <mn>2</mn> </msup> <msup> <msub> <mi>&rho;</mi> <mn>2</mn> </msub> <mn>2</mn> </msup> <mo>+</mo> <msup> <mrow> <mo>(</mo> <msub> <mi>k</mi> <mn>1</mn> </msub> <msub> <mi>&rho;</mi> <mn>1</mn> </msub> <mi>sin</mi> <msub> <mi>&beta;</mi> <mn>1</mn> </msub> <mi>tan</mi> <mi>&alpha;</mi> <mo>+</mo> <mi>h</mi> <mo>)</mo> </mrow> <mn>2</mn> </msup> </msqrt> <mo>}</mo> </mrow>其中,ρ1,ρ2,k1,k2,β1,β2为声线的位置信息,(ρ1,β1,0)为柱坐标系下声线射入沉积层的点的坐标,ρ2为反射后从沉积层反射出来的点与水听器水平面投影之间的距离,β2为反射后从沉积层反射出来的点与水听器水平面投影之间的连线与x轴的夹角,k1,k2为大于1的比例系数;cb,h,n,α,θ为海底参数,cb为沉积层中的声速,h为声源处的沉积层厚度,n为声线从水层射入沉积层的折射率,α为海底倾斜角度,θ为水平偏转角度,c0为水中的声速;H,Z,r为声源和水听器的位置信息,(0,0,H)为声源的坐标,(r,θ,Z)为水听器的坐标。
地址 150001 黑龙江省哈尔滨市南岗区南通大街145号哈尔滨工程大学科技处知识产权办公室