CN110110276A

CN110110276A - Leak position method and apparatus based on variable step recursion track

Info

Publication number: CN110110276A
Application number: CN201910202199.9A
Authority: CN
Inventors: 郑小平; 李静海; 陈茂银; 张建文; 李辉; 耿华; 邓晓娇; 朱亚征
Original assignee: Nanjing Commoner Mdt Infotech Ltd; Tsinghua University; Beijing University of Chemical Technology
Current assignee: Nanjing Commoner Mdt Infotech Ltd; Tsinghua University; Beijing University of Chemical Technology
Priority date: 2019-03-18
Filing date: 2019-03-18
Publication date: 2019-08-09
Anticipated expiration: 2039-03-18
Also published as: CN110110276B

Abstract

This application involves a kind of leak position method and apparatus based on variable step recursion track.It include: the target function model according to preset model of atmospheric diffusion building and gas concentration information association；The gas concentration information of monitoring current location obtains the target function model calculated value of current location according to gas concentration information, model of atmospheric diffusion and the target function model of current location；The gas concentration information of monitoring the next position obtains the target function model calculated value of the next position according to the gas concentration information, model of atmospheric diffusion and target function model of the next position；According to the target function model calculated value of the target function model calculated value of current location and the next position, leakage source strength and leakage source position are determined.The above method can obtain leakage source strength and leakage source position according to target function model calculated value, and realization is accurately positioned dangerization gas leakage source.

Description

Leak position method and apparatus based on variable step recursion track

Technical field

This application involves field of locating technology, in particular to a kind of leak position method based on variable step recursion track And device.

Background technique

Hazardous chemical often has diffusivity and toxic, once leak, will cause serious property loss and Casualties.After the generation of harmful influence leakage accident, harmful influence source of leaks is positioned at the first time, is at accident emergency Set, emergency management and rescue and assess accident potential hazard basis.After harmful influence leakage accident occurs, it can be leaked according to harmful influence Gas leakage situation in accident positions harmful influence source of leaks.

Traditional dangerization gas leakage source localization method is the concentration gradient and dispersal direction according to dangerization gas in region Dangerization gas leakage source is positioned, leakage accident site environment can be usually idealized.However, past in harmful influence leakage accident Toward the generation along with explosion, burning and complicated air-flow, these explosions, burning and the generation of complicated air-flow are so that site environment is more Unsatisfactoryization in number situation, the site environment of non-ideal cause that dangerization gas leakage source can not be accurately positioned.

Summary of the invention

Based on this, it is necessary in view of the above technical problems, provide a kind of leak position based on variable step recursion track Method and apparatus.

A kind of leak position method based on variable step recursion track, which comprises

According to the target function model of preset model of atmospheric diffusion building and gas concentration information association；

The gas concentration information for monitoring current location expands according to the gas concentration information of the current location, the atmosphere Model and the target function model are dissipated, the target function model calculated value of the current location is obtained；

The gas concentration information for monitoring the next position expands according to the gas concentration information of described the next position, the atmosphere Model and the target function model are dissipated, the target function model calculated value of described the next position is obtained；

It is calculated according to the target function model of the target function model calculated value of the current location and described the next position Value determines leakage source strength and leakage source position.

In one of the embodiments, the gas concentration information of the current location include monitoring time and with the monitoring Time corresponding gas concentration position；

The gas concentration information of the monitoring current location, according to the gas concentration information of the current location, described big Gas diffusion model and the target function model obtain the target function model calculated value of the current location, comprising:

The gas concentration information for monitoring the current location, according to the gas concentration information of the current location and described big Gas diffusion model, the model of atmospheric diffusion for obtaining the current location calculate concentration sequence；

Expanded according to the atmosphere of the target function model, the gas concentration information of the current location and the current location It dissipates model and calculates concentration sequence, obtain the target function model calculated value of the current location.

The gas concentration information of the monitoring current location in one of the embodiments, according to the present bit The gas concentration information and the model of atmospheric diffusion set, the model of atmospheric diffusion for obtaining the current location calculate concentration sequence Column, comprising:

The model parameter relevant to the model of atmospheric diffusion that will acquire, the monitoring time and when with the monitoring Between corresponding gas concentration position input the model of atmospheric diffusion, obtain the current location model of atmospheric diffusion calculate it is dense Degree series.

The model parameter relevant to the model of atmospheric diffusion that will acquire, institute in one of the embodiments, It states monitoring time and gas concentration corresponding with monitoring time position inputs the model of atmospheric diffusion, before further include:

According to the gas of the current location in the horizontal direction, the diffusion coefficient of vertical direction and with the monitoring time Corresponding wind speed information obtains model parameter relevant to the model of atmospheric diffusion of the current location.

The gas concentration information of the current location includes monitoring time and the monitoring in one of the embodiments, Time corresponding gas concentration size and gas concentration corresponding with monitoring time position；

It is described according to the big of the target function model, the gas concentration information of the current location and the current location Gas diffusion model calculates concentration sequence, obtains the target function model calculated value of the current location, comprising:

Expanded according to the atmosphere of the target function model, the gas concentration information of the current location and the current location It dissipates model and calculates concentration sequence, obtain first object function and adapt to degree series；Wherein, the first object function adapts to degree series With and there are first the corresponding gas concentration size of the monitoring time and gas concentration corresponding with monitoring time position Corresponding relationship；

Degree series are adapted to the first object function to be ranked up and screen, and obtain the second objective function fitness sequence Column；

Degree series are adapted to second objective function and execute preset times genetic algorithm, obtain intermediate population's fitness sequence Column；Wherein, there are the second corresponding relationships for the intermediate population's fitness sequence and second objective function adaptation degree series；

Local optimal searching is carried out to the intermediate population's fitness sequence, obtains optimal intermediate population's fitness；

According to first corresponding relationship, second corresponding relationship, the optimal intermediate population's fitness and the mesh Offer of tender exponential model obtains the target function model calculated value of the current location.

It is described in one of the embodiments, to be believed according to the target function model, the gas concentration of the current location The model of atmospheric diffusion of breath and the current location calculates concentration sequence, obtains first object function and adapts to degree series, comprising:

According to and the corresponding gas concentration size of the monitoring time and gas concentration corresponding with monitoring time position It sets, determines starting leakage source strength collection and starting leakage source position collection；

Concentration sequence is calculated according to the starting leakage source strength collection and the model of atmospheric diffusion, obtains initial atmosphere expansion It dissipates model and calculates concentration collection；

According to the starting leakage source strength collection, starting leakage source position collection, the current location gas concentration Information and the initial atmosphere diffusion model calculate concentration collection, obtain the first object function and adapt to degree series.

The target function model calculated value according to the current location and described next in one of the embodiments, The target function model calculated value of position determines leakage source strength and leakage source position, comprising:

Gas concentration position corresponding with the target function model calculated value of the current location is determined as the first leakage Gas concentration position corresponding with the target function model calculated value of described the next position is determined as the second source of leaks by source position Position；

Determine the distance between the first leakage source position and second leakage source position value；

It, will gas corresponding with the second leakage source position if the distance value is less than or equal to the first preset threshold Concentration is determined as the leakage source strength, and second leakage source position is determined as the leakage source position.

In one of the embodiments, the method also includes:

If the distance value is greater than first preset threshold, using described the next position as updated present bit It sets, and returns to the gas concentration information for executing described monitoring the next position, according to the gas concentration information of described the next position, institute The step of stating model of atmospheric diffusion and the target function model, obtaining the target function model calculated value of described the next position.

