基于探地雷达的盐碱地电导率测定方法
【专利摘要】本发明公开了一种基于探地雷达的盐碱地电导率测定方法，包括如下步骤：步骤一、使用探地雷达对土壤进行探测，获取探地雷达回波相关系数ρts；步骤二、使用公式计算土壤电导率σ，其中，ai(i＝1,2,…,n)分别为i阶系数。本发明通过在盐碱地选择不同类型的典型样地进行雷达回波探测，并将雷达回波与雷达发射机的发射子波进行相关运算，得到相关系数。然后，根据样点的实测电导率与雷达回波相关系数的对应关系，得到电导率与雷达回波相关系数关系，当根据待测盐碱地的雷达回波计算出相关系数后，就可以反演出探测点的盐碱土电导率。因此，本方法可快速的测量出探测点的电导率，进而分析出盐碱程度，具有较高的准确性。
【专利说明】
基于探地雷达的盐碱地电导率测定方法
技术领域
[0001] 本发明涉及雷达探测技术领域，特别涉及一种基于探地雷达的盐碱地电导率测定 方法。
【背景技术】
[0002] 目前，我国一些地区的盐碱化与沙化范围扩大速度正在加快，已引起各界人士的 关注。某盐碱地改良需要得到盐碱地的盐渍化程度信息，目前一般是通过测量盐碱土电导 率实现。盐碱地含盐量高，在水的溶解作用下，溶解盐离子会提高盐碱土的导电能力，因此 盐碱土的电导率比一般土壤大很多。地质雷达具有解决盐碱土分布状态和分布厚度的能 力，可通过地质雷达来了解盐碱土壤的含盐量、PH值、有机物含量、含水量等信息。
[0003] 中国专利申请号200910089204.6公开了一种土壤盐碱化的雷达遥感监测方法，基 于全极化SAR书记计算极化参数，进行裸地与植被覆盖区域的划分，生成掩膜图像，划定试 验区，利用测量数据建立含水含盐土壤介电模型，利用区域复介电常数和区域含水量，进行 盐土和碱土的区分。上述方法虽然能够获取土壤的含盐量，但过程复杂，准确性差。因此，需 要提供一种能够快速、准确的测量出土壤含盐量的方法。
【发明内容】
[0004] 本发明设计开发了一种基于探地雷达的盐碱地电导率测定方法，目的是通过雷达 回波相关系数反演出探测点的盐碱土电导率，以得到土壤盐碱化程度。
[0005] 本发明提供的技术方案为：
[0006] 基于探地雷达的盐碱地电导率测定方法，包括如下步骤：
[0007] 步骤一、使用探地雷达对土壤进行探测，获取探地雷达回波相关系数pts;
[0008] 步骤二、计算土壤电导率〇
[0009] 〇- = a0： + ￡t,pts + a2p^+ -- + aipl+-^ + anp^,
[0010] 其中，ai(i = l，2,…，n)分别为i阶系数。
[0011] 优选的是，步骤二中，n = 3。
[0012] 优选的是，至少采集四个试验点，分别测量导电率0(1)，〇(2)，…，〇(n)，并获取相 应的试验点的雷达回波相关系数Pt s(l)，Pts(2)，…，Pts(n)，代入方程
[0013] a(i) = a{) + a, pts (/) + a2pl (i) + (/), / = 1,2, ? ? ?,
[0014] 获得 a〇、ai、a2、a3l^t^。
[0015] 优选的是，米用最小二乘法计算a〇、ai、a2、a3l^t^。
[0016]优选的是，使用土壤水分温度电导率测试仪测量试验点土壤的电导率。
[0017]本发明的有益效果是：本发明提供了一种基于探地雷达的盐碱地电导率测定方 法，通过在盐碱地选择不同类型的典型样地进行雷达回波探测，并将雷达回波与雷达发射 机的发射子波进行相关运算，得到相关系数。然后，根据样点的实测电导率与雷达回波相关 系数的对应关系，得到电导率与雷达回波相关系数的拟合曲线，并求得以回波相关系数作 为自变量的电导率计算公式。得到先验公式以后，当根据待测盐碱地的雷达回波计算出相 关系数后，就可以反演出探测点的盐碱土电导率。因此，本方法可快速的测量出探测点的电 导率，进而分析出盐碱程度，具有较高的准确性。
【附图说明】
[0018] 图1为初始相位不同的同频正弦波相关系数曲线。
[0019] 图2为低电导率土壤介质回波波形图。
[0020] 图3为高电导率土壤介质回波波形图。
[0021 ]图4为反射系数曲线示意图。
[0022]图5为介电常数虚部与回波相关系数关系曲线。
[0023] 图6为测线1俯视图。
[0024] 图7为测线1雷达回波相关系数与电导率的关系示意图。
[0025] 图8为测线2俯视图。
[0026] 图9为测线2雷达回波相关系数与电导率的关系示意图。
[0027] 图10为测线1，2雷达回波相关系数与电导率关系散点图。
【具体实施方式】
[0028] 下面结合附图对本发明做进一步的详细说明，以令本领域技术人员参照说明书文 字能够据以实施。
[0029] 本发明提供了一种基于探地雷达的盐碱地电导率测定方法，分析雷达回波相关系 数与导电率的对应关系，通过获取回波相关系数得到探测点的盐碱土电导率。
[0030]考虑一个周期时间窗的某一频率信号W1表不为：
[0031] wi = Ai cos( ? it)
[0032] 令一同频率信号W2表不为：
[0033] W, - Az C〇s(〇\t -j-<3,)
[0034] 如图1所示，当物取值从0~2jt变化时，奶与《2的相关系数是相位差灼的余弦函数， 即相关系数可以用来衡量波形之间的相位差异。
[0035]盐碱土的高损耗主要是由于土壤中溶解盐离子浓度较高所导致的电导率升高，从 而使土壤的介电常数虚部增大，当介电常数虚部与实部相比不能忽略时，土壤介质的波阻 抗变成复数，进而土壤界面的反射系数r也变为复数，复反射系数会造成回波的相位偏移 (这与低损耗介质回波的无相位偏移有本质的不同），时域上体现为回波波形的畸变，如图 2、图3所示。
[0036] 考虑具有较高含水量与含盐量的土壤介质，由于电导率较高，介电常数虚部不能 忽略，令某一频率fn下土壤介电常数为e n，相对介电常数为e?，电导率为〇,对由电导率表征 的有损介质有介电常数表达式：
[0037] en=e/n-j5/〇n
[0038] 对一般损耗介质(介电常数虚部不能由低频电导率直接得到)作线性近似：
[0039] e//rn = ain8n+a〇n
[0040]其中aln，aQnS常数，aQn表征除电导率以外因素导致的介质损耗，则有：
[0041 ] ￡n-￡rn ? ￡〇 -( ￡ rn-j ￡ rn) * ￡0 -( ￡ rn_ j (EllnSn+ElOn) ) ? ￡〇
[0042]其中e〇为真空介电常数，考虑实际高电导率土壤介电常数的实部与虚部：
[0043] e/rn>〇,e//rn>〇
[0044] 土壤介质在频率fVf的波阻抗为：
[0046]其中Zo为真空波阻抗，将Zn归一化，可写为：
[0049]考虑土壤介电常数实部与虚部的符号，可得：
[0052]高损耗介质条件下电磁波穿透能力非常有限，假设地下无明显反射体，则可忽略 土壤深层产生的反射波，而仅以土壤表层反射回波参与计算。令土壤表层在频率fn的反射 系数rn:
[0054] 其中，Zo为空气波阻抗377 Q )。
[0055] 雷达回波(不考虑深层反射)包括直达波与地表反射波，因此频率为fn的雷达回波 幅度与相位可由一等效反射系数R n确定：
[0056] Rn = Adn+ r n= (Adn+ T nr)+j T ni
[0057]其中Adn为直达波系数(令雷达接收的直达波与发射子波相比仅在幅度上有衰减， 相位关系相同），图4为加入直达波系数的雷达回波等效反射系数曲线，可以看到曲线上每 一反射系数点沿逆时针的方向变化时，相位的变化是非线性的。
[0058] 对每一反射系数的相位取余弦再取反，得到取反后的频率fn雷达回波与发射子波 的相关系数(_Ptsn)，该回波相关系数与介电常数虚部的关系如图5所示，介电常数虚部与回 波相关系数呈单调上升的非线性关系（该曲线可用3次多项式拟合），考虑之前的假设一一 介电常数虚部与电导率呈线性关系：
[0059] e//rn = ain8n+a〇n (14)
[0060] 故可得电导率0与雷达回波相关系数-Ptsn呈与图5中曲线形式相近的单调关系。
[0061] 进一步，由傅式变换可将雷达发射的脉冲信号表示为通频带内各频率分量的线性 组合：
[0062] sv,(m) = /L. co^{a\t + ) + A2 cos(&2f -f ^) -f < ? > -f coslf%/1 -f (pn)
[0063] 则雷达回波可表示为各频率分量反射波的线性组合： + /'V>A cm(切y + 玖> + {，1,.、+ .厂、.M'、，十炉、）
[0064] …-' --' -4-.,. < (/I 4- / ^A.,. COS( J -f- 0.,)
[0065] 其中各频率分量匕的回波相关系数ptsn与该频率下介电常数虚部关系曲线都与图 5类似，将整体的雷达回波与雷达子波的相关系数近似表示为各频率分量相关系数的线性 组合：
[0066] Pts - Cl ? Ptsl+C2 ? Pts2+. ? ?+Cn ? Ptsn
[0067]则雷达回波与雷达发射子波的总相关系数Ptsn与电导率〇也呈近似图5的非线性关 系。
[0068]令雷达发射子波为Wt(m)，高电导率土壤介质回波为wsk(m)，取相关系数表征回波 的相位偏移：
[0070] 经盐碱地实验发现相关系数与介质电导率有较好的对应关系，令不同电导率土壤 的相关系数集合为{_Pts(l)，...，_Pt s(k)，...}(相关系数取反来获得与电导率一致的对应 关系），对应的实测土壤电导率集合为{>(1)，...， 〇(k)，...}，令
[0071] f .,v ........ i <>{? ) ：= i?,. -f a, {-p:. (k)) + a,(~p,s (k)}" -f a. (-pL, {k}}5 *
[0072] 解超定方程得到系数ao~a3的估计值，则可以根据雷达回波的相关系数pts(k)反 演土壤电导率〇(k)。
[0073]首先使用探地雷达对土壤进行探测，获取至少四个试验点处的探地雷达回波相关 系数Pts(l)，Pts(2)，…，pts(m)，并使用土壤水分温度电导率测试仪测量对应的试验点处的 电导率 〇(1)，〇(2)，…，〇(m)。其中m为大于或等于四的整数。
[0074] 然后将相关系数Pts(l)，Pts(2)，???，pts(m)和电导率〇(1)，〇(2)，???，〇(m)带入到下 式中
[0075] 〇-(/) = a{) + a,prs (i) -f a2pfs (/) + a^pl (/), i = 1,2,*??,m
[0076 ]使用最小二乘法计算ao、ai、a2、a3的解。
[0077]最终获得了相关系数与电导率的对应关系：
[0078] a = au + ci\Pts + qpl + a7>pi。
[0079] 利用上式，只需使用探地雷达对土壤进行探测，获取其中一个探测点的相关系数 后即可得到该探测点的电导率。通过电导率即可分析出该处土壤的盐碱程度。
[0080] 例如在每一次在探地雷达实验中，选取具有代表性的两条测线说明雷达回波与样 地水盐特性之间的关系。
[0081] 如图6所示，测线1样地，测线长度90m，该样地地表特征多样，测线最左侧是一片地 表开裂的干涸小水塘，由于地表开裂严重，初步估计在雨水淋溶与冲刷作用下，盐分很大一 部分被水从裂缝处带到深层的地下，因此浅表土层的含盐量应该较低;水塘右侧经过4米左 右的裸土地带后，是一片宽度2米左右的有稀疏植被的条带;条带右侧又是长度为20多米的 裸土地带;再右侧是宽度大约3米左右的具有细密植被的草地条带;草地右侧经过另外一片 裸土后，是最右侧的一片干涸的水塘，该水塘地表特征与最左侧的水塘不同，地表无裂纹， 土壤疏松，与重度盐碱地的致密土质有很大不同，疏松的土壤在雨水淋溶渗透作用下，盐分 上下输送通畅，因此初步估计含盐量也不高。
[0082] 根据雷达回波成像，可以清晰分辨干涸的水塘、裸土、稀疏植被、草地所在的区域 界限。在测线雷达经过的每一点，使用WET土壤水分温度电导率测试仪测量了土壤体积含水 量、介电常数实部与电导率。
[0083] 图7为电导率曲线与雷达回波相关系数曲线对比图，图中可以看到两条曲线的变 化趋势高度契合，体现了雷达回波相位变化与盐碱土电导率的相关关系。由于WET单点测量 范围大约数平方厘米，探地雷达波束覆盖范围大约1平方米，2种测量手段的尺度不同，因此 造成了在地表介电特征急剧变化的位置，2者的趋势虽然一致性好，但测量值绝对差异较 大，契合度相对较低。
[0084] 如图8所示，测线2样地，测线长度30m，该样地最右侧紧邻一个较大的水塘，地势从 草地右侧开始呈明显的左高右低（降水时，由于盐碱土透水性差，右侧地表的很大一部分积 水会流向右边的水塘，同时冲刷沿途的地面，造成盐度降低）。基于测线1的先验知识，根据 雷达回波成像结果，可初步分析测线15m以前（左侧）电导率较高，15米以后（右侧靠近水塘 的部分)电导率急剧降低，而含水量变化不大(虽然目测地表特征仍为典型的重度盐碱地）。 [0085]图9为WET电导率曲线与雷达回波相关系数曲线对比图，与测线1的结果类似，WET 电导率曲线与雷达回波相关系数曲线趋势的契合度仍然较好，电导率的变化趋势结果与雷 达回波成像结果也大致吻合。
[0086] 图10为测线1，2所有WET电导率与雷达回波相关系数的对应关系散点图，由研究方 案部分得到的电导率与雷达回波相关系数的对应关系，这里采用3次多项式进行曲线拟合：
[0087] y = 9.084 ? x3+20.02 ? x2+19.84 ? x+8.81
[0088] 如图10中曲线所示。利用该公式可在一定含水量条件下，根据雷达回波快速确定 当前盐碱土的电导率。
[0089] 尽管本发明的实施方案已公开如上，但其并不仅仅限于说明书和实施方式中所列 运用，它完全可以被适用于各种适合本发明的领域，对于熟悉本领域的人员而言，可容易地 实现另外的修改，因此在不背离权利要求及等同范围所限定的一般概念下，本发明并不限 于特定的细节和这里示出与描述的图例。
【主权项】
1. 基于探地雷达的盐碱地电导率测定方法，其特征在于，包括如下步骤： 步骤一、使用探地雷达对土壤进行探测，获取探地雷达回波相关系数pts; 步骤二、计算土壤电导率σ其中，ai(i = 1，2，…，η)分别为i阶系数，η为阶数。2. 根据权利要求1所述的基于探地雷达的盐碱地电导率测定方法，其特征在于，步骤二 中，η = 3〇3. 根据权利要求2所述的基于探地雷达的盐碱地电导率测定方法，其特征在于，至少采 集四个试验点，分别测量导电率σ(1)，〇(2)，…， 〇(m)，并获取相应的试验点的雷达回波相关 系数Pts(l)，Pts(2)，…，Pt s(m)，代入方程获得ao、ai、a2、a3的解。4. 根据权利要求3所述的基于探地雷达的盐碱地电导率测定方法，其特征在于，采用最 小二乘法计算ao、ai、a2、a3的解。5. 根据权利要求4所述的基于探地雷达的盐碱地电导率测定方法，其特征在于，使用土 壤水分温度电导率测试仪测量试验点土壤的电导率。
【文档编号】G01N33/24GK105928987SQ201610236199
【公开日】2016年9月7日
【申请日】2016年4月15日
【发明人】吴斌, 赵凯, 李晓峰, 姜涛, 郑兴明, 李晓洁
【申请人】中国科学院东北地理与农业生态研究所, 吉林大学