基于DTOA/DOA和牛顿迭代法的震源定位方法研究

开篇：润墨网以专业的文秘视角，为您筛选了一篇基于DTOA/DOA和牛顿迭代法的震源定位方法研究范文，如需获取更多写作素材，在线客服老师一对一协助。欢迎您的阅读与分享！

摘要：结合DTOA和DOA联合定位估计性好与牛顿迭代法收敛速度快的优点，提出了一种基于dtoa和doa初值选取和牛顿法精确迭代的震源定位方法。经过野外实验和仿真，发现该方法能有效地提高震源定位的准确性。通过对影响震源定位精度因素分析，表明减小传感器基站站盐和初至波到时差测量误差可有效提高定位精度，相同同条件下，星型布站的定位精度最高。

关键词：震源定位；时差测量；俯仰角测量；牛顿迭代法；几何精度因子

中图分类号：P315文献标识码：A文章编号：1000-0666（2013）03-0324-06

0引言

随着现代科学技术的发展，定位技术在航空、航天、交通、震源探测、海洋勘探等领域得到了广泛地应用。定位精度的要求越来越高，在硬件系统确定条件下，可以通过改进算法提高系统的定位精度（吕晶晶，姚金杰，2011；姚金杰，韩焱，2010）。本文提出了一种以牛顿迭代算法为基础，利用无源时差定位（DTOA）和测量波达方向（DOA）联合定位得到微震源的估计位置的算法，本算法充分利用DTOA和DOA联合定位算法的估计特性，能够有效地解决迭代法的初始值问题，保证算法收敛并且提高迭代算法的收敛速度。通过设置建立的冗余函数表达式和结束迭代的精度范围，作为判断迭代是否结束的条件，这样既能满足一定的精度，又能提高计算速度，本文通过对计算机仿真、数据分析和计算时间来证明该算法的有效性。

1目标定位原理

为了充分利用各种观测信息，以获取更好的定位效果，我们将观测数据中的角度信息与时间信息结合起来。本文运用的DTOA和DOA联合定位方法，主要利用主站增加对目标俯仰角的测量，形成一个包含3个子系统的冗余定位系统，各子系统分别对目标进行定位。最后对定位结果进行融合，将融合后的结果作为信源目标的初始估计值，再通过牛顿迭代算法进行精确定位。

11DTOA定位原理

无源时差定位利用主站和辅站之间接收到辐射源信号的到达时间差来定位。获得时差的信息之后，通过乘以波的传播速度就可以得到信源到主站和辅站的距离差，在三维空间定位中，得到一组双曲面方程组。在三维时差定位中至少需要4个测量站才能实现对目标的定位（陈玲，李少洪，2003），如图1所示。图1三维时差定位示意图

Fig1Schematic diagram of 3D TDOA location设震源位置坐标为P（x，y，z），基站坐标分别为Pi（xi，yi，zi）（i=0，1，2，3…n）。目标到基站的距离为ri（i=0，1，2，3…n）。目标到各分站与主站（基站P0为主站）的距离差为ri0。假设有4个基站，则根据时差定位原理有以下关系r20=（x-x0）2+（y-y0）2+（z-z0）2

r2i=（x-xi）2+（y-yi）2+（z-zi）2

r0i=r0-ri=ct0i（1）式中，t0i表示信源信号到达各分站与主站之间的时间差，c表示信号传播的速度。

1地震研究136卷第3期杨俊峰等：基于DTOA/DOA和牛顿迭代法的震源定位方法研究12DOA定位原理

测量波达方向技术在无源定位中占有极其重要的作用。随着数字信号处理技术的发展，数字相位干涉仪的测向方法由于其具有速度快、成本低、精度高的特点在工程上获得了广泛的应用。在原理上，相位干涉仪可以实现对单个脉冲测向（袁孝康，1999）。图2显示相位干涉仪（单基线）测向原理图。图2单基线相位干涉仪侧向原理图

Fig2Lateral principle graph of singlebaseline

phase interferometer假设发射信号的波长为λ，与天线视线轴的夹角为β，则到达相位干涉仪两个测向天线1、2间的相位差为=2πLsinβ1λ（2）式中，L为两个测向天线之间的基线长度，λ为入射信号的波长，β为入射信号与天线视线轴的夹角。由式（2）可以得到信号的波达角为β^=arcsin（λ12πL）（3）式中，φ是两个阵元之间输出的相位差，入射信号波长的估计值λ由信号的频率估计值f^计算得到，基线长度L则可直接测量，在上述参数确定之后，就能根据式（3）得到辐射源信号的入射角度β^。