It is described by gas corresponding with the target function model calculated value of the current location in one of the embodiments, Concentration position is determined as the first leakage source position, and gas corresponding with the target function model calculated value of described the next position is dense Degree position is determined as the second leakage source position, comprising:

If the target function model calculated value of the current location is less than or equal to the second preset threshold, will work as with described The corresponding gas concentration position of the target function model calculated value of front position is determined as first leakage source position；

If the target function model calculated value of described the next position is less than or equal to second preset threshold, will be with institute The corresponding gas concentration position of target function model calculated value for stating the next position is determined as second leakage source position.

A kind of leak position device based on variable step recursion track, described device include:

Module is constructed, for the objective function mould according to preset model of atmospheric diffusion building and gas concentration information association Type；

First determining module, for monitoring the gas concentration information of current location, the gas according to the current location is dense Information, the model of atmospheric diffusion and the target function model are spent, the target function model for obtaining the current location calculates Value；

Second determining module, for monitoring the gas concentration information of the next position, the gas according to described the next position is dense Information, the model of atmospheric diffusion and the target function model are spent, the target function model for obtaining described the next position calculates Value；

Locating module, for according to the target function model calculated value of the current location and the target of described the next position Function model calculated value determines leakage source strength and leakage source position.

A kind of computer equipment, including memory and processor, the memory are stored with computer program, the processing The step of device realizes any of the above-described the method when executing the computer program.

A kind of computer readable storage medium, is stored thereon with computer program, and the computer program is held by processor The step of any of the above-described the method is realized when row.

The above-mentioned leak position method and apparatus based on variable step recursion track, by according to preset Atmospheric Diffusion mould The target function model of type building and gas concentration information association, monitors the gas concentration information of current location, and according to current Gas concentration information, model of atmospheric diffusion and the target function model of position, the target function model for obtaining current location calculate Value, monitors the gas concentration information of the next position, and according to the gas concentration information, model of atmospheric diffusion and target of the next position Function model obtains the target function model calculated value of the next position, further according to the target function model meter of current location The positioning of source position is realized to leakage source strength and leaked to the target function model calculated value of calculation value and the next position.Above-mentioned side Method can be leaked according to the target function model calculated value of current location and the target function model calculated value of the next position Source strength and leakage source position, realization are accurately positioned dangerization gas leakage source.

Detailed description of the invention

Fig. 1 is the flow diagram of the leak position method based on variable step recursion track in one embodiment；

Fig. 2 is the flow diagram of the refinement step of step S200 in one embodiment；

Fig. 3 is the flow diagram of the refinement step of step S220 in one embodiment；

Fig. 4 is the flow diagram of the refinement step of step S221 in one embodiment；

Fig. 5 is the structural block diagram of the leak position device based on variable step recursion track in one embodiment；

Fig. 6 is the internal structure chart of computer equipment in one embodiment.

Specific embodiment

It is with reference to the accompanying drawings and embodiments, right in order to which the objects, technical solutions and advantages of the application are more clearly understood The application is further elaborated.It should be appreciated that specific embodiment described herein is only used to explain the application, not For limiting the application.

It is understood that term " first " used in this application, " second " etc. can be used to describe herein various sequences Relationship and corresponding relationship, but these sequence relations and corresponding relationship should not be limited by these terms.These terms are only used to by one Sequence relation and corresponding relationship are distinguished with another sequence relation and corresponding relationship.

In one embodiment, as shown in Figure 1, providing a kind of leak position side based on variable step recursion track Method, comprising the following steps:

Step S100, according to the target function model of preset model of atmospheric diffusion building and gas concentration information association.

In embodiments of the present invention, model of atmospheric diffusion is the gas diffusion migration shape under a kind of calculating simulation actual conditions The mathematical model of condition.Wherein, model of atmospheric diffusion can be selected according to the ambient condition locating for source of leaks, for example, being directed to The diffusion model of special weather condition and landform, enclosed type diffusion model, Hypatmism type diffusion model, mountain area model of atmospheric diffusion and Coastal atmospheric dispersion model type, is not specifically limited herein.

Wherein, gas concentration information refers to concentration of the gas released at harmful influence source of leaks in source of leaks surrounding space Information, gas concentration information herein are a kind of identifiers, are not the data of actual measurement.

Step S200 monitors the gas concentration information of current location, is expanded according to the gas concentration information of current location, atmosphere Model and target function model are dissipated, the target function model calculated value of current location is obtained.

Wherein, the gas concentration information of current location includes the monitoring time monitored, gas corresponding with monitoring time The information such as concentration and gas concentration corresponding with monitoring time position.Wherein, monitoring time, gas corresponding with monitoring time There is one-to-one relationships for bulk concentration size and gas concentration corresponding with monitoring time position.Optionally, with monitoring time Corresponding gas concentration size is monitored by ionic migration spectrometer, and gas concentration corresponding with monitoring time position is by fixed Position device navigates to.Ionic migration spectrometer monitoring is continually monitored gas concentration size in its surrounding space, when ion moves When shifting spectrometer captures the concentration information of leakage source gas for the first time, then the position that will capture leakage source gas for the first time is regarded as currently Position.

In embodiments of the present invention, target function model calculated value refers to the gas concentration according to the current location monitored The value for the target function model that information is calculated.Optionally, the gas concentration information and model of atmospheric diffusion of front position is defeated Enter the available target function model calculated value of target function model.

Step S300 monitors the gas concentration information of the next position, is expanded according to the gas concentration information of the next position, atmosphere Model and target function model are dissipated, the target function model calculated value of the next position is obtained.

Wherein, the gas concentration information of the next position and the gas concentration information of the current location in step S200 are monitored Object it is identical, the gas concentration information of the next position be also included within the monitoring time monitored at the next position, with monitoring when Between corresponding gas concentration size and gas concentration corresponding with monitoring time position.Wherein, monitoring time and monitoring time pair There is one-to-one relationships for the gas concentration size answered and gas concentration corresponding with monitoring time position.Similarly, next The calculation phase of the target function model calculated value of position and the target function model calculated value of current location in step S200 Together, details are not described herein.

Step S400 is calculated according to the target function model of the target function model calculated value of current location and the next position Value determines leakage source strength and leakage source position.

It specifically, can be according to the target function model calculated value of finally determining current location and the target of the next position Function value is back-calculated to obtain objective function calculated value corresponding with the objective function calculated value of current location and with the next position Corresponding gas concentration size and gas concentration position.Further, according to corresponding with the objective function calculated value of current location The distance between gas concentration position and gas concentration position corresponding with the objective function calculated value of the next position relationship determine Source position is leaked out, and the gas concentration at the leakage source position is sized to leakage source strength.

In above-described embodiment, target function model be may be expressed as:

Wherein, L₀Indicate leakage source position to be asked, Q₀Indicate gas leakage source intensity to be asked, t_iWhen indicating i-th It carves, L_iThe expression time is t_iWhen monitoring point position L_i=(x_i, y_i), C_Mol(L_i, t_i) indicate in time t_iWhen, position (x_i, y_i) The model of atmospheric diffusion at place calculates concentration, C_IMS(L_i, t_i) indicate in time t_iWhen, position (x_i, y_i) at ionic migration spectrometer monitoring The gas concentration size arrived.Optionally, when being positioned using global positioning system (Global Position System, GPS) (x_i, y_i) indicate longitude and latitude of the monitoring point under global coordinates system.

