CN113049062A - 一种gnss直反射信号载波干涉测量湖泊水位的装置 - Google Patents

一种gnss直反射信号载波干涉测量湖泊水位的装置 Download PDF

Info

Publication number
CN113049062A
CN113049062A CN202110270801.XA CN202110270801A CN113049062A CN 113049062 A CN113049062 A CN 113049062A CN 202110270801 A CN202110270801 A CN 202110270801A CN 113049062 A CN113049062 A CN 113049062A
Authority
CN
China
Prior art keywords
signals
module
direct
lake
water level
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202110270801.XA
Other languages
English (en)
Other versions
CN113049062B (zh
Inventor
王峰
杨东凯
张波
张国栋
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beihang University
Original Assignee
Beihang University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Beihang University filed Critical Beihang University
Priority to CN202110270801.XA priority Critical patent/CN113049062B/zh
Publication of CN113049062A publication Critical patent/CN113049062A/zh
Application granted granted Critical
Publication of CN113049062B publication Critical patent/CN113049062B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01F—MEASURING VOLUME, VOLUME FLOW, MASS FLOW OR LIQUID LEVEL; METERING BY VOLUME
    • G01F23/00—Indicating or measuring liquid level or level of fluent solid material, e.g. indicating in terms of volume or indicating by means of an alarm
    • G01F23/22—Indicating or measuring liquid level or level of fluent solid material, e.g. indicating in terms of volume or indicating by means of an alarm by measuring physical variables, other than linear dimensions, pressure or weight, dependent on the level to be measured, e.g. by difference of heat transfer of steam or water
    • G01F23/28—Indicating or measuring liquid level or level of fluent solid material, e.g. indicating in terms of volume or indicating by means of an alarm by measuring physical variables, other than linear dimensions, pressure or weight, dependent on the level to be measured, e.g. by difference of heat transfer of steam or water by measuring the variations of parameters of electromagnetic or acoustic waves applied directly to the liquid or fluent solid material
    • G01F23/284—Electromagnetic waves
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/13—Receivers
    • G01S19/33—Multimode operation in different systems which transmit time stamped messages, e.g. GPS/GLONASS
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/38—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
    • G01S19/39—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/42—Determining position
    • G01S19/43—Determining position using carrier phase measurements, e.g. kinematic positioning; using long or short baseline interferometry

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Electromagnetism (AREA)
  • Thermal Sciences (AREA)
  • Fluid Mechanics (AREA)
  • Position Fixing By Use Of Radio Waves (AREA)

Abstract

本发明公开了一种GNSS直反射信号载波干涉测量湖泊水位的装置。该装置中，双频右旋圆极化天线接收导航系统中N颗卫星的直射信号；导航系统包括GPS卫星导航系统和北斗卫星导航系统；双频左旋圆极化天线接收经待测湖泊的湖面反射的导航系统的反射信号；双模导航模块产生定位信息；四通道射频前端对直射信号和反射信号进行转换得到四路数字中频信号；多通道直/反协同处理模块根据定位信息对四路数字中频信号进行干涉处理，得到反射信号相对于直射信号的载波相位时延；湖泊水位反演模块用于根据载波相位时延和定位信息进行反演得到待测湖泊的水位。本发明应用于多个不同高度角的卫星信号，从而得到湖泊水位，测量准确度高，适应范围广。

Description

一种GNSS直反射信号载波干涉测量湖泊水位的装置
技术领域
本发明涉及湖泊水位测量领域，特别是涉及一种GNSS直反射信号载波干涉测量湖泊水位的装置。
背景技术
利用全球导航卫星系统(Global Navigation Satellites System，GNSS)反射信号进行地球表面物理参数的反演是遥感领域的新型技术之一，具有信号源广、成本低、探测设备轻等优势。该技术通过星载、机载或岸基的特殊装置接收、处理GNSS反射信号可进行地表物理参数的探测。
湖泊水位作为重要的水文参数，是内陆遥感领域长期关注的对象之一。GNSS反射信号高度计最初作为传统高度计的替代技术被提出，通过连续接收处理GNSS直射信号和经地表反射的回波信号得到反射信号相对于直射信号的路径时延，并利用GNSS卫星、接收平台以及镜面反射点之间的几何关系实现接收平台高度的反演，最后与参考高度相减得到相对高度测量。
通常的GNSS反射信号高度计都是利用伪码时延进行高度测量。基于伪码的方法最大的不足是测量精度低，难以满足湖泊水位的测量。鉴于湖泊水面粗糙度低，经湖泊表面反射的GNSS为相干信号，可以充分利于载波相位进行时延测量。目前，基于载波相位测量反射信号相对于直射信号时延的方法是单天线干涉法，通过利用反射信号对直射信号的干扰形成的振荡进行高度测量，其最大缺点是时间分辨率低，且只能应用低高度角的卫星信号。
发明内容
基于此，有必要提供一种GNSS直反射信号载波干涉测量湖泊水位的装置，以应用于多个不同高度角的卫星信号，从而得到湖泊水位。
为实现上述目的，本发明提供了如下方案：
一种GNSS直反射信号载波干涉测量湖泊水位的装置，包括：双频右旋圆极化天线、双频左旋圆极化天线、四通道射频前端、双模导航模块、多通道直/反协同处理模块和湖泊水位反演模块；
所述双频右旋圆极化天线用于接收导航系统中N颗卫星的直射信号；所述导航系统包括GPS卫星导航系统和北斗卫星导航系统；所述双频左旋圆极化天线用于接收经待测湖泊的湖面反射的所述导航系统的反射信号；所述双模导航模块与所述双频右旋圆极化天线连接；所述双模导航模块用于接收所述直射信号，并产生定位信息；所述四通道射频前端与所述双频右旋圆极化天线和所述双频左旋圆极化天线连接；所述四通道射频前端用于对所述直射信号和所述反射信号进行转换，得到四路数字中频信号；所述多通道直/反协同处理模块分别与所述四通道射频前端和所述双模导航模块连接；所述多通道直/反协同处理模块用于根据所述定位信息对四路所述数字中频信号进行干涉处理，得到所述反射信号相对于所述直射信号的载波相位时延；所述湖泊水位反演模块分别与所述多通道直/反协同处理模块和所述双模导航模块连接；所述湖泊水位反演模块用于根据所述载波相位时延和所述定位信息进行反演得到所述待测湖泊的水位。
可选的，所述定位信息包括伪随机噪声号、伪距、载波频率、方位角、高度角、纬度、经度和高度。
可选的，所述GNSS直反射信号载波干涉测量湖泊水位的装置，还包括：与所述湖泊水位反演模块和所述双模导航模块连接的湖泊水位监控模块；
所述湖泊水位监控模块用于对所述定位信息和所述待测湖泊的水位进行可视化显示，并当所述待测湖泊的水位超过预设水位时进行告警。
可选的，所述GNSS直反射信号载波干涉测量湖泊水位的装置，还包括：与所述湖泊水位监控模块连接的信息管理服务器；
所述信息管理服务器用于对所述定位信息、所述待测湖泊的水位和所述告警进行存储。
可选的，所述多通道直/反协同处理模块包括卫星选择模块和N个四通道直/反射信号干涉处理模块；
所述卫星选择模块用于根据所述伪随机噪声号、所述伪距、所述载波频率、所述方位角、所述高度角、方位角范围和高度角范围从四路所述数字中频信号中选择镜面反射点对双频左旋圆极化天线可见的卫星信号，得到四路可见数字中频信号；
所述四通道直/反射信号干涉处理模块用于对四路所述可见数字中频信号进行干涉处理，得到所述反射信号相对于所述直射信号的载波相位时延。
可选的，所述湖泊水位反演模块包括频谱估计模块、相对高度反演模块和水位计算模块；
所述频谱估计模块用于根据所述高度角，采用谱分析方法对N个所述载波相位时延进行处理，得到N个频谱；
所述相对高度反演模块用于搜索N个所述频谱的峰值，并根据N个所述峰值处的频率反演得到N个相对水位高度；
所述水位计算模块用于根据所述相对水位高度和参考面高度计算所述待测湖泊的水位。
可选的，所述四通道射频前端包括GPS射频前端和北斗射频前端；
所述GPS射频前端用于对所述GPS卫星导航系统中卫星的直射信号和所述GPS卫星导航系统中卫星的反射信号进行数字化处理，得到两路GPS数字中频信号；所述数字化处理包括下变频、滤波、增益控制以及采样量化；
所述北斗射频前端用于对所述北斗卫星导航系统中卫星的直射信号和所述北斗卫星导航系统中卫星的反射信号进行数字化处理，得到两路北斗数字中频信号；两路所述GPS数字中频信号和两路所述北斗数字中频信号构成四路所述数字中频信号。
可选的，所述四通道直/反射信号干涉处理模块，包括：载波发生器、本地码发生器、载波剥离模块、伪码剥离模块、相干积分器、非相干累加模块和相位计算模块；
所述载波发生器用于根据所述卫星选择模块输出的所述载波频率产生本地载波；
所述本地码发生器用于根据所述卫星选择模块输出的所述伪随机噪声号和所述伪距产生本地伪码；
所述载波剥离模块用于将四路所述可见数字中频信号中的直射信号分别与所述本地载波相乘，得到两路载波剥离直射信号，并将四路所述可见数字中频信号中的反射信号分别与所述本地载波相乘，得到两路载波剥离反射信号；
所述伪码剥离模块用于将两路所述载波剥离直射信号分别与所述本地伪码相乘，将两路所述载波剥离反射信号分别与所述本地伪码相乘，得到四路干涉信号；
所述相干积分器用于对所述干涉信号进行相干积分；
所述非相干累加模块用于对相干积分后的复数进行多次非相干累加，得到相关功率值；
所述相位计算模块用于根据所述相关功率值计算所述反射信号相对于所述直射信号的载波相位时延。
可选的，所述湖泊水位监控模块与所述多通道直/反协同处理模块连接；所述多通道直/反协同处理模块用于接收所述湖泊水位监控模块发送的方位角范围和高度角范围。
可选的，所述信息管理服务器与IP访问端口连接。
与现有技术相比，本发明的有益效果是：
本发明提出了一种GNSS直反射信号载波干涉测量湖泊水位的装置。该装置通过双频右旋圆极化(Right-Handed Circular Polarization，RHCP)天线和双频左旋圆极化(Left-Handed Circular Polarization，LCHP)天线两个天线对同一观测区域的GPS/北斗系统的多颗卫星的直射/反射信号进行接收、处理，从而测量湖泊水位，其中，多通道直/反协同处理模块利用N颗卫星的直射信号和反射信号进行干涉处理得到反射信号相对于直射信号的载波相位时延，湖泊水位反演模块通过N个载波相位时延构成的集合反演得到待测湖泊的水位。该装置不限于应用在低高度角的卫星信号，可以用于各种不同高度角的卫星信号，从而得到湖泊水位，适应范围更广；该装置利用载波相位时延进行湖泊水平的测量，测量精度高；该装置仅是信号接收装置，装置结构简单、成本低、功耗低；该装置利用GNSS信号作为信号源，可实施全天候观测。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
图1为本发明实施例提供的GNSS直反射信号载波干涉测量湖泊水位的装置的结构示意图；
图2为本发明实施例提供的多通道直/反协同处理模块的框图；
图3为本发明实施例提供的四通道直/反射信号干涉处理模块的框图；
图4为本发明实施例提供的湖泊水位反演模块的框图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
为使本发明的上述目的、特征和优点能够更加明显易懂，下面结合附图和具体实施方式对本发明作进一步详细的说明。
本发明利用N(N＝n+m；n为反射天线可见的GPS卫星个数，m为反射天线可见的北斗卫星个数)颗导航卫星反射信号通过干涉处理得到n+m载波相位，然后对n+m载波相位进行湖泊水位测测量。
图1为本发明实施例提供的GNSS直反射信号载波干涉测量湖泊水位的装置的结构示意图。参见图1，本实施例的GNSS直反射信号载波干涉测量湖泊水位的装置，包括：双频右旋圆极化天线、双频左旋圆极化天线、四通道射频前端、双模导航模块(导航模块)、多通道直/反协同处理模块和湖泊水位反演模块。
本实施例中，所述双频右旋圆极化天线为GPS L1/北斗B1双频右旋圆极化天线(RHCP)；所述双频左旋圆极化天线为GPS L1/北斗B1双频左旋圆极化天线。所述双频右旋圆极化天线对天安置，用于接收导航系统中N颗卫星的直射信号(射频电压信号)；所述导航系统包括GPS卫星导航系统L1和北斗卫星导航系统B1；所述双频左旋圆极化天线面向海面安置，用于接收经待测湖泊的湖面反射的所述导航系统的反射信号(射频电压信号)。
具体的，双频右旋圆极化天线用于接收GPS L1/北斗B1直射电磁信号，并将该电磁信号转换为射频电压信号Sd(t)，即：
Figure BDA0002974299830000051
其中，Adi为第i颗卫星GNSS信号幅度；Di(t)表示第i颗卫星GNSS信号调制的数据；Ci(t)为GNSS信号伪随机码；fi为第i颗卫星GNSS信号的载波频率；
Figure BDA0002974299830000061
为第i颗卫星GNSS信号初始相位，t表示时间。
双频左旋圆极化天线用于接收经湖面反射的GPS L1/北斗B1电磁信号，并将该电磁信号转换为射频电压信号Sr(t)，即：
Figure BDA0002974299830000062
其中，
Figure BDA0002974299830000063
为第i颗GNSS卫星信号在湖面的反射率；τi为第i颗卫星经湖面反射的GNSS信号相对于直射GNSS信号的时延。假设双频左旋圆极化天线的高度远小于伪码长度，则：
Figure BDA0002974299830000064
本实施例中，所述双模导航模块为GPS L1/北斗B1双模导航模块。所述双模导航模块与所述双频右旋圆极化天线连接。所述双模导航模块用于接收所述直射信号以完成导航定位，从而产生定位信息；所述定位信息包括伪随机噪声号(pseudo random noise，prn)、伪距、载波频率、高度角和方位角{prni,ρi,fi,θi,φi}，以及纬度、经度和高度{lat,lon,h}，其中，i表示第i颗卫星GNSS信号。其中，双模导航模块可以采用成熟的导航SOC芯片(Systemon Chip)。
本实施例中，所述四通道射频前端与所述双频右旋圆极化天线和所述双频左旋圆极化天线连接。所述四通道射频前端用于对所述直射信号和所述反射信号进行转换，得到四路数字中频信号。其中，四通道射频前端中两个通道连接双频右旋圆极化天线，另外两个通道连接双频左旋圆极化天线，输出四路数字中频信号给直/反信号协同处理模块完成信号处理和载波相位提取。
在实际应用中，所述四通道射频前端包括GPS射频前端和北斗射频前端，两个子射频前端由两个结构完全一样的射频通道组成，分别对GPS L1/北斗B1直射、反射信号进行数字化处理。其中，所述GPS射频前端用于对所述GPS卫星导航系统中卫星的直射信号和所述GPS卫星导航系统中卫星的反射信号进行数字化处理，得到两路GPS数字中频信号；所述数字化处理包括下变频、滤波、增益控制以及采样量化。所述北斗射频前端用于对所述北斗卫星导航系统中卫星的直射信号和所述北斗卫星导航系统中卫星的反射信号进行数字化处理，得到两路北斗数字中频信号；两路所述GPS数字中频信号和两路所述北斗数字中频信号构成四路所述数字中频信号。
具体的，四通道射频前端接收RHCP天线和LHCP天线输出的射频信号Sd(t)和Sr(t)，输出四路数字中频信号，分别对应GPS直射数字中频信号、GPS反射数字中频信号、北斗直射数字中频信号和北斗反射数字中频信号。GPS和北斗信号处理完全一样，为了描述方便，对中频信号不做GPS和北斗符合上的区分，假设直射中频信号为Sd(n)，反射中频信号为Sr(n)，即：
Figure BDA0002974299830000071
Figure BDA0002974299830000072
其中，n为自然数；Ts为采样率。
本实施例中，所述多通道直/反协同处理模块分别与所述四通道射频前端和所述双模导航模块连接；所述多通道直/反协同处理模块用于根据所述定位信息对四路所述数字中频信号进行干涉处理，得到所述反射信号相对于所述直射信号的载波相位时延。其中，所述多通道直/反协同处理模块首先接收定位信息选择直/反射信号协同处理的N路卫星信号，然后对四路所述数字中频信号进行四通道直/反干涉处理得到N个载波相位时延
Figure BDA0002974299830000073
。
实际应用中，所述多通道直/反协同处理模块包括卫星选择模块和N个四通道直/反射信号干涉处理模块(4通道直/反干涉模块)，如图2所示。所述卫星选择模块用于根据所述伪随机噪声号、所述伪距、所述载波频率、所述方位角、所述高度角、方位角范围和高度角范围从四路所述数字中频信号中选择镜面反射点对双频左旋圆极化天线可见的卫星信号，得到四路可见数字中频信号。所述四通道直/反射信号干涉处理模块用于对四路所述可见数字中频信号进行干涉处理，得到所述反射信号相对于所述直射信号的载波相位时延。
如图3所示，所述四通道直/反射信号干涉处理模块，包括：载波发生器、本地码发生器、载波剥离模块、伪码剥离模块、相干积分器、非相干累加模块和相位计算模块。所述载波发生器用于根据所述卫星选择模块输出的所述载波频率产生本地载波。所述本地码发生器用于根据所述卫星选择模块输出的所述伪随机噪声号和所述伪距产生本地伪码。所述载波剥离模块用于将四路所述可见数字中频信号中的直射信号分别与所述本地载波相乘，得到两路载波剥离直射信号，并将四路所述可见数字中频信号中的反射信号分别与所述本地载波相乘，得到两路载波剥离反射信号。所述伪码剥离模块用于将两路所述载波剥离直射信号分别与所述本地伪码相乘，将两路所述载波剥离反射信号分别与所述本地伪码相乘，得到四路干涉信号。所述相干积分器用于对所述干涉信号进行相干积分。所述相干积分器为1ms相干积分器。所述非相干累加模块用于对相干积分后的复数进行多次非相干累加，得到相关功率值。所述相位计算模块用于根据所述相关功率值计算所述反射信号相对于所述直射信号的载波相位时延。
具体的，卫星选择模块接收GPS L1/北斗B1导航模块输出的伪随机噪声号、高度角和方位角{prni,θi,φi}，以及接收高度角范围信息(截止信息)和方位角范围信息{θmask,φmask}，从而对GPS L1/北斗B1双频左旋圆极化天线视野范围内的GPS/北斗卫星进行预测，输出双频左旋圆极化天线视野范围内的卫星伪随机噪声号组成的集合{prn1,prn2,…,prnN}。预判决准则为：
θmin≤θi≤θmax,θmin,θmax∈θmask
φmin≤φi≤φmax,φmin,φmax∈φmask。
如图3所示，载波发生器根据卫星选择模块输出的载波频率fi产生本地载波，即：
Figure BDA0002974299830000081
本地码发生器根据卫星选择模块输出的伪随机噪声号和伪距信息{prni,ρi}产生本地伪码，即：
Scode_i(n)＝Ci(nTS)；
载波剥离模块将直/反射信号和本地载波进行相乘，即：
Figure BDA0002974299830000091
其中，Ads为第s颗卫星GNSS信号幅度；Ds(t)表示第s颗卫星GNSS信号调制的数据；Cs(t)为GNSS信号伪随机码；fs为第s颗卫星GNSS信号的载波频率；
Figure BDA0002974299830000092
为第s颗卫星GNSS信号初始相位。
Figure BDA0002974299830000093
其中，
Figure BDA0002974299830000094
为第s颗GNSS卫星信号湖面反射率。伪码剥离模块将直/反射信号和本地伪码进行相乘，即：
Figure BDA0002974299830000095
Figure BDA0002974299830000096
4路干涉处理模块首先利用载波剥离和伪码剥离后的直/反射信号产生4路干涉信号，即：
Figure BDA0002974299830000097
其次对4路干涉信号进行1ms相干积分，根据伪码特性，可得4路1ms相干积分后的信号，分别为第1路1ms相干积分后的信号复数相关值Icoh_i_1(m)、第2路1ms相干积分后的信号复数相关值Icoh_i_2(m)、第3路1ms相干积分后的信号复数相关值Icoh_i_3(m)和第4路1ms相干积分后的信号复数相关值Icoh_i_4(m)，即：
Figure BDA0002974299830000101
其中，Ns表示采样点总数。然后对1ms相干积分后的复数相关值进行M次非相干累加，输出4路相关功率值，分别为第1路相关功率值Ii1、第2路相关功率值Ii2、第3路相关功率值Ii3和第4路相关功率值Ii4，即：
Figure BDA0002974299830000102
相位计算模块根据M次非相干累加输出的4路相关功率值估计反射信号相对于直射信号的载波相位时延，即：
Figure BDA0002974299830000103
本实施例中，所述湖泊水位反演模块分别与所述多通道直/反协同处理模块和所述双模导航模块连接；所述湖泊水位反演模块用于根据所述载波相位时延
Figure BDA0002974299830000104
和所述定位信息进行反演得到所述待测湖泊的水位。
如图4所示，所述湖泊水位反演模块包括频谱估计模块、相对高度反演模块和水位计算模块。所述频谱估计模块用于根据所述高度角，采用谱分析方法对N个所述载波相位时延进行处理，得到N个频谱。所述相对高度反演模块用于搜索N个所述频谱的峰值，并根据N个所述峰值处的频率反演得到N个相对水位高度。所述水位计算模块用于根据所述相对水位高度和参考面高度计算所述待测湖泊的水位。
具体的，首先，根据在一定时间段内相位多通道直/反信号协同处理模块输出的N个载波相位时延的时间序列
Figure BDA0002974299830000111
得到个N相位时序的频谱。本实施例采用Lomb-Scargle谱估计方法，估计过程中时间轴由时间转化为高度角的正弦函数t＝sin[θi(t)]，具体过程表示为：
Figure BDA0002974299830000112
然后，搜索N个Lomb-Scargle谱的峰值位置{fLS_1,fLS_2,…,fLS_N}，并根据频谱峰值频率与高度关系反演得到接收天线相对于水位的高度，即：
Figure BDA0002974299830000113
其中，λ1,λ2…λN表示第1,2…,N颗GNSS卫星信号的波长。
最后，水位计算模块根据反演得到的相对水位高度{hr1,hr2,…,hrN}和参考面高度href计算水位高度，得到待测湖泊的水位，即：
Figure BDA0002974299830000114
在实际应用中，所述GNSS直反射信号载波干涉测量湖泊水位的装置，还包括：与所述湖泊水位反演模块和所述双模导航模块连接的湖泊水位监控模块；所述湖泊水位监控模块用于对所述定位信息和所述待测湖泊的水位进行可视化显示，并当所述待测湖泊的水位超过预设水位时进行告警。所述湖泊水位监控模块还可以与所述多通道直/反协同处理模块连接；所述多通道直/反协同处理模块用于接收所述湖泊水位监控模块发送的方位角范围和高度角范围。
湖泊水位监控模块提供了湖泊水位测量与用户的交互界面，接收湖泊水位反演模块的待测湖泊的水位Δh、GPS L1/北斗B1双模导航模块的定位信息{lat,lon,h}进行水位和监测站信息的显示和存储，并当水位超过预设水位时产生告警。此外根据LHCP天线的架设场景，可以设置观测视场的截止范围{θmask,φmask}。
在实际应用中，所述GNSS直反射信号载波干涉测量湖泊水位的装置，还包括：与所述湖泊水位监控模块连接的信息管理服务器；所述信息管理服务器用于对所述定位信息、所述待测湖泊的水位和所述告警进行存储。所述信息管理服务器还可以与IP访问端口连接。本实施例设置信息管理服务器，实现了对水位、定位、告警、日志以及用户信息的存储和管理，并能接收IP接口的访问。
本实施例通过双频右旋圆极化天线和双频左旋圆极化天线两个天线对同一观测区域的GPS/北斗系统的多颗卫星的直射/反射信号进行接收、处理，从而测量湖泊水位。首先在每个积分周期内，利用4通道干涉法估计N颗卫星反射信号相对于直射信号的载波相位时延，即
Figure BDA0002974299830000121
；然后利用载波相位反演得到接收平台的高度，并与参考面高度相减得到湖面高度(待测湖泊的水位)。本实施例适应范围广；利用载波相位进行湖泊水平的测量，测量精度高；装置结构简单、成本低、功耗低；利用GNSS信号作为信号源，可实施全天候观测。
本说明书中各个实施例采用递进的方式描述，每个实施例重点说明的都是与其他实施例的不同之处，各个实施例之间相同相似部分互相参见即可。
本文中应用了具体个例对本发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的方法及其核心思想；同时，对于本领域的一般技术人员，依据本发明的思想，在具体实施方式及应用范围上均会有改变之处。综上所述，本说明书内容不应理解为对本发明的限制。

Claims (10)

1.一种GNSS直反射信号载波干涉测量湖泊水位的装置，其特征在于，包括：双频右旋圆极化天线、双频左旋圆极化天线、四通道射频前端、双模导航模块、多通道直/反协同处理模块和湖泊水位反演模块；
所述双频右旋圆极化天线用于接收导航系统中N颗卫星的直射信号；所述导航系统包括GPS卫星导航系统和北斗卫星导航系统；所述双频左旋圆极化天线用于接收经待测湖泊的湖面反射的所述导航系统的反射信号；所述双模导航模块与所述双频右旋圆极化天线连接；所述双模导航模块用于接收所述直射信号，并产生定位信息；所述四通道射频前端与所述双频右旋圆极化天线和所述双频左旋圆极化天线连接；所述四通道射频前端用于对所述直射信号和所述反射信号进行转换，得到四路数字中频信号；所述多通道直/反协同处理模块分别与所述四通道射频前端和所述双模导航模块连接；所述多通道直/反协同处理模块用于根据所述定位信息对四路所述数字中频信号进行干涉处理，得到所述反射信号相对于所述直射信号的载波相位时延；所述湖泊水位反演模块分别与所述多通道直/反协同处理模块和所述双模导航模块连接；所述湖泊水位反演模块用于根据所述载波相位时延和所述定位信息进行反演得到所述待测湖泊的水位。
2.根据权利要求1所述的一种GNSS直反射信号载波干涉测量湖泊水位的装置，其特征在于，所述定位信息包括伪随机噪声号、伪距、载波频率、方位角、高度角、纬度、经度和高度。
3.根据权利要求1所述的一种GNSS直反射信号载波干涉测量湖泊水位的装置，其特征在于，还包括：与所述湖泊水位反演模块和所述双模导航模块连接的湖泊水位监控模块；
所述湖泊水位监控模块用于对所述定位信息和所述待测湖泊的水位进行可视化显示，并当所述待测湖泊的水位超过预设水位时进行告警。
4.根据权利要求3所述的一种GNSS直反射信号载波干涉测量湖泊水位的装置，其特征在于，还包括：与所述湖泊水位监控模块连接的信息管理服务器；
所述信息管理服务器用于对所述定位信息、所述待测湖泊的水位和所述告警进行存储。
5.根据权利要求2所述的一种GNSS直反射信号载波干涉测量湖泊水位的装置，其特征在于，所述多通道直/反协同处理模块包括卫星选择模块和N个四通道直/反射信号干涉处理模块；
所述卫星选择模块用于根据所述伪随机噪声号、所述伪距、所述载波频率、所述方位角、所述高度角、方位角范围和高度角范围从四路所述数字中频信号中选择镜面反射点对双频左旋圆极化天线可见的卫星信号，得到四路可见数字中频信号；
所述四通道直/反射信号干涉处理模块用于对四路所述可见数字中频信号进行干涉处理，得到所述反射信号相对于所述直射信号的载波相位时延。
6.根据权利要求2所述的一种GNSS直反射信号载波干涉测量湖泊水位的装置，其特征在于，所述湖泊水位反演模块包括频谱估计模块、相对高度反演模块和水位计算模块；
所述频谱估计模块用于根据所述高度角，采用谱分析方法对N个所述载波相位时延进行处理，得到N个频谱；
所述相对高度反演模块用于搜索N个所述频谱的峰值，并根据N个所述峰值处的频率反演得到N个相对水位高度；
所述水位计算模块用于根据所述相对水位高度和参考面高度计算所述待测湖泊的水位。
7.根据权利要求1所述的一种GNSS直反射信号载波干涉测量湖泊水位的装置，其特征在于，所述四通道射频前端包括GPS射频前端和北斗射频前端；
所述GPS射频前端用于对所述GPS卫星导航系统中卫星的直射信号和所述GPS卫星导航系统中卫星的反射信号进行数字化处理，得到两路GPS数字中频信号；所述数字化处理包括下变频、滤波、增益控制以及采样量化；
所述北斗射频前端用于对所述北斗卫星导航系统中卫星的直射信号和所述北斗卫星导航系统中卫星的反射信号进行数字化处理，得到两路北斗数字中频信号；两路所述GPS数字中频信号和两路所述北斗数字中频信号构成四路所述数字中频信号。
8.根据权利要求5所述的一种GNSS直反射信号载波干涉测量湖泊水位的装置，其特征在于，所述四通道直/反射信号干涉处理模块，包括：载波发生器、本地码发生器、载波剥离模块、伪码剥离模块、相干积分器、非相干累加模块和相位计算模块；
所述载波发生器用于根据所述卫星选择模块输出的所述载波频率产生本地载波；
所述本地码发生器用于根据所述卫星选择模块输出的所述伪随机噪声号和所述伪距产生本地伪码；
所述载波剥离模块用于将四路所述可见数字中频信号中的直射信号分别与所述本地载波相乘，得到两路载波剥离直射信号，并将四路所述可见数字中频信号中的反射信号分别与所述本地载波相乘，得到两路载波剥离反射信号；
所述伪码剥离模块用于将两路所述载波剥离直射信号分别与所述本地伪码相乘，将两路所述载波剥离反射信号分别与所述本地伪码相乘，得到四路干涉信号；
所述相干积分器用于对所述干涉信号进行相干积分；
所述非相干累加模块用于对相干积分后的复数进行多次非相干累加，得到相关功率值；
所述相位计算模块用于根据所述相关功率值计算所述反射信号相对于所述直射信号的载波相位时延。
9.根据权利要求3所述的一种GNSS直反射信号载波干涉测量湖泊水位的装置，其特征在于，所述湖泊水位监控模块与所述多通道直/反协同处理模块连接；所述多通道直/反协同处理模块用于接收所述湖泊水位监控模块发送的方位角范围和高度角范围。
10.根据权利要求4所述的一种GNSS直反射信号载波干涉测量湖泊水位的装置，其特征在于，所述信息管理服务器与IP访问端口连接。
CN202110270801.XA 2021-03-12 2021-03-12 一种gnss直反射信号载波干涉测量湖泊水位的装置 Active CN113049062B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110270801.XA CN113049062B (zh) 2021-03-12 2021-03-12 一种gnss直反射信号载波干涉测量湖泊水位的装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110270801.XA CN113049062B (zh) 2021-03-12 2021-03-12 一种gnss直反射信号载波干涉测量湖泊水位的装置

Publications (2)

Publication Number Publication Date
CN113049062A true CN113049062A (zh) 2021-06-29
CN113049062B CN113049062B (zh) 2022-04-15

Family

ID=76511986

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110270801.XA Active CN113049062B (zh) 2021-03-12 2021-03-12 一种gnss直反射信号载波干涉测量湖泊水位的装置

Country Status (1)

Country Link
CN (1) CN113049062B (zh)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113671552A (zh) * 2021-08-24 2021-11-19 北京航空航天大学 一种陆基gnss反射计的河流流速与水位反演方法
CN113671210A (zh) * 2021-08-19 2021-11-19 北京航空航天大学 陆基gnss直反射信号载波干涉测量河流流速的装置
CN114137575A (zh) * 2022-02-08 2022-03-04 浙江国遥地理信息技术有限公司 顾及卫星偏差和载噪比弧段影响的洪水探测方法
CN114355421A (zh) * 2021-12-22 2022-04-15 杭州电子科技大学 一种基于北斗卫星l4和cmc组合观测值洪水探测方法
CN115165029A (zh) * 2022-07-01 2022-10-11 中国水利水电科学研究院 基于双天线测姿的北斗水位监测装置、方法及系统
US11719732B1 (en) * 2022-07-25 2023-08-08 Divirod, Inc. Reflectometer sensor
CN116699658A (zh) * 2023-08-04 2023-09-05 山东大学 一种基于双频反射信号组合的gnss-r载波相位海面测高的方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20020130813A1 (en) * 2001-01-19 2002-09-19 Neira Manuel Martin Ocean altimetry interferometric method and device using GNSS signals
CN101608915A (zh) * 2008-06-17 2009-12-23 厦门瀛寰电子科技有限公司 全新的海浪观测装置
CN104678371A (zh) * 2015-02-11 2015-06-03 北京航空航天大学 一种基于时延修正的海面高度测量装置
CN104677463A (zh) * 2015-02-11 2015-06-03 北京航空航天大学 一种基于统计量的gnss-r多源海面测高装置
CN104678388A (zh) * 2015-02-11 2015-06-03 北京航空航天大学 一种用于提高测高精度的gnss-r正交化交错处理装置

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20020130813A1 (en) * 2001-01-19 2002-09-19 Neira Manuel Martin Ocean altimetry interferometric method and device using GNSS signals
CN101608915A (zh) * 2008-06-17 2009-12-23 厦门瀛寰电子科技有限公司 全新的海浪观测装置
CN104678371A (zh) * 2015-02-11 2015-06-03 北京航空航天大学 一种基于时延修正的海面高度测量装置
CN104677463A (zh) * 2015-02-11 2015-06-03 北京航空航天大学 一种基于统计量的gnss-r多源海面测高装置
CN104678388A (zh) * 2015-02-11 2015-06-03 北京航空航天大学 一种用于提高测高精度的gnss-r正交化交错处理装置

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
王峰 等: "机载导航卫星反射信号海面高度测量研究", 《航空学报》 *
高洪兴 等: "基于GNSS卫星反射信号的海冰厚度探测", 《电子与信息学报》 *

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113671210A (zh) * 2021-08-19 2021-11-19 北京航空航天大学 陆基gnss直反射信号载波干涉测量河流流速的装置
CN113671552A (zh) * 2021-08-24 2021-11-19 北京航空航天大学 一种陆基gnss反射计的河流流速与水位反演方法
CN114355421A (zh) * 2021-12-22 2022-04-15 杭州电子科技大学 一种基于北斗卫星l4和cmc组合观测值洪水探测方法
CN114355421B (zh) * 2021-12-22 2023-08-18 杭州电子科技大学 一种基于北斗卫星l4和cmc组合观测值洪水探测方法
CN114137575A (zh) * 2022-02-08 2022-03-04 浙江国遥地理信息技术有限公司 顾及卫星偏差和载噪比弧段影响的洪水探测方法
CN115165029A (zh) * 2022-07-01 2022-10-11 中国水利水电科学研究院 基于双天线测姿的北斗水位监测装置、方法及系统
CN115165029B (zh) * 2022-07-01 2023-03-28 中国水利水电科学研究院 基于双天线测姿的北斗水位监测装置、方法及系统
US11719732B1 (en) * 2022-07-25 2023-08-08 Divirod, Inc. Reflectometer sensor
CN116699658A (zh) * 2023-08-04 2023-09-05 山东大学 一种基于双频反射信号组合的gnss-r载波相位海面测高的方法

Also Published As

Publication number Publication date
CN113049062B (zh) 2022-04-15

Similar Documents

Publication Publication Date Title
CN113049062B (zh) 一种gnss直反射信号载波干涉测量湖泊水位的装置
CN113049777B (zh) 一种gnss直反射信号载波干涉测量土壤湿度的装置
US6549165B2 (en) Ocean altimetry interferometric method and device using GNSS signals
US8411726B2 (en) Method and apparatus for performing signal correlation at multiple resolutions to mitigate multipath interference
Yu et al. Snow depth estimation based on combination of pseudorange and carrier phase of GNSS dual-frequency signals
Zhang et al. Tide variation monitoring based improved GNSS-MR by empirical mode decomposition
Zhang et al. Feasibility of code-level altimetry using coastal BeiDou reflection (BeiDou-R) setups
Ban et al. GEO-satellite-based reflectometry for soil moisture estimation: Signal modeling and algorithm development
CN104677463A (zh) 一种基于统计量的gnss-r多源海面测高装置
CN113805208A (zh) 一种适用于导航接收机的gnss-ir测高方法
CN111521650B (zh) 土壤湿度测量中直射信号和反射信号的相位提取方法
CN116699658A (zh) 一种基于双频反射信号组合的gnss-r载波相位海面测高的方法
Wang et al. Sea-state observation using reflected BeiDou GEO signals in frequency domain
He et al. Coastal altimetry using interferometric phase from GEO satellite in quasi-zenith satellite system
Yan et al. Soil moisture estimation based on BeiDou B1 interference signal analysis
Li et al. Real-time processing of reflected GNSS signals for remote sensing: System and experiments
Connselman III et al. Backpack VLBI terminal with subscentimeter capability
Hobiger et al. Software-defined radio direct correlation GNSS reflectometry by means of GLONASS
Manandhar et al. GPS reflected signal analysis using software receiver
CN113671210B (zh) 陆基gnss直反射信号载波干涉测量河流流速的装置
CN115144884A (zh) 一种基于卫星反射信号的海面风速反演方法及芯片模组
He et al. Coastal GNSS-R phase altimetry based on the combination of L1 and L5 signals under high sea states
Garrison Anisotropy in reflected GPS measurements of ocean winds
Hobiger et al. A real-time GNSS-R system based on software-defined radio and graphics processing units
Lubeigt et al. Close-to-Ground Single Antenna GNSS-R

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant