水下导航与重力测量一体化系统Integrated system of underwater navigation and gravity measurement
技术领域Technical field
本发明涉及水下导航与重力测量一体化系统,属于组合导航和重力测量技术领域。The invention relates to an integrated system of underwater navigation and gravity measurement, and belongs to the technical field of integrated navigation and gravity measurement.
背景技术Background technique
海洋蕴藏着丰富的油气矿产资源,探索和开发海洋对国民经济发展意义重大。油气矿产资源的分布与海洋重力场密切相关,水下重力测量有助于发现它们。目前水下重力测量的困难是现有的重力测量系统难以保证水下长时间的准确定位和低成本的持续作业问题。为了解决这些问题,本发明提供了水下导航与重力测量一体化系统,它能够满足长航时、高精度和低成本的水下导航和重力测量要求。The ocean contains abundant oil and gas mineral resources, and the exploration and development of the ocean is of great significance to the development of the national economy. The distribution of oil and gas mineral resources is closely related to the ocean gravity field, and underwater gravity measurement helps to discover them. The current difficulty of underwater gravity measurement is that the existing gravity measurement system is difficult to ensure long-term accurate positioning underwater and low-cost continuous operation. In order to solve these problems, the present invention provides an integrated system of underwater navigation and gravity measurement, which can meet the requirements of long-endurance, high-precision and low-cost underwater navigation and gravity measurement.
发明内容Summary of the invention
本发明的目的是为了解决上述问题,提供一种水下导航与重力测量一体化系统,能够满足长航时、高精度和低成本的水下导航和重力测量要求。The purpose of the present invention is to solve the above-mentioned problems and provide an integrated system of underwater navigation and gravity measurement, which can meet the requirements of long-endurance, high-precision and low-cost underwater navigation and gravity measurement.
为了达到上述目的,本发明采用的技术方案是:一种水下导航与重力测量一体化系统,主要由3个激光陀螺仪、3个摆式挠性加速度计、3个重力传感器、转位与测角装置、温度传感器、温度控制装置、多普勒计程仪、深度计、2个计算机、数据采集处理平台、LCD显示和电路板组成。In order to achieve the above objective, the technical solution adopted by the present invention is: an integrated underwater navigation and gravity measurement system, which is mainly composed of 3 laser gyroscopes, 3 pendulum flexible accelerometers, 3 gravity sensors, indexing and Angle measuring device, temperature sensor, temperature control device, Doppler log, depth gauge, 2 computers, data acquisition and processing platform, LCD display and circuit board.
3个激光陀螺仪和3个摆式挠性加速度计放置在转位机构上,构成激光陀螺单轴旋转惯性测量单元。Three laser gyroscopes and three pendulum flexible accelerometers are placed on the indexing mechanism to form a single-axis rotation inertial measurement unit of the laser gyroscope.
3个重力传感器放在温度控制箱中,与激光陀螺单轴旋转惯性测量单元刚性地固定在一起,且与激光陀螺单轴旋转惯性测量单元的加速度计具有相同的姿态角。The three gravity sensors are placed in the temperature control box, rigidly fixed with the laser gyro single-axis rotating inertial measurement unit, and have the same attitude angle as the accelerometer of the laser gyro single-axis rotating inertial measurement unit.
数据采集处理平台实时处理记录激光陀螺仪、加速度计、重力传感器、转位与测角装置、温度传感器、多普勒计程仪、深度计信号。The data acquisition and processing platform processes and records the signals of laser gyroscopes, accelerometers, gravity sensors, indexing and angle measurement devices, temperature sensors, Doppler logs, and depth gauges in real time.
1个计算机接收激光陀螺仪、高精度加速度计、多普勒计程仪、深度计、转位角度信息和温度信息,以多普勒计程仪的速度信息和深度计的水深信息作为激光陀螺单轴旋转惯性测量单元的外观测量,应用扩展卡尔曼滤波,计算出水下导航与重力测量一体化系统的位置、速度和姿态。A computer receives laser gyroscope, high-precision accelerometer, Doppler log, depth gauge, indexing angle information and temperature information, and uses the speed information of the Doppler log and the water depth information of the depth gauge as the laser gyroscope The appearance measurement of the single-axis rotating inertial measurement unit applies the extended Kalman filter to calculate the position, speed and attitude of the integrated underwater navigation and gravity measurement system.
1个计算机接收重力传感器信息、温度信息,水下导航与重力测量一体化系统的位置、速度和姿态信息,计算出当地地理坐标系下的重力传感器的比力,进行重力测量改正和低通滤波处理,得到当地重力值。A computer receives the gravity sensor information, temperature information, the position, speed and attitude information of the integrated underwater navigation and gravity measurement system, calculates the specific force of the gravity sensor in the local geographic coordinate system, and performs gravity measurement correction and low-pass filtering Processing to get the local gravity value.
作为本发明的一种改进,所述的数据采集处理平台与计算机采用FPGA+DSP+ARM架构。As an improvement of the present invention, the data acquisition and processing platform and computer adopt FPGA+DSP+ARM architecture.
有益效果:Beneficial effects:
本发明提供的水下导航与重力测量一体化系统,能够满足水下运载体长航时、高精度、低成本的导航和重力测量的需求。The integrated system of underwater navigation and gravity measurement provided by the present invention can meet the requirements of long-endurance, high-precision, low-cost navigation and gravity measurement of underwater vehicles.
附图说明Description of the drawings
图1为本发明的惯性测量单元与重力传感器一体化装置结构示意图。FIG. 1 is a schematic diagram of the structure of the integrated device of the inertial measurement unit and the gravity sensor of the present invention.
图2为本发明的数据采集处理平台与计算机示意图。Figure 2 is a schematic diagram of the data acquisition and processing platform and computer of the present invention.
[根据细则91更正 17.07.2020]
图3为本发明的水下组合导航系统导航算法流程图。[Corrected according to Rule 91 17.07.2020]
Fig. 3 is a flowchart of the navigation algorithm of the underwater integrated navigation system of the present invention.
具体实施方式detailed description
下面结合具体实施例对本发明作进一步地说明。The present invention will be further described below in conjunction with specific embodiments.
如图1所示,本实施例公开了水下导航与重力测量一体化系统,主要由惯性测量单元与重力传感器一体化装置、多普勒计程仪、深度计、2个计算机、数据采集平台和LCD显示组成。As shown in Figure 1, this embodiment discloses an integrated system for underwater navigation and gravity measurement, which is mainly composed of an integrated device of an inertial measurement unit and a gravity sensor, a Doppler log, a depth gauge, two computers, and a data acquisition platform. And LCD display composition.
惯性测量单元与重力传感器一体化装置包括防护壳体1、设置在防护壳体内的3个激光陀螺仪8、3个摆式挠性加速度计9、3个重力传感器2、转位与测角装置4、温度传感器7、温度控制装置5、电路板3,6。The integrated device of the inertial measurement unit and the gravity sensor includes a protective housing 1, 3 laser gyroscopes 8 arranged in the protective housing, 3 pendulum flexible accelerometers 9, 3 gravity sensors 2, indexing and angle measuring devices 4. Temperature sensor 7, temperature control device 5, circuit board 3, 6.
三个90型二频机抖激光陀螺仪、三个高精度摆式挠性加速度计、转位与测角装置和电路板等组成激光陀螺单轴旋转惯性测量单元,转位与测角装置4由力矩电机、驱动器、光电码盘测角仪、导电滑环和机械部件组成。Three 90-type two-frequency machine-shaking laser gyroscopes, three high-precision pendulum flexible accelerometers, indexing and angle measuring devices and circuit boards, etc. constitute the laser gyroscope single-axis rotating inertial measurement unit, indexing and angle measuring device 4 It is composed of a torque motor, a driver, a photoelectric code disc goniometer, a conductive slip ring and mechanical parts.
3个重力传感器放在温度控制箱中,与激光陀螺单轴旋转惯性测量单元刚性地固定在一起,且与激光陀螺单轴旋转惯性测量单元的加速度计具有相同的姿态角。温度控制箱由温度控制板、热敏电阻、加热片组成。The three gravity sensors are placed in the temperature control box, rigidly fixed with the laser gyro single-axis rotating inertial measurement unit, and have the same attitude angle as the accelerometer of the laser gyro single-axis rotating inertial measurement unit. The temperature control box is composed of a temperature control board, a thermistor, and a heating plate.
数据采集处理平台与计算机采用FPGA+DSP+ARM架构,如图2所示,FPGA和ARM负责数据采集和处理,DSP作为计算机负责导航与重力计算。The data acquisition and processing platform and computer adopt FPGA+DSP+ARM architecture, as shown in Figure 2, FPGA and ARM are responsible for data acquisition and processing, and DSP as the computer is responsible for navigation and gravity calculation.
数据采集处理平台实时处理记录激光陀螺仪、加速度计、重力传感器、转位与测角装置、温度传感器、多普勒计程仪、深度计信号。The data acquisition and processing platform processes and records the signals of laser gyroscopes, accelerometers, gravity sensors, indexing and angle measurement devices, temperature sensors, Doppler logs, and depth gauges in real time.
计算机执行组合导航计算和重力实时数据处理任务,1个计算机接收激光陀螺仪、高精度加速度计、多普勒计程仪、深度计、转位角度信息和温度信息,以多普勒计程仪的速度信息和深度计的水深信息作为激光陀螺单轴旋转惯性测量单元的外观测量,应用扩展卡 尔曼滤波,计算出水下导航与重力测量一体化系统的位置、速度和姿态,其算法流程如图3所示。The computer performs integrated navigation calculations and real-time gravity data processing tasks. A computer receives laser gyroscope, high-precision accelerometer, Doppler log, depth gauge, indexing angle information and temperature information, and uses Doppler log The speed information and the depth information of the depth gauge are used as the appearance measurement of the single-axis rotating inertial measurement unit of the laser gyro. The extended Kalman filter is applied to calculate the position, speed and attitude of the integrated underwater navigation and gravity measurement system. The algorithm flow is shown in the figure. 3 shown.
水下组合导航系统滤波状态方程为The filter state equation of the underwater integrated navigation system is
其中,X
k是系统状态向量,Φ
k+1/k是状态转移矩阵,Γ
k+1是系统的噪声转换矩阵,W
K是噪声矩阵。状态向量为:
Among them, X k is the system state vector, Φ k+1/k is the state transition matrix, Γ k+1 is the noise transition matrix of the system, and W K is the noise matrix. The state vector is:
其中,
分别为捷联惯导数学平台的3个误差角,δV
E,δV
N,δV
UP分别为东向、北向和天向速度误差,δL,δλ,δh分别为经度、纬度和高度误差,δG
x,δG
y和δG
z分别为陀螺仪X,Y,Z轴的零位漂移,δA
x,δA
y,δA
z分别为加速度计X,Y,Z轴的零位偏置。W噪声矩阵由陀螺仪和加速度计的白噪声组成:
in, These are the three error angles of the strapdown inertial navigation mathematics platform, δV E , δV N , δV UP are the east, north and sky velocity errors, δL, δλ, and δh are the longitude, latitude, and altitude errors, respectively, and δG x , ΔG y and δG z are the zero drifts of the X, Y, and Z axes of the gyroscope, and δA x , δA y and δA z are the zero offsets of the X, Y, and Z axes of the accelerometer, respectively. The W noise matrix consists of the white noise of the gyroscope and accelerometer:
W=[w
gx,w
gy,w
gz,w
ax,w
ay,w
az]
T (2)
W=[w gx ,w gy ,w gz ,w ax ,w ay ,w az ] T (2)
状态转移矩阵:State transition matrix:
Φ
k+1/k≈E
n+F
k·Δt (3)
Φ k+1/k ≈E n +F k ·Δt (3)
其中,E
n是单位矩阵,F
k是组合系统误差方程矩阵,F=[f
i,j],i,j=1,...15,f
i,j非零项为:f
1,9=n
N;f
1,13=c11;f
1,14=c12;f
1,15=c13;f
2,7=-f
1,8=n
h;f
2,13=c21;f
2,14=c22;f
2,15=c23;f
3,8=-f
2,9=n
E;f
3,13=c31;f
3,14=c32;f
3,15=c33;
f
6,3=1;f
7,9=-f
9,7=-ω
N,
f
8,9=-f
9,8=ω
E,
f
8,10=c21;f
8,11=c22;f
8,12=c23;
Wherein, E n is a unit matrix, F k is the combination of systematic error matrix equation, F = [f i, j ], i, j = 1, ... 15, f i, j nonzero entries as: f 1,9 = N N ; f 1,13 = c11; f 1,14 = c12; f 1,15 = c13; f 2,7 = -f 1,8 = n h ; f 2,13 = c21; f 2,14 = C22; f 2,15 = c23; f 3,8 = -f 2,9 = n E ; f 3,13 = c31; f 3,14 = c32; f 3,15 = c33; f 6,3 = 1; f 7,9 = -f 9,7 = -ω N , f 8,9 =-f 9,8 =ω E , f 8,10 = c21; f 8,11 = c22; f 8,12 = c23;
f
9,10=c31,f
9,11=c32,f
9,12=c33;f
7,10=c11;f
7,11=c12;f
7,12=c13。
f 9,10 = c31, f 9,11 = c32, f 9,12 = c33; f 7,10 = c11; f 7,11 = c12; f 7,12 = c13.
c
ij为姿态矩阵元素,U为地球旋转角速率,n
E,n
N,n
h分别为东北天方向上的加速度计比力。
c ij is the element of the attitude matrix, U is the angular rate of the earth's rotation, n E , n N , and n h are the accelerometer specific forces in the northeast sky direction respectively.
水下组合导航系统的观测方程为The observation equation of the underwater integrated navigation system is
Z
k+1=H
k+1X
k+1+V
k+1 (4)
Z k+1 = H k+1 X k+1 +V k+1 (4)
其中Z
k+1是观测向量,H
k+1是观测矩阵,V
k+1是观测噪声矩阵,具体如下:
Among them, Z k+1 is the observation vector, H k+1 is the observation matrix, and V k+1 is the observation noise matrix. The details are as follows:
其中,V
E,I和V
N,I分别是激光陀螺单轴旋转捷联惯性导航系统计算得到的东向和北向速度,V
E,DVL和V
N,BVL是多普勒计程仪输出的东向和北向,h
I是激光陀螺单轴旋转捷联惯性导航系统计算得到的深度,h
D是深度计输出的深度。观测矩阵H
k+1=[h
i,j]i=1,...5,j=1,...15中的非零项为:h
1,1=h
2,2=h
3,3=1,h
1,8=-V
N,h
2,8=V
E。
Among them, V E, I and V N, I are the east and north velocities calculated by the laser gyro single-axis rotation strapdown inertial navigation system , respectively, and V E, DVL and V N, BVL are output from the Doppler log. For east and north directions, h I is the depth calculated by the single-axis rotation strapdown inertial navigation system of the laser gyro, and h D is the depth output by the depth gauge. The non-zero items in the observation matrix H k+1 =[h i,j ]i=1,...5,j=1,...15 are: h 1,1 = h 2, 2 = h 3, 3 = 1, h 1,8 = -V N , h 2,8 = V E.
1个计算机接收重力传感器信息、温度信息,水下导航与重力测量一体化系统的位置、速度和姿态信息,计算出当地地理坐标系下的重力传感器的比力,对重力传感器比力进行重力测量改正,如厄特弗斯改正,零点漂移改正、滞后效应改正等,对改正后的数据进行Fir低通滤波处理,得到当地重力值。A computer receives the gravity sensor information, temperature information, the position, speed and attitude information of the integrated underwater navigation and gravity measurement system, calculates the gravity sensor's specific force in the local geographic coordinate system, and conducts the gravity measurement of the gravity sensor's specific force Correction, such as Uttforth correction, zero drift correction, hysteresis effect correction, etc., perform Fir low-pass filtering on the corrected data to obtain the local gravity value.
本发明说明书中未作详细描述的内容属于本领域专业技术人员的公知技术,应当指出:对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和等同替换,这些对本发明权利要求进行改进和等同替换后的技术方案,均落入本发明的保护范围。The content not described in detail in the specification of the present invention belongs to the well-known technology of those skilled in the art. It should be pointed out that for those of ordinary skill in the art, several improvements can be made without departing from the principle of the present invention. And equivalent replacements, these technical solutions after improvements and equivalent replacements to the claims of the present invention all fall into the protection scope of the present invention.