The above-mentioned leak position method based on variable step recursion track, by being constructed according to preset model of atmospheric diffusion With the target function model of gas concentration information association, the gas concentration information of current location is monitored, and according to current location Gas concentration information, model of atmospheric diffusion and target function model obtain the target function model calculated value of current location, monitoring The gas concentration information of the next position, and according to the gas concentration information, model of atmospheric diffusion and objective function mould of the next position Type obtains the target function model calculated value of the next position, further according to the target function model calculated value of current location and The positioning of source position is realized to leakage source strength and leaked to the target function model calculated value of the next position.The above method can Leakage source strength is obtained according to the target function model calculated value of the target function model calculated value of current location and the next position With leakage source position, realization is accurately positioned dangerization gas leakage source.

In one of the embodiments, as shown in Fig. 2, the flow diagram of the refinement step for step S200, wherein when The gas concentration information of front position includes monitoring time and gas concentration corresponding with monitoring time position；Step S200, monitoring The gas concentration information of current location, according to gas concentration information, model of atmospheric diffusion and the target function model of current location, The target function model calculated value of current location is obtained, specifically includes the following steps:

Step S210 monitors the gas concentration information of current location, according to the gas concentration information and atmosphere of current location Diffusion model, the model of atmospheric diffusion for obtaining current location calculate concentration sequence.

Wherein, model of atmospheric diffusion calculates concentration sequence and refers to according to the gas concentration information of known current location and big The sequence about model of atmospheric diffusion that gas diffusion model is calculated.Optionally, the gas concentration information of current location is defeated Enter model of atmospheric diffusion, available model of atmospheric diffusion calculates concentration sequence.

Step S220, according to target function model, the gas concentration information of current location and the Atmospheric Diffusion of current location Model calculates concentration sequence, obtains the target function model calculated value of current location.

In embodiments of the present invention, by the gas concentration information of the current location monitored and the Atmospheric Diffusion of current location Model calculates concentration sequence inputting target function model, and what is obtained is the sequence comprising several function model calculated values, right The sequence is ranked up and searches for, available one optimal function model calculated value, by this optimal function model meter Calculation value is as above-mentioned target function model calculated value.

In above-described embodiment, it is one of the committed step for positioning source of leaks that target function model calculated value, which is calculated, is obtained Target function model calculated value is taken to provide basis for positioning source of leaks.

Step S210 in one of the embodiments, monitors the gas concentration information of current location, according to current location Gas concentration information and model of atmospheric diffusion, the model of atmospheric diffusion for obtaining current location calculate concentration sequence, comprising:

The relevant model parameter of will acquire and model of atmospheric diffusion, monitoring time and gas corresponding with monitoring time Concentration position inputs model of atmospheric diffusion, and the model of atmospheric diffusion for obtaining current location calculates concentration sequence.

Optionally, model of atmospheric diffusion can roll into a ball model for Gauss cigarette.Gauss cigarette group model may be expressed as:

Wherein, C (x, y, t) indicates that the Atmospheric models when time is t at position (x, y) calculate concentration, and v is wind speed, σ_x、σ_y、 And σ_zRespectively gas x (trunnion axis), y (trunnion axis), z (vertical direction) direction diffusion coefficient.Wherein wind speed v includes wind Speed and the information such as direction, measured by ultrasonic wind velocity indicator.

In above-described embodiment, selection Gauss cigarette group's model, will model relevant to Gauss cigarette group's model as objective for implementation Parameter (for example, gas in the horizontal direction, the diffusion coefficient of vertical direction and wind speed information corresponding with monitoring time etc.), prison It surveys the time and gas concentration corresponding with monitoring time position input Gauss cigarette rolls into a ball model, obtain Gauss cigarette group model and calculate concentration Sequence, it only includes Q in concentration sequence that Gauss cigarette group model at this time, which calculates,₀One unknown parameter, wherein will only include unknown ginseng Measure Q₀Gauss cigarette group model calculate concentration sequence as model of atmospheric diffusion calculating concentration sequence, model of atmospheric diffusion is calculated Concentration sequence is as the leakage source strength determined in objective function and the basis of leakage source position.

The model parameter relevant to model of atmospheric diffusion that will acquire in one of the embodiments, monitoring time and Gas concentration corresponding with monitoring time position inputs model of atmospheric diffusion, before further include:

According to the gas of current location in the horizontal direction, the diffusion coefficient of vertical direction and wind corresponding with monitoring time Fast information obtains model parameter relevant to the model of atmospheric diffusion of current location.Wherein, wind speed letter corresponding with monitoring time It ceases and refers under the currently monitored time, the speed of the wind at gas concentration corresponding with monitoring time position and direction.

In above-described embodiment, model of atmospheric diffusion and current gas in the horizontal direction with the diffusion Chengdu phase of vertical direction It closes, and related to monitoring point wind direction at that time and wind speed, therefore, in calculating process, needs first to obtain and Atmospheric Diffusion mould The relevant model parameter of type.

In one of the embodiments, as shown in figure 3, the flow diagram of the refinement step for step S220, wherein when The gas concentration information of front position includes monitoring time, gas concentration size corresponding with monitoring time and corresponding with monitoring time Gas concentration position；Step S220, according to the big of target function model, the gas concentration information of current location and current location Gas diffusion model calculates concentration sequence, obtains the target function model calculated value of current location, specifically includes the following steps:

Step S221, according to target function model, the gas concentration information of current location and the Atmospheric Diffusion of current location Model calculates concentration sequence, obtains first object function and adapts to degree series.

Wherein, first object function adapts to degree series with gas concentration size corresponding with monitoring time and and monitoring time There are the first corresponding relationships for corresponding gas concentration position.

Specifically, the gas concentration information of the current location monitored and the model of atmospheric diffusion of current location are calculated dense Degree series input target function model, and available first object function adapts to degree series, therefore first object function fitness Sequence with current location gas concentration information (monitoring time, gas concentration size corresponding with monitoring time and with monitoring when Between corresponding gas concentration position) between there is an one-to-one relationships, this one-to-one relationship can be claimed Make the first corresponding relationship.For example, the gas concentration information sequence of current location can be expressed as A_a1、A_a2......A_an, the first mesh Scalar functions, which adapt to degree series, can be expressed as T_t1、T_t2......T_tn, wherein an is equal to tn, A_a1With T_t1It is corresponding, A_a2With T_t2It is right It answers, A_anWith T_tnIt is corresponding.

Step S222 adapts to degree series to first object function and is ranked up and screens, and obtains the adaptation of the second objective function Degree series.

First object function adaptation degree series are ranked up, after rejecting extreme value, obtain the second objective function fitness sequence Column, it is corresponding that the second objective function at this time adapts to the element that degree series are adapted to first object function in degree series.Therefore, Two objective functions adapt to degree series, and with the gas concentration information of current location, (monitoring time, gas corresponding with monitoring time are dense Spend size and gas concentration corresponding with monitoring time position) between there is also a corresponding relationships.For example, obtain Two objective functions, which adapt to degree series, can be expressed as T_tt1、T_tt2......T_ttn, wherein ttn is less than or equal to tn, T_tt1、 T_tt2......T_ttnIn each element with first object function adapt to degree series T_t1、T_t2......T_tnIn only one Element is corresponding, and first object function adapts to degree series T_t1、T_t2......T_tnIn each element also gas with current location Bulk concentration information sequence A_a1、A_a2......A_anIn only one element it is corresponding, therefore, T_tt1、T_tt2......T_ttnIn Gas concentration information sequence A of each element with current location_a1、A_a2......A_anIn only one element it is corresponding.

