专利名称：大面积土壤湿度测量方法
技术领域：
本发明涉及农业环境监测领域，尤其是一种对大面积裸露土壤湿度进行测量的方 法。
背景技术：
土壤湿度即土壤含水量，是土壤的重要物理性质。土壤湿度是土壤_植被_大气 系统能量与物质交换及土地资源和环境管理的重要参数之一，监测地表土壤湿度的时空变 化特性对研究自然及生态环境问题有重要的意义。目前土壤湿度的测量方法主要有下面几种，它们有着各自的优点和缺点。烘干称重法，直观准确，但其操作复杂。电容法(电容式土壤湿度传感器及其制作 方法，专利号200510041258)、电阻法(土壤湿度分层测量系统，专利号200710018822)、电 导率测试法(一种土壤湿度测试仪，专利号200910098049)等利用土壤中电容、电阻或电导 率与与含水量相关的关系，将两个电极埋入土壤中，测定土壤块的电容、电阻或电导率，求 出土壤湿度。但这些方法都是接触式方法，需要将探头实地插入土壤中，费时费力。射线法包括中子散射法、Y射线法、X射线法等(土壤湿度和雪中储水量测量仪， 申请号89103760)。这些方法价格昂贵，且存在潜在的辐射危害，近年来已经在发达国家遭 到弃用。时域反射仪法(Topp G C,Davis J L. Measurement of soil water content using Time Domain Reflectometry(TDR) :A field evaluation. Soil Sci. Soc. Am. J. 1997， 61(6) 1576-1585)目前已成为研究土壤水分的基本仪器设备。其通过测量高频电磁脉冲 在土壤中的传播速度求得土壤的介电常数，从而计算出土壤的含水量。但其是取样式测量， 测管周围的土层会受到严重扰动，也会导致较大的测量误差。此外，TDR电路复杂，设备昂
虫
贝ο专利(土壤湿度估计的系统和方法，专利号200910098049)提出利用自然界的闪 电等信号进行土壤湿度的测量，则其易受到闪电的随机性影响，效果难以保证。应用GPS信号的反射信号测量土壤湿度是最近几年发展的一项新技术。GPS信号 位于L波段，可以直接穿透云雾，完全不受云盖的影响而到达待测地面。利用现有的GPS系 统，可以全天候进行试验，快速获得大量数据，并在短时间内对大片区域进行探测。此外，由 于其只需要接收机而不需要发射机，会显著的降低系统费用、复杂度和功耗。GPS信号测量 还是一种非接触式的土壤湿度测量方式，可以大大地减少劳动强度。国外研究者曾利用装备有两套天线系统的接收机在飞机上进行了土壤湿度遥感 实验(参见丁金才主编的《GPS气象学及其应用》，气象出版社2009年)，武汉大学电波传 播实验室曾利用两个天线和两个接收机，通过分别接收直达信号和经由地面反射的反射信 号，然后根据电磁反射模型和土壤湿度的半经验介电常数模型来反演得到土壤湿度，具体 实施细节可见文章(严颂华张训械，基于GNSS-R信号的土壤湿度反演研究，电波科学学报 2010年第1期)。该方法可利用直达信号与反射信号的反射功率的差异，具有一定的实用价值。但这些方法需要依靠复杂的接收机设计。主要问题是需要两个天线和两个接收 机。由于这种方法主要利用的是接收信号功率的比值，两个天线及接收机必须具有很好的 一致性，即具有相同的增益才不至于影响反演的结果。而完全一致的天线和接收机生产会 增加生产和调试的成本。一种替换的方法，即进行通道校正的方法需要额外的校准过程，给 应用带来麻烦。(严颂华张训械，基于GNSS-R信号的土壤湿度反演研究，电波科学学报2010 年第1期)。可见，本领域亟待研究简单易实施的大面积土壤湿度测量方法。
发明内容
针对现有技术的局限性，本发明的目的是利用GPS直达信号及由地面土壤反射回 来的信号，提供一种对于大面积裸露土壤进行湿度遥测的方法，应用于农业生产。本发明提供的技术方案是一种大面积土壤湿度测量方法，设置一个垂直极化接收 天线、一个GPS接收机，以及进行信号处理的计算机，测量过程包括以下步骤，步骤1，利用垂直极化天线同时接收来自空中的直达信号和来自地面的反射信号， 形成干涉信号；所得干涉信号由GPS接收机接收后输出至进行信号处理的计算机；计算机 连续采样获取干涉信号的功率序列。步骤2，计算机对干涉信号的功率序列进行一阶和二阶微分运算，并根据运算结果 给出干涉信号的所有峰值点和谷值点；然后通过峰值点的最小值和谷值点的最大值搜索拐占.
^ \\\ 步骤3，根据步骤2搜索所得拐点获取对应的入射角值；利用入射角值和介电常数 的关系求解此时对应的介电常数；利用介电常数与土壤湿度的半经验模型计算土壤湿度。本发明提出了一种利用直达波和反射波的干涉技术以及拐点的搜索技术，满足了 节省接收机成本的基本要求。本发明的优势在于其出色的实用性能完全不需要复杂的接 收机改造，不需要另外获取直达通道的相干功率，节省了硬件设计与制造的费用；利用了 GPS卫星的自身移动获取不同入射角时的幅度序列，天空中多颗卫星可利用，信息利用率较 高，具有良好的精度和稳健性；利用了一维信号序列处理，运算量较小，能满足实时性要求。 实施本发明所提供技术方案，在提高探测性能的同时，大幅降低了设备的研制成本和维护 费用。
图1为本发明实施例的系统结构示意图；图2为本发明实施例的干涉信号及其一阶与二阶微分运算结果示意图；图3为本发明实施例的干涉信号的峰值、谷值和拐点判断示意图；图4为本发明实施例的干涉信号处理流程图。
具体实施例方式针对现有方法的局限性，本发明的目的是利用GPS直达信号及由地面土壤反射回 来的信号的干涉现象，提供一种对于大面积裸露土壤进行湿度遥测的方法，应用于农业生产。该方法利用垂直极化天线接收直达信号与土壤反射信号形成的干涉信号；利用干涉信 号的特征，即干涉信号功率序列的相位拐点对应着一个特定的信号入射角；而根据此入射 角来求取土壤表面的介电常数，从而进一步反演土壤湿度。本发明考虑到，如果一个天线能够同时既接收到直达波信号，又接收到反射波信 号，那么，就可以避开双天线的要求。但此时反射功率的比值也消失了，必须利用另外的技 术。当天线既接收到直达信号，又接收到反射信号时，会发生干涉现象。由于直达信号与反 射信号来自于同一个信号源，所以他们有着相同的频率。但反射信号来自于地面的反射， 所以相比直达信号有着功率的衰减和固定的距离差引起的相移。这样的两列电波(直达 信号与反射信号)满足相干叠加的条件，会出现明显的干涉条纹。由于两列波的相位差和 功率差与土壤介电常数和入射角有关。所以，干涉条纹并不均勻。根据仿真和实测观察该 图形可知，合成信号的图形实际上是一个调制信号，其包络为指数函数，根据调制信号的规 律，对应于入射角Ytl时信号整体会发生相位突变，包络在这里出现一个拐点。该拐点(入 射角Ytl)的位置与土壤介电常数ε有关。如果我们能够找到这个拐点及其对应的信号 入射角Ytl，就能够求解得到对应的土壤介电常数ε，然后根据介电常数与土壤湿度之间 的对应关系如 Topp 经验公式(参见 Topp G C, Davis J L, Annan A P. Electromagnetic determination of soil water content measurement in coaxial transmission lines. Water Resources Research. 1980，(2) :574 583)，来求取土壤、湿度。目前GPS卫星星座由21颗工作卫星和3颗在轨备用卫星组成。GPS距离地表 20200km,这种卫星的分布使得在全球任何地方、任何时间地面接收机都可观测到4颗以上 的卫星。位于地平线以上的卫星颗数随着时间和地点的不同而不同，最少可见到4颗，最多 可见到11颗。这样就保证了接收信号的连续性。GPS卫星发射的信号为右旋圆极化波，经过近似光滑的地表反射后，变为左旋圆极 化波。在左旋圆极化波和右旋圆极化波中，都存在水平极化和垂直极化的成分。在菲涅尔 假设下，直达波经过地面散射后的极化特性改变过程可以用地面反射电场表示，极化电场 可以表示为极化矩阵与入射信号的乘积 其中，圮、￡；分别为入射复数场的右旋和左旋成分，玛、拉分别为反射复数场 的右旋和左旋部分，极化矩阵可表示为 式中Γ\，Γ ,分别为在光滑表面电波水平极化和垂直极化分量的菲涅尔反射系 数，可由Fresnel公式给出如下 其中，γ为电波射线入射角，ε为地表相对介电常数。在本发明中，我们设置接收天线波束指向水平方向，提取天线中的垂直极化分量，装置设置如附图1所示，实施例设置有一个垂直极化天线1 (以下简称天线)、一个GPS接 收机2，以及进行信号处理的计算机3。垂直极化天线1与GPS接收机2连接，GPS接收机 2与计算机3连接。垂直极化天线1作为接收天线使用，其上不仅接收到来自空中的GPS 卫星4的右旋直达信号，而且接收到来自地面5反射的左旋反射信号。直达信号和反射信 号都来自于GPS卫星同时发出的信号，因此他们具有相同的初相。若设直达信号水平分量 Esh，则反射信号的水平分量的功率为rh · Esh ；若设直达信号垂直分量Esv，则反射信号的垂 直分量的功率为Γν·Ε3ν，记为Ε 。由于反射信号来源于地面的反射，根据电波反射的镜像 原理，其路程差相比直达信号要多，为此引起相位相比直达信号的滞后。由于GPS发射信号 源距离GPS接收机2数千米以外，当GPS信号通过两条路径，到达GPS接收机2时，可以视 为平面波而非球面波，这样反射信号具有下面的几个特点1、与直达信号具有相同的频率； 2、相比直达信号，相位上多了由于路程差而增加的相位差；3、功率上增加了由于反射而引 起的衰减量。根据附图1的图示，若设天线假设的高度为h，直达信号的垂直极化分量为Esv，则 反射信号的垂直分量E 相比直达信号的垂直分量Esv，存在的相位延迟为φ = (2hcos γ/λ)Χ23 = 43ih*cosy/A(4)其中λ为19cm，是GPS Ll信号的波长，、为信号的入射角。同时反射信号的垂 直分量E 的功率相比直达信号的衰减可以利用反射系数1\来度量，即反射信号的垂直分 量E 的功率大小为Γ vEsv，因此在t时刻天线上接收到的垂直极化信号功率为两者的叠加 (干涉信号)Sv (t) = Esv (t) +Erv (t) = Esv (t) + Γ vEsv (t) exp (j Φ t) = Esv (t) (1+ Γ vexp (j Φ t)) (5)由于我们需要获取的是干涉信号在不同角度时的功率。一种方式是卫星不动，控 制接收机天线在不同角度旋转以改变入射角，显然这个控制机构成本是较高的。实际上我 们可以设定接收机天线静止而利用GPS卫星自身的运动。GPS卫星在空中沿着固定的轨道 飞行，则在一天之中GPS信号会从不同角度照射到地面，且该角度变化是随时刻t连续变化 的。因此，我们能够获取到干涉信号在不同角度的功率。实施例的GPS接收机2与计算机3通过USB接口连接。当天线维持不动时，干涉 信号被GPS接收机2连续采样后，通过USB接口进入计算机3，成为一个随角度变化的功率 序列S(i)，其中i为角度变化值。接下来就是在功率序列Si中寻找拐点。这里首先利用五点滑动平均法对功率序 列s(i)进行平滑，以滤除各种冲击干扰和毛刺，然后对S⑴进行一阶和二阶微分处理，通 过对一阶和二阶处理的结果判断来寻找功率序列的最大值和最小值；然后在最大值和最小 值中寻找拐点；提取该拐点对应的入射角Y O值，根据公式cos2 Y0- ε2-ε +sin2 Y0 = O来求得ε。最后利用半经验模型计算土壤湿度。可见，本发明的优势在于其出色的实用性能1、只采用一个天线和接收机就可以 实现土壤湿度反演，不需要另外获取直达通道的相干功率，节省了硬件设计与制造的费用； 2、利用了 GPS卫星的自身移动获取不同入射角时的功率序列，天空中多颗GPS卫星数据可 利用，结果可互相校验，信息利用率较高，具有良好的精度和稳健性；3、利用了一维信号序 列处理，运算量较小，能满足实时性要求。在提高探测性能的同时，大幅降低了设备的研制成本和维护费用。本发明的关键在于通过接收干涉信号，将土壤湿度估计问题转化为功率曲线中的 拐点搜索问题，并利用介电常数与土壤湿度之间的关系得到土壤湿度。以下分步骤对本发明的一个实施例，即天线及接收装置为地基时进行描述，参见 图4 步骤1、干涉信号的接收如附图1所示，在高度h (具体实施时可设Im < h < IOm)处设置一颗垂直放置的 微带天线作为垂直极化天线1，其接收波束指向水平方向，则该天线接收到的是GPS卫星4 发来的直达信号中垂直分量与经由地面而来的反射信号中的垂直分量的干涉信号。为了增 加算法的稳定性，从所有GPS卫星信号中按照高信噪比和仰角分布适度的原则选择信噪比 高，仰角处于30-60度范围的GPS卫星作为使用的卫星，Sv(t) 二 max匁⑴,j为接收到的所有N颗GPS卫星的PRN号·
H其中S((t)为利用垂直极化天线接收到的第j颗卫星的干涉信号利用GPS接收机2采集其功率回波信号作为待用的功率序列S⑴。功率序列S⑴ 为天线上接收到的两路垂直极化信号功率的叠加(干涉信号)S(i) = Ev(i) + rv(i)Ev(i)exp(j^ (i)) = Ev(i) (1+ Γ v(i) exp (j4 Ji h ·cos y (i)/ λ)) (6)其中Θχρ(」Φ⑴)表示eJ4>a)其中Ev(i)，rv(i)，φ⑴表示在第i个入射角时直达信号功率、反射系数大小和 反射信号相移的大小。实际接收的干涉信号s(i)的图形如附图2所示。其中实线6为干 涉信号，虚线7为对干涉信号进行一阶微分运算的结果，实线8为对干涉信号进行二阶微分 运算的结果，圆圈9为一阶微分运算结果中的过零值。incident angle ()表示入射角。步骤2、干涉信号中拐点的搜索。理论分析和仿真实验表明，实际接收到的干涉信号功率序列S(i)会存在扰动。扰 动的原因包括接收机中噪声的干扰，土壤地表面的不平等等，使得接收到的包络并非如理 想的那样恰好服从一个平滑的分布规律，若简单的在存在干扰的情况下进入相位突变点的 搜索，会引起错误的估计，因此进行平滑预处理的效果更好。实施例采用了一些低运算量的 处理方式，即首先对其进行平滑，然后通过一阶和二阶微分处理获得干涉信号的峰值和谷 值，最后再进行搜索拐点，由于维数为一维，该搜索的运算量不大，可满足实时性要求。具体步骤如下1)对干涉信号的样本S(i)进行五点平滑滤波。以去除功率信号中的噪声干扰和
毛刺。实施例采用图像处理中使用的五点平滑滤波法来去掉毛刺。基本公式为 1 <+2S](i)=-YS(i) (7)
5 ,-2其中S1⑴表示滤波处理后的结果。2)对S1 (i)以一阶微分算子进行计算。其一阶微分运算的结果为S2(i) = S1 (i+1)-S1(I) ,i = l, ...N (8)
7
3)对S2(i)再进行一次微分运算，即对S1G)进行二阶微分算子计算。结果为S3(i) = s2(i+l)-s2(i) ,i = l, ...N (9)则一阶微分计算结果S2 (i)为0的地方对应着干涉信号S1 (i)的峰值和谷值点处， 而峰值和谷值的确定可以通过判断s3(i)的正负号来确定，当S3(i) < 0时，i点对应干涉 信号S1⑴的波峰值，即最大值；当S3⑴>0时，i点对应干涉信号S1G)的波谷值，即最 小值。附图2给出了入射角incident angle ()为1_50度时干涉信号及其一阶与二阶微分
运算结果。求取拐点。拐点是干涉信号S1G)的所有峰值点中的最小点和所有谷值点中的最 大点。若两者所处的角度相同，则任取一个作为入射角Yci即可。若两者的所处的角度有 细微差别，则取其平均值作为入射角Υο。拐点的搜索图示见附图3。其中10为峰值点，11 为谷值点，12为判断出的拐点。4)步骤3、根据入射角值和介电常数的关系求得土壤湿度。由公式cos2 Y。· ε 2- ε +sin2 Yq = 0来求得对应的介电常数ε步骤4、根据介电常数与土壤湿度的关系式求得土壤湿度。虽然Dobson等人在1. 4GHz到18GHz范围内，经过对多种土壤实验，提出了针对含 水混合物的介电常数模型。在这个模型里，介电常数被描述为 其中，es> ea> ε fw和ε bw分别为土壤固体物、空气、自由水和束缚水的介电常数。 Va, Vfw, Vbw分别是空气、自由水和束缚水的体积比。在这一方法中，通过对土壤结构的详尽 描述来计算自由水和束缚水。但在实际应用中，很多要求土壤物理参数都无法准确获得，应 用性较差。因此，实施例采用简单实用的Topp模型。该模型建立了土壤体积含水量(wv， cm3 - cm-3)与介电常数ε之间的关系，即介电常数与土壤湿度的半经验模型，可以用一个三 次多项经验公式表示wv =-5· 3Χ1(Γ2+2· 92Χ1(Γ2 ε-5· 5Χ1(Γ4 ε 2+4· 3Χ1(Γ6 ε 3(11)于是可根据介电常数ε反演计算土壤湿度。另外，本发明的用途除了利用GPS信号测量土壤湿度之外，还可进一步应用于接 收两种以上的不同频率不同极化的信号，例如俄罗斯GL0NASS系统或我国的北斗系统的信号。本发明以实施例，即地基土壤湿度测量为参考来描述。然而，本发明的实施还可用 于塔基，山基等固定平台或车载等移动平台。
权利要求
一种大面积土壤湿度测量方法，其特征在于设置一个垂直极化天线、一个GPS接收机，以及进行信号处理的计算机，测量过程包括以下步骤，步骤1，利用垂直极化天线同时接收来自空中的直达信号和来自地面的反射信号，形成干涉信号；所得干涉信号由GPS接收机接收后输出至进行信号处理的计算机；计算机连续采样获取干涉信号的功率序列；步骤2，计算机对干涉信号的功率序列进行一阶和二阶微分运算，并根据运算结果给出干涉信号的所有峰值点和谷值点；然后通过峰值点的最小值和谷值点的最大值搜索拐点；步骤3，根据步骤2搜索所得拐点获取对应的入射角值；利用入射角值和介电常数的关系求解此时对应的介电常数；利用介电常数与土壤湿度的半经验模型计算土壤湿度。
全文摘要
本发明提供一种大面积土壤湿度测量方法，设置一个垂直极化接收天线、一个GPS接收机，以及进行信号处理的计算机，测量过程包括利用GPS信号及其反射信号的干涉效应来进行处理。这种技术方案利用广泛存在的GPS信号作为信号源，可以全天候大面积测量且成本低廉，并利用干涉效应更进一步简化了接收装置。
文档编号G01N22/04GK101900692SQ201010207399
公开日2010年12月1日 申请日期2010年6月18日 优先权日2010年6月18日
发明者严颂华, 龚健雅 申请人:武汉大学