准确地获取农田多源异构信息是实施精准农业的基础^［1-3］。目前，基于无线传感器网络（Wireless Sensor Networks，WSNs）的农业环境和作物生长监测已成为农业生产监测的重要方向^［1,2］。由于监测面积较小、供电便利等原因，WSNs技术在温室大棚、智慧养殖等设施农业场景下应用效果较好^［4,5］。与设施小气候环境^［6-8］不同，在农田场景数据采集中，由于WSNs监测需要覆盖整个耕作区域，覆盖面积大，在单跳模式下，传感器节点与sink节点距离较远导致数据采集过程中消耗大量能量。而且农田中物联网监测采集时间往往会持续数月，而高频次、高能耗的数据传输会严重地制约传感器网络稳定工作的时长^［9,10］。

为实现低能耗的数据传输路径对数据进行采集，Mittal等^［11］利用分簇算法以及双跳通信实现负载均衡和能量最小化，延长远程数据采集监测时间^［12］，但未将数据收集质量考虑进去，若将该方法应用到农田数据采集中，会降低数据的精确性^［13］。Verma等^［14,15］采用群体智能算法，利用剩余能量和节点部署密度得到适应度函数，实现对信道的优化选择，并通过网络节点分簇和多个数据汇聚的方式实现小区域数据采集^［16］，但该方法若应用在农田中，可能由于监测的农田面积较大导致低效的数据采集。樊志平等^［17］提出了一种基于能量有效负载均衡的多路径路由策略（Energy-Efficient Load Balancing Multipath Routing Scheme，EMR），该算法生成多条源节点与目的节点间的数据传输链路，通过链路总跳数和传输能耗计算权值的方式，选择能耗低、时延小的数据采集方式。但EMR算法未考虑每一跳的传输距离，存在两中继转发节点间通信距离较远的情况，并且EMR算法是一种平面型路由算法，全网节点都需按时与sink节点进行数据上传，若利用该算法对农田环境及作物数据进行采集，则会导致网络能耗开销较大，部分节点快速死亡，无法精确地、长时间地进行数据采集，且会造成采集的数据冗余度较高。Nouh和Geta^［18］提出了一种贪婪外围无状态路由算法（Greedy Perimeter Stateless Routing-Algorithm，GPSR-A）该算法在保证数据传输距离最短的同时，网络传输时延较低。但根据邻居节点相对位置的单一路由选择方式，未考虑转发节点的剩余能量，导致部分节点高频率的参与数据转发而加速能量衰竭，最终影响网络的监测质量与生命周期。孙想等^［19］提出了一种能量高效的农田WSNs拓扑关联路由算法（Topolopy Based Cluster-head Switching-Energy Approximation，TBCS-EA），该算法以节点位置及能量对于簇首数据汇集阶段的能耗影响，通过动态构建网络拓扑结构来均衡节点间能耗差异，并通过能量逼近的簇首轮换方式，降低网络算法开销。此算法在小规模网络中性能较好，但由于农田网络监测数据采集过程中存在规模大、周期长等特点，TBCS-EA无法满足针对农田的网络监测数据采集。

大规模农田目前并无明确的定义，但多认为大规模农田应适宜于农业机械化作业生产的需求，一般认为地块面积应在70,000 m²以上。大规模农田网络监测数据采集是一种能量供给困难的WSNs应用，对网络传输的能耗与时延都有着严苛的要求。现有农田WSNs研究多采用增加中继转发节点的方式，构建网络数据汇集中的多跳数据传输链路，以多次短距离数据传输降低远端节点在数据汇集中的能量损耗^［20］。但网络结构，尤其是多跳节点的路径选择也在很大程度上决定了网络能耗、时效、可靠性等性能。由此，本研究提出了一种基于粒子群（Particle Swarm Optimization，PSO）和模拟退火（Simulated Annealing，SA）协同优化的农田物联网混合多跳路由算法（Particle Swarm Optimization and Simulated Annealing Multi-hop Routing Algorithm for Hybrid Networks，PSMR），通过选择具有较高剩余能量和较低节点度的优势节点，动态划分网络结构，利用多特征值加权的方式筛选与sink节点间传输能耗低和跳次均衡的数据传输链路，最终实现WSNs对农田高效、可靠的多源数据采集。

