1 引言

全球性海平面的上升可能会对沿海地区的社会活动产生重大的影响，因此，监测和了解海平面变化显得非常重要[1]。GNSS-R技术是众多测量技术中的一种。该技术由欧州航天局科学家Martin-Neira教授于1993年首次提出[2]，之后迅速成为国内外一个研究的热点。

随着改革的深化和科学技术的快速发展，农业档案管理工作水平都应伴随时代的脉搏，与农业发展形势的需要相适应。农业工作涉及面广、管理环节多、工作要求高，档案管理工作要在农业开发建设中提供更好的服务。

作为一种新型的海面高度测量技术，其具备众多的优点，包括测量范围全球覆盖、接收机成本低、拥有大量信号源等[2-3]。目前，在国内外的研究中，该技术实现测高目的研究有多种方法，主要有码测高[4-6]、相位测高[7-13]和基于信号的信噪比(SNR)观测值测高[14-16]等。其中，码测高和相位测高一般除了需要设置一个接收直射信号的右旋圆极化(RHCP)的天线，还需要设置一个专门接收反射信号的左旋圆极化(LHCP)的天线和一个专门研制的GNSS-R接收机。基于SNR观测值测高仅采用一个RHCP天线和一个普通的测量型接收机，同时对直射信号和反射信号进行接收。而且，该种测高方法所达到的测量精度比码测高方法所达到的测量精度要高；在海面粗糙度较大的情况下，其测高的性能比相位测高的性能要好[16]。

在过去的时间内，大部分的测高研究都是基于GPS系统所开展的，而作为我国独立自主研究的北斗卫星导航系统，由于其组网时间较晚，在GNSS-R领域的相关研究较少。但是，随着北斗系统的发展和其独特的混合星座的特点，该系统在该领域拥有重大的发展前景，例如在风速反演[17]，土壤湿度反演[18]，海面波高反演[19]和海面测高[20-21]等方面。

基于上述原因，本文依据基于SNR观测值测高的方法，首次采用北斗系统3个频段(B1、 B2、B3)信号成功地实现对潮位高度进行反演。而此研究中所分析的数据来源于本团队于2014年9月5—31日，在中国浙江省舟山市的大洋山岛上所收集的。由于实验期间，岛上环境复杂和天气多变，在反演潮位高度过程中，我们对干扰源和海面粗糙度作了一定的分析。

2 单天线GNSS-R海面高度反演原理

2.1 北斗系统简介

依据BDS-SIS-ICD-2.1所述，北斗卫星导航系统简称北斗系统，英文缩写为BDS，其空间星座由5颗地球静止轨道(GEO)卫星、27颗中圆地球轨道(MEO)卫星和3颗倾斜地球同步轨道(IGSO)卫星组成。其信号结构中，B1、B2信号由I、Q两个支路的“测距码+导航电文”正交调制在载波上构成，分别如式(1)和式(2)所示：

 

 

(1)

 

 

(2)

式中，上角标j表示卫星编号；AB1I表示B1I信号振幅；AB2I表示B2I信号振幅；AB1Q表示B1Q信号振幅；AB2Q表示B2Q信号振幅；CB1I表示B1I信号测距码；CB2I表示B2I信号测距码；CB1Q表示B1Q信号测距码；CB2Q表示B2Q信号测距码；DB1I表示调制在B1I测距码上的数据码；DB2I表示调制在 B2I 测距码上的数据码；DB1Q表示调制在B1Q测距码上的数据码；DB2Q表示调制在B2Q测距码上的数据码；f1表示B1信号载波频率；f2表示B2信号载波频率；φB1I表示B1I信号载波初相；φB2I表示B2I信号载波初相；φB1Q表示B1Q信号载波初相；φB2Q表示B2Q信号载波初相。

关于B3信号，BDS-SIS-ICD-2.1并没有给出相关详细的介绍。

2.2 基本反演原理

由于天线放置的位置周围存在较为丰富的植被和人造建筑，如图6所示。因此，一方面为了保证被跟踪的卫星反射信号落在海面上，其方位角需限制在180°～260°，其仰角需限制为8°～25°。另一方面，周围的环境(植被和人造建筑)将对最终的反演结果造成影响。所以在数据分析上需要使用本文在数据处理过程中提到的判断原则，避免结果出现较大偏差。基于上述原因，我们所选择跟踪的卫星及其相应的频段、类型等参数如表3所示。需要说明的是，本次实验并没有选择GEO卫星进行反演，这是因为GEO卫星的信号在岸基观测中，其仰角不变或者变化非常小。