Step S223 adapts to degree series to the second objective function and executes preset times genetic algorithm, it is suitable to obtain intermediate population Answer degree series.

Wherein, there are the second corresponding relationships for intermediate population's fitness sequence and the second objective function adaptation degree series.

Optionally, in this step after performing 60 generation genetic algorithm iteration, an intermediate population is obtained, in practical behaviour In work to execute genetic algorithm preset times and be not specifically limited.Wherein, the intermediate population's fitness sequence that obtains and the The element that two objective functions adapt in degree series is corresponding, and there is also a corresponding relationships, and the second objective function adapts to Degree series are corresponding with the element in first object function adaptation degree series, and first object function adapts to degree series and current location Gas concentration information (monitoring time and the corresponding gas concentration size of monitoring time and gas corresponding with monitoring time are dense Spend position) between there is also a corresponding relationships.Therefore, the gas concentration information of intermediate population's fitness and current location (the monitoring time and corresponding gas concentration size of monitoring time and gas concentration corresponding with monitoring time position) there is also One one-to-one relationship.For example, obtained intermediate population's fitness sequence can be expressed as M_m1、M_m2......M_mn, wherein Mn is less than or equal to ttn, and ttn is less than or equal to tn, M_m1、M_m2......M_mnIn each element with the second objective function fitness Sequence T_tt1、T_tt2......T_ttnIn only one element it is corresponding, T_tt1、T_tt2......T_ttnIn each element with One objective function adapts to degree series T_t1、T_t2......T_tnIn only one element it is corresponding, it is also dense with the gas of current location Spend information sequence A_a1、A_a2......A_anIn only one element it is corresponding, therefore, M_m1、M_m2......M_mnIn every unitary Element with current location gas concentration information sequence A_a1、A_a2......A_anIn only one element it is corresponding.

Step S224 carries out local optimal searching to intermediate population's fitness sequence, obtains optimal intermediate population's fitness.

In embodiments of the present invention, intermediate population's fitness sequence is ranked up, chooses intermediate population's fitness and is in The individual execution pattern searching algorithm of preset ratio carries out local optimal searching, optimal intermediate population's fitness is obtained, by step Explanation in S223 is it is found that optimal intermediate population's fitness at this time and current location gas concentration information sequence A_a1、 A_a2......A_anIn only one element it is corresponding, that is, correspond to unique one group gas concentration size corresponding with monitoring time With gas concentration corresponding with monitoring time position.Optionally, the individual that intermediate population's fitness is in preset ratio can root It, can be in preceding 20% intermediate population's fitness according to specific implementation situation setting.

Step S225, according to the first corresponding relationship, the second corresponding relationship, optimal intermediate population's fitness and objective function mould Type obtains target function model calculated value.

