发明名称 |
基于射线声学的三维倾斜海底参数快速测量方法 |
摘要 |
本发明提供的是一种基于射线声学的三维倾斜海底参数快速测量方法。声源发出的声信号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>ρ</mi> <mn>1</mn> </msub> <mn>2</mn> </msup> <mrow> <msup> <msub> <mi>ρ</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>ρ</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>ρ</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>ρ</mi> <mn>1</mn> </msub> <mi>sin</mi> <msub> <mi>β</mi> <mn>1</mn> </msub> <mi>tan</mi> <mi>α</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>ρ</mi> <mn>2</mn> </msub> <mn>2</mn> </msup> <mrow> <msup> <msub> <mi>ρ</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>ρ</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>ρ</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>ρ</mi> <mn>1</mn> </msub> <mi>sin</mi> <msub> <mi>β</mi> <mn>1</mn> </msub> <mi>tan</mi> <mi>α</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>ρ</mi> <mn>1</mn> </msub> <mi>sin</mi> <msub> <mi>β</mi> <mn>1</mn> </msub> <mi>tan</mi> <mi>α</mi> <mo>-</mo> <mrow> <mo>(</mo> <msub> <mi>k</mi> <mn>1</mn> </msub> <msub> <mi>ρ</mi> <mn>1</mn> </msub> <mi>sin</mi> <msub> <mi>β</mi> <mn>1</mn> </msub> <mi>tan</mi> <mi>α</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>ρ</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>ρ</mi> <mn>1</mn> </msub> <mi>sin</mi> <msub> <mi>β</mi> <mn>1</mn> </msub> <mi>tan</mi> <mi>α</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>ρ</mi> <mn>2</mn> </msub> <mi>sin</mi> <msub> <mi>β</mi> <mn>2</mn> </msub> <mi>tan</mi> <mi>α</mi> <mo>-</mo> <mrow> <mo>(</mo> <msub> <mi>k</mi> <mn>1</mn> </msub> <msub> <mi>ρ</mi> <mn>1</mn> </msub> <mi>sin</mi> <msub> <mi>β</mi> <mn>1</mn> </msub> <mi>tan</mi> <mi>α</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>ρ</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>ρ</mi> <mn>1</mn> </msub> <mi>sin</mi> <msub> <mi>β</mi> <mn>1</mn> </msub> <mi>tan</mi> <mi>α</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>α</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>ρ</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>ρ</mi> <mn>1</mn> </msub> <mn>2</mn> </msup> </msqrt> <mrow> <mo>(</mo> <mi>r</mi> <mi>cos</mi> <mi>θ</mi> <mo>+</mo> <msub> <mi>ρ</mi> <mn>2</mn> </msub> <mi>cos</mi> <msub> <mi>β</mi> <mn>2</mn> </msub> <mo>-</mo> <msub> <mi>ρ</mi> <mn>1</mn> </msub> <mi>cos</mi> <msub> <mi>β</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>ρ</mi> <mn>1</mn> </msub> <msub> <mi>ρ</mi> <mn>2</mn> </msub> <mi>sin</mi> <mrow> <mo>(</mo> <msub> <mi>β</mi> <mn>2</mn> </msub> <mo>-</mo> <msub> <mi>β</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>ρ</mi> <mn>1</mn> <mn>2</mn> </msubsup> <mo>+</mo> <msup> <mi>H</mi> <mn>2</mn> </msup> </msqrt> <mo>+</mo> <msqrt> <msubsup> <mi>ρ</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>ρ</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>ρ</mi> <mn>1</mn> </msub> <mi>sin</mi> <msub> <mi>β</mi> <mn>1</mn> </msub> <mi>tan</mi> <mi>α</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>ρ</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>ρ</mi> <mn>1</mn> </msub> <mi>sin</mi> <msub> <mi>β</mi> <mn>1</mn> </msub> <mi>tan</mi> <mi>α</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号哈尔滨工程大学科技处知识产权办公室 |