CN106951958A - A kind of mixing artificial bee colony algorithm of inverting the earth parameter - Google Patents
A kind of mixing artificial bee colony algorithm of inverting the earth parameter Download PDFInfo
- Publication number
- CN106951958A CN106951958A CN201710207227.7A CN201710207227A CN106951958A CN 106951958 A CN106951958 A CN 106951958A CN 201710207227 A CN201710207227 A CN 201710207227A CN 106951958 A CN106951958 A CN 106951958A
- Authority
- CN
- China
- Prior art keywords
- search
- colony algorithm
- bee colony
- artificial bee
- formula
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/004—Artificial life, i.e. computing arrangements simulating life
- G06N3/006—Artificial life, i.e. computing arrangements simulating life based on simulated virtual individual or collective life forms, e.g. social simulations or particle swarm optimisation [PSO]
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- General Health & Medical Sciences (AREA)
- Biomedical Technology (AREA)
- Biophysics (AREA)
- Computational Linguistics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Evolutionary Computation (AREA)
- Artificial Intelligence (AREA)
- Molecular Biology (AREA)
- Computing Systems (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Software Systems (AREA)
- Health & Medical Sciences (AREA)
- Medicines Containing Material From Animals Or Micro-Organisms (AREA)
Abstract
The embodiment of the invention discloses a kind of mixing artificial bee colony algorithm of inverting the earth parameter, for solve to there is no in the prior art it is a kind of can efficiently, quick obtaining globally optimal solution and can effectively cover earth's surface to up to a hundred kilometers of underground depth ground resistivity measure and inversion method technical problem.The embodiment of the present invention includes:Obtain the object function for carrying out comprehensive inversion to the earth electric resistance of soil parameter by four-electrode method and magnetotelluric method;Globally optimal solution is scanned for the object function by using the artificial bee colony algorithm of embedded conjugate gradient algorithms progress Local Search, and obtains globally optimal solution.
Description
Technical field
The present invention relates to the mixing people worker bee of the earth soil resistivity fields of measurement, more particularly to a kind of inverting the earth parameter
Group's algorithm.
Background technology
During the debugging of DC transmission system or monopole failure, straight-flow system is connect with greatly returning to the method for operation to direct current
Earth polar, which is flowed into, may up to count kA earth current.The continuous action of this earth current is to being grounded the height of operation in power system
Piezoelectric forces transformer, long distance oil-gas pipeline, earthquake monitoring stations and networks etc. cause serious adverse effect.In order to study earth current
The regularity of distribution in ground, domestic scholars propose a variety of models, such as power network direct current distributed simulation model and the distribution of pipeline direct current
The ground electromagnetic field appraising model of model, earthquake monitoring stations and networks, example shows that these methods achieve better application effect.
On the basis of DC current distribution theory model is further improved, accurate rational ground resistivity how is chosen
Parameter model is the key for improving the model solution degree of accuracy.In reality, resistivity distribution shows higher complexity, and this causes me
The power transformer of ground connection operation faces totally different D.C. magnetic biasing risk in the AC system of state different regions, while D.C. magnetic biasing
The improvement of harm also becomes sufficiently complex.And there is no at present one kind can efficiently, quick obtaining globally optimal solution and can be effective
Covering earth's surface to up to a hundred kilometers of underground depth ground resistivity measurement and inversion method.
The content of the invention
The embodiments of the invention provide a kind of mixing artificial bee colony algorithm of inverting the earth parameter, solve in the prior art
There is no it is a kind of can efficiently, quick obtaining globally optimal solution and can effectively cover earth's surface to the big of up to a hundred kilometers of underground depth
Earth resistivity measures the technical problem with inversion method.
A kind of mixing artificial bee colony algorithm of inverting the earth parameter provided in an embodiment of the present invention, including:
Obtain the object function for carrying out comprehensive inversion to the earth electric resistance of soil parameter by four-electrode method and magnetotelluric method;
The artificial bee colony algorithm for carrying out Local Search by using embedded conjugate gradient algorithms is carried out to the object function
Globally optimal solution is searched for, and obtains globally optimal solution.
Alternatively, described obtain carries out comprehensive inversion by four-electrode method and magnetotelluric method to the earth electric resistance of soil parameter
Object function includes:
Obtained by formula one by four-electrode method and magnetotelluric method to the earth electric resistance of soil parameter progress comprehensive inversion
Object function, the formula one is specially:
MinF=ω1F1+ω2F2;
Wherein, F is the object function of comprehensive inversion, F1And F2The respectively inverting target letter of four-electrode method and magnetotelluric method
Number, ω1And ω2Respectively F1And F2Weight;
The F1Tried to achieve by formula two, the formula two is specially:
The F2Tried to achieve by formula three, the formula three is specially:
Wherein, m1、di、ρaAnd ρMRespectively the number of four-electrode method measurement data, pole span, the inverting value of apparent resistivity and
Apparent resistivity measured value, m2、fiThe respectively number of magnetotelluric method measurement data, characteristic frequency.
Alternatively, the ω1For 1, the ω2For 1.
Alternatively, it is described to carry out the artificial bee colony algorithm of Local Search to the mesh by using embedded conjugate gradient algorithms
Scalar functions scan for globally optimal solution, and obtain globally optimal solution and include:
The search of global solution is carried out to the object function by artificial bee colony algorithm, and in the artificial bee colony algorithm
Observe embedded conjugate gradient algorithms in the link of honeybee progress Local Search and carry out Local Search.
Alternatively, the search for carrying out global solution to the object function by artificial bee colony algorithm, and in the people
Embedded conjugate gradient algorithms progress Local Search is specifically included in the link of the observation honeybee progress Local Search of work ant colony algorithm:
S1, carry out using artificial bee colony algorithm searching the activity of obtaining, and set the number N=100 solved in population;
S2, initial population S is randomly generated, row information of going forward side by side exchange and local selection course, all people's work honeybee retain
The greedy selection course in preferable nectar source, search bee searches for the random selection process in new nectar source, and observation honeybee is according to related to nectar source
Probable value PiSelect nectar source, the probable value PiTried to achieve by formula four, the formula four is specially:
Wherein, SiIt is that (i=1,2 ..., N) is the vector of a D dimension, D is soil model ginseng for the solution in initial population S
Several numbers, F (Si) it is solution SiTarget function value, i=1,2 ..., N.
S3, observation honeybee are according to the probable value PiA position candidate is produced from the position in green molasses source:
Sij1=Sij+φij(Sij-Soj)
Wherein, o is the nectar source numbering different from i, and j is randomly selected subscript (1≤j≤D), φijFor between [0,1]
Random number, for controlling the generation of nectar source position in S neighborhoods, position candidate Sij1Represent green molasses source position SijWith in neighborhood with
One nectar source S of machineojBetween relativity;
S4, observation the honeybee embedded conjugate gradient algorithms in Local Search link, and take Local Search times NLimit=3, if
Nectar source position SiBy limited number of time NLimitAfter the cyclic search of gathering honey honeybee and observation honeybee, it is impossible to be modified, that is, abandon the position
Put, and gathering honey honeybee is changed into search bee, and green molasses source is replaced in one nectar source of random search, search bee is by formula five by searching at random
Rope determines new nectar source, and the formula five is specially:
Wherein,For in kth generation circulationBound.
As can be seen from the above technical solutions, the embodiment of the present invention has advantages below:
The embodiments of the invention provide a kind of mixing artificial bee colony algorithm of inverting the earth parameter, including:Obtain and pass through four
Pole method and magnetotelluric method carry out the object function of comprehensive inversion to the earth electric resistance of soil parameter;By using embedded conjugate gradient
The artificial bee colony algorithm that algorithm carries out Local Search scans for globally optimal solution to the object function, and obtains global optimum
The mixing artificial bee colony algorithm of embedded Local conjugation gradient algorithm search is by simulating bee colony gathering honey row in solution, the embodiment of the present invention
To realize the solution of optimization problem, i.e., nectar source is determined by gathering honey honeybee, and nectar source is plucked and nectar source information is remembered, then with sight
Examine honeybee share nectar source information;Observation honeybee makes a choice by certain selection strategy in neighbouring nectar source;It is abandoned at nectar source
Gathering honey honeybee is changed into search bee and the new nectar source of random search, wherein, the search of Local conjugation gradient algorithm is embedded in, is substantially increased
Local search ability, meets the requirement of search globally optimal solution, and greatly accelerates the solution of soil inversion problem, and formation more has
Ground resistivity measurement and inversion method of the covering earth's surface of effect to up to a hundred kilometers of underground depth.
Brief description of the drawings
In order to illustrate more clearly about the embodiment of the present invention or technical scheme of the prior art, below will be to embodiment or existing
The accompanying drawing to be used needed for having technology description is briefly described, it should be apparent that, drawings in the following description are only this
Some embodiments of invention, for those of ordinary skill in the art, without having to pay creative labor, may be used also
To obtain other accompanying drawings according to these accompanying drawings.
Fig. 1 illustrates for a kind of flow of mixing artificial bee colony algorithm of inverting the earth parameter provided in an embodiment of the present invention
Figure;
Fig. 2 is specific for the mixing artificial bee colony algorithm of embedded Local conjugation gradient algorithm search provided in an embodiment of the present invention
Schematic flow sheet;
Fig. 3 is conjugate gradient method algorithm flow chart provided in an embodiment of the present invention.
Embodiment
The embodiments of the invention provide a kind of mixing artificial bee colony algorithm of inverting the earth parameter, for solving prior art
In there is no it is a kind of can efficiently, quick obtaining globally optimal solution and can effectively cover earth's surface to up to a hundred kilometers of underground depth
Ground resistivity measures the technical problem with inversion method.
To enable goal of the invention of the invention, feature, advantage more obvious and understandable, below in conjunction with the present invention
Accompanying drawing in embodiment, the technical scheme in the embodiment of the present invention is clearly and completely described, it is clear that disclosed below
Embodiment be only a part of embodiment of the invention, and not all embodiment.Based on the embodiment in the present invention, this area
All other embodiment that those of ordinary skill is obtained under the premise of creative work is not made, belongs to protection of the present invention
Scope.
Referring to Fig. 1, a kind of mixing artificial bee colony algorithm of inverting the earth parameter provided in an embodiment of the present invention includes:
101st, the target letter for carrying out comprehensive inversion to the earth electric resistance of soil parameter by four-electrode method and magnetotelluric method is obtained
Number;
First, obtained and the earth electric resistance of soil parameter is carried out comprehensively instead by four-electrode method and magnetotelluric method by formula one
The object function drilled, the formula one is specially:
MinF=ω1F1+ω2F2;
Wherein, F is the object function of comprehensive inversion, F1And F2The respectively inverting target letter of four-electrode method and magnetotelluric method
Number, ω1And ω2Respectively F1And F2Weight, 1 can be taken respectively;
The F1Tried to achieve by formula two, the formula two is specially:
The F2Tried to achieve by formula three, the formula three is specially:
Wherein, m1、di、ρaAnd ρMRespectively the number of four-electrode method measurement data, pole span, the inverting value of apparent resistivity and
Apparent resistivity measured value, m2、fiThe respectively number of magnetotelluric method measurement data, characteristic frequency.
102nd, the artificial bee colony algorithm of Local Search is carried out to the object function by using embedded conjugate gradient algorithms
Globally optimal solution is scanned for, and obtains globally optimal solution.
Then, the search of global solution is carried out to the object function by artificial bee colony algorithm, and in the artificial bee colony
Embedded conjugate gradient algorithms carry out Local Search in the link of the observation honeybee progress Local Search of algorithm, as shown in Fig. 2 being insertion
The mixing artificial bee colony algorithm flow chart of Local conjugation gradient algorithm search, it is comprised the following steps that:
S1, carry out using artificial bee colony algorithm searching the activity of obtaining, and set the number N=100 solved in population;
S2, initial population S is randomly generated, row information of going forward side by side exchange and local selection course, all people's work honeybee retain
The greedy selection course in preferable nectar source, search bee searches for the random selection process in new nectar source, and observation honeybee is according to related to nectar source
Probable value PiSelect nectar source, the probable value PiTried to achieve by formula four, the formula four is specially:
Wherein, SiIt is that (i=1,2 ..., N) is the vector of a D dimension, D is soil model ginseng for the solution in initial population S
Several numbers, F (Si) it is solution SiTarget function value, i=1,2 ..., N.
S3, observation honeybee are according to the probable value PiA position candidate is produced from the position in green molasses source:
Sij1=Sij+φij(Sij-Soj)
Wherein, o is the nectar source numbering different from i, and j is randomly selected subscript (1≤j≤D), φijFor between [0,1]
Random number, for controlling the generation of nectar source position in S neighborhoods, position candidate Sij1Represent green molasses source position SijWith in neighborhood with
One nectar source S of machineojBetween relativity;
S4, observation the honeybee embedded conjugate gradient algorithms in Local Search link, and take Local Search times NLimit=3, if
Nectar source position SiBy limited number of time NLimitAfter the cyclic search of gathering honey honeybee and observation honeybee, it is impossible to be modified, that is, abandon the position
Put, and gathering honey honeybee is changed into search bee, and green molasses source is replaced in one nectar source of random search, search bee is by formula five by searching at random
Rope determines new nectar source, and the formula five is specially:
Wherein,For in kth generation circulationBound, NLimitCan value be 3.As shown in figure 3, being conjugation
Gradient method algorithm flow chart.
It is apparent to those skilled in the art that, for convenience and simplicity of description, the system of foregoing description,
The specific work process of device and unit, may be referred to the corresponding process in preceding method embodiment, will not be repeated here.
The above, the above embodiments are merely illustrative of the technical solutions of the present invention, rather than its limitations;Although with reference to before
Embodiment is stated the present invention is described in detail, it will be understood by those within the art that:It still can be to preceding
State the technical scheme described in each embodiment to modify, or equivalent is carried out to which part technical characteristic;And these
Modification is replaced, and the essence of appropriate technical solution is departed from the spirit and scope of various embodiments of the present invention technical scheme.
Claims (5)
1. a kind of mixing artificial bee colony algorithm of inverting the earth parameter, it is characterised in that including:
Obtain the object function for carrying out comprehensive inversion to the earth electric resistance of soil parameter by four-electrode method and magnetotelluric method;
The artificial bee colony algorithm for carrying out Local Search by using embedded conjugate gradient algorithms is scanned for the object function
Globally optimal solution, and obtain globally optimal solution.
2. the mixing artificial bee colony algorithm of inverting the earth parameter according to claim 1, it is characterised in that the acquisition is led to
Cross four-electrode method and magnetotelluric method includes to the object function of the earth electric resistance of soil parameter progress comprehensive inversion:
The target for carrying out comprehensive inversion to the earth electric resistance of soil parameter by four-electrode method and magnetotelluric method is obtained by formula one
Function, the formula one is specially:
MinF=ω1F1+ω2F2;
Wherein, F is the object function of comprehensive inversion, F1And F2The respectively inversion objective function of four-electrode method and magnetotelluric method,
ω1And ω2Respectively F1And F2Weight;
The F1Tried to achieve by formula two, the formula two is specially:
The F2Tried to achieve by formula three, the formula three is specially:
Wherein, m1、di、ρaAnd ρMRespectively the number of four-electrode method measurement data, pole span, the inverting value of apparent resistivity and apparent are electric
Resistance rate measured value, m2、fiThe respectively number of magnetotelluric method measurement data, characteristic frequency.
3. the mixing artificial bee colony algorithm of inverting the earth parameter according to claim 2, it is characterised in that the ω1For 1,
The ω2For 1.
4. the mixing artificial bee colony algorithm of inverting the earth parameter according to claim 2, it is characterised in that described by adopting
The artificial bee colony algorithm that Local Search is carried out with embedded conjugate gradient algorithms scans for globally optimal solution to the object function,
And obtain globally optimal solution and include:
The search of global solution is carried out to the object function by artificial bee colony algorithm, and in the observation of the artificial bee colony algorithm
Embedded conjugate gradient algorithms carry out Local Search in the link of honeybee progress Local Search.
5. the mixing artificial bee colony algorithm of inverting the earth parameter according to claim 4, it is characterised in that described to pass through people
Work ant colony algorithm carries out the search of global solution to the object function, and carries out part in the observation honeybee of the artificial bee colony algorithm
Embedded conjugate gradient algorithms progress Local Search is specifically included in the link of search:
S1, carry out using artificial bee colony algorithm searching the activity of obtaining, and set the number N=100 solved in population;
S2, initial population S is randomly generated, row information of going forward side by side exchange and local selection course, all people's work honeybee retains preferable
The greedy selection course in nectar source, search bee searches for the random selection process in new nectar source, and observation honeybee is according to the probability related to nectar source
Value PiSelect nectar source, the probable value PiTried to achieve by formula four, the formula four is specially:
Wherein, SiIt is that (i=1,2 ..., N) is the vector of a D dimension, D is soil model parameter for the solution in initial population S
Number, F (Si) is solution SiTarget function value, i=1,2 ..., N.
S3, observation honeybee are according to the probable value PiA position candidate is produced from the position in green molasses source:
Sij1=Sij+φij(Sij-Soj)
Wherein, o is the nectar source numbering different from i, and j is randomly selected subscript (1≤j≤D), φijTo be random between [0,1]
Number, for controlling the generation of nectar source position in S neighborhoods, position candidate Sij1Represent green molasses source position SijWith it is random in neighborhood
One nectar source SojBetween relativity;
S4, observation the honeybee embedded conjugate gradient algorithms in Local Search link, and take Local Search times NLimit=3, if nectar source
Position SiBy limited number of time NLimitAfter the cyclic search of gathering honey honeybee and observation honeybee, it is impossible to be modified, that is, abandon the position, and
Gathering honey honeybee is changed into search bee, and green molasses source is replaced in one nectar source of random search, and search bee is true by random search by formula five
Determine new nectar source, the formula five is specially:
Wherein,For in kth generation circulationBound.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710207227.7A CN106951958A (en) | 2017-03-31 | 2017-03-31 | A kind of mixing artificial bee colony algorithm of inverting the earth parameter |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710207227.7A CN106951958A (en) | 2017-03-31 | 2017-03-31 | A kind of mixing artificial bee colony algorithm of inverting the earth parameter |
Publications (1)
Publication Number | Publication Date |
---|---|
CN106951958A true CN106951958A (en) | 2017-07-14 |
Family
ID=59474941
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710207227.7A Pending CN106951958A (en) | 2017-03-31 | 2017-03-31 | A kind of mixing artificial bee colony algorithm of inverting the earth parameter |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106951958A (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111126591A (en) * | 2019-10-11 | 2020-05-08 | 重庆大学 | Magnetotelluric deep neural network inversion method based on space constraint technology |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104123455A (en) * | 2014-07-22 | 2014-10-29 | 中国地质科学院矿产资源研究所 | Three-dimensional dip inversion method for nonlinear conjugate gradient of magnetotelluric field |
CN105510193A (en) * | 2015-11-30 | 2016-04-20 | 浙江工业大学 | Spectral extinction particle size distribution inversion method based on hybrid intelligent optimization |
CN105976048A (en) * | 2016-04-28 | 2016-09-28 | 苏州泛能电力科技有限公司 | Power transmission network extension planning method based on improved artificial bee colony algorithm |
CN106526329A (en) * | 2016-11-10 | 2017-03-22 | 广东电网有限责任公司电力科学研究院 | Method and device for measuring resistivity from land surface to deep earth |
-
2017
- 2017-03-31 CN CN201710207227.7A patent/CN106951958A/en active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104123455A (en) * | 2014-07-22 | 2014-10-29 | 中国地质科学院矿产资源研究所 | Three-dimensional dip inversion method for nonlinear conjugate gradient of magnetotelluric field |
CN105510193A (en) * | 2015-11-30 | 2016-04-20 | 浙江工业大学 | Spectral extinction particle size distribution inversion method based on hybrid intelligent optimization |
CN105976048A (en) * | 2016-04-28 | 2016-09-28 | 苏州泛能电力科技有限公司 | Power transmission network extension planning method based on improved artificial bee colony algorithm |
CN106526329A (en) * | 2016-11-10 | 2017-03-22 | 广东电网有限责任公司电力科学研究院 | Method and device for measuring resistivity from land surface to deep earth |
Non-Patent Citations (1)
Title |
---|
陈红安: "融合经典搜索策略的微粒群算法及其应用", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111126591A (en) * | 2019-10-11 | 2020-05-08 | 重庆大学 | Magnetotelluric deep neural network inversion method based on space constraint technology |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104466959B (en) | Power system critical circuits discrimination method and system | |
CN106789190B (en) | Vulnerability assessment and route optimization method for power communication network | |
CN103942461B (en) | Water quality parameter Forecasting Methodology based on online online-sequential extreme learning machine | |
CN104601596B (en) | Data-privacy guard method in a kind of Classification Data Mining system | |
CN105471637B (en) | A kind of complex network node importance appraisal procedure and system | |
Zhou et al. | A proof of concept study for criminal network analysis with interactive strategies | |
CN106909989A (en) | A kind of grid disturbance Forecasting Methodology and device | |
CN104268626B (en) | Power cable service life estimation method and system | |
CN106253270A (en) | Electric system vulnerable line identifying method and system | |
CN104122592B (en) | Time shift magnetotelluric signal acquisition and inversion method | |
CN107491694A (en) | Method for quantitative evaluation SCADA system information security fragility | |
CN107819790A (en) | The recognition methods of attack message and device | |
CN104268629A (en) | Complex network community detecting method based on prior information and network inherent information | |
CN108090677A (en) | A kind of critical infrastructures reliability assessment method | |
CN106099891A (en) | Marine wind electric field leading-out terminal sea cable is carried out the analysis method that shunt reactor is distributed rationally | |
CN106951958A (en) | A kind of mixing artificial bee colony algorithm of inverting the earth parameter | |
CN106951615B (en) | Electric power transportation network Topology Structure Design method based on close female algorithm | |
CN106680664A (en) | Power transmission line fault point location comprehensive optimization method | |
CN105515859B (en) | The method and system of community's detection are carried out to symbolic network based on similarity of paths | |
CN108566642A (en) | A kind of two-dimentional union feature authentication method based on machine learning | |
CN104021083B (en) | Method of testing and device | |
CN106789906A (en) | Betting data analysis method and device | |
CN116663052A (en) | Power data privacy protection method, system, equipment and medium under multiparty collaboration | |
CN106645862A (en) | Intelligent identification method of transformer station ground network for maximum step potential difference and maximum contact potential difference | |
CN104168188B (en) | Internet AS deduction and router division method based on SVI |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20170714 |