Wherein, after step S224 obtains optimal intermediate population's fitness, because during optimal intermediate population's fitness is contained in Between in Population adaptation degree series, and intermediate population's fitness sequence with the second objective function adapts to degree series there are the second corresponding passes System then can adapt to find out and optimal intermediate population's fitness in degree series from the second objective function according to the second corresponding relationship Element in corresponding second objective function fitness is (wherein, for the ease of subsequent descriptions, by the second objective function fitness sequence Element in column in the second objective function fitness corresponding with optimal intermediate population's fitness is referred to as the second object element), Further, because the second objective function adaptation degree series are corresponding with the element in first object function adaptation degree series, then It can adapt to find out with the element in first object function fitness (wherein, in order to just in degree series from first object function In subsequent descriptions, first object function is adapted into element corresponding with first object element in degree series and is referred to as first object member Element), further because first object function adapt to degree series with gas concentration size corresponding with monitoring time and with monitoring There are the first corresponding relationships between time corresponding gas concentration position, then according to the first corresponding relationship, can determine and first The corresponding current location gas concentration information of object element (monitoring time, gas concentration size corresponding with monitoring time and with The corresponding gas concentration position of monitoring time), it will current location gas concentration information input mesh corresponding with first object element Offer of tender exponential model obtains determining target function model calculated value.

In above-described embodiment, concentration sequence is calculated according to target function model, the concentration of gas and model of atmospheric diffusion Column obtain first object function and adapt to degree series, and adapt to execute after degree series pre-process to first object function default The genetic algorithm of number obtains intermediate population's fitness sequence, further, obtains optimal intermediate Population adaptation with local optimal searching Degree.According to the first corresponding relationship, the second corresponding relationship, optimal intermediate population's fitness and target function model obtain specifically with The corresponding target function model meter of the corresponding gas concentration size of monitoring time and gas concentration corresponding with monitoring time position Calculation value.

In one of the embodiments, as shown in figure 4, being the flow diagram of the refinement step of step S221, according to mesh Offer of tender exponential model, the gas concentration information of current location and model of atmospheric diffusion calculate concentration sequence, obtain first object function Degree series are adapted to, specifically includes the following steps:

Step S2211, according to and the corresponding gas concentration size of monitoring time and gas concentration corresponding with monitoring time Position determines starting leakage source strength collection and starting leakage source position collection.

Wherein, to leakage source position L=(x, y) and leakage source strength Q₀Real coding is carried out, and carries out random initializtion, Available starting leakage source strength collection and starting leakage source position collection.Optionally, to leakage source position L=(x, y) and leakage Source strength Q₀Carry out real coding and the step for random initializtion, be according to gas concentration size corresponding with monitoring time and Gas concentration corresponding with monitoring time position carries out real coding and random initializtion, at this time gas corresponding with monitoring time After bulk concentration size and gas concentration corresponding with monitoring time position refer to the gas concentration information for monitoring current location, with Any corresponding gas concentration size of monitoring time and gas concentration corresponding with any monitoring time position.

Step S2212 calculates concentration sequence according to starting leakage source strength collection and model of atmospheric diffusion, obtains initial atmosphere Diffusion model calculates concentration collection.

Wherein, it only includes unknown parameter leakage source strength Q that model of atmospheric diffusion, which calculates concentration sequence,₀, right in step S2211 Leak source strength Q₀It is initialized, by the Q after initialization₀It inputs model of atmospheric diffusion and calculates concentration sequence, can obtain not Initial atmosphere diffusion model containing unknown parameter calculates concentration collection.

Step S2213 believes according to the gas concentration of starting leakage source strength collection, starting leakage source position collection, current location Breath and initial atmosphere diffusion model calculate concentration collection, obtain first object function and adapt to degree series.

Source of leaks strength set, starting leakage source position collection, the gas concentration information of current location and initial atmosphere are spread Model calculates in concentration collection input target function model formula (1), obtains first object function and adapts to degree series.

In above-described embodiment, obtain first object function adapt to degree series be determine objective function in leakage source strength with It leaks source position and basis is provided.

Step S400 in one of the embodiments, according to the target function model calculated value and next bit of current location The target function model calculated value set determines leakage source strength and leakage source position, comprising:

Gas concentration position corresponding with the target function model calculated value of current location is determined as the first source of leaks position It sets, gas concentration position corresponding with the target function model calculated value of the next position is determined as the second leakage source position；Really The distance between fixed first leakage source position and the second leakage source position value；If distance value is less than or equal to the first preset threshold, Gas concentration corresponding with the second leakage source position is then sized to leakage source strength, the second leakage source position is determined as Leak source position.

Optionally, the method for determination of the next position are as follows:

Wherein, Δ L (i) is the traveling step-length of i-th of time step, κ_iFor a random generation number (κ_i∈ [0,1]), α is most small step It is long, and ρ ∈ (1,3] it is step-length regulatory factor；Δ Θ (i) is the direction of travel of i-th of time step, ∈_iFor a random generation number (∈_i∈ [0,1]), θ_kFor direction regulatory factor.

Since initial position, the step-length and direction generated by formula (3) and (4) is gradually carried forward gas search, works as i= When 1, direction regulatory factor θ is taken₁For upwind side, when i is more than or equal to 2, direction regulatory factor θ_iDirection is to be referred to by current location To the first leakage source position.Hereafter each traveling time step, it is random to generate number κ_iAnd ∈_iIt is all different, so that in traveling process Step-length and direction realize dynamic adjust.

Wherein it is determined that leakage source strength and the foundation of leakage source position are, when the first leakage source position and the second source of leaks When the distance between position is less than or equal to a preset threshold value, illustrate that the leakage source position released counter twice is very close, Show that source of leaks estimation result is sufficiently stable, is then sized to leak by gas concentration corresponding with the second leakage source position Source strength is determined as the second leakage source position to leak source position.

Optionally, if above-mentioned distance value is greater than the first preset threshold, using described the next position as updated current Position, and the gas concentration information for executing monitoring the next position is returned, according to the gas concentration information, Atmospheric Diffusion of the next position Model and target function model, the step of obtaining the target function model calculated value of the next position.

In embodiments of the present invention, if the distance between the first leakage source position and the second leakage source position are greater than first in advance If when threshold value, illustrating the leakage source position wide apart released counter twice, showing that source of leaks estimation result is stable not enough.This When, then using the gas concentration information of described the next position as the gas concentration information of updated current location, the next position The next position target function model calculated value of the target function model calculated value as updated current location, next bit Determining the second leakage source strength and the second leakage source position are set as updated first leakage source strength and the first source of leaks Position.Further, and by the method for determination of above-mentioned the next position, the gas concentration information of the next position is reacquired, again Be calculated and determined the second leakage source position and second leakage source position, until first leakage source position and second leakage source position it Between distance value be less than or equal to the first preset threshold, then will be sized to let out with the second corresponding gas concentration in leakage source position Drain-source intensity is determined as the second leakage source position to leak source position.

It, can be according to the target function model calculated value of current location and the objective function of the next position in above-described embodiment Model calculation value, obtains leakage source strength and leakage source position, and realization is accurately positioned dangerization gas leakage source.

In one of the embodiments, by gas concentration position corresponding with the target function model calculated value of current location It is determined as the first leakage source position, gas concentration position corresponding with the target function model calculated value of the next position is determined as Second leakage source position, comprising:

If the target function model calculated value of current location be less than or equal to the second preset threshold, by with current location The corresponding gas concentration position of target function model calculated value is determined as the first leakage source position.

In embodiments of the present invention, the target function model calculated value of current location is less than or equal to the second preset threshold One preset condition of convergence indicates that model of atmospheric diffusion calculates concentration C_Mol(L_i, t_i) and ionic migration spectrometer monitoring concentration C_IMS (L_i, t_i) difference quadratic sum be less than or equal to setting the limits of error, quadratic sum be less than or equal to the limits of error show anti-knot fruit very Close to the actual information in actual leakage source, then it will meet gas concentration corresponding to the condition of convergence and be sized to the first leakage Source strength will meet gas concentration position corresponding to the condition of convergence and be determined as the first leakage source position.

Optionally, the target function model calculated value of current location is greater than the second preset threshold, then returns to execution monitoring and work as The step of front position gas concentration information, until the target function model calculated value of current location is less than or equal to the second default threshold Value, then be determined as the first leakage source position for gas concentration position corresponding with the target function model calculated value of current location. Specifically, if the target function model calculated value of current location is greater than the second preset threshold, the gas for updating current location is dense Spend information, wherein the gas concentration information for updating current location refers to the current location for resurveying or updating for calculating Gas concentration information acquisition time and position, and by the gas concentration information input of new current location preset Atmospheric Diffusion Model obtains updated model of atmospheric diffusion and calculates concentration sequence, and by the gas concentration information of updated current location Concentration sequence inputting target function model is calculated with updated model of atmospheric diffusion, obtains updated target function model meter Calculation value.The size relation of more updated target function model calculated value and the second preset threshold, until updated target Function model calculated value is less than or equal to the second preset threshold, then will be corresponding with the target function model calculated value of current location Gas concentration is sized to the first leakage source strength, and gas corresponding with the target function model calculated value of current location is dense Degree position is determined as the first leakage source position.

It optionally, will be under if the target function model calculated value of the next position is less than or equal to the second preset threshold The corresponding gas concentration position of the target function model calculated value of one position is determined as the second leakage source position.

In embodiments of the present invention, gas concentration position corresponding with the target function model calculated value of the next position is true The mode for being set to the second leakage source position is identical as the mode that current location is determined as the first leakage source position, no longer superfluous herein It states.

In above-described embodiment, the can be obtained according to the target function model calculated value of current location and the default condition of convergence One leakage source strength and the first leakage source position, and the target function model calculated value in current location is unsatisfactory for default convergence item The parameter of model involved in localization method is constantly corrected in the case where part, obtains the first leakage source strength and the first source of leaks position It sets.

In a specific embodiment, specific step is as follows for leak position:

Since initial position, the step-length and direction generated with step S300 Chinese style (3) and formula (4) is gradually transported forward It is dynamic.At this point, taking θ_iFor upwind side (upwind side can be determined by ultrasonic aerovane), robot location L_i= (x_i, y_i) can be updated by formula (5) and formula (6), and advance by the position updated:

x_i+1=x_i+ΔL(i)·cos(ΔΘ(i)) (5)

y_i+1=y_i+ΔL(i)·sin(ΔΘ(i)) (6)

Once ionic migration spectrometer is in position L_i=(x_i, y_i), i=k₁Place measures concentration data, before robot stops immediately Into, and n moment, the concentration data input type (3) that this n moment is measured, the source of leaks position estimated are stopped at this It sets

By robot present position L_i=(x_i, y_i), (i >=k₁) it is directed toward the leakage source position estimatedDirection as formula (2) the new direction regulatory factor in, i.e.,It is passed by formula (3) and formula (4) It pushes away and generates step-length and direction, and further march to moment i >=k by the position that formula (5) and formula (6) update₂, in-position L_i= (x_i, y_i), (i=k₂):

If in L_i=(x_i, y_i), (i=k₂) there is concentration data output, then it is stopped again at this n moment, and by this n The position input type (3) for the concentration data and robot that moment measures, the leakage source position estimated

By robot present position L_i=(x_i, y_i), (i >=k₂) it is directed toward the leakage source position estimatedDirection as formula (2) the new direction regulatory factor in, i.e.,It is passed by formula (1) and formula (2) It pushes away and generates step-length and direction, and further march to moment i=k by the position that formula (5) and formula (6) update₃, in-position L_i= (x_i, y_i), (i=k₃):

It repeats the above steps, obtains preceding i=k_mThe source of leaks at a moment estimates positionIf adjacent two It is secondary to estimate source of leaks position deviationThen show source of leaks estimation result It is sufficiently stable.Wherein ε is the limits of error of a setting.

In one embodiment, as shown in figure 5, providing a kind of leak position dress based on variable step recursion track It sets, comprising: building module 501, the first determining module 502, the second determining module 503 and locating module 504, in which:

Module 501 is constructed, for the target letter according to preset model of atmospheric diffusion building and gas concentration information association Exponential model；

First determining module 502, for monitoring the gas concentration information of current location, according to the gas concentration of current location Information, model of atmospheric diffusion and target function model obtain the target function model calculated value of current location；

Second determining module 503, for monitoring the gas concentration information of the next position, according to the gas concentration of the next position Information, model of atmospheric diffusion and target function model obtain the target function model calculated value of the next position；

Locating module 504, for according to the target function model calculated value of current location and the objective function of the next position Model calculation value determines leakage source strength and leakage source position.

The first determining module 502 is also used to execute the gas concentration letter of monitoring current location in one of the embodiments, Breath, according to the gas concentration information and model of atmospheric diffusion of current location, the model of atmospheric diffusion calculating for obtaining current location is dense Degree series；It is calculated according to the model of atmospheric diffusion of target function model, the gas concentration information of current location and current location dense Degree series obtain the target function model calculated value of current location.

The first determining module 502 is also used to execute will acquire and model of atmospheric diffusion in one of the embodiments, Relevant model parameter, monitoring time and gas concentration corresponding with monitoring time position input model of atmospheric diffusion, are worked as The model of atmospheric diffusion of front position calculates concentration sequence.

The first determining module 502 is also used to execute the gas according to current location in level in one of the embodiments, Direction, the diffusion coefficient of vertical direction and wind speed information corresponding with monitoring time, obtain the Atmospheric Diffusion with current location The relevant model parameter of model.

The first determining module 502 is also used to execute according to target function model, current location in one of the embodiments, Gas concentration information and current location model of atmospheric diffusion calculate concentration sequence, obtain first object function fitness sequence Column；Wherein, it is corresponding with gas concentration size corresponding with monitoring time and with monitoring time to adapt to degree series for first object function There are the first corresponding relationships for gas concentration position；Degree series are adapted to first object function to be ranked up and screen, and obtain the Two objective functions adapt to degree series；Degree series are adapted to the second objective function and execute preset times genetic algorithm, obtain among Group adapts to degree series；Wherein, there are the second corresponding relationships for intermediate population's fitness sequence and the second objective function adaptation degree series； Local optimal searching is carried out to intermediate population's fitness sequence, obtains optimal intermediate population's fitness；According to the first corresponding relationship, second Corresponding relationship, optimal intermediate population's fitness and target function model, obtain the target function model calculated value of current location.

The first determining module 502 is also used to execute according to gas corresponding with monitoring time in one of the embodiments, Concentration and gas concentration corresponding with monitoring time position, determine starting leakage source strength collection and starting leakage source position Collection；Concentration sequence is calculated according to starting leakage source strength collection and model of atmospheric diffusion, it is dense to obtain the calculating of initial atmosphere diffusion model Degree collection；Expanded according to starting leakage source strength collection, starting leakage source position collection, the gas concentration information of current location and initial atmosphere It dissipates model and calculates concentration collection, obtain first object function and adapt to degree series.

Locating module 504 is also used to execute by the target function model meter with current location in one of the embodiments, The corresponding gas concentration position of calculation value is determined as the first leakage source position, by the target function model calculated value pair with the next position The gas concentration position answered is determined as the second leakage source position；It determines between the first leakage source position and the second leakage source position Distance value；It, will gas concentration size corresponding with the second leakage source position if distance value is less than or equal to the first preset threshold It is determined as leaking source strength, is determined as the second leakage source position to leak source position.

If locating module 504 is also used to execute distance value greater than the first preset threshold in one of the embodiments, incite somebody to action The next position returns to the gas concentration information for executing monitoring the next position as updated current location, according to next bit Gas concentration information, model of atmospheric diffusion and the target function model set obtain the target function model calculated value of the next position The step of.

If the target function model that locating module 504 is also used to execute current location in one of the embodiments, calculates Value is less than or equal to the second preset threshold, then by gas concentration position corresponding with the target function model calculated value of current location It is determined as the first leakage source position；If the target function model calculated value of the next position is less than or equal to the second preset threshold, Gas concentration position corresponding with the target function model calculated value of the next position is determined as the second leakage source position.

Specific restriction about the leak position device based on variable step recursion track may refer to above for base In the restriction of the leak position method of variable step recursion track, details are not described herein.It is above-mentioned based on variable step recursion track Modules in leak position device can be realized fully or partially through software, hardware and combinations thereof.Above-mentioned each module It can be embedded in the form of hardware or independently of in the processor in computer equipment, calculating can also be stored in a software form In memory in machine equipment, the corresponding operation of the above modules is executed in order to which processor calls.

In one embodiment, a kind of computer equipment is provided, which can be server, internal junction Composition can be as shown in Figure 6.The computer equipment include by system bus connect processor, memory, network interface and Database.Wherein, the processor of the computer equipment is for providing calculating and control ability.The memory packet of the computer equipment Include non-volatile memory medium, built-in storage.The non-volatile memory medium is stored with operating system, computer program and data Library.The built-in storage provides environment for the operation of operating system and computer program in non-volatile memory medium.The calculating The database of machine equipment is used for memory leak source location data.The network interface of the computer equipment is used for logical with external terminal Cross network connection communication.To realize a kind of source of leaks based on variable step recursion track when the computer program is executed by processor Localization method.

It will be understood by those skilled in the art that structure shown in Fig. 6, only part relevant to application scheme is tied The block diagram of structure does not constitute the restriction for the computer equipment being applied thereon to application scheme, specific computer equipment It may include perhaps combining certain components or with different component layouts than more or fewer components as shown in the figure.

In one embodiment, a kind of computer equipment, including memory and processor are provided, memory is stored with meter Calculation machine program, processor perform the steps of when executing computer program

Monitor current location gas concentration information, according to the gas concentration information of current location, model of atmospheric diffusion and Target function model obtains the target function model calculated value of current location；

Monitor the next position gas concentration information, according to the gas concentration information of the next position, model of atmospheric diffusion and Target function model obtains the target function model calculated value of the next position；

According to the target function model calculated value of the target function model calculated value of current location and the next position, determination is let out Drain-source intensity and leakage source position.

Monitoring current location is also performed the steps of when processor executes computer program in one of the embodiments, Gas concentration information, according to the gas concentration information and model of atmospheric diffusion of current location, the atmosphere for obtaining current location expands It dissipates model and calculates concentration sequence；Expanded according to the atmosphere of target function model, the gas concentration information of current location and current location It dissipates model and calculates concentration sequence, obtain the target function model calculated value of current location.

In one of the embodiments, processor execute computer program when also perform the steps of will acquire with The relevant model parameter of model of atmospheric diffusion, monitoring time and gas concentration corresponding with monitoring time position input Atmospheric Diffusion Model, the model of atmospheric diffusion for obtaining current location calculate concentration sequence.

It also performs the steps of when processor executes computer program in one of the embodiments, according to current location Gas in the horizontal direction, the diffusion coefficient of vertical direction and wind speed information corresponding with monitoring time, obtain and present bit The relevant model parameter of the model of atmospheric diffusion set.

It also performs the steps of when processor executes computer program in one of the embodiments, according to objective function The model of atmospheric diffusion of model, the gas concentration information of current location and current location calculates concentration sequence, obtains first object Function adapts to degree series；Wherein, first object function adapt to degree series with gas concentration size corresponding with monitoring time and with There are the first corresponding relationships for the corresponding gas concentration position of monitoring time；To first object function adapt to degree series be ranked up and Screening obtains the second objective function and adapts to degree series；Degree series are adapted to the second objective function and execute preset times genetic algorithm, Obtain intermediate population's fitness sequence；Wherein, intermediate population's fitness sequence and the second objective function adapt to degree series there are the Two corresponding relationships；Local optimal searching is carried out to intermediate population's fitness sequence, obtains optimal intermediate population's fitness；According to first pair It should be related to, the second corresponding relationship, optimal intermediate population's fitness and target function model, obtain the objective function mould of current location Type calculated value.

When also performed the steps of when processor executes computer program in one of the embodiments, according to monitoring Between corresponding gas concentration size and gas concentration corresponding with monitoring time position, determine starting leakage source strength collection and initial Leak source position collection；Concentration sequence is calculated according to starting leakage source strength collection and model of atmospheric diffusion, obtains initial atmosphere diffusion Model calculates concentration collection；According to starting leakage source strength collection, starting leakage source position collection, current location gas concentration information and Initial atmosphere diffusion model calculates concentration collection, obtains first object function and adapts to degree series.

In one of the embodiments, processor execute computer program when also perform the steps of by with current location The corresponding gas concentration position of target function model calculated value be determined as the first leakage source position, by the target with the next position The corresponding gas concentration position of function model calculated value is determined as the second leakage source position；Determine the first leakage source position and second Leak the distance between source position value；It, will be with the second leakage source position pair if distance value is less than or equal to the first preset threshold The gas concentration answered is sized to leakage source strength, is determined as the second leakage source position to leak source position.

Distance value is also performed the steps of when if processor executing computer program in one of the embodiments, to be greater than First preset threshold, then using the next position as updated current location, and the gas for returning to execution monitoring the next position is dense It spends information and the mesh of the next position is obtained according to the gas concentration information, model of atmospheric diffusion and target function model of the next position The step of scalar functions model calculation value.

Current location is also performed the steps of when if processor executing computer program in one of the embodiments, Target function model calculated value is less than or equal to the second preset threshold, then by the target function model calculated value pair with current location The gas concentration position answered is determined as the first leakage source position；If the target function model calculated value of the next position is less than or equal to Gas concentration position corresponding with the target function model calculated value of the next position is then determined as second and let out by the second preset threshold Leakage location.

In one embodiment, a kind of computer readable storage medium is provided, computer program is stored thereon with, is calculated Machine program performs the steps of when being executed by processor

Monitoring present bit is also performed the steps of when computer program is executed by processor in one of the embodiments, The gas concentration information set obtains the atmosphere of current location according to the gas concentration information and model of atmospheric diffusion of current location Diffusion model calculates concentration sequence；According to target function model, the gas concentration information of current location and the atmosphere of current location Diffusion model calculates concentration sequence, obtains the target function model calculated value of current location.

Also performing the steps of when computer program is executed by processor in one of the embodiments, will acquire And the relevant model parameter of model of atmospheric diffusion, monitoring time and gas concentration corresponding with monitoring time position input atmosphere expand Model is dissipated, the model of atmospheric diffusion for obtaining current location calculates concentration sequence.

It also performs the steps of when computer program is executed by processor in one of the embodiments, according to present bit The gas set in the horizontal direction, the diffusion coefficient of vertical direction and wind speed information corresponding with monitoring time, obtain with currently The relevant model parameter of the model of atmospheric diffusion of position.

It also performs the steps of when computer program is executed by processor in one of the embodiments, according to target letter The model of atmospheric diffusion of exponential model, the gas concentration information of current location and current location calculates concentration sequence, obtains the first mesh Scalar functions adapt to degree series；Wherein, first object function adapt to degree series with gas concentration size corresponding with monitoring time and There are the first corresponding relationships for gas concentration corresponding with monitoring time position；Degree series are adapted to first object function to be ranked up And screening, it obtains the second objective function and adapts to degree series；Degree series are adapted to the second objective function and execute preset times heredity calculation Method obtains intermediate population's fitness sequence；Wherein, intermediate population's fitness sequence and the second objective function adapt to degree series presence Second corresponding relationship；Local optimal searching is carried out to intermediate population's fitness sequence, obtains optimal intermediate population's fitness；According to first Corresponding relationship, the second corresponding relationship, optimal intermediate population's fitness and target function model, obtain the objective function of current location Model calculation value.

Basis and monitoring are also performed the steps of when computer program is executed by processor in one of the embodiments, Time corresponding gas concentration size and gas concentration corresponding with monitoring time position, determine starting leakage source strength collection and just Begin leakage source position collection；Concentration sequence is calculated according to starting leakage source strength collection and model of atmospheric diffusion, obtains initial atmosphere expansion It dissipates model and calculates concentration collection；According to starting leakage source strength collection, starting leakage source position collection, current location gas concentration information Concentration collection is calculated with initial atmosphere diffusion model, first object function is obtained and adapts to degree series.

Also performed the steps of when computer program is executed by processor in one of the embodiments, by with present bit The corresponding gas concentration position of the target function model calculated value set is determined as the first leakage source position, by the mesh with the next position The corresponding gas concentration position of scalar functions model calculation value is determined as the second leakage source position；Determine the first leakage source position and The distance between two leakage source positions value；It, will be with the second leakage source position if distance value is less than or equal to the first preset threshold Corresponding gas concentration is sized to leakage source strength, is determined as the second leakage source position to leak source position.

If it is big also to perform the steps of distance value when computer program is executed by processor in one of the embodiments, In the first preset threshold, then using the next position as updated current location, and the gas for executing monitoring the next position is returned to Concentration information obtains the next position according to the gas concentration information, model of atmospheric diffusion and target function model of the next position The step of target function model calculated value.

If also performing the steps of current location when computer program is executed by processor in one of the embodiments, Target function model calculated value be less than or equal to the second preset threshold, then by the target function model calculated value with current location Corresponding gas concentration position is determined as the first leakage source position；If the target function model calculated value of the next position is less than or waits In the second preset threshold, then gas concentration position corresponding with the target function model calculated value of the next position is determined as second Leak source position.

For the positioning for realizing harmful influence source of leaks, need to acquire enough monitoring data, using based on monitoring data and danger The target function model that change product model of atmospheric diffusion is constituted calculates leakage source position.In terms of acquiring harmful influence monitoring information, The application using ionic migration spectrometer as monitoring means, can the harmful influence gas to extremely low concentration carry out mobile monitoring, utilize Input data of the gas concentration information for the current location that ion mobility spectrometry measures as target function model.It is being based on variable step It is dense with ionic migration spectrometer measures at a certain position harmful influence of a certain moment in terms of the leak position method of recursion track What the matching degree that degree calculates concentration with harmful influence model of atmospheric diffusion was monitored as objective function, foundation based on ionic migration spectrometer Leak position model, using the parameter for constantly correcting location model based on genetic algorithm and pattern search hybrid optimization algorithm, Obtained optimal source of leaks parameter is estimated leakage source strength and leakage source position.

Those of ordinary skill in the art will appreciate that realizing all or part of the process in above-described embodiment method, being can be with Relevant hardware is instructed to complete by computer program, the computer program can be stored in a non-volatile computer In read/write memory medium, the computer program is when being executed, it may include such as the process of the embodiment of above-mentioned each method.Wherein, To any reference of memory, storage, database or other media used in each embodiment provided herein, Including non-volatile and/or volatile memory.Nonvolatile memory may include read-only memory (ROM), programming ROM (PROM), electrically programmable ROM (EPROM), electrically erasable ROM (EEPROM) or flash memory.Volatile memory may include Random access memory (RAM) or external cache.By way of illustration and not limitation, RAM is available in many forms, Such as static state RAM (SRAM), dynamic ram (DRAM), synchronous dram (SDRAM), double data rate sdram (DDRSDRAM), enhancing Type SDRAM (ESDRAM), synchronization link (Synchlink) DRAM (SLDRAM), memory bus (Rambus) direct RAM (RDRAM), direct memory bus dynamic ram (DRDRAM) and memory bus dynamic ram (RDRAM) etc..

Each technical characteristic of above embodiments can be combined arbitrarily, for simplicity of description, not to above-described embodiment In each technical characteristic it is all possible combination be all described, as long as however, the combination of these technical characteristics be not present lance Shield all should be considered as described in this specification.

The several embodiments of the application above described embodiment only expresses, the description thereof is more specific and detailed, but simultaneously It cannot therefore be construed as limiting the scope of the patent.It should be pointed out that coming for those of ordinary skill in the art It says, without departing from the concept of this application, various modifications and improvements can be made, these belong to the protection of the application Range.Therefore, the scope of protection shall be subject to the appended claims for the application patent.

Claims

1. a kind of leak position method based on variable step recursion track, which is characterized in that the described method includes:

The gas concentration information for monitoring current location, according to the gas concentration information of the current location, the Atmospheric Diffusion mould Type and the target function model obtain the target function model calculated value of the current location；

The gas concentration information for monitoring the next position, according to the gas concentration information of described the next position, the Atmospheric Diffusion mould Type and the target function model obtain the target function model calculated value of described the next position；

According to the target function model calculated value of the target function model calculated value of the current location and described the next position, really Surely leakage source strength and leakage source position.

2. the method according to claim 1, wherein when the gas concentration information of the current location includes monitoring Between and gas concentration corresponding with monitoring time position；

The gas concentration information of the monitoring current location, expands according to the gas concentration information of the current location, the atmosphere Model and the target function model are dissipated, the target function model calculated value of the current location is obtained, comprising:

The gas concentration information for monitoring the current location expands according to the gas concentration information of the current location and the atmosphere Model is dissipated, the model of atmospheric diffusion for obtaining the current location calculates concentration sequence；

According to the Atmospheric Diffusion mould of the target function model, the gas concentration information of the current location and the current location Type calculates concentration sequence, obtains the target function model calculated value of the current location.

3. according to the method described in claim 2, it is characterized in that, the gas concentration information of the monitoring current location, According to the gas concentration information and the model of atmospheric diffusion of the current location, the Atmospheric Diffusion mould of the current location is obtained Type calculates concentration sequence, comprising:

The model parameter relevant to the model of atmospheric diffusion that will acquire, the monitoring time and with the monitoring time pair The gas concentration position answered inputs the model of atmospheric diffusion, and the model of atmospheric diffusion for obtaining the current location calculates concentration sequence Column.

4. according to the method described in claim 3, it is characterized in that, described will acquire related to the model of atmospheric diffusion Model parameter, the monitoring time and gas concentration corresponding with monitoring time position input the Atmospheric Diffusion mould Type, before further include:

According to the gas of the current location in the horizontal direction, the diffusion coefficient of vertical direction and corresponding with the monitoring time Wind speed information, obtain model parameter relevant to the model of atmospheric diffusion of the current location.

5. according to the method described in claim 2, it is characterized in that, when the gas concentration information of the current location includes monitoring Between and the corresponding gas concentration size of the monitoring time and gas concentration corresponding with monitoring time position；

It is described to be expanded according to the atmosphere of the target function model, the gas concentration information of the current location and the current location It dissipates model and calculates concentration sequence, obtain the target function model calculated value of the current location, comprising:

According to the Atmospheric Diffusion mould of the target function model, the gas concentration information of the current location and the current location Type calculates concentration sequence, obtains first object function and adapts to degree series；Wherein, the first object function adapt to degree series with The corresponding gas concentration size of monitoring time and gas concentration position corresponding with the monitoring time are corresponding there are first Relationship；

Degree series are adapted to the first object function to be ranked up and screen, and are obtained the second objective function and are adapted to degree series；

Degree series are adapted to second objective function and execute preset times genetic algorithm, obtain intermediate population's fitness sequence； Wherein, there are the second corresponding relationships for the intermediate population's fitness sequence and second objective function adaptation degree series；

According to first corresponding relationship, second corresponding relationship, the optimal intermediate population's fitness and the target letter Exponential model obtains the target function model calculated value of the current location.

6. according to the method described in claim 5, it is characterized in that, described according to the target function model, the present bit The model of atmospheric diffusion of the gas concentration information and the current location set calculates concentration sequence, obtains the adaptation of first object function Degree series, comprising:

According to and the corresponding gas concentration size of the monitoring time and gas concentration corresponding with monitoring time position, really Determine starting leakage source strength collection and starting leakage source position collection；

Concentration sequence is calculated according to the starting leakage source strength collection and the model of atmospheric diffusion, obtains initial atmosphere diffusion mould Type calculates concentration collection；

According to the starting leakage source strength collection, starting leakage source position collection, the current location gas concentration information Concentration collection is calculated with the initial atmosphere diffusion model, the first object function is obtained and adapts to degree series.

7. the method according to claim 1, wherein the target function model meter according to the current location The target function model calculated value of calculation value and described the next position determines leakage source strength and leakage source position, comprising:

Gas concentration position corresponding with the target function model calculated value of the current location is determined as the first source of leaks position It sets, gas concentration position corresponding with the target function model calculated value of described the next position is determined as the second source of leaks position It sets；

It, will gas concentration corresponding with the second leakage source position if the distance value is less than or equal to the first preset threshold It is sized to the leakage source strength, second leakage source position is determined as the leakage source position.

8. the method according to the description of claim 7 is characterized in that the method also includes:

If the distance value is greater than first preset threshold, using described the next position as updated current location, and The gas concentration information for executing described monitoring the next position is returned to, according to the gas concentration information of described the next position, described big Gas diffusion model and the target function model, the step of obtaining the target function model calculated value of described the next position.

9. the method according to the description of claim 7 is characterized in that the target function model meter by with the current location The corresponding gas concentration position of calculation value is determined as the first leakage source position, will calculate with the target function model of described the next position It is worth corresponding gas concentration position and is determined as the second leakage source position, comprising:

If the target function model calculated value of the current location is less than or equal to the second preset threshold, will be with the present bit The corresponding gas concentration position of the target function model calculated value set is determined as first leakage source position；

If the target function model calculated value of described the next position be less than or equal to second preset threshold, will with it is described under The corresponding gas concentration position of the target function model calculated value of one position is determined as second leakage source position.

10. a kind of leak position device based on variable step recursion track, which is characterized in that described device includes:

Module is constructed, for the target function model according to preset model of atmospheric diffusion building and gas concentration information association；

First determining module is believed for monitoring the gas concentration information of current location according to the gas concentration of the current location Breath, the model of atmospheric diffusion and the target function model, obtain the target function model calculated value of the current location；

Second determining module is believed for monitoring the gas concentration information of the next position according to the gas concentration of described the next position Breath, the model of atmospheric diffusion and the target function model, obtain the target function model calculated value of described the next position；

Locating module, for according to the target function model calculated value of the current location and the objective function of described the next position Model calculation value determines leakage source strength and leakage source position.

11. a kind of computer equipment, including memory and processor, the memory are stored with computer program, feature exists In the step of processor realizes any one of claims 1 to 9 the method when executing the computer program.

12. a kind of computer readable storage medium, is stored thereon with computer program, which is characterized in that the computer program The step of method described in any one of claims 1 to 9 is realized when being executed by processor.