CN113671552B - 一种陆基gnss反射计的河流流速与水位反演方法 - Google Patents
一种陆基gnss反射计的河流流速与水位反演方法 Download PDFInfo
- Publication number
- CN113671552B CN113671552B CN202110973188.8A CN202110973188A CN113671552B CN 113671552 B CN113671552 B CN 113671552B CN 202110973188 A CN202110973188 A CN 202110973188A CN 113671552 B CN113671552 B CN 113671552B
- Authority
- CN
- China
- Prior art keywords
- signal
- river
- reflected signal
- water level
- doppler frequency
- 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.)
- Active
Links
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 72
- 238000000034 method Methods 0.000 title claims abstract description 65
- 230000001427 coherent effect Effects 0.000 claims abstract description 48
- 238000009825 accumulation Methods 0.000 claims abstract description 26
- 238000001228 spectrum Methods 0.000 claims description 31
- 238000002310 reflectometry Methods 0.000 claims description 23
- 238000000819 phase cycle Methods 0.000 claims description 12
- 238000005070 sampling Methods 0.000 claims description 8
- 230000000737 periodic effect Effects 0.000 claims description 6
- 230000001186 cumulative effect Effects 0.000 claims 2
- 238000001514 detection method Methods 0.000 abstract description 11
- 238000012544 monitoring process Methods 0.000 description 15
- 238000005259 measurement Methods 0.000 description 9
- 238000012545 processing Methods 0.000 description 8
- 238000004364 calculation method Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 4
- 238000012937 correction Methods 0.000 description 3
- 238000000605 extraction Methods 0.000 description 3
- 238000006243 chemical reaction Methods 0.000 description 2
- 230000003595 spectral effect Effects 0.000 description 2
- 101100412102 Haemophilus influenzae (strain ATCC 51907 / DSM 11121 / KW20 / Rd) rec2 gene Proteins 0.000 description 1
- 101100356020 Haemophilus influenzae (strain ATCC 51907 / DSM 11121 / KW20 / Rd) recA gene Proteins 0.000 description 1
- 238000005311 autocorrelation function Methods 0.000 description 1
- 230000000295 complement effect Effects 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 230000001934 delay Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
Classifications
-
- 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/52—Determining velocity
-
- 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/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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/30—Assessment of water resources
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)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明涉及一种陆基GNSS反射计的河流流速与水位反演方法，包括：获取陆基场景中GNSS信号在河流表面的反射信号；根据反射信号利用高斯拟合法估计反射信号多普勒频率；根据反射信号多普勒频率利用反射信号双基几何关系反演河流流速；获取反射信号的相干信号和非相干信号；根据相干信号和非相干信号确定累加比；当累加比小于设定阈值时，则利用模糊场景进行水位反演，确定接收机相对于河流表面的高度；当累加比大于或者等于设定阈值时，则利用相干场景进行水位反演，确定接收机相对于河流表面的高度；根据接收机相对于河流表面的高度和接收机相对于大地的高程确定河流水位。本发明能够实现对河流流速和水位进行低成本和高精度的探测。
Description
技术领域
本发明涉及河流测量领域，特别是涉及一种陆基GNSS反射计的河流流速与水位反演方法。
背景技术
GNSS反射信号(Global Navigation Satellite System-Reflection)，利用现有的导航卫星信号通过协同接收、处理直达和经地表反射的导航卫星信号实现地表物理参数的反演和监测，是一种双基机会源探测手段，具有成本低、功耗低、重量轻的优点。其L波段在降雨条件下信号衰减小，适合于降雨等恶劣条件的监测。另外，陆基探测可利用成熟的多通道导航射频、相关器芯片完成信号处理实现快速处理和监测。
我国在河流监测方面目前常用手段为基于接触式缆道无线测流的水位站观测。该方式在财力、物力方面均消耗巨大，且很难满足野外刮风、降雨、洪水汛期等复杂环境的需求，使得我国偏远、气候恶劣地区的水文站点稀疏或者缺失，获得的河流水文信息极其有限。遥感技术的发展使得流速摄像法、甚高频(UHF)地波雷达、主被动微波雷达以及星载多光谱等非接触式河流监测手段涌现。尽管上述非接触式监测手段在一定程度上补充了较大尺度及欠发达地区河流水文信息的观测，但存在如下不足：1)均是对河流边界、水位、流速等单一参数进行针对性监测(例如甚高频地波雷达仅用于流速监测)；2)均采用河流边界、水位、流速等参数和实测河流流量建立经验关系实现河流监测；3)星载高度计、SAR、多光谱等遥感手段由于分辨率和大气层影响，在狭窄河流区域、降雨等恶劣条件受到制约，且时间采样率低、实时性差、成本和功耗高。随着日益增长的河流水文信息需求，迫切需要一种便捷、成本低廉、且支持多参数一体化实时监测的方法。
发明内容
本发明的目的是提供一种陆基GNSS反射计的河流流速与水位反演方法，以实现对河流流速和水位进行低成本和高精度的探测。
为实现上述目的，本发明提供了如下方案：
一种陆基GNSS反射计的河流流速与水位反演方法，包括：
获取陆基场景中GNSS信号在河流表面的反射信号；
根据所述反射信号利用高斯拟合法估计反射信号多普勒频率；所述反射信号多普勒频谱为估计的河流流速引起的GNSS反射信号多普勒频率；
根据所述反射信号多普勒频率利用反射信号双基几何关系反演河流流速；
获取所述反射信号的相干信号和非相干信号；
根据所述相干信号和所述非相干信号确定累加比；
判断所述累加比是否小于设定阈值，得到第一判断结果；
当所述第一判断结果表示所述累加比小于设定阈值时，则利用模糊场景进行水位反演，确定接收机相对于河流表面的高度；
当所述第一判断结果表示所述累加比大于或者等于设定阈值时，则利用相干场景进行水位反演，确定接收机相对于河流表面的高度；
根据所述接收机相对于河流表面的高度和接收机相对于大地的高程确定河流水位。
可选的，所述获取陆基场景中GNSS信号在河流表面的反射信号之后，还包括：
获取接收机中接收的反射功率和直射功率；
根据所述反射功率和所述直射功率确定反射率；
判断所述反射率是否大于反射率阈值，得到第二判断结果；
若所述第二判断结果表示所述反射率大于反射率阈值，则根据所述反射信号利用高斯拟合法估计反射信号多普勒频率；
若所述第二判断结果表示所述反射率小于或者等于反射率阈值，则返回步骤“获取接收机中接收的反射功率和直射功率”。
可选的，所述根据所述反射信号利用高斯拟合法估计反射信号多普勒频率，具体包括：
根据所述反射信号利用Welch谱估计法进行估计，得到反射信号相对于直射信号多普勒频率；
根据所述反射信号相对于直射信号多普勒频率利用高斯拟合法估计反射信号多普勒频率。
可选的，所述根据所述反射信号相对于直射信号多普勒频率利用高斯拟合法估计反射信号多普勒频率，具体包括：
根据如下公式估计反射信号多普勒频率：
其中，为估计的河流流速引起的GNSS反射信号多普勒频率，ψ(f)为反射信号相对于直射信号多普勒频率，/>为估计高斯分布标准差，f为多普勒频率，μ为多普勒偏移频谱在高斯分布下的均值，σ为多普勒偏移频谱在高斯分布下的标准差。
可选的，所述河流流速的反演公式为：
其中，为河流流速，/>为估计的河流流速引起的GNSS反射信号多普勒频率，φ为河流流向和GNSS卫星方位向之间的夹角，θ为GNSS卫星高度角，fc为GNSS卫星信号载波频率，c为光速。
可选的，所述当所述第一判断结果表示所述累加比小于设定阈值时，则利用模糊场景进行水位反演，确定接收机相对于河流表面的高度，具体包括：
采用非相干叠加法提取所述反射信号的相干成分；
根据所述相干成分进行非相干叠加得到所述反射信号的时延波形；
跟踪所述反射信号的时延波形得到反射信号相对于直射信号的时延；
根据所述反射信号相对于直射信号的时延利用双基几何关系确定接收机相对于河流表面的高度；
利用通道交叉法对所述接收机相对于河流表面的高度进行校正，确定最终的接收机相对于河流表面的高度。
可选的，所述当所述第一判断结果表示所述累加比大于或者等于设定阈值时，则利用相干场景进行水位反演，确定接收机相对于河流表面的高度，具体包括：
对所述反射信号进行载波相位估计，得到相位时序；
对所述相位时序进行相位补偿，得到补偿后的相位时序；
对补偿后的相位时序利用非周期采样的长序列谱估计法进行谱估计，确定接收机相对于河流表面的高度。
可选的，所述河流水位的计算公式如下：
Hr＝Hrec-hrec
其中，Hr为河流水位，Hrec为接收机相对大地参考坐标的高程，hrec为接收机相对于河流表面的高度。
根据本发明提供的具体实施例，本发明公开了以下技术效果：
本发明提供的陆基GNSS反射计的河流流速与水位反演方法，通过将GNSS反射信号技术应用到陆河边界与河流流速反演上，信号源丰富，成本低，通过累加比选择相干场景或者模糊场景进行行水位反演，提供反演精度。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
图1为本发明提供的陆基GNSS反射计的河流流速与水位反演方法流程图；
图2为本发明提供的陆基GNSS反射计的河流流速与水位反演方法在实际应用中的流程图；
图3为本发明提供的河流流速和水位反演时序图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
本发明的目的是提供一种陆基GNSS反射计的河流流速与水位反演方法，以实现对河流流速和水位进行低成本和高精度的探测。
为使本发明的上述目的、特征和优点能够更加明显易懂，下面结合附图和具体实施方式对本发明作进一步详细的说明。
如图1所示，本发明提供的一种陆基GNSS反射计的河流流速与水位反演方法，包括：
步骤101：获取陆基场景中GNSS信号在河流表面的反射信号。
在实际应用中，步骤101，还包括：
获取接收机中接收的反射功率和直射功率。根据所述反射功率和所述直射功率确定反射率。判断所述反射率是否大于反射率阈值，得到第二判断结果。若所述第二判断结果表示所述反射率大于反射率阈值，则根据所述反射信号利用高斯拟合法估计反射信号多普勒频率；若所述第二判断结果表示所述反射率小于或者等于反射率阈值，则返回步骤“获取接收机中接收的反射功率和直射功率”。
步骤102：根据所述反射信号利用高斯拟合法估计反射信号多普勒频率；所述反射信号多普勒频谱为估计的河流流速引起的GNSS反射信号多普勒频率。
其中步骤102，具体包括：
根据所述反射信号利用Welch谱估计法进行估计，得到反射信号相对于直射信号多普勒频率。
根据所述反射信号相对于直射信号多普勒频率利用高斯拟合法估计反射信号多普勒频率。其中，所述根据所述反射信号相对于直射信号多普勒频率利用高斯拟合法估计反射信号多普勒频率，具体包括：根据如下公式估计反射信号多普勒频率：
其中，为估计的河流流速引起的GNSS反射信号多普勒频率，ψ(f)为反射信号相对于直射信号多普勒频率，/>为估计高斯分布标准差，f为多普勒频率，μ为多普勒偏移频谱在高斯分布下的均值，σ为多普勒偏移频谱在高斯分布下的标准差。
步骤103：根据所述反射信号多普勒频率利用反射信号双基几何关系反演河流流速。其中，所述河流流速的反演公式为：
其中，为河流流速，/>为估计的河流流速引起的GNSS反射信号多普勒频率，φ为河流流向和GNSS卫星方位向之间的夹角，θ为GNSS卫星高度角，fc为GNSS卫星信号载波频率，c为光速。
步骤104：获取所述反射信号的相干信号和非相干信号。
步骤105：根据所述相干信号和所述非相干信号确定累加比。
步骤106：判断所述累加比是否小于设定阈值，得到第一判断结果；当所述第一判断结果表示所述累加比小于设定阈值时，则执行步骤107；当所述第一判断结果表示所述累加比大于或者等于设定阈值时，则执行步骤108。
步骤107：利用模糊场景进行水位反演，确定接收机相对于河流表面的高度；其中，步骤107，具体包括：
采用非相干叠加法提取所述反射信号的相干成分。
根据所述相干成分进行非相干叠加得到所述反射信号的时延波形。
跟踪所述反射信号的时延波形得到反射信号相对于直射信号的时延。
根据所述反射信号相对于直射信号的时延利用双基几何关系确定接收机相对于河流表面的高度。
利用通道交叉法对所述接收机相对于河流表面的高度进行校正，确定最终的接收机相对于河流表面的高度。
步骤108：利用相干场景进行水位反演，确定接收机相对于河流表面的高度；其中，步骤108具体包括：
对所述反射信号进行载波相位估计，得到相位时序。
对所述相位时序进行相位补偿，得到补偿后的相位时序。
对补偿后的相位时序利用非周期采样的长序列谱估计法进行谱估计，确定接收机相对于河流表面的高度。
步骤109：根据所述接收机相对于河流表面的高度和接收机相对于大地的高程确定河流水位。其中，所述河流水位的计算公式如下：
Hr＝Hrec-hrec
其中，Hr为河流水位，Hrec为接收机相对大地参考坐标的高程，hrec为接收机相对于河流表面的高度。
如图2-3所示，本发明还提供了一种陆基GNSS反射计的河流流速与水位反演方法在实际应用中的具体流程，步骤如下：
在信号层面，河流流速主要被调制在GNSS反射信号频率上。在陆基场景，由于GNSS信号到达河流表面的多普勒频率和直达信号的相同，因此反射信号相对于直达信号的多普勒频率偏差仅由河流流动引起，即Δf＝fflow，其中，△f为反射信号到达天线时的频率与直射信号到达天线时候的差值，fflow为河流流动引起的频率差值。
由于Welch法对谱估计数据进行分段时，不同数据段可存在重叠，较一般周期图法有更好的方差性能，且数据段的窗函数有多种选择可更好的控制旁瓣泄露，因此，通过Welch谱估计得到GNSS反射信号相对于直射信号的多普勒频率：
基于welch谱，通过数据统计计算，得到反射信号的多普勒频率。此处为welch谱公式。其中，w(n)为窗函数；M为窗函数长度；S为数据分段个数；ψi(f)为第i个序列基于welch谱公式的频率函数，f为频率，n为在各自序列中的点序数，i为序列数，j为虚数，ψ(f)为多序列计算后得到的多普勒函数，yi(n)为分段数据，具体为：
yi(n)＝Yrmax((i-1)K+n)
其中，Yimax(n)表示反射信号功率峰值对应的复数相干值序列；K表示分段之间的重叠的数据点数。
如图2所示，在判断符合场景的反射率后，对峰值频率求解。分析岸基场景的GNSS海面反射信号的谱特性时发现GNSS反射信号的谱近似呈高斯函数分布。其中反射率为反射功率与直射功率之比，两个值可以从接收机获取。基于此，通过高斯拟合的方法估计GNSS反射信号的频谱特征值：
其中，为谱估计的峰值频率，即估计的河流流速引起的GNSS反射信号多普勒频率。
接下来实现流速反演，GNSS反射信号相对于直达信号的多普勒频率由河流流速在GNSS反射信号径向引起，可根据GNSS反射信号双基几何关系得到：
其中，vflow为河流径流流速；φ为河流流向和GNSS卫星方位向之间的夹角；θ为GNSS卫星高度角；fc为GNSS卫星信号载波频率；c为光速。根据上式以及估计的GNSS反射信号相对于直达信号的多普勒频率反演河流流速：
其中，为河流径流流速反演值。
下一步进行相干指示计算，相对于伪码，载波测距的精度高，水位反演精度高。但是载波测量需要GNSS反射信号相位连续可估计。由于河流粗糙度的影响，GNSS反射信号既包含了相干信号，也包含了非相干信号。相干信号载波相位连续可估计，而非相干成分的相位服从均匀随机分布。因此，当GNSS反射信号主成分为相干信号，而非相干成分可忽略时，GNSS反射信号的载波相位连续可估计，在此场景可利用载波相位测量水位；当非相干成分不可忽略时，GNSS反射信号的载波相位呈现随机性，在此场景只能利用伪码测量水位。基于此，采用相干/模糊场景判决的伪码/载波协同的水位反演策略。当场景判决为相干时，采用载波相位测量；而为非相干时，采用伪码测量。本课题采用上述式定义的相干/非相干累加比进行相干和模糊场景的判决，即
其中，CIR为相干/非相干累加比，Yi(τ)为反射信号，其中分子为反射信号在一段时间内(时间序列i内)相干累加信号，分母为反射信号在一段时间内(时间序列i内)非相干累加信号，Tc为判决阈值。通过理论研究和试验数据分析确定上述阈值。当相干/非相干累加比大于阈值时，GNSS反射信号的载波相位连续可估计，利用载波相位反演水位，而小于阈值为模糊场景，载波相位不可估计，利用伪码反演水位。通过伪码和载波协同保证水位测量的连续性和高精度。
另外，模糊场景下由于非相干成分的影响使载波相位不可估计，因此只能采用伪码进行水位反演，如图3(a)所示，图3(a)为利用模糊场景进行水位反演的示意图。
GNSS反射信号的相干成分和非相干信号的时延波形分别服从Λ(τ)和∫Λ(τ-τ(ρ))dρ，τ是相干信号下，直反信号时延差，ρ是反射非相干信号的微分向量，τ(ρ)是在ρ下的非相干信号下的时延偏差，其中，Λ(τ)为直/反信号的自相关函数。非相干成分使得GNSS反射信号的时延波形并非和直达信号的完全一样，通过跟踪时延波形峰值点将得到有偏差的反射相对于直达信号时延，使得反演的水位存在一定偏差。针对海面高度测量，已通过跟踪其他时延波形特征值(例如：时延波形前沿分位点)进行反射相对于直达信号时延的估计。但是该类特征点均是基于非相干散射场景提出的，并不适合相干和非相干散射混合模糊场景。基于此，通过如图3所示的首先相干累加的方法提取GNSS反射信号的相干成分，然后非相干累加得到GNSS反射信号时延波形。由于射频通道存在通道时延，因此对射频通道交替转换周期的前半周期和后半周期输出的时延波形分别进行非相干累加。
如图3(b)所示，图3(b)为利用相干场景进行水位反演的示意图，由于时延波形采样率较低，无法满足高度反演精度的要求，因此采用三次样条法对上述得到的时延波形进行插值得到光滑的高采样时延波形。通过跟踪光滑时延波形得到GNSS反射信号相对于直射信号的时延然后通过双基几何关系得到接收机相对于河流表面的高度hrec：
由于直达和反射信号通道时延不一致，导致了通过时延波形估计的反射信号相当于直达信号的时延存在一个固定的通道偏差。为了消除该偏差，本发明采用通道交叉校正。假设真实的反射相对于直达信号的时延为τspec，射频通道1相对于通道2的偏差为Δτ。当通道交替转换前半周期，反射和直达信号分别通过通道1和2进行放大和下变频，估计时延为其中，Δτ为固定通道偏差，而在后半周期，反射和直达信号分别通过通道2和1进行放大和下变频，估计时延为/>利用前半周期和后半周期时延波形反演的高度分别为：
从上式可以得到通道不一致校正的方法：
其中，hrec1为利用前半周期时延波形反演的高度，hrec2为利用后半周期时延波形反演的高度。
载波相位采用四象鉴相器进行估计，假设反射信号功率峰值对应的复数相关值为Yrmax(i)＝Irmax(i)+jQrmax(i)，则四象鉴相器输出的相位时序为：
其中，Irmax(i)为反射信号I路最大值，Qrmax(i)为反射信号Q路最大值。
GNSS反射信号的载波相位由两部分组成，一部分为河流流动导致的多普勒频率产生的相位，一部分是反射信号相对于直达信号时延产生的相位。在利用载波相位进行水位反演前，必须对相位进行补偿消除多普勒频移产生的相位。首先对估计相位进行相位解缠得到连续的线性变化的相位：
其中，nT为观测期间估计相位经历的周期数。通过估计的反射信号多普勒频率对解缠后的相位进行补偿，即：
其中，Tint为相干积分时间，为估计的信号载波相位，/>为估计的解缠前的载波相位。
通过长期观测避免载波相位测距时的整周模糊度问题。载波相位和观测高度、卫星高度角存在关系：
其中，λ为波长。
GNSS卫星运动过程中卫星高度角是时间的函数。以卫星高度角为自变量对相位求导：
其中，为估计的载波相位频率。
由反射相对于直达信号时延产生的相位的变化周期和高度成正比例关系。因此，本发明通过补偿后的相位进行谱估计得到接收机相对于河流表面的高度。由于卫星高度角呈非周期性采样特征，因此本发明拟采用非周期采样的LS(Lomb-Scargle)法对时序进行谱估计。相比于利用短时序复数相关值估计河流流速，由于GNSS卫星动态性导致的大范围高度角变化需长时间，因此该处为长序列谱估计。
最终，假设接收机相对大地参考坐标的高程为Hrec，则河流水位为：
Hr＝Hrec-hrec
GNSS反射信号(Global Navigation Satellite System-Reflection)通过协同接收、处理直达和经地表反射的导航卫星信号实现地表物理参数的反演和监测，具备以下优势：1)利用现有的导航卫星信号，是一种双基机会源探测手段，成本低、功耗低、重量轻；2)L波段在降雨条件下信号衰减小，适合于降雨等恶劣条件的监测；3)陆基探测可利用成熟的多通道导航射频、相关器芯片完成信号处理实现快速处理和监测；4)协同多信号特征量观测实现多参数一体化探测(信号时延测水位；信号功率测海冰)。本发明瞄准河流监测，另辟蹊径通过信号处理、信息提取和河流监测应用三个层面的研究，旨在为河流监测提供一种快速、低成本、精准的水位、流速探测手段。
针对河流流速反演，本技术使用短时序谱估计反演河流流速的方法，利用相干/模糊场景判决的伪码/载波协同的河流水位反演方法，实现高精度水位反演。通过陆基接收机接收处理直、反信号，实现河流流速与水位高度反演方法，为河流流速与水位探测提出一种新型的低成本、一体化的探测手段。
在模糊场景下通过相干成分时延波形提取消除非相干成分，提高时延波形跟踪精度及水位反演精度；在相干场景下通过短、长时序两步谱估计法消除流速和水位在载波相位调制上的耦合，利用载波相位高精度反演水位。
本发明的优点在于：
一、本发明提出了一种陆基GNSS反射计的河流流速反演方法，将GNSS反射信号技术应用在陆河边界与河流流速反演上，信号源丰富，成本低。
二、本发明利用了短时序谱估计，通过峰值求解，实现河流流速反演。
三、本发明通过对直达和反射信号的射频通道在一定周期内进行交替转换，校正射频链路增益。
四、本发明建立粗糙度因子与相干/非相干累加比之间的粗糙度矫正，并采用阈值判决得到陆表指示函数，实现陆河边界提取与河流宽度获取。
本说明书中各个实施例采用递进的方式描述，每个实施例重点说明的都是与其他实施例的不同之处，各个实施例之间相同相似部分互相参见即可。
本文中应用了具体个例对本发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的方法及其核心思想；同时，对于本领域的一般技术人员，依据本发明的思想，在具体实施方式及应用范围上均会有改变之处。综上所述，本说明书内容不应理解为对本发明的限制。
Claims (8)
1.一种陆基GNSS反射计的河流流速与水位反演方法，其特征在于，包括：
获取陆基场景中GNSS信号在河流表面的反射信号；
根据所述反射信号利用高斯拟合法估计反射信号多普勒频率；所述反射信号多普勒频谱为估计的河流流速引起的GNSS反射信号多普勒频率；
根据所述反射信号多普勒频率利用反射信号双基几何关系反演河流流速；
获取所述反射信号的相干信号和非相干信号；
根据所述相干信号和所述非相干信号确定累加比；
CIR为相干/非相干累加比，Yi(τ)为反射信号，其中分子为反射信号在一段时间内相干累加信号，分母为反射信号在一段时间内非相干累加信号；
判断所述累加比是否小于设定阈值，得到第一判断结果；
当所述第一判断结果表示所述累加比小于设定阈值时，则利用模糊场景进行水位反演，确定接收机相对于河流表面的高度；
当所述第一判断结果表示所述累加比大于或者等于设定阈值时，则利用相干场景进行水位反演，确定接收机相对于河流表面的高度；
根据所述接收机相对于河流表面的高度和接收机相对于大地的高程确定河流水位。
2.根据权利要求1所述的陆基GNSS反射计的河流流速与水位反演方法，其特征在于，所述获取陆基场景中GNSS信号在河流表面的反射信号之后，还包括：
获取接收机中接收的反射功率和直射功率；
根据所述反射功率和所述直射功率确定反射率；
判断所述反射率是否大于反射率阈值，得到第二判断结果；
若所述第二判断结果表示所述反射率大于反射率阈值，则根据所述反射信号利用高斯拟合法估计反射信号多普勒频率；
若所述第二判断结果表示所述反射率小于或者等于反射率阈值，则返回步骤“获取接收机中接收的反射功率和直射功率”。
3.根据权利要求1所述的陆基GNSS反射计的河流流速与水位反演方法，其特征在于，所述根据所述反射信号利用高斯拟合法估计反射信号多普勒频率，具体包括：
根据所述反射信号利用Welch谱估计法进行估计，得到反射信号相对于直射信号多普勒频率；
根据所述反射信号相对于直射信号多普勒频率利用高斯拟合法估计反射信号多普勒频率。
4.根据权利要求3所述的陆基GNSS反射计的河流流速与水位反演方法，其特征在于，所述根据所述反射信号相对于直射信号多普勒频率利用高斯拟合法估计反射信号多普勒频率，具体包括：
根据如下公式估计反射信号多普勒频率：
其中，为估计的河流流速引起的GNSS反射信号多普勒频率，ψ(f)为反射信号相对于直射信号多普勒频率，/>为估计高斯分布标准差，f为多普勒频率，μ为多普勒偏移频谱在高斯分布下的均值，σ为多普勒偏移频谱在高斯分布下的标准差。
5.根据权利要求1所述的陆基GNSS反射计的河流流速与水位反演方法，其特征在于，所述河流流速的反演公式为：
其中，为河流流速，/>为估计的河流流速引起的GNSS反射信号多普勒频率，Φ为河流流向和GNSS卫星方位向之间的夹角，θ为GNSS卫星高度角，fc为GNSS卫星信号载波频率，c为光速。
6.根据权利要求1所述的陆基GNSS反射计的河流流速与水位反演方法，其特征在于，所述当所述第一判断结果表示所述累加比小于设定阈值时，则利用模糊场景进行水位反演，确定接收机相对于河流表面的高度，具体包括：
采用非相干叠加法提取所述反射信号的相干成分；
根据所述相干成分进行非相干叠加得到所述反射信号的时延波形；
跟踪所述反射信号的时延波形得到反射信号相对于直射信号的时延；
根据所述反射信号相对于直射信号的时延利用双基几何关系确定接收机相对于河流表面的高度；
利用通道交叉法对所述接收机相对于河流表面的高度进行校正，确定最终的接收机相对于河流表面的高度。
7.根据权利要求1所述的陆基GNSS反射计的河流流速与水位反演方法，其特征在于，所述当所述第一判断结果表示所述累加比大于或者等于设定阈值时，则利用相干场景进行水位反演，确定接收机相对于河流表面的高度，具体包括：
对所述反射信号进行载波相位估计，得到相位时序；
对所述相位时序进行相位补偿，得到补偿后的相位时序；
对补偿后的相位时序利用非周期采样的长序列谱估计法进行谱估计，确定接收机相对于河流表面的高度。
8.根据权利要求1所述的陆基GNSS反射计的河流流速与水位反演方法，其特征在于，所述河流水位的计算公式如下：
Hr＝Hrec-hrec
其中，Hr为河流水位，Hrec为接收机相对大地参考坐标的高程，hrec为接收机相对于河流表面的高度。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110973188.8A CN113671552B (zh) | 2021-08-24 | 2021-08-24 | 一种陆基gnss反射计的河流流速与水位反演方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110973188.8A CN113671552B (zh) | 2021-08-24 | 2021-08-24 | 一种陆基gnss反射计的河流流速与水位反演方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113671552A CN113671552A (zh) | 2021-11-19 |
| CN113671552B true CN113671552B (zh) | 2024-02-06 |
Family
ID=78545467
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110973188.8A Active CN113671552B (zh) | 2021-08-24 | 2021-08-24 | 一种陆基gnss反射计的河流流速与水位反演方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113671552B (zh) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115290160A (zh) * | 2022-08-03 | 2022-11-04 | 哈尔滨工程大学 | 一种基于北斗水位反演的无人机动态水位监测系统及方法 |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5585799A (en) * | 1995-06-08 | 1996-12-17 | The United States Of America As Represented By The Secretary Of The Army | Microwave doppler radar system for detection and kinematic measurements of river ice |
| KR20120060415A (ko) * | 2010-12-02 | 2012-06-12 | 아이에스테크놀로지 주식회사 | 신호의 믹싱기법을 이용한 코히어런트 도플러 유속분포 측정 방법 및 장치 |
| CN104677463A (zh) * | 2015-02-11 | 2015-06-03 | 北京航空航天大学 | 一种基于统计量的gnss-r多源海面测高装置 |
| CN104865581A (zh) * | 2015-03-19 | 2015-08-26 | 中国科学院空间科学与应用研究中心 | 一种基于gnss-r技术的河流流速测量方法与系统 |
| WO2017067033A1 (zh) * | 2015-10-19 | 2017-04-27 | 叶雷 | 一种gnss-r地表探测装置和方法 |
| CN111580080A (zh) * | 2020-05-29 | 2020-08-25 | 大连海事大学 | 一种基于船载gnss-r的有效波高反演方法及系统 |
Family Cites Families (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6549165B2 (en) * | 2001-01-19 | 2003-04-15 | Agence Spatiale Europeenne | Ocean altimetry interferometric method and device using GNSS signals |
| US11726200B2 (en) * | 2020-01-21 | 2023-08-15 | The Regents Of The University Of Colorado, A Body Corporate | Adaptive hybrid tracking algorithms for radio signal parameters estimations |
| AU2020103449A4 (en) * | 2020-11-16 | 2021-01-28 | China University Of Mining And Technology | Method for monitoring the water level of reservoir by using GNSS triple-frequency phase combination data |
| CN113049062B (zh) * | 2021-03-12 | 2022-04-15 | 北京航空航天大学 | 一种gnss直反射信号载波干涉测量湖泊水位的装置 |
-
2021
- 2021-08-24 CN CN202110973188.8A patent/CN113671552B/zh active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5585799A (en) * | 1995-06-08 | 1996-12-17 | The United States Of America As Represented By The Secretary Of The Army | Microwave doppler radar system for detection and kinematic measurements of river ice |
| KR20120060415A (ko) * | 2010-12-02 | 2012-06-12 | 아이에스테크놀로지 주식회사 | 신호의 믹싱기법을 이용한 코히어런트 도플러 유속분포 측정 방법 및 장치 |
| CN104677463A (zh) * | 2015-02-11 | 2015-06-03 | 北京航空航天大学 | 一种基于统计量的gnss-r多源海面测高装置 |
| CN104865581A (zh) * | 2015-03-19 | 2015-08-26 | 中国科学院空间科学与应用研究中心 | 一种基于gnss-r技术的河流流速测量方法与系统 |
| WO2017067033A1 (zh) * | 2015-10-19 | 2017-04-27 | 叶雷 | 一种gnss-r地表探测装置和方法 |
| CN111580080A (zh) * | 2020-05-29 | 2020-08-25 | 大连海事大学 | 一种基于船载gnss-r的有效波高反演方法及系统 |
Non-Patent Citations (3)
| Title |
|---|
| A Random Model and Simulation for Generating GNSS Ocean Reflected Signals;Bowen Li et al;《IEEE GEOSCIENCE AND REMOTE SENSING LETTERS》;第16卷(第7期);第1-5页 * |
| GNSS 多径信号模型及测高方法;李惟;《北京航空航天大学学报》;第第44卷卷(第第6期期);第1-7页 * |
| 基于GNSS-R的海面风场监测系统研究;路勇 等;《武汉大学学报(信息科学版)》;第34卷(第04期);第470-473+503页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113671552A (zh) | 2021-11-19 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Lehner et al. | Wind and wave measurements using complex ERS-2 SAR wave mode data | |
| US6525688B2 (en) | Location-determination method and apparatus | |
| CA2193814C (en) | Apparatus and method for mitigating range-doppler ambiguities in pulse-doppler radars | |
| US6191730B1 (en) | Two-channel fast-sequencing high-dynamics GPS navigation receiver | |
| EP2067054B1 (en) | Mobile-unit positioning device | |
| US20020145557A1 (en) | Location-determination method and apparatus | |
| CN109507668B (zh) | 一种基于导航卫星信号的双基成像方法 | |
| MXPA01007775A (es) | Un navegador de ins/gps adaptable profundamente integrado con rastreo de codigo de intervalo extendido. | |
| CN105277955A (zh) | 一种用于gps卫星导航的基带信号处理方法 | |
| Yu et al. | GNSS-based model-free sea surface height estimation in unknown sea state scenarios | |
| Hobiger et al. | GLONASS-R: GNSS reflectometry with a frequency division multiple access-based satellite navigation system | |
| CN113671552B (zh) | 一种陆基gnss反射计的河流流速与水位反演方法 | |
| CN102141611B (zh) | 斜视合成孔径雷达多普勒模糊数快速测定方法 | |
| KR20170004900A (ko) | 서브캐리어 트래킹 모호성을 해결하기 위해 비모호 판별기를 계산하는 gnss 수신기 | |
| CN113419232B (zh) | 一种陆基gnss反射计的河流边界与宽度探测方法 | |
| CN116699658A (zh) | 一种基于双频反射信号组合的gnss-r载波相位海面测高的方法 | |
| US3195136A (en) | Time and distance of closest approach detector of moving object | |
| Zuffada et al. | Coherence time and statistical properties of the GPS signal scattered off the ocean surface and their impact on the accuracy of remote sensing of sea surface topography and winds | |
| Egido et al. | Soil moisture monitorization using GNSS reflected signals | |
| Wang | A landslide monitoring technique based on dual-receiver and phase difference measurements | |
| Zhang et al. | Multi-angle fusion of SS-BiSAR images using Compass-2/Beidou-2 Satellites as opportunity illuminators | |
| Liu et al. | Doppler ambiguity resolving in compressed azimuth time and range frequency domain | |
| Hobiger et al. | A real-time GNSS-R system based on software-defined radio and graphics processing units | |
| Azmani et al. | Soil moisture estimation using land-reflected gps l2c bi-static radar measurements | |
| Jales | GNSS-Reflectometry: Techniques for scatterometric remote sensing |
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 |