2算法121DTOA和DOA联合定位算法求初始值利用相位干涉仪测向方法，各辅站测得的目标的俯仰角分别为sinα1=（z-z1）r1，（4）sinα2=（z-z2）r2，（5）sinα3=（z-z3）r3（6）式中，αi（i=1，2，3）表示各辅站测得的俯仰角；z表示目标辐射源的高度。

将式（1）和式（4）组成方程组（其中，i=1，2），称为子系统Ⅰ，同理，将式（1）和式（5）组成方程组（其中，i=2，3），称为子系统Ⅱ，式（1）和式（6）组成方程组（其中，i=1，3），称为子系统Ⅲ，这三个子系统均可完成对目标的空间三维定位，我们选择这三个子系统组成一个冗余定位系统，则该系统的定位方程组为：sinα1（z-z1）r1

r20=（x-x0）2+（y-y0）2+（z-z0）2

r2i=（x-xi）2+（y-yi）2+（z-zi）2

r0i=r0-ri=ct0i，（7）其中，i=1，2，可改写为z=（r0-r01）sinα1+z1

（x0-x1）x+（y0-y1）y+（z0-z1）z=k1+r0r01

（x0-x2）x+（y0-y2）y+（z0-z2）z=k2+r0r02（8）其中，ki=112r20i+（x20+y20+z20）-（x2i+y2i+z2i）

（i=1，2），可改矩阵形式为AX=F（9）其中，A=01011

x0-x11y0-y11z0-z1

x0-x21y0-y21z0-z2，A是可逆矩阵，可以求得X的最小二乘解为X^=（ATA）-1ATF（10）在消除定位模糊以后，3个子系统一共会有3个定位解，如何利用这3个定位解来获取更加精确的目标真实位置是测向测时差联合定位所要解决的关键问题。可以采用简化加权最小二乘（SWLS）的方法对目标位置进行估计（吕明，郭士民，2007；Carter，1987）。

假设3个子系统所得到的目标位置分别为X^i（i=1，2，3），定位协方差矩阵分别为Pi（i=1，2，3），各协方差矩阵之间互不相关。SWLS融合后的目标位置为（Petre，Li，2006；毛永毅，白菊蓉，2006）X^=x^y^z^T=∑31i=1p-1i-1∑31i=1p-1iX^i（11）22牛顿迭代定位算法

牛顿法是一种使用导数的算法，每一步迭代方向都是沿着当前点函数值下降的方向。因此，需设定合适的初值才能保证算法是否收敛及其收敛速度。根据牛顿迭代算法的条件，将所建立目标定位模型转换为（柳辉，2007）fi（x，y，z）=（x-x0）2+（y-y0）2+（z-z0）2-（x-xi）2+（y-yi）2+（z-zi）2-ct0i，（12）其雅克比矩阵为f′（x，y，z）=f1/x1f1/y1f1/z

（13）当雅克比矩阵为非奇异阵，则震源坐标用牛顿迭代法表示为x（k+1）

y（k+1）

z（k+1）=x（k）

y（k）

z（k）-fi（x，y，z）-1f1（x（k），y（k），z（k））

f2（x（k），y（k），z（k））

f3（x（k），y（k），z（k））

（14）2计算机仿真及分析

在野外条件下，通过实验完成震源定位。首先测算震动波在试验场地介质中的传播速度，再根据加速度传感器接收到的震动波数据，结合震源定位的原理及算法，得出震源的位置（邹翔宇，徐翊峰，2010）。

震动波在该实验场地介质中的传播速度测算具体方法是：在距离待爆炸位置3 m、5 m、7 m、10 m、15 m的直线上分别依次布设5个加速度传感器A、B、C、D、E，埋深均为05 m，在同一位置爆炸3次，通过数据采集仪采集震动波数据，并将数据储存在计算机里，利用MATLAB对信号进行分析处理，计算出两个加速度传感器之间接收到震动信号的时间差，根据两个传感器距离震源的距离差计算震动波的传播速度v=Δd/Δt。通过计算测得的震动波传播速度如表1所示。表1震动波传播速度测算

Tab1Measurement of vibratioon wave propagation velocity

测试次数1速度/m·s-1113857213727313799平均值13794

通过加速度传感器接收到的震动波分析处理可大致得到震动波的传播规律，在距震源约7 m处的震动信号衰减为弹性波。以此为依据在待测震源周围布设加速度传感器，传感器基站的位置分别为（6，3，-05），（-4，6，-04），（-6，-5，-042），（7，-4，-028），（0，5，-038）。数据采集仪的采样率为10 kHz，采样时间为6553 8 s，采用手动触发方式。在同一地点爆炸3次，通过测算的震动波传播速度，传感器基站间接收到震动波的时差及震动定位算法以解算出震源位置。加速度传感器接收到的局部放大信号如图3所示。图3局部震动信号