为对农田网络高效、低能耗进行定量分析与研究，首先需要建立农田网络监测模型。模型构建主要分为两方面，一方面需要定义网络模型对网络节点的拓扑组织关系，另一方面需要定义网络的能耗模型以便进行能耗分析与仿真。

由于大规模农田场景WSNs覆盖范围大，传感器节点数目多，成簇结构网络中“簇首-Sink节点”间单跳的数据模式易导致远端簇首节点快速死亡^［21］。因此，本研究采用一种下层成簇、簇首多跳的混合型网络结构，如图1所示。

本研究提出的PSMR算法采用了上下分层的混合型网络结构，该结构中下层采用成簇机制，提高组网效率，而且相邻区域节点成簇为网络数据融合提供基础；上层采用平面多跳结构，有利于均衡不同区域能耗。为便于分析，可将大规模农田WSNs中参与数据汇集的节点抽象化为一个无向赋权图G（V，E）。其中V代表的是WSNs当前轮次的簇首集合，V=｛C₁，C₂，…，C_n｝，每一个节点的通信半径为d_o，e是相邻节点C_i与C_j间的数据传输路径，e=（C_i，C_j）∈ E，即E为相邻节点间传输路径集合。G中的一条数据传输链路P（C₁，sink）是C₁到sink节点间边的有序组合序列，P（C₁，sink）=（（C₁，C₂），（C₂，C₃），…，（C_i，sink））。而路由优化的目的就是寻找满足簇首与sink节点间低传输能耗与时延要求的数据传输链路P。路由算法需要解决的问题即是为每一个节点找到一条到sink节点的路径，使得网络整体上能耗性能最优。结合参考文献［21，22］中的仿真模型，对网络模型作以下假设。

（1）N个传感器节点随机布置在M×M m²的正方形农田中。

（2）传感器节点具有环境感知、数据融合及传输功能。

（3）传感器和sink节点在部署后，在各自初始位置保持静止状态。

（4）传感器节点能量有限且无法补充。

（5）所有的传感器节点都有簇内传输与簇间传输两种传输模式，但只有当选择簇首后才可以进行较大功率的簇间数据传输。

根据节点间数据传输距离的差异性，本研究采用了两种不同的能耗模型：当节点数据传输距离小于距离阈值d_o时，采用自由空间衰减模型计算节点能耗，其路径损耗指数为2；当节点数据传输距离大于等于距离阈值d_o时，采用多路径衰减模型计算节点能耗，其路径损耗指数为4。然后由此节点向距离为d的另一节点发送k b数据时能耗E_tr的计算方式为^［23］：

其中，E_ele为节点处理1 b数据所需能量；ε_fa和ε_mf分别为两种模型中功率放大器的能耗系数；距离阈值d_o由能耗放大系数决定，是ε_fa与ε_mf比的算术平方根。

传感器接收k b数据的能耗E_re为：

混合结构网络的路由构建主要分为两部分，首先是底层的成簇，网络节点根据节点位置、能耗等参数确定权重，竞争成为簇首，未成为簇首节点选择邻近的簇首加入。完成成簇后，簇首再以平面路由方式组网连接，形成各自到sink节点的路由路径。

由于簇首承担着局部数据的融合与上传任务，簇首的能耗远大于普通节点，因此本研究提出的PSMR算法采用剩余能量和节点度计算权值的方式筛选簇首，使节点密集区域的高能节点有更高的概率担任簇首，簇首竞争权值如式（3）所示。

其中，r是网络当前轮次；α与β为权重因子；T（i，r）是节点i在第r次的竞争权值；N（i，r）是节点i在第r次的节点度；E（i，r）是节点i在第r次的剩余能量。当节点i的竞争权值T大于全部邻居节点时成为簇首。当网络结构划分结束后，普通节点按照传输距离选择最近的簇首加入，并开始初始化簇首在与sink节点在数据汇集中的传输链路。

在簇首与sink节点的初始数据传输链路构建中，当前簇首根据邻居节点的通信能耗、剩余能量、与sink节点间距离作为选择下一跳节点进行数据转发的依据，具体公式如下。

