主权项 |
等梯度声线跟踪的声速剖面自适应分层方法,其特征在于,包括以下步骤:(1)设置曲线拟合点数n和阈值f;假设声波经历的水柱分为n个等梯度层,假设声波在每层的传播速度是等梯度的;在水深为z<sub>i</sub>的i层,分别用C<sub>i</sub>和θ<sub>i</sub>表示声波的传播速度和入射角度;第i层内的声速梯度g<sub>i</sub>表示为:g<sub>i</sub>=(C<sub>i+1</sub>‑C<sub>i</sub>)/Δz<sub>i</sub>式中,Δz<sub>i</sub>为第i层的水层厚度,C<sub>i+1</sub>为第i+1层声速;声波的传播满足Snell法则,在声速常梯度变化的情况下,波束在第i层内的实际传播轨迹为一连续的、对应一定曲率半径R<sub>i</sub>的弧段,R<sub>i</sub>的表达为:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>R</mi><mi>i</mi></msub><mo>=</mo><mfrac><mrow><mo>-</mo><mn>1</mn></mrow><mrow><msub><mi>pg</mi><mi>i</mi></msub></mrow></mfrac></mrow>]]></math><img file="FDA0000851845110000011.GIF" wi="205" he="147" /></maths>式中p为Snell系数;层i内声线的水平位移为:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>y</mi><mi>i</mi></msub><mo>=</mo><mfrac><mrow><msup><mrow><mo>[</mo><mn>1</mn><mo>-</mo><msup><mrow><mo>(</mo><msub><mi>pC</mi><mi>i</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>]</mo></mrow><mrow><mn>1</mn><mo>/</mo><mn>2</mn></mrow></msup><mo>-</mo><msup><mrow><mo>[</mo><mn>1</mn><mo>-</mo><msup><mrow><mo>(</mo><mi>p</mi><mo>(</mo><mrow><msub><mi>C</mi><mi>i</mi></msub><mo>+</mo><msub><mi>g</mi><mi>i</mi></msub><msub><mi>Δz</mi><mi>i</mi></msub></mrow><mo>)</mo><mo>)</mo></mrow><mn>2</mn></msup><mo>]</mo></mrow><mrow><mn>1</mn><mo>/</mo><mn>2</mn></mrow></msup></mrow><mrow><msub><mi>pg</mi><mi>i</mi></msub></mrow></mfrac></mrow>]]></math><img file="FDA0000851845110000012.GIF" wi="901" he="166" /></maths>式中,θ<sub>i</sub>为第i层声线的入射角,θ<sub>i+1</sub>为第i+1层声线的入射角;波束在该层经历的弧段长度为:S<sub>i</sub>=R<sub>i</sub>(θ<sub>i</sub>‑θ<sub>i+1</sub>)则经历该段的时间为:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>t</mi><mi>i</mi></msub><mo>=</mo><mfrac><mrow><mi>arcsin</mi><mo>[</mo><mi>p</mi><mrow><mo>(</mo><msub><mi>C</mi><mi>i</mi></msub><mo>+</mo><msub><mi>g</mi><mi>i</mi></msub><msub><mi>Δz</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>]</mo><mo>-</mo><mi>arcsin</mi><mrow><mo>(</mo><msub><mi>pC</mi><mi>i</mi></msub><mo>)</mo></mrow></mrow><mrow><msubsup><mi>pg</mi><mi>i</mi><mn>2</mn></msubsup><msub><mi>Δz</mi><mi>i</mi></msub></mrow></mfrac><mi>ln</mi><mo>[</mo><mn>1</mn><mo>+</mo><mfrac><mrow><msub><mi>g</mi><mi>i</mi></msub><msub><mi>Δz</mi><mi>i</mi></msub></mrow><msub><mi>C</mi><mi>i</mi></msub></mfrac><mo>]</mo></mrow>]]></math><img file="FDA0000851845110000013.GIF" wi="1133" he="150" /></maths>(2)已知第i层的水平位移y<sub>i</sub>和经历该段位移的时间t<sub>i</sub>,得出第i层的平均声速C<sub>i</sub>(i=1,2,…,n),从顶层开始选取连续的n个声速值C<sub>i</sub>(i=1,2,…,n)进行曲线拟合,得到拟合曲线C=f(x),曲线是一个关于声速C<sub>i</sub>(i=1,2,…,n)与水平位移x有关的变化曲线;(3)拟合的曲线C求其最大曲率,曲率计算公式为:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mi>R</mi><mo>=</mo><mfrac><mrow><mo>|</mo><msup><mi>C</mi><mrow><mo>′</mo><mo>′</mo></mrow></msup><mo>|</mo></mrow><msup><mrow><mo>(</mo><mn>1</mn><mo>+</mo><msup><mi>C</mi><mrow><mo>′</mo><mn>2</mn></mrow></msup><mo>)</mo></mrow><mrow><mn>3</mn><mo>/</mo><mn>2</mn></mrow></msup></mfrac></mrow>]]></math><img file="FDA0000851845110000014.GIF" wi="316" he="164" /></maths>其中,C'和C”代表声速曲线的一阶和二阶导数。(4)将曲率半径R与阈值f进行对比,若小于阈值f则将n层声速剖面合并为一层;若大于阈值f则释放第一层,释放n个点中的第一个点,下移一个点继续选取n个点声速值进行拟合,直至整个声速剖面数据遍历完成。 |