板涧河调蓄水库主要建筑物，包括大坝、溢洪道、导流泄洪洞、连通洞及补水泵站。大坝为2级建筑物，溢洪道和泄洪洞为3级建筑物，补水泵站为2级建筑物。洪水标准为50年一遇设计，1 000年一遇校核。

  

图1 单天线GNSS-R反演模型Fig.1 The GNSS-R altimetry with a single antenna

 

(3)

[8] Lowe S T， Zuffada C， Chao Yi， et al. 5-cm-Precision aircraft ocean altimetry using GPS reflections[J]. Geophysical Research Letters， 2002， 29(10): 13-1-13-4.

φ=2πλ-1ρ，

(4)

式中，λ表示GNSS信号的载波波长；ρ表示反射信号相对于直射信号的路径延迟。根据图1所示，该路径延迟的表达如式(5)所示：

ρ=2hsin θ，

(5)

式中，h表示天线和反射面之间的垂直高度(简称反射高度)；θ表示卫星仰角。

结合式(3)至式(5)，再对原始的SNR数据SNRraw(单位：dB)进行线性尺度变换及提取出变换后的数据SNRlinear中的多径振荡分量δSNRlinear(单位：V)后，将可以构造出该多径振荡分量δSNRlinear与h的模型[14]，如式(6)所示：

δSNRlinear=Acos (4πhλ-1sin θ)，

(6)

式中，(单位：V)。由式(6)可计算出反射面的高度h。

2.3 数据处理流程

在获得所选择的卫星的SNRraw后，首先对其进行相应的线性尺度变换从而获得SNRlinear。为了去除SNRlinear中的直射信号分量，需要对SNRlinear进行相应的低阶多项式消除，从而得到δSNRlinear。该低阶多项式的阶数选择视实际环境而定，在本实验中该阶数设为2阶则可以实现去除直射信号分量的目的。最后，在忽略式(6)中A值的微小变化及通过对δSNRlinear进行L-S周期图法(Lomb Scargle Periodogram，LSP)[23]分析后，就可以计算出式(6)中的反射面高度h。需要说明的是，此处采用LSP分析而不是传统的傅里叶变换，是因为LSP不需要要求数据是均匀采样的[14]。

出于实际环境的考虑，即直射天线周围可能存在植被和人工建筑，会一定程度地干扰最终的反演结果。特别地，在受风速影响下的海面粗糙度过大时，来自海面的反射信号的相干分量的能量将大大降低，而来自植被和人工建筑的干扰并不受风速影响，甚至在暴雨天气中反而有所增强，从而导致最终的反演结果偏差过大。为了避免出现这种情况，我们在从LSP分析中提取出天线到反射面的垂直高度h前，需要设定一个判决原则。而该判决原则需依据实验中干扰源的分布和特性所设定的。

首先，在所得到的LSP分析中提取出最高峰值及该峰值所对应的反射高度。以图2为例，我们分别定义为Rate(b)和b。然后，通过式(7)进行计算。

  

图2 LSP分析例子Fig.2 An example of LSP analysis

 

门限值

(7)

当通过式(7)计算后所得到的值高于一个设定的门限值时，我们认定反射高度的结果满足判决原则，反之则不满足。一般而言，该门限值过低时，将造成漏警率过大；过高时，将造成虚警率过大；而基于我们大量数据分析所得，该门限值定为6时最佳。

数据处理整体的流程图如图3所示。

总而言之，路吉阿诺斯说得上是古希腊文明的一个“流氓鬼”，他那嘲讽的棘刺并未因了一千八百年之久远而钝挫，有着不以古而遂湮灭的生命伟力。他以洞悉一切的眼力、识见对所处时代的一切“虚妄”作出淋漓尽致的指斥与绝不妥协的批判，最大限度地驱除人们思想上的痼疾，重现真实。因此，路吉阿诺斯是他那个时代少有的唯理者，清醒者与批判者，被周作人誉称为古希腊的“殿军人物”。

  

图3 数据处理流程图Fig.3 The main procedure of data processing using SNR method