其中，W( i, j )为链路选择权重；a、b、c是选择中继转发簇首的权重因子；i是数据发送节点；j是节点i的邻居节点；dist是两节点间的传输距离；E( j )是第j个节点的剩余能量；E_o是网络节点的初始能量。该公式表明当选择下一中转节点进行数据转发时，优先选择剩余能量多，与节点i之间通信距离更接近于d_o且离sink节点最近的相邻簇首，以此构造出初始数据发送簇首与sink节点间的数据传输链路。网络初始数据转发树如图2所示，远端簇头通过多跳转发传输至sink节点，避免了远距离传输高额能耗，可有效改善网络的能耗均衡性能。

WSNs应用场景不同，数据传输的要求也不同。WSNs在大规模农田中部署时，收集不同位置的作物生长及环境信息，并通过对这些信息分析后，对不同位置的作物采取不同措施，因此，需要环境监测精度高、数据传输可靠性强、节点适应性强。根据农田中不同作物的生长周期不同，WSNs在大规模农田进行多源异构数据监测和收集中属于一种能耗低的数据传输路径，使其能够在大规模农田中实现长时间的数据收集。本研究通过节点剩余能量、节点度加权选择簇首，采用成簇结构实现异构网络高效动态组网，通过簇首间多跳数据结构解决簇首远距离传输能耗过高问题，利用PSO与SA对数据传输路径的适应函数进行协同优化，提高算法的收敛速度，实现sink节点加速采集簇首中聚合的数据。

通过对大规模农田WSNs数据传输过程进行分析，发现大规模农田WSNs对数据汇集阶段的传输能耗、质量、时延都有着较为严苛的需求，因此本研究采用多参数加权求和的方式构建了数据传输链路优化的适应度函数，如公式（5）所示。

其中，fitness为适应度函数；w₁、w₂、w₃为权重因子；HopCount_min为数据传输链路的最小跳数；Energy_max为同一条数据传输链路的簇首平均最大能量值；dist_min是链路平均最小的单次数据转发距离，如公式（6）所示。

其中，dist()用于计算节点i与数据转发节点NextNode()内距离最短的节点间距离。每一次数据转发的路径长度，决定了转发节点的传输能耗。公式（6）使网络每一次数据转发的距离尽可能短且相近，使参与数据汇集的簇首能耗保持在较低的同一水平。

HopCount_min具体计算如公（7）所示。

其中，NodeCount( )用于统计该链路数据上传中转发簇首的数目，目的是使参与一次数据汇集的中间转发节点尽可能的少，减少数据转发过程中的传输时延，为农田WSNs决策提供高时效性的数据支撑。

由于转发节点承担着网络额外的能量消耗，高能量节点进行数据转发，在保证数据传输质量的同时，可以有效均衡网络负载。Energy_max具体计算如公式（8）所示。

其中，E( )是簇首的剩余能量值。

以网络传输距离、拟消耗传输能量和传输跳数构造的公式（5）适应度函数，本节提出一种基于PSO和SA协同优化方式，实现大规模农田基于WSNs的低能耗、高质量数据传输链路进行优质的、高效的数据采集。

PSO以其实现简单、收敛速度快、适用性强的特点，非常适合应用于节点计算能力薄弱、数据传输时效性要求较高的农田WSNs数据传输链路优化阶段。通过随机定义粒子位置和速度，计算每个粒子的适应值。由PSO算法率计算粒子的应适度函数值，P_i粒子的pbest值可以表示为P_i = （p_i1，p_i2，…，p_ij）。最佳粒子是当前情况群体中所经历的历史最优位置向量Pg = （ p_g1，p_g2，…，p_gD）。粒子的速度和位置可以用下面的等式表示。

其中，x_i，j和V_i，j表示在t迭代中当前速度和位置；p_i，j表示在第i个粒子的j维分量更新前的历史最佳位置；p_gj表示种群历史最佳位置；φ为权重值；c₁、c₂为学习因子；r₁、r₂为均匀随机数。

但PSO在优化过程中需要进行大量的迭代运算，加速了网络的能量损耗，且PSO易陷入局部最优解而浪费大量的计算资源。针对这一问题，本研究引入具有突跳能力的SA算法用于对PSO进一步处理，以此避免陷入局部最优节解。利用PSO得到粒子的应适度函数值P_i以及最优粒子的适应度为P_g，通过公式（12）得到新的适应度。

