CN105510081A - Sewage sampling vehicle - Google Patents

Abstract

A sewage sampling vehicle comprises a sampling vehicle body and a navigator mounted on the sampling vehicle body, wherein the navigator specifically comprises a signal module, a processing module and a generating module. An optimized path algorithm is adopted, various cost factors during sampling are considered, the optimization effect is good, the solution efficiency is high, the performance is stable, the global searching capability is enhanced, the sampling running cost can be saved to the largest extent, and a good energy saving effect can be realized.

Description

Sewage sampling vehicle

Technical Field

The invention relates to the field of sewage sampling, in particular to a sewage sampling vehicle.

Background

The industrial sewage is an important pollution source in the modernization, and in order to control the index of sewage discharge, the sampling is required to be carried out on each sampling point regularly.

Since each sampling point of sewage sampling work is scattered in a region with a long distance, the sewage sampling vehicle is an important tool for sewage sampling work. How to select a path which can save the running cost of the sampling vehicle to the maximum extent according to different destinations and the required time for reaching each destination is an urgent problem to be solved.

Disclosure of Invention

In order to solve the problems, the invention provides a sewage sampling vehicle.

The purpose of the invention is realized by adopting the following technical scheme:

a sewage sampling vehicle is used for sampling sewage of a plurality of remote destinations and comprises a sampling vehicle and a navigator arranged on the sampling vehicle, and is characterized in that the navigator specifically comprises a signal module, a processing module and a generating module;

the signal module is used for receiving a plurality of sampling destinations of the current round and predicted required time for reaching each destination, which are input by a user;

the processing module is used for selecting an optimal path according to the sampling destination of the round and the geographical environment information input in advance, and specifically comprises the following steps:

a simulation module:

wherein minS is the lowest cost in the sampling process; m is the total number of the current sampling vehicles; u is the number of destinations; b₀Carbon emission cost per unit distance; omega₀Is the carbon emission coefficient; phi (₀Fuel consumption per unit distance when no load exists; f. of_ijDistance between destination i (i ═ 1,2, …, U) to destination j (j ═ 1,2, …, U); c is the load capacity of the sampling vehicle; h is the maximum load capacity of the sampling vehicle; phi (^*Fuel consumption per unit distance when fully loaded;

T₁in order for the sampling vehicle to arrive at the loss factor in advance,for cost penalty in arriving at destination i earlier at time G, T₂In order to sample the vehicle late loss factor,for cost loss when the destination i is reached with a delay to the time O, an early arrival loss coefficient and a late arrival loss coefficient are used to take the reference point condition, T, of the sample car to each destination₁And T₂A coefficient set artificially;

an opportunity module: assuming a total of R nodes, γ_ij(t) represents the intensity of the tracker between node i and node j at time t, γ_ij(0) If the sampling vehicle selects a transfer direction during movement according to the intensity of the tracker, the probability that the sampling vehicle K (K is 1, 2.. multidot.m) is transferred from the node i to the node j is as follows:

wherein, g ∈ A_k；A_k＝{0,1,...,R-1}-B_kRepresenting the set of points the sampling vehicle k is allowed to select next, dynamically changing over time, B_k(k is 1,2,.. multidot.m) is a taboo table of the kth sampling vehicle and is used for recording points sampled by the sampling vehicle k;the heuristic factor represents the expected degree of the t time from the node i to the node j, and is generally takenψ is an information heuristic factor, μ is an expected heuristic factor, α (i, j) is a time degree of the next destination;relative importance of time degree for next destination;

an update module: introducing an optimization variable X_ij(t) satisfying X_ij(t+1)＝σX(t)[1-X_ij(t)]And obtaining an optimized tracker updating rule by taking sigma as a control variable:

γ_ij(t+1)＝(1-ζ)γ_ij(t)+Δγ_ij(t)+чX_ij(t)

wherein,

F_kis the length of the path taken by the kth sampling vehicle in the current cycle, I is a constant of the intensity of the tracker,showing the intensity of the tracker left on the path (i, j) by the kth sampling vehicle in the cycle, wherein zeta is the global volatility factor of the tracker and is zeta ∈ [0, 1%]And ζ is a parameter that is dynamically adjusted according to the following formula:wherein ζ_minIs a minimum value set manually; delta gamma_ij(t) represents the sum of the intensities of the tracers left on the path (i, j) by all the sampling vehicles in the cycle, ч is an adjustable coefficient;

an initial module: initializing parameters and adjusting tracking elements of each path by setting the iteration number DD to be 0; produce a range of [0,1]If p < given constant p₀The next node j is selected according to the following equation: wherein l ∈ A_k(ii) a Otherwise, selecting the next node j according to the probability formula in the opportunity module, and adding the j into the array B_kRepeating the steps until all the node tasks are completed to obtain an initial set S of the simulation algorithm_i；

An optimal solution module: generating a new set of feasible solutions S from the current initial set_jThe target value variation Δ S is S_j-S_iIf Δ S < 0, then accept the new feasible solution S_jIs the optimal solution; otherwise the effect of the deviation is taken into account: r ═ exp (- Δ S/N (t)), where N is the amount that changes over time, and S is accepted if r > 1_jIs the optimal solution, otherwise does not accept the new feasible solution, and the optimal solution is still S_i；

A judging module: after finding out the optimal solution, judging whether the new path has an overload phenomenon, if so, regenerating a feasible solution, and if not, accepting the new feasible solution as the optimal solution; when the current optimal solution is smaller than a certain specific value, updating the tracker; if the current round lists B_kIf there is no data update, then a [0,1 ] is generated]Random number of range u, if e₁+e₂+,...,e_i-1＜u＜e₁+e₂+,...,e_iThen the probability of selection is e_iThe candidate sampling vehicle is used as the next target node;

a generation module: outputting the calculated optimal path by making the iteration number DD equal to DD +1, and if DD is less than DD_maxClearing B according to a formula N (t +1) ═ N (t) v according to a tracker updating rule_kList, where v ∈ [0,1]Returning to the initial module to regenerate the random number p; if DD is DD ═ DD_maxAnd outputting the optimal solution as the optimal path.

Drawings

The invention is further illustrated by means of the attached drawings, but the embodiments in the drawings do not constitute any limitation to the invention, and for a person skilled in the art, other drawings can be obtained on the basis of the following drawings without inventive effort.

Fig. 1 is a block diagram of the present invention.

Reference numerals: a signal module-1; an analog module-3; opportunity module-5; an update module-7; an initial module-9; an optimal solution module-11; a judgment module-13; and generating a module-15.

Detailed Description

The invention is further described with reference to the following examples.

The sewage sampling vehicle shown in fig. 1 is used for sampling sewage at a plurality of remote destinations, and comprises a sampling vehicle and a navigator installed on the sampling vehicle, wherein the navigator specifically comprises a signal module 1, a processing module and a generating module 15;

the system comprises a signal module 1, a data processing module and a data processing module, wherein the signal module is used for receiving a plurality of sampling destinations of the current round and predicted required time for reaching each destination, which are input by a user;

the processing module is used for selecting an optimal path according to the sampling destination of the round and the geographical environment information input in advance, and specifically comprises the following steps:

and (3) an analog module:

wherein minS is the lowest cost in the sampling process; m is the total number of the current sampling vehicles; u is the number of destinations; b₀Carbon emission cost per unit distance; omega₀Is the carbon emission coefficient; phi (₀Fuel consumption per unit distance when no load exists; f. of_ijDistance between destination i (i ═ 1,2, …, U) to destination j (j ═ 1,2, …, U); c is the load capacity of the sampling vehicle; h is the maximum load capacity of the sampling vehicle; phi (^*Fuel consumption per unit distance when fully loaded;

T₁in order for the sampling vehicle to arrive at the loss factor in advance,for cost penalty in arriving at destination i earlier at time G, T₂In order to sample the vehicle late loss factor,for cost loss when the destination i is reached with a delay to the time O, an early arrival loss coefficient and a late arrival loss coefficient are used to take the reference point condition, T, of the sample car to each destination₁And T₂A coefficient set artificially;

the opportunity module 5: assuming a total of R nodes, γ_ij(t) represents the intensity of the tracker between node i and node j at time t, γ_ij(0) If the sampling vehicle selects a transfer direction during movement according to the intensity of the tracker, the probability that the sampling vehicle K (K is 1, 2.. multidot.m) is transferred from the node i to the node j is as follows:

wherein, g ∈ A_k；A_k＝{0,1,...,R-1}-B_kRepresenting the set of points the sampling vehicle k is allowed to select next, dynamically changing over time, B_k(k is 1,2,.. multidot.m) is a taboo table of the kth sampling vehicle and is used for recording points sampled by the sampling vehicle k;the heuristic factor represents the expected degree of the t time from the node i to the node j, and is generally takenψ is an information heuristic factor, μ is an expected heuristic factor, α (i, j) is a time degree of the next destination;relative importance of time degree for next destination;

the updating module 7: introducing an optimization variable X_ij(t) satisfying X_ij(t+1)＝σX(t)[1-X_ij(t)]And obtaining an optimized tracker updating rule by taking sigma as a control variable:

γ_ij(t+1)＝(1-ζ)γ_ij(t)+Δγ_ij(t)+чX_ij(t)

wherein,

F_kis the length of the path taken by the kth sampling vehicle in the current cycle, I is a constant of the intensity of the tracker,showing the intensity of the tracker left on the path (i, j) by the kth sampling vehicle in the cycle, wherein zeta is the global volatility factor of the tracker and is zeta ∈ [0, 1%]And ζ is a parameter that is dynamically adjusted according to the following formula:wherein ζ_minIs a minimum value set manually; delta gamma_ij(t) represents the sum of the intensities of the tracers left on the path (i, j) by all the sampling vehicles in the cycle, ч is an adjustable coefficient;

an initial module 9: initializing parameters and adjusting tracking elements of each path by setting the iteration number DD to be 0; produce a range of [0,1]If p < given constant p₀The next node j is selected according to the following equation: wherein l ∈ A_k(ii) a Otherwise, selecting the next node j according to the probability formula in the opportunity module 5, and adding j into the node jArray B_kRepeating the steps until all the node tasks are completed to obtain an initial set S of the simulation algorithm_i；

The optimal solution module 11: generating a new set of feasible solutions S from the current initial set_jThe target value variation Δ S is S_j-S_iIf Δ S < 0, then accept the new feasible solution S_jIs the optimal solution; otherwise the effect of the deviation is taken into account: r ═ exp (- Δ S/N (t)), where N is the amount that changes over time, and S is accepted if r > 1_jIs the optimal solution, otherwise does not accept the new feasible solution, and the optimal solution is still S_i；

The judging module 13: after finding out the optimal solution, judging whether the new path has an overload phenomenon, if so, regenerating a feasible solution, and if not, accepting the new feasible solution as the optimal solution; when the current optimal solution is smaller than a certain specific value, updating the tracker; if the current round lists B_kIf there is no data update, then a [0,1 ] is generated]Random number of range u, if e₁+e₂+,...,e_i-1＜u＜e₁+e₂+,...,e_iThen the probability of selection is e_iThe candidate sampling vehicle is used as the next target node;

the generation module 15: outputting the calculated optimal path by making the iteration number DD equal to DD +1, and if DD is less than DD_maxClearing B according to a formula N (t +1) ═ N (t) v according to a tracker updating rule_kList, where v ∈ [0,1]Returning to the initial module 9, and regenerating the random number p; if DD is DD ═ DD_maxAnd outputting the optimal solution as the optimal path.

The invention adopts an optimized path algorithm, considers various cost factors in the sampling process, has good optimization effect, high solving efficiency and stable performance, enhances the global searching capability, can save the sampling operation cost to the maximum extent and can play a good energy-saving effect.

Finally, it should be noted that the above embodiments are only used for illustrating the technical solutions of the present invention, and not for limiting the protection scope of the present invention, although the present invention is described in detail with reference to the preferred embodiments, it should be understood by those skilled in the art that modifications or equivalent substitutions can be made on the technical solutions of the present invention without departing from the spirit and scope of the technical solutions of the present invention.

Claims (1)

1. A sewage sampling vehicle is used for sampling sewage of a plurality of remote destinations and comprises a sampling vehicle and a navigator arranged on the sampling vehicle, and is characterized in that the navigator specifically comprises a signal module, a processing module and a generating module;

the signal module is used for receiving a plurality of sampling destinations of the current round and predicted required time for reaching each destination, which are input by a user;

the processing module is used for selecting an optimal path according to the sampling destination of the round and the geographical environment information input in advance, and specifically comprises the following steps:

a simulation module:

$\begin{array}{c}\min S=\underset{m=1}{\overset{m}{\&Sigma;}}\underset{i=0}{\overset{U}{\&Sigma;}}\underset{i=0}{\overset{U}{\&Sigma;}}{b}_0{\mathrm{\&omega;}}_0{\mathrm{\&Phi;}}_0{f}_{ij}{y}_{ijk}+\underset{m=1}{\overset{m}{\&Sigma;}}\underset{i=0}{\overset{U}{\&Sigma;}}\underset{i=0}{\overset{U}{\&Sigma;}}{b}_0{\mathrm{\&omega;}}_0\frac{{\mathrm{\&Phi;}}^{*}-{\mathrm{\&Phi;}}_0}{H}{c}_i{f}_{ij}{y}_{ijk}+T_1\underset{i=0}{\overset{U}{\&Sigma;}}\left(G_i-t_i\right)\\ +T_2\underset{i=0}{\overset{U}{\&Sigma;}}\left(t_i-O_i\right)\end{array}$

wherein minS is the lowest cost in the sampling process; m is the total number of the current sampling vehicles; u is the number of destinations; b₀Carbon emission cost per unit distance; omega₀Is the carbon emission coefficient; phi (₀Fuel consumption per unit distance when no load exists; f. of_ijDistance between destination i (i ═ 1,2, …, U) to destination j (j ═ 1,2, …, U); c is the load capacity of the sampling vehicle; h is the maximum load capacity of the sampling vehicle; phi (^*Fuel consumption per unit distance when fully loaded;

T₁in order for the sampling vehicle to arrive at the loss factor in advance,for cost penalty in arriving at destination i earlier at time G, T₂In order to sample the vehicle late loss factor,for cost loss when the destination i is reached with a delay to the time O, an early arrival loss coefficient and a late arrival loss coefficient are used to take the reference point condition, T, of the sample car to each destination₁And T₂A coefficient set artificially;

an opportunity module: assuming a total of R nodes, γ_ij(t) represents the intensity of the tracker between node i and node j at time t, γ_ij(0) If the sampling vehicle selects a transfer direction during movement according to the intensity of the tracker, the probability that the sampling vehicle K (K is 1, 2.. multidot.m) is transferred from the node i to the node j is as follows:

wherein, g ∈ A_k；A_k＝{0，1，...，R-1}-B_kRepresenting the set of points the sampling vehicle k is allowed to select next, dynamically changing over time, B_k(k is 1,2,.. multidot.m) is a taboo table of the kth sampling vehicle and is used for recording points sampled by the sampling vehicle k;the heuristic factor represents the expected degree of the t time from the node i to the node j, and is generally takenψ is an information heuristic factor, μ is an expected heuristic factor, α (i, j) is a time degree of the next destination;relative importance of time degree for next destination;

an update module: introducing an optimization variable X_ij(t) satisfying X_ij(t+1)＝σX(t)[1-X_ij(t)]And obtaining an optimized tracker updating rule by taking sigma as a control variable:

γ_ij(t+1)＝(1-ζ)γ_ij(t)+Δγ_ij(t)+чX_ij(t)

wherein, ${\mathrm{\&Delta;\&gamma;}}_{ij}\left(t\right)={\&Sigma;}_{k=1}^m{\mathrm{\&Delta;\&gamma;}}_{ij}^k\left(t\right),$

F_kis the length of the path taken by the kth sampling vehicle in the current cycle, I is a constant of the intensity of the tracker,showing the intensity of the tracker left on the path (i, j) by the kth sampling vehicle in the cycle, wherein zeta is the global volatility factor of the tracker and is zeta ∈ [0, 1%]And ζ is a parameter that is dynamically adjusted according to the following formula:wherein ζ_minIs a minimum value set manually; delta gamma_ij(t) represents the sum of the intensities of the tracers left on the path (i, j) by all the sampling vehicles in the cycle, ч is an adjustable coefficient;

an initial module: initializing parameters and adjusting tracking elements of each path by setting the iteration number DD to be 0; produce a range of [0,1]If p < given constant p₀The next node j is selected according to the following equation: wherein l ∈ A_k(ii) a Otherwise, selecting the next node j according to the probability formula in the opportunity module, and adding the j into the array B_kRepeating the steps until all the node tasks are completed to obtain an initial set S of the simulation algorithm_i；

An optimal solution module: generating a new set of feasible solutions S from the current initial set_jThe target value variation Δ S is S_j-S_iIf Δ S < 0, then accept the new feasible solution S_jIs the optimal solution; otherwise the effect of the deviation is taken into account: r ═ exp (- Δ S/N (t)), where N is the amount that changes over time, and S is accepted if r > 1_jIs the optimal solution, otherwise does not accept the new feasible solution, and the optimal solution is still S_i；

A judging module: after finding out the optimal solution, judging whether the new path has an overload phenomenon, if so, regenerating a feasible solution, and if not, accepting the new feasible solution as the optimal solution; when the current optimal solution is smaller than a certain specific value, updating the tracker; if the current round lists B_kIf there is no data update, then a [0,1 ] is generated]Random number of range u, if e₁+e₂+，...，e_i-1＜u＜e₁+e₂+，...，e_iThen the probability of selection is e_iThe candidate sampling vehicle is used as the next target node;

a generation module: outputting the calculated optimal path by making the iteration number DD equal to DD +1, and if DD is less than DD_maxClearing B according to a formula N (t +1) ═ N (t) v according to a tracker updating rule_kList of υ ∈ [0,1 [ ]]Returning to the initial module to regenerate the random number p; if DD is DD ═ DD_maxAnd outputting the optimal solution as the optimal path.