3 实验场景和数据

本次实验于2014年9月5—31日，在中国浙江省舟山市的大洋山岛上开展，其相关的实验场景如图4所示。实验中，天线位置为30.578 657 434°N，122.064 522 466°E，高29.171 1 m，天线型号为华信HX-BS781A型RHCP天线，相关的物理特性如表1所示；接收机型号为和芯星通UR370型接收机，如图5所示，其相关的物理特性如表2所示。

我国不属于判例法国家，这就决定了我国不能像美国司法实践那样以判例形式确立作者享有的立法未否定的利益。然而，我国著作权法又不存在与法国或意大利类似的赋予作者对作品使用行为的全面控制权条款。因此，只有借助于兜底条款的适用来解决这一困境。一方面，兜底条款的设置目的在于弥补作品使用行为列举的不足；另一方面，兜底条款的存在也给法官提供了可适用的法律依据。比如对于那种与有名权利控制的行为本质上具有同质性的被诉行为，且对该行为的控制对作者而言是一项重要的财产利益应给予法律的保护，在这种情形下，法院可适用兜底条款来对此类权益进行保护。

  

图4 潮位反演实验场景Fig.4 Construction of the Dayangshan port experiment

 

表1 华信HX-BS781A天线物理特性

 

Tab.1 The parameters of the antenna HX-BS781A

  

天线尺寸接头形式工作温度存储温度湿度173.4mm×62.6mmTNC阴头-40～85℃-55～85℃95%不冷凝

  

图5 和芯星通UR370型接收机Fig.5 The GNSS receiver of type unicoremum UR370

 

表2 和芯星通UR370型接收机物理特性

 

Tab.2 The parameters of the receiver UR370

  

尺寸工作温度存储温度湿度100mm×60mm×11.4mm-40～85℃-55～95℃95%不冷凝

有别于传统GNSS-R测高反演的装置，单天线GNSS-R反演过程仅需采用一个水平放置的右旋圆极化(RHCP)的天线，同时对直射信号和来自海面的反射信号进行接收，如图1所示。RHCP天线会在最大程度上接收来自直射信号的能量的同时抑制反射信号。但是，来自反射信号的能量并没有完全被抑制[14]。基于此，RHCP天线所采集的信号的信噪比(Signal-to-Noise Ratio， SNR)数据与直射信号、反射信号的能量关系可如式(3)[22]所示：

  

图6 天线周围环境Fig.6 The surrounding environment for the experimental location

 

表3 跟踪卫星参数

 

Tab.3 Three signals analyzed in the experiment

  

频段载波波长PRN号卫星类型B10.1920cm9、10IGSO12、14MEOB20.2483cm9、10IGSO12、14MEOB30.2363cm9、10IGSO12、14MEO

考虑到北斗系统是混合星座，需对该系统3种类型的星座就本文所研究的SNR测高作讨论。首先是本次实验没有使用的GEO类型的星座，该类型的星座在岸基实验中，其观测的仰角是不变的或是变化极小，因此并不适合SNR岸基测高。但是这并不能否定其在GNSS-R测高领域的研究价值，在我们早期的研究中，GEO在其他测高方法上，例如码测高上表现出的性能比IGSO和MEO优越[20]。而对于IGSO类型的卫星和MEO类型的卫星，虽然两者卫星类型存在较大的差异，例如IGSO的轨道高度(35 786 km)高于MEO的轨道高度(21 528 km)，但通过观察图9到图13，我们发现在两种类型的卫星在各自频段上的SNR大致相同上，同时两者的反演精度也大致相同。所以我们有理由认为在使用北斗系统时，可以同时采用IGSO和MEO进行SNR测高。

4 反演结果分析

依据所选择的卫星，我们绘制了这些卫星于实验时间内，在海面的镜面反射点路径，如图8所示。作为后续分析需要，我们分别选择这4颗卫星在3个频段上具有代表性的数据处理流程的例子作展示(图9至图12)。

  

图7 潮位对比数据和风速对比数据位置示意图(UTC 2014/9/5-2014/9/31)Fig.7 The locations of tide data observation and wind speed data observation(UTC 2014/9/5-2014/9/31)

  

图8 北斗9(红线)、10(蓝线)、12(绿线)、14(粉线)号卫星镜面反射点路径(UTC 2014/9/5-2014/9/31)Fig.8 Specular point motion paths of GPS 9(red line)，10(bule line)，12(green line)，14(pink line)

  

图9 北斗系统9号卫星3个频段信号分析过程(自上而下各行依次表示北斗系统3个频段的原始SNR，经线性尺度化后的SNR，多径振荡分量和LSP分析)Fig.9 The main procedure of BDS PRN9’s SNR observation on three frequencies (for each of BDS three frequencies， raw SNR， linear-scale SNR， detrended SNR and the LSP analysis on detrended SNR are showed from top to bottom of the figure， respectively)

  

图10 北斗系统10号卫星3个频段信号分析过程(自上而下各行依次表示北斗系统3个频段的原始SNR，经线性尺度化后的SNR，多径振荡分量和LSP分析)Fig.10 The main procedure of BDS PRN10’s SNR observation on three frequencies (for each of BDS three frequencies， raw SNR， linear-scale SNR， detrended SNR and LSP analysis on detrended SNR are showed from top to bottom of the figure， respectively)

  

图11 北斗系统12号卫星3个频段信号分析过程(自上而下各行依次表示北斗系统3个频段的原始SNR，经线性尺度化后的SNR，多径振荡分量和LSP分析)Fig.11 The main procedure of BDS PRN12’s SNR observation on three frequencies (for each of BDS three frequencies， raw SNR， linear-scale SNR， detrended SNR and LSP analysis on detrended SNR are showed from top to bottom of the figure， respectively)

  

图12 北斗系统14号卫星3个频段信号分析过程(自上而下各行依次表示北斗系统3个频段的原始SNR，经线性尺度化后的SNR，多径振荡分量和LSP分析)Fig.12 The main procedure of BDS PRN14’s SNR observation on three frequencies (for each of BDS three frequencies， raw SNR， linear-scale SNR， detrended SNR and LSP analysis on detrended SNR are showed from top to bottom of the figure， respectively)

由于LSP分析所提取的是天线到海面的高度，在与实际潮位信息进行对比分析前，需要进行转换[14，24]，如式(8)所示：

htide(t)=[h(t0)+Htide(t0)]-h(t)，

(8)

式中，htide(t)为t时刻时反演的潮位高度；h(t0)为LSP在某一参考时刻t0所提取的天线到海面的高度；Htide(t0)为该时刻t0实际潮位高度；h(t)为LSP在t时刻时所提取的天线到海面的高度。在本文中，我们选用了2014年9月5日的所有反演结果取均值代入到h(t0)，以及选择与反演结果所对应时刻的实际潮位高度取均值代入到Htide(t0)。图13显示的是这4颗卫星在3个频段上的潮位反演结果。图14、图15和表4对该结果做了相应的定量分析。

1.CMV视网膜脉络膜炎的诊断和治疗：典型的症状包括飞蚊症、漂浮物、盲点或外周视野缺损，患者常表现为快速视力下降，眼底检查表现为“番茄炒鸡蛋样”改变，沿血管分布的浓厚的黄白色视网膜损伤，伴或不伴视网膜内出血；确诊有赖于眼底镜检查。

 

表4 北斗系统3个频段定量分析结果

 

Tab.4 The precise of BDS-R altimetry relative to tides

  

频段标准差相关系数B133.82cm95.98%B238.48cm94.68%B334.56cm96.08%

如图9至图13所示，北斗卫星的B1频段的SNR最小；B3频段的SNR略高于B1；作为SNR最高的B2频段，其SNR高于B1频段约为5 dB。而依据表4，B1频段的反演精度略高于B3频段，B3频段反演结果的相关系数略高于B1频段；B1、B3频段的反演精度和相关系数明显高于B2频段，这是因为B2频段仍处于调试中(BDS-SIS-ICD-2.0)，其性能可能不如B1和B3频段。尽管各频段的反演结果得精度有所差异，但总体上看，北斗系统的各个频段都能成功地反演潮位高度。

