CN114047554A - Earth resistivity model modeling method and device, computer equipment and storage medium - Google Patents
Earth resistivity model modeling method and device, computer equipment and storage medium Download PDFInfo
- Publication number
- CN114047554A CN114047554A CN202111304750.4A CN202111304750A CN114047554A CN 114047554 A CN114047554 A CN 114047554A CN 202111304750 A CN202111304750 A CN 202111304750A CN 114047554 A CN114047554 A CN 114047554A
- Authority
- CN
- China
- Prior art keywords
- target function
- inversion
- earth
- inversion target
- weight parameter
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 161
- 238000003860 storage Methods 0.000 title claims abstract description 15
- 230000006870 function Effects 0.000 claims abstract description 181
- 239000002689 soil Substances 0.000 claims abstract description 46
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 22
- 238000005259 measurement Methods 0.000 claims description 59
- 238000004590 computer program Methods 0.000 claims description 19
- 239000013598 vector Substances 0.000 claims description 12
- 230000003247 decreasing effect Effects 0.000 claims description 10
- 238000005457 optimization Methods 0.000 claims description 5
- 230000007423 decrease Effects 0.000 claims description 4
- 230000035772 mutation Effects 0.000 claims description 3
- 238000010586 diagram Methods 0.000 description 14
- 230000005684 electric field Effects 0.000 description 9
- 230000008569 process Effects 0.000 description 6
- 238000004364 calculation method Methods 0.000 description 5
- 238000004891 communication Methods 0.000 description 5
- 238000011835 investigation Methods 0.000 description 5
- 230000007704 transition Effects 0.000 description 5
- 238000001514 detection method Methods 0.000 description 4
- 238000009826 distribution Methods 0.000 description 4
- 238000012360 testing method Methods 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 3
- 238000009933 burial Methods 0.000 description 2
- 238000010291 electrical method Methods 0.000 description 2
- 239000000523 sample Substances 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000010438 heat treatment Methods 0.000 description 1
- 239000004973 liquid crystal related substance Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000007935 neutral effect Effects 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 230000007480 spreading Effects 0.000 description 1
- 238000003892 spreading Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
- G01V3/38—Processing data, e.g. for analysis, for interpretation, for correction
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R27/00—Arrangements for measuring resistance, reactance, impedance, or electric characteristics derived therefrom
- G01R27/02—Measuring real or complex resistance, reactance, impedance, or other two-pole characteristics derived therefrom, e.g. time constant
- G01R27/20—Measuring earth resistance; Measuring contact resistance, e.g. of earth connections, e.g. plates
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/30—Assessment of water resources
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Geophysics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Remote Sensing (AREA)
- Geology (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Environmental & Geological Engineering (AREA)
- Geophysics And Detection Of Objects (AREA)
- Measurement Of Resistance Or Impedance (AREA)
Abstract
The application relates to a method and a device for modeling a georesistivity model, computer equipment and a storage medium. The method comprises the following steps: respectively adopting a quadrupole method to obtain first earth resistivity data of a shallow layer, adopting a controllable source audio frequency earth electromagnetic method to measure second earth resistivity data of a middle layer and adopting an earth electromagnetic method to measure third earth resistivity data of a deep layer; inverting the first ground resistivity data to obtain a first inversion target function, inverting the second ground resistivity data to obtain a second inversion target function, and inverting the third ground resistivity data to obtain a third inversion target function; and acquiring a comprehensive inversion target function according to the first, second and third inversion target functions, inverting the comprehensive inversion target function by adopting a differential evolution algorithm, and taking the comprehensive inversion target function with the inverted soil parameters as a ground resistivity model. By adopting the method, the uniform modeling can be accurately carried out on the earth resistivity of a wide area covering tens of kilometers from the earth surface to the underground.
Description
Technical Field
The present application relates to the field of power system grounding technologies, and in particular, to a method and an apparatus for modeling earth resistivity, a computer device, and a storage medium.
Background
With the development of the grounding technology field of the power system, it is a feasible idea to extend the size of the grounding electrode towards the vertical direction, so that when the deep layer of the electrode site has a low earth resistivity, the current can be led to the underground deep. However, the earth electrode has a large burying depth and a wide range of the ground layer, so that earth resistivity exploration is very complicated. In order to better research the distribution rule of the earth resistivity, an earth resistivity modeling method is proposed. However, the traditional earth resistivity modeling method does not consider the distribution rule of earth resistivity at different depths, and cannot accurately and uniformly model the earth resistivity of a wide area covering tens of kilometers from the earth surface to the underground.
Disclosure of Invention
Based on this, it is necessary to provide a method, an apparatus, a computer device and a storage medium capable of accurately modeling a wide-area earth resistivity model covering tens of kilometers from the surface to the underground.
A method of modeling a georesistivity model, the method comprising:
measuring first earth resistivity data of a shallow layer by adopting a quadrupole method;
measuring second earth resistivity data of the middle layer by adopting a controllable source audio magnetotelluric method;
measuring third earth resistivity data of a deep layer by adopting an earth electromagnetic method; the shallow layer, the middle layer and the deep layer are divided according to depth;
inverting the first ground resistivity data to obtain a first inversion target function;
inverting the second ground resistivity data to obtain a second inversion target function;
inverting the third ground resistivity data to obtain a third inversion target function;
acquiring a comprehensive inversion target function according to the first inversion target function, the second inversion target function and the third inversion target function;
and performing soil parameter inversion on the comprehensive inversion target function by adopting a differential evolution algorithm, and taking the comprehensive inversion target function with the inverted soil parameters as a ground resistivity model.
In one embodiment, the step of performing soil parameter inversion on the synthetic inversion target function by using a differential evolution algorithm comprises the following steps:
initializing a population by setting an inversion initial value of a comprehensive inversion target function, wherein the population refers to a population formed by taking each soil parameter in the comprehensive inversion target function as an individual;
iteratively executing variation operation on individuals in the current generation population to generate variation individuals, performing cross operation on the current generation population and the variation individuals to generate experimental individuals, selecting excellent individuals between the experimental individuals and the individuals in the current generation population to form a next generation population, and outputting the soil parameters in the latest generation population as the soil parameters after inversion until termination conditions are met.
In one embodiment, the step of performing mutation operations on individuals in the contemporary population to generate variant individuals comprises:
and randomly selecting three different individuals in the population, and after the vectors of any two individuals are weighted differentially, superposing the vectors of the other two individuals with the vectors of the other one individual to obtain a variant individual.
In one embodiment, the step of obtaining the synthetic inversion objective function according to the first inversion objective function, the second inversion objective function, and the third inversion objective function includes:
configuring a first weight parameter of a first inversion target function, a second weight parameter of a second inversion target function and a third weight parameter of a third inversion target function;
and carrying out weighted summation on the first inversion target function, the second inversion target function and the third inversion target function to obtain a comprehensive inversion target function.
In one embodiment, the value of the first weight parameter decreases with increasing depth of sounding; the value of the second weight parameter is increased and then decreased along with the increase of the sounding depth; and when the value of the second weight parameter is reduced, the value of the third weight parameter is increased along with the increase of the sounding depth and is kept unchanged when the sounding depth reaches the maximum sounding depth.
In one embodiment, the earth resistivity model modeling method further comprises:
configuring an initial value of a first weight parameter to be 1; configuring an initial value of the second weight parameter to 0; the initial value of the third weight parameter is configured to be 0.
When the depth measurement is gradually increased to the first threshold value, the first weight parameter is reduced, the second weight parameter is increased, and when the depth measurement is increased to the first threshold value, the first weight parameter is reduced to 0, and the second weight parameter is increased to 1;
when the depth measurement is gradually increased from the first threshold value to the second threshold value, the second weight parameter is reduced, and the third weight parameter is increased, and when the depth measurement is increased to the second threshold value, the second weight parameter is reduced to 0, and the third weight parameter is increased to 1.
A georesistivity model modeling apparatus, the apparatus comprising:
the first earth resistivity acquisition module is used for measuring first earth resistivity data of a shallow layer by adopting a quadrupole method;
the second earth resistivity acquisition module is used for measuring second earth resistivity data of the middle layer by adopting a controllable source audio frequency earth electromagnetic method;
the third earth resistivity acquisition module is used for measuring deep third earth resistivity data by adopting an earth electromagnetic method; the shallow layer, the middle layer and the deep layer are divided according to depth;
the first inversion function establishing module is used for inverting the first ground resistivity data to obtain a first inversion target function;
the second inversion function establishing module is used for inverting the second ground resistivity data to obtain a second inversion target function;
a third inversion function establishing module, configured to invert the third ground resistivity data to obtain a third inversion target function;
a comprehensive inversion function establishing module, configured to obtain a comprehensive inversion target function according to the first inversion target function, the second inversion target function, and the third inversion target function;
and the optimization module is used for performing soil parameter inversion on the comprehensive inversion target function by adopting a differential evolution algorithm and taking the comprehensive inversion target function with the inverted soil parameters as a ground resistivity model.
A computer device comprising a memory storing a computer program and a processor implementing the steps of the above-described earth-resistivity model modeling method when the processor executes the computer program.
A computer-readable storage medium, having stored thereon a computer program which, when being executed by a processor, carries out the steps of the above-mentioned earth-resistivity model modelling method.
According to the method and the device for modeling the earth resistivity, the computer equipment and the storage medium, the earth resistivity data of different depths are obtained by respectively measuring the earth resistivity of different depths by adopting a quadrupole method, a controllable source audio frequency earth electromagnetic method and an earth electromagnetic method; then, mixing an inversion target function obtained by inverting the ground resistivity data obtained based on a quadrupole method, an inversion target function obtained by inverting the ground resistivity data obtained based on a controllable source audio ground method and an inversion target function obtained by inverting the ground resistivity obtained based on a ground electromagnetic method to obtain a comprehensive inversion target function; and finally, performing soil parameter inversion on the comprehensive inversion target function by adopting a differential evolution algorithm, and taking the comprehensive inversion target function with the inverted soil parameters as a ground resistivity model, thereby realizing uniform modeling on the wide-area ground resistivity covering from the ground surface to tens of kilometers underground.
Drawings
FIG. 1 is a diagram of an environment in which a modeling method for a georesistivity model may be applied in one embodiment;
FIG. 2 is a schematic flow diagram of a method for modeling a georesistivity model in one embodiment;
FIG. 3 is a schematic diagram of a quadrupole measurement;
FIG. 4 is a schematic diagram of the measurement of a controllable source audio magnetotelluric method;
FIG. 5 is a schematic view of a geoelectromagnetic method of measurement;
FIG. 6 is a flowchart illustrating the steps of a method for modeling a georesistivity model in one embodiment;
FIG. 7 is a schematic flow chart diagram illustrating steps of a method for modeling a georesistivity model in another embodiment;
FIG. 8 is a schematic diagram of measuring the pole address of an extra-high voltage DC grounding pole;
FIG. 9 is a schematic diagram of ground resistance measurement of a ground electrode;
FIG. 10 is a block diagram of the structure of the earth resistivity model modeling apparatus in one embodiment;
FIG. 11 is a block diagram of the structure of the synthetic inversion function building block in one embodiment;
FIG. 12 is a block diagram of the structure of an optimization module in one embodiment;
FIG. 13 is a diagram illustrating an internal structure of a computer device according to an embodiment.
Detailed Description
In order to make the objects, technical solutions and advantages of the present application more apparent, the present application is described in further detail below with reference to the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are merely illustrative of the present application and are not intended to limit the present application.
The modeling method of the earth resistivity model provided by the application can be applied to the application environment shown in FIG. 1. Wherein the terminal 102 communicates with the server 104 via a network. The user may send, via the mobile terminal 102, electric field and electromagnetic data that characterize earth resistivity obtained based on a quadrupole method, electric field and electromagnetic data that characterize earth resistivity obtained based on a controlled source audio earth method, and electric field and electromagnetic data that characterize earth resistivity obtained based on an earth electromagnetic method to the server 104. The server 104 performs the steps of the earth resistivity modeling method based on the electric field and electromagnetic data described above to model earth resistivity. The terminal 102 may be, but not limited to, various personal computers, notebook computers, smart phones, tablet computers, and portable wearable devices, and the server 104 may be implemented by an independent server or a server cluster formed by a plurality of servers.
In one embodiment, as shown in fig. 2, a method for modeling earth resistivity is provided, which is illustrated by applying the method to the server in fig. 1, and includes the following steps:
and S20, measuring the first earth resistivity data of the shallow layer by a quadrupole method.
And S40, measuring second earth resistivity data of the middle layer by adopting a controllable source audio frequency earth electromagnetic method.
And S60, measuring the third earth resistivity data of the deep layer by adopting an earth electromagnetic method.
S80, inverting the first ground resistivity data to obtain a first inversion objective function.
And S100, inverting the second ground resistivity data to obtain a second inversion target function.
And S120, inverting the third ground resistivity data to obtain a third inversion target function.
And S140, acquiring a comprehensive inversion target function according to the first inversion target function, the second inversion target function and the third inversion target function.
And S160, performing soil parameter inversion on the comprehensive inversion target function by adopting a differential evolution algorithm, and taking the comprehensive inversion target function with the inverted soil parameters as a ground resistivity model.
Wherein the shallow layer, the middle layer and the deep layer are obtained by dividing according to depth.
The quadrupole method is a measuring method for deducing underground resistivity distribution by utilizing space potential difference formed by current backflow, has higher sensitivity to shallow earth resistivity, and is therefore suitable for measuring shallow earth resistivity data.
Specifically, the measurement scheme of the quadrupole method is shown in FIG. 3, P1And P4Is a current pole, P2And P3The measuring center point is taken as a center, a current pole and a voltage pole are respectively arranged on two sides of the measuring center point, d is the distance between the electrodes, and b is the depth of the measuring pole.
The specific measuring process for measuring the shallow first earth resistivity data by adopting the quadrupole method comprises the following steps:
survey the measurement locations to determine a measurement center point.
Keeping the position of the measuring center point unchanged, wiring in the north-south direction (or arbitrarily selecting a certain direction) for polar distance diThe polar distance is selected according to the principle of 1-2-3-5-7.
Keeping the position of the measurement center point constant, the same measurements were made in an east-west (or perpendicular to the last measurement direction) wiring manner and a polar pitch arrangement, and the difference in apparent resistivity was observed in both directions.
Measurement m is completed according to the above steps1After the group polar distance is arranged, according to a first inversion target function F1Performing one-dimensional inversion of soil parameters:
where ρ isa(di) And ρm(di) The measurements are set i, the forward evolution of the apparent resistivity at depth di and the measured value. The calculation method of the forward value of apparent resistivity of the quadrupole method is disclosed by other documents and is not described in detail herein.
In addition, when the pole pitch exceeds 100m, the quadrupole method causes inaccurate measurement data due to the wire mutual inductance effect, and meanwhile, the quadrupole method is difficult to penetrate when encountering a high-resistance layer, so that the maximum measurement pole pitch of the quadrupole method is generally set to be about 100m, namely the maximum depth of measurement of the quadrupole method is 100m, and the earth resistivity cannot be measured when the depth exceeds 100 m.
The controllable source audio frequency geoelectromagnetic method is a frequency domain sounding method developed for overcoming the defect of poor signal strength of the geoelectromagnetic method, and has the characteristics of high signal-to-noise ratio and strong anti-interference capability. However, the controllable source audio frequency geoelectromagnetic method is greatly influenced by frequency, and is easily interfered when high-frequency shallow data is measured, so that the controllable source audio frequency geoelectromagnetic method is not as accurate as a quadrupole method for the shallow data. Generally, the depth measurement of the controllable source audio frequency geoelectromagnetic method can reach 2.5km-3km, so in the embodiment of the invention, a quadrupole method is adopted to measure the first georesistivity data of the shallow layer, and a controllable source audio frequency geoelectromagnetic method is adopted to measure the second georesistivity data of the middle layer.
In particular, FIG. 4 illustrates a measurement scheme of the controlled source audio magnetotelluric method.
The measuring step of measuring the second earth resistivity data of the middle layer by adopting the controllable source audio frequency earth electromagnetic method comprises the following steps:
the transmitter and the emission sources a and B are arranged and the electrodes and poles of the controlled source audio magnetotelluric method are arranged at the measurement point.
And determining the lowest measuring frequency of the pole addresses, thereby determining the maximum measuring depth.
Measuring the apparent resistivity corresponding to the measurement point with frequency fiVarying data, extracting electric field E in frequency range in receiving systemxAnd a magnetic field HyCalculating the frequency fiThe corresponding apparent resistivity:
wherein Z is the apparent resistivity, T is the period, fiThe frequency of the ith measurement.
Take out m2A characteristic point according to a second inversion target function F2Performing one-dimensional inversion of soil parameters:
where ρ isa(fi) And ρm(fi) Frequency f of the ith characteristic pointiThe calculation method of the forward value and the measured value of the corresponding apparent resistivity and the forward value of the apparent resistivity of the controllable source audio magnetotelluric method are disclosed by other documents and are not repeated herein.
And determining detection by the inversion result and the measurement data together, and judging whether to repeat the measurement of the steps or not through the detection.
And moving the measuring points and repeating the steps until the measurement of all the measuring points is completed.
Magnetotelluric methods theoretically can reach a probe depth of up to 100 kilometers, and the probe depth of the method depends on the lowest frequency that can be received by the measuring field device. However, since the method needs to measure the data of the natural magnetic field and the electric field corresponding to the low frequency, the shallow layer and the middle layer georesistivity data measured by the method can be easily submerged in the measurement site, and therefore, in the embodiment of the invention, only the third georesistivity data of the deep layer is measured by the geoelectromagnetic method.
Specifically, the arrangement of the geoelectromagnetic method is shown in fig. 5. The step of measuring deep earth resistivity data by using an earth electromagnetic method comprises the following steps:
arranging electrodes and magnetic poles at the measuring center point, respectively arranging non-polarized electrodes in four directions of east, south, west and north, taking the receiver as the center, and measuring the potential difference V in the horizontal direction x and the vertical direction yxAnd VyAnd arranges the magnetic poles in both x and y directions of the second and fourth quadrants as shown in the drawing to measure two quadrant components HxAnd HyRespectively, the components of the magnetic field. The electric field in the x and y directions measured by the non-polarized electrodes is calculated as follows:
wherein E isxElectric field in the x direction, EyElectricity in the y directionThe field, D, is the distance between two non-polarized electrodes in the same direction.
Monitoring electric field and magnetic field data in a period of time, extracting electric field and magnetic field in frequency range, and calculating frequency fiThe corresponding apparent resistivity:
wherein Z is1(fi) Is the apparent resistivity in the xy direction, and Z2(fi) Is the apparent resistivity in the yx direction.
Take out m3A characteristic point according to a third inversion target function F3Performing one-dimensional inversion of soil parameters:
where ρ isa(fi) And ρm(fi) Respectively being the ith characteristic point frequency fiA forward value and a measured value of the corresponding apparent resistivity. The calculation method of the forward value of the apparent resistivity of the magnetotelluric method is disclosed by other documents and is not described in detail herein.
And determining the depth measurement by the inversion result and the measurement data together, and judging whether the requirements are met or not through the depth measurement. If not, the above geoelectromagnetic method measuring steps are repeated.
Adding the first inversion target function, the second inversion target function and the third inversion target function to obtain a comprehensive inversion target function, wherein a comprehensive inversion target function minF can be expressed as:
minF＝F1+F2+F3
in one embodiment, as shown in fig. 6, the step of obtaining the synthetic inversion objective function according to the first inversion objective function, the second inversion objective function and the third inversion objective function includes:
and S141, configuring a first weight parameter of the first inversion target function, a second weight parameter of the second inversion target function and a third weight parameter of the third inversion target function according to the depth measurement.
The first weight parameter, the second weight parameter and the third weight parameter are configured according to the sounding, so that the value of the earth resistivity model for the resistivity at different depths is more biased to a method for measuring the depth more accurately.
And S142, carrying out weighted summation on the first inversion target function, the second inversion target function and the third inversion target function to obtain a comprehensive inversion target function.
Specifically, the first weight parameter, the second weight parameter and the third weight parameter can be configured by the following formula to obtain the comprehensive inversion objective function minF:
min F＝ω1F1+ω2F2+ω3F3
wherein, ω is1、ω2And ω3Respectively representing a first weight parameter, a second weight parameter and a third weight parameter.
The first weight parameter, the second weight parameter, and the third weight parameter are all related to depth sounding. Specifically, the first weight parameter is related to the pole pitch, and the second weight parameter and the third weight parameter are related to the frequency.
In one embodiment, the value of the first weight parameter decreases with increasing depth of investigation, the value of the second weight parameter increases first and then decreases with increasing depth of investigation, and the value of the third weight parameter increases with increasing depth of investigation, and remains unchanged when the depth of investigation reaches the maximum depth of investigation.
In one embodiment, the earth resistivity modeling method further comprises:
s130, configure the initial value of the first weight parameter as 1, configure the initial value of the second weight parameter as 0, and configure the initial value of the third weight parameter as 0.
In one embodiment, the step S141 of configuring the values of the first weight parameter, the second weight parameter and the third weight parameter according to the depth measurement includes:
when the depth measurement is gradually increased to the first threshold value, the first weight parameter is reduced, the second weight parameter is increased, and when the depth measurement is increased to the first threshold value, the first weight parameter is reduced to 0, and the second weight parameter is increased to 1;
when the depth measurement is gradually increased from the first threshold value to the second threshold value, the second weight parameter is reduced, and the third weight parameter is increased, and when the depth measurement is increased to the second threshold value, the second weight parameter is reduced to 0, and the third weight parameter is increased to 1.
The typical depth measurement range of the quadrupole method is 0-100m, the typical depth measurement range of the controllable source audio frequency geoelectromagnetic method is 10-2500 m, and the typical depth measurement range of the geoelectromagnetic method is 2.5-100 km. Therefore, typical depth measurement ranges of the quadrupole method, the controllable source audio frequency magnetotelluric method and the magnetotelluric method are overlapped to a certain extent, the overlapped part is called a transition layer, and for example, the transition layer between the shallow layer and the middle layer is 10 m-100 m. By setting the first threshold and the second threshold, the shallow transition layer result in the earth resistivity model is more biased to a quadrupole method with higher shallow measurement precision, the middle transition layer result is more biased to a controllable source audio earth electromagnetic method with higher middle measurement precision, and the deep transition layer result is more biased to an earth electromagnetic method with higher deep measurement precision, so that a more accurate wide-area earth resistivity model is finally obtained.
Specifically, in the shallow layer, the earth resistivity model mainly refers to the measurement result of the quadrupole method, the model output result is related to a first weight parameter and a second weight parameter, the initial value of the first weight parameter is set to 0, and the initial value of the second weight parameter is set to 1. With the increase of the sounding depth, the first weight parameter is gradually reduced, and the second weight parameter is gradually increased. When the sounding depth is increased to a first threshold value for distinguishing the shallow layer from the middle layer, the first weight parameter is reduced to 0, the second weight parameter is increased to 1, at the moment, the second weight parameter is gradually reduced along with the continuous increase of the sounding depth, the third weight parameter is increased from 0, the earth resistivity model mainly refers to the measurement result of the controllable source audio frequency earth electromagnetic method, and the model output result is related to the second weight parameter and the third weight parameter. When the sounding depth is increased to a second threshold value for distinguishing the middle layer from the deep layer, the second weight parameter is reduced to 0, the third weight parameter is increased to 1, at the moment, the third weight parameter is kept unchanged along with the continuous increase of the sounding depth, and the earth resistivity model mainly refers to the measurement result of the earth electromagnetic method. The shallow layer, the middle layer and the deep layer can be divided according to the measuring depth, and the first threshold value and the second threshold value are set according to the actual measuring environment.
In one embodiment, as shown in fig. 7, the step S160 of performing soil parameter inversion on the synthetic inversion target function by using a differential evolution algorithm includes:
s161, initializing a population by setting an inversion initial value of the comprehensive inversion target function.
Wherein, the population refers to the population formed by taking each soil parameter in the comprehensive inversion target function as an individual. The inversion initial value comprises a first weight parameter, a second weight parameter, a third weight parameter, a maximum inversion layer number, a depth, inversion accuracy and other parameters needing to be preset with the initial value.
Specifically, N D-dimensional vectors are generated in the search space and cover the entire search space.
Where G represents an evolution algebra, the initial population G ═ 0, and the individual s can be represented as:
and S162, iteratively executing variation operation on individuals in the contemporary population to generate variation individuals, performing cross operation on the contemporary population and the variation individuals to generate experimental individuals, selecting excellent individuals among the experimental individuals and the individuals in the contemporary population to form a next generation population, and outputting the soil parameters in the latest generation population as the soil parameters after inversion until a termination condition is met.
Specifically, the step of performing mutation operation on individuals in the contemporary population to generate variant individuals may be:
and randomly selecting three different individuals in the population, and obtaining the variant individual by overlapping the vector of any two individuals with the vector of the remaining individual after the vector of any two individuals is weighted differentially.
Specifically, the process of generating variant individuals can be represented as follows:
wherein r is1,r1,r1E (1, 2.. eta., N) and unlike i, k is a scaling factor.
And then performing cross operation on the current generation population and the variant individuals to generate experimental individuals.
The above-described interleaving operation can be expressed as shown in the following formula:
wherein,
for the variant, rand (j) is [0,1 ]]Random numbers are uniformly distributed among the random numbers; CR ∈ [0,1 ]]Is the cross probability; rnbr (i) is a random integer between {1, 2., D }.
Selecting excellent individuals from the experimental individuals and the individuals in the current generation population to form the next generation population.
And if the judgment result shows that the termination condition is met, outputting the soil parameters in the latest generation of population as the soil parameters after inversion.
If the algorithm reaches the preset inversion accuracy, the termination condition is considered to be met, and a feasible solution is output; and if the algorithm exceeds the preset operation amount and does not reach the preset inversion accuracy, determining that a termination condition is met, and outputting the current optimal solution. The preset inversion accuracies include an inversion accuracy F1 of the first inverted target function, an inversion accuracy F2 of the second inverted target function, an inversion accuracy F3 of the third inverted target function, and an inversion accuracy F of the synthetic inverted target function.
The preset operation amount comprises the iteration number, the running time and the like of the algorithm.
In one embodiment, F is1、F2And F3Both set to 6% and F set to 18%. When F is present1、F2And F3When the sum of the values is less than 6% and the F is less than 18%, the differential evolution algorithm is considered to meet the termination condition, and a feasible solution is output.
If the termination condition is not satisfied, S162 is iteratively executed.
Because the quadrupole method is an electrical method, the controllable source audio frequency magnetotelluric method and the magnetotelluric method are electromagnetic methods, the inversion theories of the electrical method and the electromagnetic methods cannot be applied across methods, and the normalization of different data to a uniform size is the basis of joint inversion.
The differential evolution algorithm has the advantages of fast convergence, few control parameters, simple setting, good robustness and the like. Therefore, the soil parameter inversion is carried out on the comprehensive inversion target function by adopting the differential evolution algorithm, so that the efficiency of the modeling process of the earth resistivity model is higher, namely the convergence of the algorithm can be realized through fewer iteration times.
It should be understood that, although the steps in the flowcharts of fig. 2 or 6 or 7 are shown in sequence as indicated by the arrows, the steps are not necessarily performed in sequence as indicated by the arrows. The steps are not performed in the exact order shown and described, and may be performed in other orders, unless explicitly stated otherwise. Moreover, at least a portion of the steps in fig. 2, 6 or 7 may include multiple steps or multiple stages, which are not necessarily performed at the same time, but may be performed at different times, which are not necessarily performed in sequence, but may be performed alternately or alternately with other steps or at least a portion of the steps or stages in other steps.
In addition, the earth resistivity model modeling method provided in the embodiment of the application can be applied to earth resistivity modeling of the pole address of the extra-high voltage direct current grounding pole.
The direct current grounding electrode is used as an important component in an extra-high voltage direct current transmission system, plays a role in clamping neutral point potential in system operation, and provides ground channels for rated current in monopolar earth operation or unbalanced current in bipolar asymmetric operation. When the direct current system is operated in a single-pole ground, strong rated system current flows into the ground through the grounding electrode to raise the ground potential, and meanwhile, the current is scattered in the ground to cause heating of surrounding soil, so that the stable operation of the grounding electrode is seriously influenced when the temperature is too high. Therefore, the selection of the grounding electrode address is very important in the DC engineering design stage. Due to the strict requirements of the critical performance indicators of the dc grounding electrode, it is necessary to have a sufficiently large current spreading area of the grounding electrode. The current land resource is in shortage, and the size of most conventional direct current grounding electrodes is large, so that the problems of difficult site selection and difficult land acquisition of the direct current grounding electrodes are highlighted day by day.
Therefore, it is a feasible idea to extend the size of the ground electrode in the vertical direction in practical engineering. When the electrode site deep layer ground resistivity is lower, the current can be led to the underground deep, certain advantages are achieved in the aspects of reducing the ground resistance and the step voltage, meanwhile, the occupied area is reduced, and the method is an effective ground electrode arrangement scheme. The burial depth of the direct current grounding electrode can be in the range of hundreds to kilometers at present, and the electrical and thermal characteristics of the grounding electrode can be determined to a great extent by the electrode site earth resistivity in the range of 10 times the burial depth.
However, at present, the traditional earth resistivity modeling method does not consider the distribution rule of earth resistivity at different depths, and cannot accurately and uniformly model the earth resistivity of a wide area covering tens of kilometers from the earth surface to the underground. Therefore, the embodiment of the application provides a combined measuring method and a ground resistivity modeling method which comprehensively use the geodetic method, the controllable source audio geodetic method and the geodetic electromagnetic method, so that the wide-area ground resistivity covering tens of kilometers from the earth surface to the underground can be more accurately and uniformly modeled.
Specifically, the details of the embodiment of the present invention are explained in detail with an extra-high voltage dc ground electrode address.
The ground electrode address requires detection of earth resistivity of at least 10km depth. According to the three methods with different advantageous detection depths, the earth resistivity test is divided into three sections with different depths, and different methods are adopted for testing. The quadrupole method is used for measuring earth resistivity data of 0m-100m, the controlled source audio frequency earth electromagnetic method is used for measuring earth resistivity data of 0m-3km, and the earth electromagnetic method is used for measuring earth resistivity data of 0-68 km. A schematic view of the pole site measurement is shown in fig. 8. Due to the limitation of field measurement conditions, the wiring of the transverse, longitudinal and quadrupole method adopts a three-longitudinal-one-transverse method, namely S1-S4 in the figure, and the maximum pole pitch is 100 m. The controllable source audio frequency geoelectromagnetic method adopts a straight line method arrangement form, namely L2-L6 in the figure. L2 is 70m away from L3, the distance between every two of L3-L6 is 50m, the length of each measuring line is 700m, and the measuring point interval in each line is 12.5 m. Essentially covering the ground electrode address area. It is difficult to arrange all the geoelectromagnetic measuring points in a straight line according to the field disturbance situation. The geoelectromagnetic method has four measuring points which are M1-M4 in the figure.
The method in the embodiment of the application is used for obtaining the comprehensive inversion target function based on the shallow first earth resistivity data, the middle second earth resistivity data and the deep third earth resistivity data which are obtained by the quadrupole method, the controllable source audio frequency earth electromagnetic method and the earth electromagnetic method through measurement, and performing soil parameter inversion on the comprehensive inversion target function by adopting a differential evolution algorithm. Setting the initial value of the number of the large ground layers in the differential evolution algorithm to be 8, and if the minF is met<18％，F1<6％，F2<6% and F3<6% "outputs feasible solution, if not, minF<18％，F1<6％，F2<6% and F3<6% ", the soil layer number is added with 1 to continue inversion until the preset calculation is metThe process stop conditions. According to the inversion idea, inversion is carried out on a single measuring point.
The first weight parameter, the second weight parameter and the third weight parameter are configured according to the sounding, so that the earth resistivity model is more biased to data of a quadrupole method for the depth of 0-100m, the audio earth electromagnetic method for the controllable source is more biased to the depth of 100m-3000m, the earth electromagnetic method is more biased below the depth of 3km, and finally the wide-area earth resistivity model of the deep well grounding electrode address is obtained, as shown in table 1.
TABLE 1 Wide-area earth resistivity model for deep-well earth electrode site
In order to verify the accuracy of the earth resistivity measuring method and the modeling method provided by the embodiment of the application, an earth resistance test is carried out on the direct current earth electrode. And carrying out grounding resistance test work of the grounding electrode by adopting a three-pole method. The measurement uses a 220V alternating current power supply, the rectifier converts alternating current into direct current, and the maximum output is 7A. The voltage is measured at the tail end of the deep well grounding electrode drainage cable, the measurement schematic diagram is shown in fig. 9, and the measurement result is shown in table 2.
TABLE 2 measurement of grounding resistance of deep well grounding electrode
The earth resistivity model obtained in the embodiment of the application is used for modeling and calculating the earth resistance of the earth electrode in CDEGS simulation software, the calculation result is 0.138 omega, and compared with an actual measurement value, the error is only 0.036 omega, so that the reliability of the earth resistivity model obtained by the earth resistivity modeling method in the embodiment of the application is proved.
In one embodiment, as shown in fig. 10, there is provided a ground resistivity model modeling apparatus including: the device comprises a measuring module, a comprehensive inversion function establishing module and an optimizing module, wherein:
the first earth resistivity acquisition module 20 is configured to measure shallow first earth resistivity data by a quadrupole method;
the second earth resistivity acquisition module 40 is used for measuring second earth resistivity data of the middle layer by adopting a controllable source audio frequency earth electromagnetic method;
a third earth resistivity obtaining module 60, configured to measure deep third earth resistivity data by using an earth electromagnetic method; the shallow layer, the middle layer and the deep layer are divided according to depth;
a first inversion function establishing module 80, configured to invert the first ground resistivity data to obtain a first inversion target function;
a second inversion function establishing module 100, configured to invert the second ground resistivity data to obtain a second inversion target function;
a third inversion function establishing module 120, configured to invert the third ground resistivity data to obtain a third inversion target function;
a comprehensive inversion function establishing module 140, configured to obtain a comprehensive inversion target function according to the first inversion target function, the second inversion target function, and the third inversion target function;
and the optimization module 160 is configured to perform soil parameter inversion on the comprehensive inversion target function by using a differential evolution algorithm, and use the comprehensive inversion target function with the inverted soil parameters as a ground resistivity model.
In one embodiment, as shown in FIG. 11, the synthetic inversion function building block includes:
the configuration unit 141 is configured to configure a first weight parameter of the first inversion target function, a second weight parameter of the second inversion target function, and a third weight parameter of the third inversion target function according to the depth measurement.
And a superposition unit 142, configured to perform weighted summation on the first inversion target function, the second inversion target function, and the third inversion target function, so as to obtain a comprehensive inversion target function.
In one embodiment, the synthetic inversion function building module further comprises:
a weight initial value setting unit 130, configured to configure the initial value of the first weight parameter as 1, the initial value of the second weight parameter as 0, and the initial value of the third weight parameter as 0.
In one embodiment, the configuration unit 141 includes:
the weight updating unit is used for decreasing the first weight parameter and increasing the second weight parameter when the depth measurement is gradually increased to be before the first threshold value, and when the depth measurement is increased to be before the first threshold value, the first weight parameter is decreased to 0 and the second weight parameter is increased to 1;
when the depth of sounding gradually increases within a range from a first threshold value to a second threshold value, the second weight parameter is decreased, and the third weight parameter is increased, and when the depth of sounding increases to the second threshold value, the second weight parameter is decreased to 0, and the third weight parameter increases to 1.
In one embodiment, as shown in FIG. 12, the optimization module 160 includes:
an initializing unit 161, configured to initialize a population by setting an inversion initial value of the synthetic inversion target function.
And the iterative execution unit 162 is configured to iteratively execute the steps of performing a variation operation on individuals in the current generation population to generate variation individuals, performing a crossover operation on the current generation population and the variation individuals to generate experimental individuals, and selecting excellent individuals from the experimental individuals and the individuals in the current generation population to form a next generation population until a termination condition is met, and outputting a soil parameter in the latest generation population as an inverted soil parameter.
In one embodiment, the variant individual generating unit includes:
and the differential weighting unit is used for randomly selecting three different individuals in the population, and obtaining the variant individual by superposing the vector of any two individuals with the vector of the remaining individual after differential weighting.
For specific definition of the earth resistivity model modeling device, reference may be made to the above definition of the earth resistivity model modeling method, which is not described herein again. The respective modules in the earth resistivity model modeling apparatus described above may be implemented in whole or in part by software, hardware, and a combination thereof. The modules can be embedded in a hardware form or independent from a processor in the computer device, and can also be stored in a memory in the computer device in a software form, so that the processor can call and execute operations corresponding to the modules.
In one embodiment, a computer device is provided, which may be a terminal, and its internal structure diagram may be as shown in fig. 13. The computer device includes a processor, a memory, a communication interface, a display screen, and an input device connected by a system bus. Wherein the processor of the computer device is configured to provide computing and control capabilities. The memory of the computer device comprises a nonvolatile storage medium and an internal memory. The non-volatile storage medium stores an operating system and a computer program. The internal memory provides an environment for the operation of an operating system and computer programs in the non-volatile storage medium. The communication interface of the computer device is used for carrying out wired or wireless communication with an external terminal, and the wireless communication can be realized through WIFI, an operator network, NFC (near field communication) or other technologies. The computer program is executed by a processor to implement a method of modeling a georesistivity model. The display screen of the computer equipment can be a liquid crystal display screen or an electronic ink display screen, and the input device of the computer equipment can be a touch layer covered on the display screen, a key, a track ball or a touch pad arranged on the shell of the computer equipment, an external keyboard, a touch pad or a mouse and the like.
Those skilled in the art will appreciate that the architecture shown in fig. 13 is merely a block diagram of some of the structures associated with the disclosed aspects and is not intended to limit the computing devices to which the disclosed aspects apply, as particular computing devices may include more or less components than those shown, or may combine certain components, or have a different arrangement of components.
In one embodiment, a computer device is provided, comprising a memory and a processor, the memory having a computer program stored therein, the processor implementing the following steps when executing the computer program:
and S20, measuring the first earth resistivity data of the shallow layer by a quadrupole method.
And S40, measuring second earth resistivity data of the middle layer by adopting a controllable source audio frequency earth electromagnetic method.
And S60, measuring the third earth resistivity data of the deep layer by adopting an earth electromagnetic method.
S80, inverting the first ground resistivity data to obtain a first inversion objective function.
And S100, inverting the second ground resistivity data to obtain a second inversion target function.
And S120, inverting the third ground resistivity data to obtain a third inversion target function.
And S140, acquiring a comprehensive inversion target function according to the first inversion target function, the second inversion target function and the third inversion target function.
And S160, performing soil parameter inversion on the comprehensive inversion target function by adopting a differential evolution algorithm, and taking the comprehensive inversion target function with the inverted soil parameters as a ground resistivity model.
Wherein the shallow layer, the middle layer and the deep layer are obtained by dividing according to depth.
In one embodiment, the processor, when executing the computer program, further performs the steps of:
s130, configure the initial value of the first weight parameter as 1, configure the initial value of the second weight parameter as 0, and configure the initial value of the third weight parameter as 0.
In one embodiment, the processor, when executing the computer program, further performs the steps of:
s141, configuring a first weight parameter of the first inversion target function, a second weight parameter of the second inversion target function and a third weight parameter of the third inversion target function.
And S142, carrying out weighted summation on the first inversion target function, the second inversion target function and the third inversion target function to obtain a comprehensive inversion target function.
In one embodiment, the processor, when executing the computer program, further performs the steps of:
s161, initializing a population by setting an inversion initial value of the comprehensive inversion target function.
And S162, iteratively executing variation operation on individuals in the contemporary population to generate variation individuals, performing cross operation on the contemporary population and the variation individuals to generate experimental individuals, selecting excellent individuals among the experimental individuals and the individuals in the contemporary population to form a next generation population, and outputting the soil parameters in the latest generation population as the soil parameters after inversion until a termination condition is met.
If it is determined that the termination condition is not satisfied, S162 is performed.
In one embodiment, a computer-readable storage medium is provided, having a computer program stored thereon, which when executed by a processor, performs the steps of:
and S20, measuring the first earth resistivity data of the shallow layer by a quadrupole method.
And S40, measuring second earth resistivity data of the middle layer by adopting a controllable source audio frequency earth electromagnetic method.
And S60, measuring the third earth resistivity data of the deep layer by adopting an earth electromagnetic method.
S80, inverting the first ground resistivity data to obtain a first inversion objective function.
And S100, inverting the second ground resistivity data to obtain a second inversion target function.
And S120, inverting the third ground resistivity data to obtain a third inversion target function.
And S140, acquiring a comprehensive inversion target function according to the first inversion target function, the second inversion target function and the third inversion target function.
And S160, performing soil parameter inversion on the comprehensive inversion target function by adopting a differential evolution algorithm, and taking the comprehensive inversion target function with the inverted soil parameters as a ground resistivity model.
Wherein the shallow layer, the middle layer and the deep layer are obtained by dividing according to depth.
In one embodiment, the computer program when executed by the processor further performs the steps of:
s141, configuring a first weight parameter of the first inversion target function, a second weight parameter of the second inversion target function and a third weight parameter of the third inversion target function.
And S142, carrying out weighted summation on the first inversion target function, the second inversion target function and the third inversion target function to obtain a comprehensive inversion target function.
In one embodiment, the computer program when executed by the processor further performs the steps of:
s130, configure the initial value of the first weight parameter as 1, configure the initial value of the second weight parameter as 0, and configure the initial value of the third weight parameter as 0.
In one embodiment, the computer program when executed by the processor further performs the steps of:
s161, initializing a population by setting an inversion initial value of the comprehensive inversion target function.
And S162, iteratively executing variation operation on individuals in the contemporary population to generate variation individuals, performing cross operation on the contemporary population and the variation individuals to generate experimental individuals, selecting excellent individuals among the experimental individuals and the individuals in the contemporary population to form a next generation population, and outputting the soil parameters in the latest generation population as the soil parameters after inversion until a termination condition is met.
If it is determined that the termination condition is not satisfied, S162 is performed.
It will be understood by those skilled in the art that all or part of the processes of the methods of the embodiments described above can be implemented by hardware instructions of a computer program, which can be stored in a non-volatile computer-readable storage medium, and when executed, can include the processes of the embodiments of the methods described above. Any reference to memory, storage, database or other medium used in the embodiments provided herein can include at least one of non-volatile and volatile memory. Non-volatile Memory may include Read-Only Memory (ROM), magnetic tape, floppy disk, flash Memory, optical storage, or the like. Volatile Memory can include Random Access Memory (RAM) or external cache Memory. By way of illustration and not limitation, RAM can take many forms, such as Static Random Access Memory (SRAM) or Dynamic Random Access Memory (DRAM), among others.
The technical features of the above embodiments can be arbitrarily combined, and for the sake of brevity, all possible combinations of the technical features in the above embodiments are not described, but should be considered as the scope of the present specification as long as there is no contradiction between the combinations of the technical features.
The above-mentioned embodiments only express several embodiments of the present application, and the description thereof is more specific and detailed, but not construed as limiting the scope of the invention. It should be noted that, for a person skilled in the art, several variations and modifications can be made without departing from the concept of the present application, which falls within the scope of protection of the present application. Therefore, the protection scope of the present patent shall be subject to the appended claims.
Claims (10)
1. A method of modeling a georesistivity model, the method comprising:
measuring first earth resistivity data of a shallow layer by adopting a quadrupole method;
measuring second earth resistivity data of the middle layer by adopting a controllable source audio magnetotelluric method;
measuring third earth resistivity data of a deep layer by adopting an earth electromagnetic method; the shallow layer, the middle layer and the deep layer are divided according to depth;
inverting the first earth resistivity data to obtain a first inversion target function;
inverting the second earth resistivity data to obtain a second inversion target function;
inverting the third ground resistivity data to obtain a third inversion target function;
acquiring a comprehensive inversion target function according to the first inversion target function, the second inversion target function and the third inversion target function;
and performing soil parameter inversion on the comprehensive inversion target function by adopting a differential evolution algorithm, and taking the comprehensive inversion target function with the inverted soil parameters as a ground resistivity model.
2. The method of claim 1, wherein the step of performing soil parameter inversion on the synthetic inversion objective function using a differential evolution algorithm comprises:
initializing a population by setting an inversion initial value of the comprehensive inversion target function, wherein the population is a population formed by taking each soil parameter in the comprehensive inversion target function as an individual;
iteratively executing variation operation on individuals in the current generation population to generate variation individuals, performing cross operation on the current generation population and the variation individuals to generate experimental individuals, selecting excellent individuals between the experimental individuals and the individuals in the current generation population to form a next generation population, and outputting soil parameters in the latest generation population as the soil parameters after inversion until termination conditions are met.
3. The method of claim 2, wherein the step of performing mutation operations on individuals in the contemporary population to generate variant individuals comprises:
and randomly selecting three different individuals in the population, and obtaining the variant individual by overlapping the vector of any two individuals with the vector of the remaining individual after the vector of any two individuals is weighted differentially.
4. The method of any of claims 1-3, wherein obtaining a synthetic inverted objective function from the first, second, and third inverted objective functions comprises:
configuring a first weight parameter of the first inversion target function, a second weight parameter of the second inversion target function and a third weight parameter of the third inversion target function according to depth measurement;
and carrying out weighted summation on the first inversion target function, the second inversion target function and the third inversion target function to obtain the comprehensive inversion target function.
5. The method of claim 4, wherein the value of the first weight parameter decreases with increasing depth of sounding;
the value of the second weight parameter is increased and then decreased along with the increase of sounding depth;
and when the value of the second weight parameter is reduced, the value of the third weight parameter is increased along with the increase of the sounding depth and is kept unchanged when the sounding depth reaches the maximum sounding depth.
6. The method of claim 4, further comprising:
configuring an initial value of the first weight parameter to 1; configuring an initial value of the second weight parameter to 0; configuring an initial value of the third weight parameter to be 0.
7. The method of claim 6, wherein the step of configuring the values of the first, second and third weight parameters according to depth of observation comprises:
when the depth measurement is gradually increased to a first threshold value, decreasing the first weight parameter and increasing the second weight parameter, and when the depth measurement is increased to the first threshold value, decreasing the first weight parameter to 0 and increasing the second weight parameter to 1;
when the depth of sounding gradually increases within a range from a first threshold value to a second threshold value, the second weight parameter is decreased, and the third weight parameter is increased, and when the depth of sounding increases to the second threshold value, the second weight parameter is decreased to 0, and the third weight parameter increases to 1.
8. A georesistivity model modeling apparatus, the apparatus comprising:
the first earth resistivity acquisition module is used for measuring first earth resistivity data of a shallow layer by adopting a quadrupole method;
the second earth resistivity acquisition module is used for measuring second earth resistivity data of the middle layer by adopting a controllable source audio frequency earth electromagnetic method;
the third earth resistivity acquisition module is used for measuring deep third earth resistivity data by adopting an earth electromagnetic method; the shallow layer, the middle layer and the deep layer are divided according to depth;
the first inversion function establishing module is used for inverting the first ground resistivity data to obtain a first inversion target function;
the second inversion function establishing module is used for inverting the second ground resistivity data to obtain a second inversion target function;
a third inversion function establishing module, configured to invert the third ground resistivity data to obtain a third inversion target function;
a comprehensive inversion function establishing module, configured to obtain a comprehensive inversion target function according to the first inversion target function, the second inversion target function, and the third inversion target function;
and the optimization module is used for performing soil parameter inversion on the comprehensive inversion target function by adopting a differential evolution algorithm and taking the comprehensive inversion target function with the inverted soil parameters as a ground resistivity model.
9. A computer device comprising a memory and a processor, the memory storing a computer program, characterized in that the processor, when executing the computer program, implements the steps of the method of any of claims 1 to 7.
10. A computer-readable storage medium, on which a computer program is stored, which, when being executed by a processor, carries out the steps of the method of any one of claims 1 to 7.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111304750.4A CN114047554B (en) | 2021-11-05 | 2021-11-05 | Earth resistivity model modeling method, apparatus, computer device and storage medium |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111304750.4A CN114047554B (en) | 2021-11-05 | 2021-11-05 | Earth resistivity model modeling method, apparatus, computer device and storage medium |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114047554A true CN114047554A (en) | 2022-02-15 |
| CN114047554B CN114047554B (en) | 2024-04-02 |
Family
ID=80207322
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111304750.4A Active CN114047554B (en) | 2021-11-05 | 2021-11-05 | Earth resistivity model modeling method, apparatus, computer device and storage medium |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114047554B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117270072A (en) * | 2023-09-19 | 2023-12-22 | 云南大学 | Gravity magnetic potential field imaging inversion method and system based on improved differential evolution algorithm |
Citations (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20080059075A1 (en) * | 2006-09-04 | 2008-03-06 | Daniele Colombo | Methods and apparatus for geophysical exploration via joint inversion |
| US20090083006A1 (en) * | 2007-09-20 | 2009-03-26 | Randall Mackie | Methods and apparatus for three-dimensional inversion of electromagnetic data |
| CN102879828A (en) * | 2012-10-11 | 2013-01-16 | 甘肃铁道综合工程勘察院有限公司 | Audio magnetotelluric optimized collection method |
| CN103869173A (en) * | 2014-02-26 | 2014-06-18 | 国家电网公司 | Method for measuring earth resistivity distribution from earth surface to underground tens of kilometers |
| CN106094045A (en) * | 2016-06-16 | 2016-11-09 | 华北电力大学 | A kind of method utilizing mt 3-d inversion data to set up horizontal soil model |
| CN106443189A (en) * | 2016-12-05 | 2017-02-22 | 国家电网公司 | Method and system for three-dimensional detection on earth electrode field and surrounding soil resistivity |
| CN106526329A (en) * | 2016-11-10 | 2017-03-22 | 广东电网有限责任公司电力科学研究院 | Method and device for measuring resistivity from land surface to deep earth |
| CN107630697A (en) * | 2017-09-26 | 2018-01-26 | 长江大学 | Based on the formation resistivity joint inversion method for boring electromagnetic wave resistivity logging |
| CN110007357A (en) * | 2019-05-16 | 2019-07-12 | 核工业航测遥感中心 | A kind of aviation TEM and aviation MT joint inversion method |
| CN111965712A (en) * | 2020-10-21 | 2020-11-20 | 国网江西省电力有限公司电力科学研究院 | Method for correcting static effect of controllable source audio magnetotelluric method |
| CN112130215A (en) * | 2019-06-24 | 2020-12-25 | 中国石油天然气集团有限公司 | Electromagnetic exploration data processing method and device |
| CN113297526A (en) * | 2021-06-03 | 2021-08-24 | 国家电网有限公司 | Horizontal layered soil structure joint inversion method based on Wenner quadrupole and magnetotelluric data |
-
2021
- 2021-11-05 CN CN202111304750.4A patent/CN114047554B/en active Active
Patent Citations (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20080059075A1 (en) * | 2006-09-04 | 2008-03-06 | Daniele Colombo | Methods and apparatus for geophysical exploration via joint inversion |
| US20090083006A1 (en) * | 2007-09-20 | 2009-03-26 | Randall Mackie | Methods and apparatus for three-dimensional inversion of electromagnetic data |
| CN102879828A (en) * | 2012-10-11 | 2013-01-16 | 甘肃铁道综合工程勘察院有限公司 | Audio magnetotelluric optimized collection method |
| CN103869173A (en) * | 2014-02-26 | 2014-06-18 | 国家电网公司 | Method for measuring earth resistivity distribution from earth surface to underground tens of kilometers |
| CN106094045A (en) * | 2016-06-16 | 2016-11-09 | 华北电力大学 | A kind of method utilizing mt 3-d inversion data to set up horizontal soil model |
| CN106526329A (en) * | 2016-11-10 | 2017-03-22 | 广东电网有限责任公司电力科学研究院 | Method and device for measuring resistivity from land surface to deep earth |
| CN106443189A (en) * | 2016-12-05 | 2017-02-22 | 国家电网公司 | Method and system for three-dimensional detection on earth electrode field and surrounding soil resistivity |
| CN107630697A (en) * | 2017-09-26 | 2018-01-26 | 长江大学 | Based on the formation resistivity joint inversion method for boring electromagnetic wave resistivity logging |
| CN110007357A (en) * | 2019-05-16 | 2019-07-12 | 核工业航测遥感中心 | A kind of aviation TEM and aviation MT joint inversion method |
| CN112130215A (en) * | 2019-06-24 | 2020-12-25 | 中国石油天然气集团有限公司 | Electromagnetic exploration data processing method and device |
| CN111965712A (en) * | 2020-10-21 | 2020-11-20 | 国网江西省电力有限公司电力科学研究院 | Method for correcting static effect of controllable source audio magnetotelluric method |
| CN113297526A (en) * | 2021-06-03 | 2021-08-24 | 国家电网有限公司 | Horizontal layered soil structure joint inversion method based on Wenner quadrupole and magnetotelluric data |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117270072A (en) * | 2023-09-19 | 2023-12-22 | 云南大学 | Gravity magnetic potential field imaging inversion method and system based on improved differential evolution algorithm |
| CN117270072B (en) * | 2023-09-19 | 2024-04-19 | 云南大学 | Gravity magnetic potential field imaging inversion method and system based on improved differential evolution algorithm |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114047554B (en) | 2024-04-02 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN112949134B (en) | Earth-well transient electromagnetic inversion method based on non-structural finite element method | |
| CN110346654B (en) | Electromagnetic spectrum map construction method based on common kriging interpolation | |
| CN107766666B (en) | Three-dimensional time domain electromagnetic anomalous diffusion simulation method based on fractional order differential method | |
| CN106526329B (en) | It is a kind of measurement earth's surface to deep layer ground resistivity method and device | |
| CN113051779B (en) | Numerical simulation method of three-dimensional direct-current resistivity method | |
| CN115238550A (en) | Self-adaptive unstructured grid landslide rainfall geoelectric field numerical simulation calculation method | |
| CN114047554B (en) | Earth resistivity model modeling method, apparatus, computer device and storage medium | |
| CN109343134A (en) | A kind of Transient Electromagnetic Method in Mine data analysis interpretation method and system | |
| CN108388750B (en) | Novel radial grounding resistance calculation method | |
| CN113406707A (en) | Magnetotelluric multi-scale and multi-time-period detection method | |
| Jing et al. | Research on earth resistivity measuring and modeling of HVDC deep-well grounding electrode sites | |
| CN113297526B (en) | Horizontal layered soil structure joint inversion method based on Wenner quadrupole and magnetotelluric data | |
| CN113204905B (en) | Numerical simulation method for finite element by contact induced polarization method | |
| CN112882124B (en) | Method for delineating beneficial parts of three-dimensional uranium ore formation | |
| CN110820814B (en) | Pile foundation detection device and method | |
| CN112462432B (en) | High-precision detection method and device for fire area cavity of coal field | |
| Gomes et al. | Embedded system to grounding grid diagnosis of energized substations | |
| CN113238285A (en) | Resistivity calculation method and system for geophysical charging method exploration and terminal | |
| CN114488327A (en) | Horizontal magnetic field and vertical magnetic field combined measurement method in well based on ground base point | |
| CN113868919A (en) | Simplified method for electromagnetic wave logging while drilling 3D simulation | |
| CN113221320A (en) | Method for rapidly calculating geomagnetic storm induction ground electric field based on apparent resistivity | |
| Brignone et al. | An algorithm for the semi-analytical computation of fields emitted in layered ground by HVDC electrodes | |
| Jing et al. | The research on new type of Earth resistivity exploring method for HVDC deep-well Earth electrode | |
| CN114152989B (en) | Method for calculating coordinates of pseudo-section diagram of two-dimensional asymmetric quadrupole electrical method | |
| CN112632206B (en) | Lightning feature analysis method and device, computer equipment and storage medium |
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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |