等梯度声线跟踪的声速剖面自适应分层方法
【技术领域】
[0001] 本发明设及一种水下声学定位方法，尤其是设及等梯度声线跟踪的声速剖面自适 应分层方法。
【背景技术】
[0002] 在水下声学定位中，由于声波在海水中的传播轨迹是弯曲的，要获得高精度的定 位结果，大多采用分层等梯度声线跟踪方法进行定位。在深海定位中，声线传播距离长，声 速剖面层数多，采用该方法运算量大，会明显降低定位的计算效率。
【发明内容】
[0003] 针对上述技术问题，本发明提供一种计算简单、计算效率高的水下定位方法，具体 的技术方案为：
[0004] 等梯度声线跟踪的声速剖面自适应分层方法，包括W下步骤：
[0005] (1)设置曲线拟合点数η和阔值f;
[0006] 假设声波经历的水柱分为η个等梯度层，假设声波在每层的传播速度是等梯度的； 在水深为Zi的i层，分别用Ci和0i表示声波的传播速度和入射角度;第i层内的声速梯度gi表 示为：
[0007] gi= (Ci+i-Ci)/Δzi
[000引式中，AZi为第i层的水层厚度，Ci+i为第i+1层声速；
[0009] 声波的传播满足Snell法则，在声速常梯度变化的情况下，波束在第i层内的实际 传播轨迹为一连续的、对应一定曲率半径Ri的弧段，虹的表达为：
[0010]
[0011] 式中P为Snell系数；
[0012] 层i内声线的水平位移为：
[0013]
[0014] 式中，θι为第i层声线的入射角，θι+ι为第i+1层声线的入射角；
[0015] 波束在该层经历的弧段长度为：
[0016] Si=Ri(目广目i+i)
[0017] 则经历该段的时间为：
[001 引
[0019] (2)已知第i层的水平位移yi和经历该段位移的时间ti，得出第i层的平均声速Ci(i =1，2^-，11)，从顶层开始选取连续的11个声速值(：1。= 1，2，-，，11)进行曲线拟合，得到拟合 曲线C=f(x)，曲线是一个关于声速。。=1，2，一，11)与水平位移^有关的变化曲线；
[0020] (3)拟合的曲线C求其最大曲率，曲率计算公式为：
[0021]
[0022] 其中，C'和C'代表声速曲线的一阶和二阶导数。
[0023] (4)将曲率半径R与阔值f进行对比，若小于阔值f则将η层声速剖面合并为一层;若 大于阔值f则释放第一层，释放η个点中的第一个点，下移一个点继续选取η个点声速值进行 拟合，直至整个声速剖面数据遍历完成。
[0024]本发明提供的等梯度声线跟踪的声速剖面自适应分层方法，将声速剖面数据中变 化很小的声速层进行过滤，适当抽稀，保留声速梯度变化较大的数据层，从而减少迭代中的 运算量。声速剖面经过自适应分层后，用新的声速剖面进行分层等梯度声线跟踪，从而完成 水下定位解算。
【附图说明】
[0025]图1为本发明的流程图；
[0026]图2为实施例中所用到的中国南海声速剖面数据；
[0027]图3为实施例中阔值为1(Τ?时声速剖面自适应分层的结果；
[0028]图4为实施例中阔值为时声速剖面自适应分层结果。
【具体实施方式】
[0029]下面结合附图和具体实施例对本发明进行详细说明。
[0030] 对中国南海500米水深声速剖面进行自适应分层，并用自适应分层后的结果进行 定位：
[0031] 1、定位数据
[003^ 实验选取拟合点数为7,进行S次曲线拟合，阔值设置为ιοΛιοΛιοΛ... 共20个阔值。
[0033] 模拟四个浮标布设成正方形形状，边长为1600m，给浮标的Ξ维坐标分别加上2米 的余弦波动。海底模拟4个应答器A、B、C、D，坐标分别为(0,0，-500)、（0,40，-500)、（40,0，-500)、（40，40，-500)，单位均为m，总共模拟100个历元。水域的声速剖面图如图2,为中国南 海声速剖面数据，采样间隔为Im。误差模拟如下:换能器定位中误差x、y、z方向各10cm，应答 器时延偏差为8cm的系统误差，测量时间造成10cm的测距中误差。
[0034] 2、声速剖面自适应分层结果
[003引共取10-1，10-2，10-3，...，10-^，20个阔值，其中最大阔值为10-咐可将原始500层声 速剖面数据分为145层，如图3所示;最小阔值1(Γ2咐可将500层声速剖面数据分为450层，如 图4所示。
[0036] 3、利用自适应分层后的数据进行定位
[0037] 实验中的电脑配置为:CPU型号为双核化ntiumE2200，CPU主频为2.20細Ζ，内存大 小为3G。
[0038] (1)原始声速剖面定位结果
[0042]
[0043]从上表结果中可W看出，当拟合点数为7时，对于不同的阔值，对声速剖面有不同 程度的抽稀，层数从145层到450层，其中，阔值从10^1到ΚΓ"分层结果一致，均为145层。定位 的时间从11s到34s，定位的效率较原始结果有了不同幅度的提升。定位的精度与原始定位 误差接近，不同阔值定位误差差异很小。对于不同的阔值，其定位RMS值均在0.36m左右，相 差不大，与未经过自适应分层的定位结果相当。
[0044]本发明提出的声速剖面自适应分层方法能够对声速剖面进行合理的抽稀，减少声 速剖面的层数。通过用自适应分层后的数据进行定位可W看出，定位结果均能在保证定位 精度的前提下提高定位效率。
【主权项】
1.等梯度声线跟踪的声速剖面自适应分层方法，其特征在于，包括以下步骤： (1) 设置曲线拟合点数η和阈值f; 假设声波经历的水柱分为η个等梯度层，假设声波在每层的传播速度是等梯度的；在水 深为21的1层，分别用(：1和01表示声波的传播速度和入射角度;第i层内的声速梯度8 1表示 为： gi= (Ci+i-Ci)/ Δ zi 式中，A zi为第i层的水层厚度，Ci+1为第i+1层声速； 声波的传播满足Snell法则，在声速常梯度变化的情况下，波束在第i层内的实际传播 轨迹为一连续的、对应一定曲率半径Ri的弧段，Ri的表达为：式中p为Snell系数； 层i内声线的水平位移为：式中，为第i层声线的入射角，θ1+1为第i+Ι层声线的入射角； 波束在该层经历的弧段长度为： Si = Ri(0i-0i+1) 则经历该段的时间为：(2) 已知第i层的水平位移yi和经历该段位移的时间^，得出第i层的平均声速α(? = 1， 2,…，η)，从顶层开始选取连续的η个声速值匕^二^广^"进行曲线拟合^导到拟合曲线 C = f(x)，曲线是一个关于声速Ci(i = l，2,···，!!)与水平位移X有关的变化曲线； (3) 拟合的曲线C求其最大曲率，曲率计算公式为：其中，C '和C"代表声速曲线的一阶和二阶导数。 (4) 将曲率半径R与阈值f进行对比，若小于阈值f则将η层声速剖面合并为一层;若大于 阈值f则释放第一层，释放η个点中的第一个点，下移一个点继续选取η个点声速值进行拟 合，直至整个声速剖面数据遍历完成。
【专利摘要】本发明涉及一种水下声学定位方法，尤其是涉及等梯度声线跟踪的声速剖面自适应分层方法，包括设置曲线拟合点数n，阈值f；从顶层开始选取连续的n个声速值进行曲线拟合；对于拟合的曲线求取最大曲率半径；将曲率半径与阈值f进行对比，若小于阈值f则将n层声速剖面合并为一层；若大于阈值f则释放第一层，从下一层继续拟合曲线，直至整个声速剖面数据遍历完成。本发明提供的等梯度声线跟踪的声速剖面自适应分层方法，将声速剖面数据中变化很小的声速层进行过滤，适当抽稀，保留声速梯度变化较大的数据层，从而减少迭代中的运算量。声速剖面经过自适应分层后，用新的声速剖面进行分层等梯度声线跟踪，从而完成水下定位解算。
【IPC分类】G01S1/72
【公开号】CN105445694
【申请号】CN201510809031
【发明人】王振杰, 李圣雪, 王毅, 聂志喜, 吴会胜
【申请人】中国石油大学(华东)
【公开日】2016年3月30日
【申请日】2015年11月19日