Fig3Local vibration signal

利用时窗能量比法判断初至波到时，对传感器1接收到的震动信号进行预处理及时窗能量比，如图4所示。

以同样的方法判断其它各传感器接收到的初至波到时，结合各传感器基站位置及定位算法解算出的震源位置如表2所示。对运用本文算法所得定位结果进行误差分析如表3所示。

在定位过程中，由于传感器基站的布设及埋深较浅，使得Z方向的误差较大，加之震动波在传播过程中引入的误差，传感器基站站址测量误差，初至波到时估计误差等因数影响导致震源位置与实际位置之间有偏差。

从表3可知，本文算法定位结果的均方根误差值远远小于DTOA和DOA联合定位、LSNewton法及ChanNewton法定位结果的均方根误差值，定位结果比较理想。图4判断初至波时刻图

Fig4Determine the time of the first arrival表3定位结果误差分析

Tab3Error analysis of positioning results

实验1X轴误差/m1Y轴误差/m1Z轴误差/m1RMSE/m110000 610002 110042 410042 5210000 810002 810064 410064 5310000 710002 710062 010062 13震源定位精度分析

通过对震源定位进行精度分析，可以了解影响定位精度的因素以及这些因素对定位精度产生的影响，因此对震源定位进行精度分析是很必要的。下文对在不同的布站形式、基站数目、初至波到达时间差估计误差、及在不同站址测量误差情况下的定位精度变化情况进行仿真分析，仿真结果可供工程参考。图5 不同误差情况下星型布站GDOP分布图

（a）a=003 m，b=05 ms；（b）a=001 m，b=50 us

Fig5GDOP distribution in starmodel station

position under different error conditions设震源范围是x方向0～30 m，y方向是-5～+5 m，z方向-1 m，对4个传感器基站为星型布局时进行仿真分析。图5a中存在a=003 m的传感器基站站址测量误差，b=05 ms的初至波到达时间差测量误差。图5b中存在a=001 m的传感器基站站址测量误差，b=50 us的初至波到达时间差测量误差。由图5b可以看出，减小传感器基站站址和初至波达到时差测量误差可有效地提高定位精度。

以4个传感器基站为例，在相同的传感器基站站址测量误差a=001 m和初至波达到时差测量误差b=50 us的情况下，对不同布设的定位精度进行仿真，如图6所示。从图中4种布站下的GDOP图中可知，星型布站的定位精度最高；平行四边形布站的定位次之，其在X轴方向的定位精度要高于Y轴；菱形布站和倒三角布站的定位精度较差。图6不同布站形式GDOP图

（a）星型布站；（b）平行四边形布站；（c）菱形布站；（d）倒三角布站

Fig6GDOP in different forms of station arrangement

（a）starmodel station position；（b）parallelogrammodel station position；（c）dimodmodel

station position；（d）upsidedown trianglemodel station position4结语

本文结合DTOA和DOA联合定位法估计性好和牛顿迭代法收敛速度快的优点，提出了DTOA、DOA和牛顿迭代法相结合的地下震源定位算法。避免了DTOA和DOA联合定位法受测量误差影响较大和牛顿法对初始值选择敏感的问题。经过在野外条件下实验及仿真，实验结果表明本文算法有效地提高了震源定位的精度。在此基础上对震源系统中存在不同站址测量误差、不同初至波到达时差估计误差、不同布站形式情况下的定位精度进行了仿真分析，结果表明减小传感器基站站址测量误差和初至波达到时差测量误差可有效地提高定位精度，星型布站的定位精度要优于平行四边形布站、菱形布站和倒三角布站的定位精度。参考文献：

陈玲，李少洪2003无源侧向测时差定位算法研究[J].电子信息学报，（6）：771-776

柳辉2007解非线性方程组的牛顿迭代法及其应用[J].重庆工学院学报，21（8）：95-98

吕晶晶，姚金杰2011基于最小二乘法和牛顿迭代法的空中目标定位[J].微电子学于计算机，28（9）：110-115

吕明，郭士民2007基于数据融合的时差定位处理算法的应用[J].仪器仪学报，28（4）：100-102

毛永毅，白菊蓉2006空间四站时差定位中的模糊及无线研究[J].电讯技术，（6）：53-57

姚金杰，韩焱2010基于粒子群和牛顿迭代法的目标定位方法研究[J].计算机应用研究，27（5）：1 700-1 713

袁孝康1999相位干涉仪侧向定位研究[J].上海航天，（3）：1-7

邹翔宇，徐翊峰2010基于无线传感器网络的矿震震源定位[J].煤矿安全，（10）：53-57