本发明涉及一种监测水库水位的方法,尤其适用于防洪抗旱中使用的利用水库大坝变形监测系统的gnss数据进行水库水位监测的方法。
背景技术:
水是人类赖以生存的宝贵自然资源,同时也是灾害源。合理地存贮、利用、监测水资源变化是经济、生态可持续发展的重要保证。水库作为重要的水利基础设施,在水运交通、水力发电、防洪抗旱、生态环境等方面发挥着至关重要的作用,其安全与否直接影响到下游人民的生命财产安全。目前,传统的水库水位监测一般是采用人工水尺测量,或是超声波测量。人工方法存在一定的局限性:不能实现连续观测,观测时空分辨率不高,特别是在暴雨汛期,观测受限且存在安全隐患。超声波及其他手段通常需要建立一整套监测系统,监测成本较高。国内许多库区和大坝都建立了gnss形变监测系统,水库大坝变形监测gnss观测系统的测站通常设置在大坝外观或临近水域的高边坡,这为我们利用gnss-interferometricreflectometry(gnss-ir)测量水库水面高度的变化提供了可能。已有的gnss监测系统在进行大坝形变监测的同时还可以成为一种无需投入额外监测成本的水库水位监测手段。
技术实现要素:
针对上述技术的不足之处,提供一种步骤简单,检测效果好,能够解决目前现有水库水位监测方式存在的问题和不足,经济性好的利用gnss三频相位组合数据监测水库水位的方法。
为实现上述技术目的,本发明的一种利用gnss三频相位组合数据监测水库水位的方法,在水坝的高边坡位置设置有接收gnss信号的接收机及天线,其于步骤为:通过架设的水库大坝高坡临水面的gnss变形监测接收机数据,利用三频相位组合获取干涉信号,提取三个载波信号,通过消除几何距离和电离层延迟的组合方式,获取干涉时序信号,利用小波分析进行谱分析获取干涉信号的主频率,将时序信号与卫星高度角一一对应,形成以卫星高度角为变量的时序信号,并对该非等间隔信号进行频谱分析,获取主频率,建立主频率和已知水库水位的模型,并通过模型获取实时水库水位。
具体的:
a首先从卫星观测值中提取三频相位观测数据及卫星高度角;
b利用除几何距离因素和电离层影响的方式将三频相位观测数据组成三频相位组合观测值,具体利用公式:
c将以时间为变量的三频相位组合观测值m1,2,5(t)序列替换为时间对应的卫星高度角的正弦为变量的序列m1,2,5(sinθ);
d采用小波变换对三频相位组合观测值m1,2,5(sinθ)进行频谱分析,,通过频谱分析获取三频相位组合观测值m1,2,5(sinθ)的主频率;
e利用三频相位组合观测值m1,2,5(sinθ)的主频率与已知的水库水位高度建立线性模型,利用此模型获取接收机天线相位中心到水库水面的垂直高度;
f利用垂直高度与标准水面高度差值获取水库水位高度。
将水面看作是镜面反射,反射信号来自于镜像点sp,反射信号与直射信号的叠加信号,即干涉信号相对于直射信号的路径延迟为与水面距离接收机天线相位中心的垂直高度h相关;
gnss接收机接收到的信号表示为直射信号与反射信号的叠加:
其中,
利用三频相位组合观测值获取干涉信号β,三频相位组合观测值m1,2,5(t)的组合用以消除信号传播的几何距离以及电离层的影响,令
得到三频相位组合观测值表达式:
式中,
导航卫星信号经过电离层,相位延迟为:
m1,2,5r只与卫星高度角、信号波长和天线高度相关;
λiψi(t)=di+δi+δ,i=1,2,5
di为卫星与接收机天线相位中心之间的几何距离,δi为电离层延迟,δ为钟差和噪声;
由于相位组合的原则为无几何距离和消除电离层,三频相位组合为:
m1,2,5(t)=m1,2,5r=κ1β1(t)+κ2β2(t)+κ5β5(t)
线性相位组合的有效波长为:
根据频率与水面与接收机天线相位中心之间的垂直高度之间的对应关系,建立线性模型:
h=kf+b
h为水面与接收机天线相位中心之间的垂直高度,f为时间序列的主频,k、b为模型系数,利用实际观测值进行k、b为模型系数的确定。
有益效果:
本技术利用已在水库大坝边坡上建设且运行的变形监测系统的观测数据,能够在不增加额外设施和建设成本的基础上,实现对于水库水位的自动监测。利用三频信号高精度的特性,使得利用三频信号进行水库水位反演精度更高,可用数据更多;由于数据及设备是已有的,且大坝上均有安装相关的设备,因此无需额外硬件搭建费用,本方法的具有较高的可用性及经济性;无需人工监测且数据量大,能够实现近实时的水库水位监测。
附图说明
图1是本发明利用大坝变形监测gnss三频相位组合数据进行水库水位监测方法示意图;图2是本发明使用的干涉相位模型。
具体实施方式
下面结合附图对本发明的实施例做进一步说明。
传统大地接收机天线所接收的信号是由卫星直射信号与水面或地面反射的信号叠加产生的干涉信号,图1为信号干涉的几何图形,因此将水面看作是镜面反射,首先从卫星观测值中提取三频相位观测数据及卫星高度角;反射信号来自于镜像点sp,反射信号与直射信号的叠加信号(即干涉信号)相对于直射信号的路径延迟为与水面距离接收机天线相位中心的垂直高度h相关。
图2为干涉信号的相位关系,横坐标为信号的同相分量,纵坐标为信号的正交分量。
gnss接收机接收到的信号为为直射信号与反射信号的叠加:
其中,
βi(t)是由于来自水面反射信号干扰而产生的干涉相位。
然后利用三频相位组合观测值获取干涉信号β,具体的:三频相位组合观测值m1,2,5(t)的组合用以消除信号传播的几何距离以及电离层的影响,令
其中,
u=κ1n1+κ2n2+κ5n5
导航卫星信号经过电离层,相位延迟为:
其中,tec为传播路径的电子含量,f为卫星信号的频率。延迟大小与信号频率呈反比。三频相位组合观测值m1,2,5(t)为:m1,2,5(t)=m1,2,5d+m1,2,5r,
其中,
m1,2,5r只与卫星高度角、信号波长和天线高度相关;
λiψi(t)=di+δi+δ,i=1,2,5
di为卫星与接收机天线相位中心之间的几何距离,δi为电离层延迟,δ为钟差和噪声;
由于相位组合的原则为无几何距离和消除电离层,三频相位组合表示为:
m1,2,5(t)=m1,2,5r=κ1β1(t)+κ2β2(t)+κ5β5(t)
线性相位组合的有效波长为
然后选取适合的干涉信号β:
在低卫星高度角的情况下,反射信号较强,因此在选择干涉信号时,一般选取卫星高度角为10~30度范围内的数据。并根据观测接收机与水库的方位及镜像点的位置,确定方位角的取值范围。
然后获取干涉信号β的主频率
对获取的干涉信号转换为以卫星高度角的正弦为变量的时间序列。序列变为非等间隔的序列,利用傅里叶变换获取主频率不再适用。选择小波变换进行时间序列的主频率获取。
根据频率与天线高度(反射面与接收机天线相位中心之间的垂直高度)之间的对应关系,建立线性模型:
h=kf+b
h为天线高度,f为时间序列的主频,k、b为模型系数,利用实际观测值进行k、b为模型系数的确定。
获取水库水位变化的检测信息
t1时刻根据所得的干涉信号的时间序列频率,代入线性模型,求得天线高度h1,在t2计算天线高度h2,不同时刻的高度相减就可以获取水库水位的变化,如果得知水库水位基准高度h0,对任意时刻的获取的天线高度hi与基准相减,便可以实时获取水库水位的高度。
PHP网站源码坑梓至尊标王深圳模板推广坑梓网站seo优化光明企业网站设计坑梓关键词排名包年推广南联网站改版荷坳关键词按天扣费同乐百度竞价包年推广吉祥百搜词包大鹏标王大浪推广网站沙井建站龙岗网站优化推广永湖网站优化软件宝安网站搜索优化吉祥建网站福永网站建设设计荷坳SEO按天扣费福田企业网站建设坑梓网络推广松岗模板推广东莞网站排名优化松岗外贸网站设计西乡模板制作松岗seo网站推广大鹏seo石岩seo网站推广坪山推广网站惠州网站建设设计罗湖网站优化按天扣费歼20紧急升空逼退外机英媒称团队夜以继日筹划王妃复出草木蔓发 春山在望成都发生巨响 当地回应60岁老人炒菠菜未焯水致肾病恶化男子涉嫌走私被判11年却一天牢没坐劳斯莱斯右转逼停直行车网传落水者说“没让你救”系谣言广东通报13岁男孩性侵女童不予立案贵州小伙回应在美国卖三蹦子火了淀粉肠小王子日销售额涨超10倍有个姐真把千机伞做出来了近3万元金手镯仅含足金十克呼北高速交通事故已致14人死亡杨洋拄拐现身医院国产伟哥去年销售近13亿男子给前妻转账 现任妻子起诉要回新基金只募集到26元还是员工自购男孩疑遭霸凌 家长讨说法被踢出群充个话费竟沦为间接洗钱工具新的一天从800个哈欠开始单亲妈妈陷入热恋 14岁儿子报警#春分立蛋大挑战#中国投资客涌入日本东京买房两大学生合买彩票中奖一人不认账新加坡主帅:唯一目标击败中国队月嫂回应掌掴婴儿是在赶虫子19岁小伙救下5人后溺亡 多方发声清明节放假3天调休1天张家界的山上“长”满了韩国人?开封王婆为何火了主播靠辱骂母亲走红被批捕封号代拍被何赛飞拿着魔杖追着打阿根廷将发行1万与2万面值的纸币库克现身上海为江西彩礼“减负”的“试婚人”因自嘲式简历走红的教授更新简介殡仪馆花卉高于市场价3倍还重复用网友称在豆瓣酱里吃出老鼠头315晚会后胖东来又人满为患了网友建议重庆地铁不准乘客携带菜筐特朗普谈“凯特王妃P图照”罗斯否认插足凯特王妃婚姻青海通报栏杆断裂小学生跌落住进ICU恒大被罚41.75亿到底怎么缴湖南一县政协主席疑涉刑案被控制茶百道就改标签日期致歉王树国3次鞠躬告别西交大师生张立群任西安交通大学校长杨倩无缘巴黎奥运