由于实际环境的限制，与反演结果作对比分析的潮位数据和风速数据选择在与实验场地相临近的区域，如图7所示。其中潮位对比数据区域(上海芦潮港)位于实验场地西北方向约30 km处(http://www.chinaports.com/)，而风速对比数据区域(舟山嵊泗)位于实验场地东北方向约35 km处(http://lishi.tianqi.com/)。潮位、风速数据均直接通过网上所显示的记录获得。在进行反演结果与潮位数据作对比时，我们对潮位数据进行了三阶样条插值。

  

图13 北斗系统3个频段潮位结果Fig.13 The altimetry results of BDS SNR observation on three frequencies

  

图14 北斗系统3个频段潮位偏差Fig.14 The residual between the altimetry and tide of BDS SNR observation on three frequencies

由于北斗系统仍处于发展阶段，在本实验开展时，北斗系统MEO卫星数量相对较少，且因实验场地的限制，只能在某个范围内对卫星反射信号进行分析，故实验中各个反演结果的时间间隔较大(1天中只有约3个反演结果)，未能很好地反演潮位变化趋势。随着北斗卫星的增加，并且开发BeiDou/GPS融合的方法，可以得到更加密集的测量结果，时间分辨率也会进一步增加。

为了找出造成该错误结果的原因，我们分析了相邻近区域的风速数据，如图17所示。我们发现在该段时间内，其风速比平时有很明显地增强。而风速的过大将造成海面粗糙度增大，最终使来自海面的反射信号的相干分量大大降低，这一点与图16中属于海面反射高度范围内的幅值过小相对应。同时，另一个不能忽视的因素是在该时间段内持续发生了大暴雨，这在一定程度上增强了天线周围的植被和土壤对反演结果的干扰。该干扰表现在图16的LSP分析中，相应的反射高度范围(10 m以内)的幅值与平常时(图9到图12)相比有明显的增强。

  

图15 反演高度与潮位高度的相关系数图Fig.15 Correlation coefficent between altimetry and tide

  

图16 不满足判决原则的LSP分析Fig.16 The LSP analysis of BDS during the duration of no altimetry results

  

图17 相邻区域日平均风速图Fig.17 The wind speed data of neighboring areas during the experiment

5 结论

本文依据SNR测高原理，开展了近1个月的相关的实验，成功地通过北斗系统3个频段的多个卫星信号实现潮位高度反演。其中，B1频段反演结果的STD为33.82 cm，相关系数为95.98%；B2频段反演结果的STD为38.48 cm，相关系数为94.68%；B3频段的反演结果的STD为34.56 cm，相关系数为96.08%。在我们的实验中，就反演性能而言，B1和B3较为接近，而B2较差，这是因为B2频段仍处于调试中(BDS-SIS-ICD-2.0)，其性能可能不如B1和B3频段。

总之，将法律认同作为实现“一带一路”倡议的基础性工作，这要求我国在充分研究和尊重现有国际商事规则的基础上，提出既适合中国国情，又适应世界发展潮流的现代商事规则。只有在充分考虑国际社会可接受程度的基础上，完善我国的商事法律制度并对外推广才能真正在推动法律趋同的过程中获得法律认同。因此，在“一带一路”倡议中向世界推行的商事规则，不可能等到我国商事规则完善之后才进行，我国应当双管齐下，在现代商法理念的指导下构建与完善国内商事法律制度的同时，积极探寻符合世界商事活动要求的商事规则，从而确保我国“一带一路”倡议的准确方向和最终目标的实现，提升我国在当代世界商事规则形成上的话语权。

在图13中，北斗系统3个频段在22—25日之间没有出现相应的反演结果。图16显示了该时间段内4颗卫星在3个频段上具有代表性的LSP分析。从图中我们可以看出幅值较高处所对应的反射高度在10 m以内，这对应于实际环境中的植被区域，如图6所示。而且，在该时间段内，4颗卫星在3个频段上的所有LSP分析都未能满足判决原则。这说明，我们所提出的判决原则能够有效地排除错误的反演结果。

本文展示的是一种简易且经济的潮位高度测量技术，该技术仅需单台普通的测量型天线和接收机，可以依赖已布设在我国海岸附近的GNSS天线对海面进行较高时间分辨率的潮位高度反演，降低了潮位观测成本，实现高空间分辨率的岸基潮位观测。

在接下来的工作中，我们将继续针对该方法进行分析，同时，拟在其他平台环境下开展相关实验。

致谢:感谢欧州航天局Martin Neira 教授、Fran Fabra博士，北京航空航天大学杨东凯教授、李伟强博士在本研究中提供的指导和帮助。

参考文献：

[1] Bindoff N L， Willebrand J， Artale V， et al. Observations: Oceanic Climate Change and Sea Level Climate Change 2007: The Physical Science Basis[M]. Cambridge: Cambridge University Press， 2007: 385-433.

[2] Martin-Neira M. A passive reflectometry and interferometry system (PARIS): application to ocean altimetry[J]. ESA Journal， 1993， 17(4): 331-355.

这是国内首次关于北斗3个频段信号的SNR测高研究，对相关的研究具有一定的参考价值和借鉴意义。在分析过程中，我们采用了一个简单而有效的判决原则以实现高精度自动反演过程。同时，我们引入相邻区域的风速数据作对比分析，分析了极端风速情况下对该方法的影响。

[3] Auber J C， Bibault A， Rigal J M. Characterization of multipath on land and sea at GPS frequencies[C]//Proceedings of the 7th International Technical Meeting of the Satellite Division of the Institute of Navigation. Salt Lake City， UT:Institute of Navigation， 1994: 1155-1171.

[4] Germain O， Ruffini G. A revisit to the GNSS-R code range precision[Z]. eprint arXiv: physics/0606180， 2006.

[5] Stosius R， Beyerle G， Helm A， et al. Simulation of space-borne tsunami detection using GNSS-Reflectometry applied to tsunamis in the Indian Ocean[J]. Natural Hazards and Earth System Sciences， 2010， 10(6): 1359-1372.

△求购《飘泊诗人李商隐》（第二版），胡良清、刘敬堂合著。长江文艺出版社出版。要求内页完整，封面无破损。价格可谈。联系人：良清（QQ：550623123）

[6] Carreno-Luengo H， Camps A， Ramos-Pérez I， et al. Experimental Evaluation of GNSS-Reflectometry Altimetric precision using the P(Y) and C/A signals[J]. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing， 2014， 7(5): 1493-1500.

每个人对同一事物的主动程度是不尽相同的，这与他们的兴趣、爱好直接相关。同时，也与“主动性”的醒悟程度紧密关联。具体表现是：第一种程度是不需要人告诉就能出色地完成学习任务；次之，就是别人告诉一次就能完成学习；再次之，就是经过多次提醒才去完成（惰性开始出现）；更次之，就是只有在形式所迫时才能把事情做好，这种人是在磨洋工（惰性向懒惰转化）；最次等的是即使有人追着他，告诉他怎么去做，他也不会把事情做好（完全的懒惰）。

将我院在2016年4月至2017年4月收治的100例高血压患者纳入本次实验,按照随机抽取的方式将其分为联合组(50例)与对照组(50例)。联合组中男性30例,女性20例;年龄54—79岁,平均(67.1±4.1)岁;病程1—14年,平均(6.8±2.3)年。对照组中男性28例,女性22例;年龄52—80岁,平均(66.8±4.0)岁;病程1—16年,平均(7.0±2.2)年。采取SPSS19.0统计学软件分析联合组、对照组患者以上基线资料的差异,得出P>0.05,证实本次实验可行。

[7] Martin-Neira M， Colmenarejo P， Ruffini G， et al. Altimetry precision of 1 cm over a pond using the wide-lane carrier phase of GPS reflected signals[J]. Canadian Journal of Remote Sensing， 2002， 28(3): 394-403.

式中，Pd表示直射信号的能量；Pr表示反射信号的能量；φ表示直射信号和反射信号之间的相位差。该相位差与反射信号相对于直射信号的路径延迟ρ的关系如式(4)所示：

前一阶段滴滴顺风车恶性事件闹得沸沸扬扬，同时也把滴滴打车客服人员推向风口浪尖。究竟是客服外包引起的服务质量下降，还是客服人员墨守成规不作为，进而延误了警方介入的最佳时机，留待法律界人士与公众讨论。但不可否认，这起案件从一个侧面说明了客服人员在公司中的重要地位。

[9] Ruffini G， Soulat F， Caparrini M， et al. The GNSS-r eddy experiment Ⅰ: altimetry from low altitude aircraft[Z]. eprint arXiv: physics/0310092， 2003.

[10] Helm A. Ground-based GPS altimetry with the L1 Open GPS receiver using carrier phase delay observations of reflected GPS signals[R]. Potsdam: Geoforschungszentrum Potsdam， 2008.

[11] Semmling A M， Beyerle G， Stosius R， et al. Detection of Arctic Ocean tides using interferometric GNSS-R signals[J]. Geophysical Research Letters， 2011， 38(4): L04103.

[12] Fabra F， Cardellach E， Rius A， et al. Phase altimetry with dual polarization GNSS-R over sea ice[J]. IEEE Transactions on Geoscience and Remote Sensing， 2012， 50(6): 2112-2121.

[13] Zhang Yun， Liu Fengling， Gu Qiming， et al. Study of accurate ocean-altimetry with GNSS-R[C]//Proceedings of 2014 IEEE International Geoscience and Remote Sensing Symposium. Melbourne， VIC， Australia: IEEE， 2014: 1410-1413.

[14] Larson K M， Löfgren J S， Haas R. Coastal sea level measurements using a single geodetic GPS receiver[J]. Advances in Space Research， 2013， 51(8): 1301-1310.

1061 Clinical characteristics of patients with cryptogenic stroke and patent foramen ovale

[15] Larson K M， Ray R D， Nievinski F G， et al. The accidental tide gauge: a GPS reflection case study from Kachemak Bay， Alaska[J]. IEEE Geoscience and Remote Sensing Letters， 2013， 10(5): 1200-1204.

[16] Löfgren J S， Haas R. Sea level measurements using multi-frequency GPS and GLONASS observations[J]. EURASIP Journal on Advances in Signal Processing， 2014， 2014(1): 50.

[17] Li Weiqiang， Fabra F， Yang Dongkai， et al. Initial results of typhoon wind speed observation using coastal GNSS-R of BeiDou GEO satellite[J]. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing， 2016， 9(10): 4720-4729.

[18] 邹文博， 张波， 洪学宝， 等. 利用北斗GEO卫星反射信号反演土壤湿度[J]. 测绘学报， 2016， 45(2): 199-204.

Zou Wenbo， Zhang Bo， Hong Xuebao， et al. Soil moisture retrieval using reflected signals of BeiDou GEO satellites[J]. Acta Geodaetica et Cartographica Sinica， 2016， 45(2): 199-204.

12月18日上午10时，庆祝改革开放40周年大会在北京人民大会堂举行，中共中央总书记、国家主席、中央军委主席习近平出席大会并发表重要讲话。六师一〇六团组织机关、事业单位全体人员，及社区、医院、学校、连队“两委”人员在团机关集中观看庆祝改革开放40周年大会直播。

[19] Wang F， Yang D K， Li W Q， et al. A new retrieval method of significant wave height based on statistics of scattered BeiDou GEO Signals[C]//Proceedings of the 28th International Technical Meeting of The Satellite Division of the Institute of Navigatio. Tampa， Florida: Tampa Convention Center， 2015.

[20] Zhang Yun， Tian Luman， Meng Wanting， et al. Feasibility of code-level altimetry using coastal BeiDou reflection (Beidou-R) setups[J]. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing， 2015， 8(8): 4130-4140.

[21] Zhang Yun， Li Binbin， Tian Luman， et al. Phase altimetry using reflected signals from BeiDou GEO satellites[J]. IEEE Geoscience and Remote Sensing Letters， 2016， 13(10): 1410-1414.

[22] Bilich A， Larson K M， Axelrad P. Modeling GPS phase multipath with SNR: case study from the Salar de Uyuni， Boliva[J]. Journal of Geophysical Research: Solid Earth， 2008， 113(B4): B04401.

[23] Levin M P. Numerical recipes in Fortran 90: the art of parallel scientific computing[J]. IEEE Concurrency， 1998， 6(4): 79.

[24] Larson K M， Nievinski F G. GPS snow sensing: results from the EarthScope Plate Boundary Observatory[J]. GPS Solutions， 2013， 17(1): 41-52.

1.实现源头数据采集的信息化。围绕采油队源头数据采集，建立单油井、单水井的运行参数、投入耗材、产出油量计量的数据自动采集网络，为生产管理提供全面、可靠、准确的没有人工干预的一性资料。针对联合站仪器仪表应用广泛的实际，强化仪器仪表的自动化数据采集。根据生产实际需求，实现大罐计量、原油盘库、管道防泄露在线检测的数据采集、数据远传。