其中，f ( )是算法构建的适应度函数；S( )用于更新当前温度t下粒子i的适应度。

当粒子的当前适应度比以前更差时，更新将以一定的概率l_i被接受，如公式（13）所示。

其中，Fit-1表示粒子i的先前适应度；Fit代表粒子i的当前适应度；Te代表退火温度，该温度控制过程以优化搜索最优值的方向。

其中，Te₀代表初始温度，当迭代次数逐渐增加时，退火温度逐渐减小。基于退火模拟更新位置步骤如下。

①根据公式（12）计算粒子i的适应值。

②如果新的适应值比先前的适应值好，那么粒子i的位置将被更新，否则跳到步骤③。

③生成一个介于0和1之间的随机数，表示为r。

④用公式（13）计算粒子i的概率。

⑤如果l_i大于r，那么粒子i的位置将被更新，否则将被拒绝。

⑥重复步骤①至⑤，直到所有粒子都被计算一次。

⑦通过公式（14）更新退火温度。

PSMR是在成簇型路由算法的基础上，对上层簇首与sink节点间的数据传输路径进行优化。在开始阶段，簇首竞争及网络结构换分完成后，按照公式（4）生成多条簇首与sink节点间的数据传输链路作为初始搜索粒子，并以网络传输距离、拟消耗传输能量、传输跳数构造适应度函数，筛选符合大规模农田WSNs监测应用中传输能耗低、质量数据传输高的链路，具体的筛选过程如下。

a. 将节点与下一跳候选节点的剩余能量及两点间的传输距离作为约束，构建源节点与sink节点间多条多跳的数据传输链路，作为初始粒子。

b. 对粒子的位置和速度进行初始化，并分别计算出其相应的个体极值与全局极值，并将适应度和当前位置分别存放在代表个体信息的P_i和全局信息的P_g内。

c. 确定当前温度下粒子的适应值P_i，并与全局极值进行比较及替换操作。

d. 通过公式（6）和公式（7）对粒子的速度和位置信息进行更新。

e. 判断目标函数值是否达到最优值或算法达到最大迭代次数，若未达到算法停止要求，则重复步骤b到e；否则输出当前全局极值，并将当前粒子作为源节点与sink节点间通信的数据传输链路。

簇首按照构建的数据传输链路将局部监测数据逐步的向sink节点汇集，并当网络完成一轮数据采集与上传后，重新开始网络结构与数据传输链路的构建过程，当监测区域内的节点出现大面积死亡、网络陷入瘫痪状态或网络运行到预先设定的最大轮次时，停止数据采集与传输工作，并向用户发送监测结束消息。

PSMR的流程如图3所示。

为考察PSMR算法性能，本研究选择WSNs路由算法最常见的网络生命周期、节点剩余能量标准差和网络每轮能量消耗3个指标，在Matlab平台上进行了仿真试验。仿真场景参数如表1所示^［24,25］。在性能对比方面，首先将PSMR算法与PSO和SA算法进行了纵向对比，并进一步将PSMR与EMR、GPSR-A算法进行横向性能对比分析。

在大规模农田WSNs监测中，评价网络性能的核心指标参数为网络生命周期，当节点死亡数目超过总节点数目的30%后，网络出现大面积空洞，采集的数据无法完全代表监测区域的环境情况，因此本研究将网络开始工作到节点死亡数目达到总结点数目30%的这一段时间定义为网络的生命周期。

如前文所述，PSMR采用PSO与SA结合方法对路由路径进行优化，由此本研究首先对比PSMR与PSO、SA的网络生命周期，从图4中可以看出存活节点数在保持一段时间不变后均呈现逐渐下降的趋势。其中SA算法的保持时间最短，约在200轮时起出现首个死亡节点，PSO和PSMR出现首个死亡节点的时间约在500轮。SA算法的生命周期最短为600轮左右，PSO算法的生命周期约在1000轮左右，PSMR算法的生命周期大致相同为1260轮左右。从图4中可以看出，PSO算法收敛速度较快，因此前期性能与PSMR相当，而SA算法收敛速度慢，前期的性能较差，造成较早的节点死亡。PSO算法易陷入局部极值，因此在中后期出现存活节点数量的迅速下降，而SA算法存活节点数量呈缓慢下降趋势。图4的结果也印证了3.2小节的讨论分析。

如图5所示，EMR、GPSR-A和PSMR 3种算法网络生命周期均呈现出保持稳定不变到逐渐下降的趋势。其中EMR的网络生命周期最短，为800轮左右；PSMR和GPSR-A算法的网络生命周期大致相同，为1260轮左右。以网络生命周期作为评价性能依据，PSMR较EMR提升了57%，并且PSMR的第一个节点死亡轮数较GPSR-A算法更晚，在网络生命周期内，节点存活率高于GPSR-A算法。可见，基于PSMR的数据采集质量更优，更好地满足了大规模农田WSNs长时间的监测需求。

对比EMR、GPSR-A和PSMR 3种算法的网络剩余能量标准差变化情况，结果如图6所示。3种算法的剩余能量标准差总体上呈现逐渐上升再下降的趋势。其中PSMR的网络剩余能量标准差最小，平均约为0.025 J；其次分别是EMR和GPSR-A算法，网络剩余能量标准差平均值分别约为0.04和0.065 J。PSMR通过动态划分网络结构的方式起到了均衡簇首和普通节点能耗差异的效果，并且多特征值的数据传输链路优化保证了簇首在数据转发中的低能量损耗，因此PSMR的能量均衡性较好。EMR通过节点剩余能量和通信能耗选择数据转发节点，缩小了节点的能耗差异，但该种方式的路径选择导致每一跳的传输能耗都存在较多差异，且易导致数据转发次数过多；由于离sink节点近的簇首频繁参与数据转发，能量开销较大，因此网络能量标准差有所下降。GPSR-A算法中全网节点都需要构建与sink节点间的数据传输链路，并以一种贪婪选择的方式选取转发节点，该种方式可以有效的缩减数据转发次数，但节点间的数据转发距离差异较大，且平面网络结构下距离sink节点较近的节点因承担着更重的数据转发任务容易早衰，因此网络剩余能量标准差最大。

对比EMR、GPSR-A和PSMR 3种算法的每轮能量消耗情况，结果如图7所示，可以看出PSMR、EMR由于周期性的簇首轮换及构建数据传输链路更新导致每轮能量消耗存在着较大的波动。PSMR考虑了传输跳数与间距，每轮能量消耗较EMR更小。而EMR由于有过多的节点参与了数据转发，因此网络的能量消耗较大。GPSR-A算法每一跳都保持了较长的通信距离，在网络初期节点能耗较大，而随着部分节点死亡和数据传输链路的重新构建，网络轮能量消耗速度随着网络周期轮数的增加呈逐渐降低的趋势。结合图5可以看出，在800轮之后，网络存活节点数量已出现明显下降

通过与EMR、GPSR-A算法的性能对比，可见PSMR的网络生命周期更长、节点剩余能量标准差更低，每轮网络能量消耗与GPSR-A相当，低于EMR，延缓了节点死亡的速度，说明根据分簇的方式动态构建混合型网络结构与多特征值选取数据传输链路的优化方法在提高节点能量利用率、均衡网络负载上效果较好，进而延长了网络生命周期。

本研究分析了传统成簇路由算法在农田WSNs应用中存在的问题，设计了一种动态成簇的混合型网络结构。该网络异构以网络传输距离、拟消耗传输能量和传输跳数为基础构造路由优化函数，并引入PSO和SA的混合算法形成数据传输链路优化策略PSMR。算法仿真试验结果表明，PSMR算法降低了簇首在远距离数据传输中的能量损耗，网络生命周期可达到1260轮，同时实现了0.025 J的网络剩余能量标准差，能够均衡了节点间的能耗差异，延长网络生命周期。

传统路由算法研究多针对10,000 m²及以下的地块面积，本研究中讨论的是90,000 m²及以上的地块大小，方法更加适用于大规模农田应用场景。下一步的研究工作考虑在基础设施薄弱、能量约束严格的大规模农田中部署基于PSMR链路优化策略的WSNs监测网络，探索实际环境中的物联网优化性能。