CN114047554B - Earth resistivity model modeling method, apparatus, computer device and storage medium - Google Patents
Earth resistivity model modeling method, apparatus, computer device and storage medium Download PDFInfo
- Publication number
- CN114047554B CN114047554B CN202111304750.4A CN202111304750A CN114047554B CN 114047554 B CN114047554 B CN 114047554B CN 202111304750 A CN202111304750 A CN 202111304750A CN 114047554 B CN114047554 B CN 114047554B
- Authority
- CN
- China
- Prior art keywords
- objective function
- inversion
- weight parameter
- earth
- inversion objective
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 159
- 230000006870 function Effects 0.000 claims abstract description 182
- 239000002689 soil Substances 0.000 claims abstract description 46
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 25
- 238000005259 measurement Methods 0.000 claims description 59
- 238000004590 computer program Methods 0.000 claims description 20
- 239000013598 vector Substances 0.000 claims description 12
- 230000003247 decreasing effect Effects 0.000 claims description 10
- 230000007423 decrease Effects 0.000 claims description 9
- 230000007704 transition Effects 0.000 claims description 9
- 230000035772 mutation Effects 0.000 claims description 8
- 238000005457 optimization Methods 0.000 claims description 6
- 238000010586 diagram Methods 0.000 description 11
- 230000005684 electric field Effects 0.000 description 11
- 238000001514 detection method Methods 0.000 description 5
- 238000012360 testing method Methods 0.000 description 5
- 238000004891 communication Methods 0.000 description 4
- 230000008569 process Effects 0.000 description 3
- 238000009933 burial Methods 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000000691 measurement method Methods 0.000 description 2
- 230000003068 static effect Effects 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 1
- 238000012512 characterization method Methods 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
- 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
- 230000005404 monopole Effects 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
- 230000035945 sensitivity Effects 0.000 description 1
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)
- Measurement Of Resistance Or Impedance (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The application relates to a method, a device, computer equipment and a storage medium for modeling a ground resistivity model. The method comprises the following steps: the method comprises the steps of respectively obtaining first earth resistivity data of a shallow layer by a four-pole method, measuring second earth resistivity data of a middle layer by a controllable source audio earth electromagnetic method and measuring third earth resistivity data of a deep layer by the earth electromagnetic method; inverting the first earth resistivity data to obtain a first inversion objective function, inverting the second earth resistivity data to obtain a second inversion objective function, and inverting the third earth resistivity data to obtain a third inversion objective function; and obtaining a comprehensive inversion objective function according to the first, second and third inversion objective functions, inverting the comprehensive inversion objective function by adopting a differential evolution algorithm, and taking the comprehensive inversion objective function with the inverted soil parameters as a ground resistivity model. The method can accurately perform unified modeling on wide area earth resistivity covering tens of kilometers from the surface to the underground.
Description
Technical Field
The application relates to the technical field of grounding of power systems, in particular to a method and a device for modeling earth resistivity, computer equipment and a storage medium.
Background
With the development of the field of grounding technology of the power system, extending the size of the grounding electrode towards the vertical direction is a feasible idea, and when the deep ground resistivity of the electrode address is low, current can be led to the deep ground. However, the earth electrode has a large buried depth and a wide range of stratum, so that the earth resistivity is very complex to explore. In order to better study the distribution rule of the earth resistivity, an earth resistivity modeling method is proposed. However, the conventional earth resistivity modeling method does not consider the distribution rule of earth resistivity at different depths, and cannot accurately model the wide area earth resistivity covering tens of kilometers from the surface to the ground in a unified manner.
Disclosure of Invention
In view of the foregoing, it is desirable to provide a method, apparatus, computer device, and storage medium capable of accurately modeling a wide area earth resistivity model covering tens of kilometers from the surface to the subsurface.
A method of modeling a earth resistivity model, the method comprising:
Measuring first earth resistivity data of the shallow layer by using a quadrupole method;
measuring second earth resistivity data of the middle layer by adopting a controllable source audio frequency earth electromagnetic method;
measuring third earth resistivity data of the 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 objective function;
inverting the second earth resistivity data to obtain a second inversion objective function;
inverting the third earth resistivity data to obtain a third reflection objective function;
acquiring a comprehensive inversion objective function according to the first inversion objective function, the second inversion objective function and the third inversion objective function;
and carrying out soil parameter inversion on the comprehensive inversion objective function by adopting a differential evolution algorithm, and taking the comprehensive inversion objective function with the inverted soil parameters as a ground resistivity model.
In one embodiment, the step of inverting the soil parameters of the synthetic inversion objective function using a differential evolution algorithm includes:
initializing a population by setting an inversion initial value of the comprehensive inversion objective function, wherein the population refers to the population formed by taking each soil parameter in the comprehensive inversion objective function as an individual;
And iteratively executing the steps of performing mutation operation on individuals in the current population to generate variant individuals, performing cross operation on the current population and the variant individuals to generate experimental individuals, selecting excellent individuals between the experimental individuals and the individuals in the current population to form a next-generation population, and outputting soil parameters in the latest-generation population as soil parameters after inversion until termination conditions are met.
In one embodiment, the step of performing a mutation operation on individuals in the current population to generate mutated individuals includes:
three different individuals in the population are randomly selected, and the vectors of any two individuals are differentially weighted and then are overlapped with the vector of the rest one individual to obtain a variant individual.
In one embodiment, the step of obtaining the synthetic inversion objective function from the first inversion objective function, the second inversion objective function, and the third inversion objective function comprises:
configuring a first weight parameter of a first inversion objective function, a second weight parameter of a second inversion objective function and a third weight parameter of a third inversion objective function;
and carrying out weighted summation on the first inversion objective function, the second inversion objective function and the third inversion objective function to obtain the comprehensive inversion objective function.
In one embodiment, the value of the first weight parameter decreases as the sounding increases; the value of the second weight parameter increases and decreases with the increase of the sounding; 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, and the sounding is kept unchanged when the sounding reaches the maximum sounding.
In one embodiment, the earth resistivity model modeling method further includes:
configuring an initial value of a first weight parameter to be 1; configuring an initial value of the second weight parameter to be 0; the initial value of the third weight parameter is configured to be 0.
When the sounding gradually increases to the first threshold, the first weight parameter is reduced, the second weight parameter is increased, and when the sounding increases to the first threshold, the first weight parameter is reduced to 0, and the second weight parameter is increased to 1;
and when the sounding is gradually increased within the range from the first threshold value to the second threshold value, decreasing the second weight parameter, and increasing the third weight parameter, and when the sounding is increased to the second threshold value, decreasing the second weight parameter to 0, and increasing the third weight parameter to 1.
A device for modeling a earth resistivity model, the device comprising:
the first earth resistivity acquisition module is used for measuring first earth resistivity data of the shallow layer by using 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 building module is used for inverting the first earth resistivity data to obtain a first inversion objective function;
the second inversion function building module is used for inverting the second earth resistivity data to obtain a second inversion objective function;
the third reflection function building module is used for inverting the third earth resistivity data to obtain a third reflection target function;
the comprehensive inversion function building module is used for obtaining a comprehensive inversion objective function according to the first inversion objective function, the second inversion objective function and the third inversion objective function;
and the optimization module is used for inverting the soil parameters of the comprehensive inversion objective function by adopting a differential evolution algorithm, and taking the comprehensive inversion objective 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 computer program is executed.
A computer readable storage medium having stored thereon a computer program which when executed by a processor implements the steps of the above-described earth resistivity model modeling method.
The earth resistivity modeling method, the earth resistivity modeling device, the computer equipment and the storage medium are used for respectively measuring the earth resistivities with different depths by adopting a quadrupole method, a controllable source audio earth electromagnetic method and an earth electromagnetic method to obtain earth resistivity data with different depths; then, mixing an inversion objective function obtained by inversion of the earth resistivity data obtained based on a four-pole method, an inversion objective function obtained by inversion of the earth resistivity data obtained based on a controllable source audio earth method and an inversion objective function obtained by inversion of the earth resistivity obtained based on a magnetotelluric method, and obtaining a comprehensive inversion objective function; and finally, carrying out soil parameter inversion on the comprehensive inversion objective function by adopting a differential evolution algorithm, and taking the comprehensive inversion objective function with the inverted soil parameters as a ground resistivity model, thereby realizing unified modeling of wide area ground resistivity covering tens of kilometers from the earth surface to the underground.
Drawings
FIG. 1 is an application environment diagram of a method for modeling a earth resistivity model in one embodiment;
FIG. 2 is a flow diagram of a method for modeling a resistivity-to-earth model in one embodiment;
FIG. 3 is a schematic illustration of a quadrupole measurement;
FIG. 4 is a schematic illustration of a measurement of a controllable source audio magnetotelluric method;
FIG. 5 is a schematic illustration of magnetotelluric measurements;
FIG. 6 is a flow chart diagram of steps of a method for modeling a resistivity model of the earth in one embodiment;
FIG. 7 is a flowchart illustrating steps of a method for modeling a resistivity model of the earth in another embodiment;
FIG. 8 is a schematic diagram of extra-high voltage DC grounding electrode address measurement;
FIG. 9 is a schematic diagram of a ground electrode ground resistance measurement;
FIG. 10 is a block diagram of a device for modeling a earth resistivity model in one embodiment;
FIG. 11 is a block diagram of the architecture of the synthetic inversion function building block in one embodiment;
FIG. 12 is a block diagram of the architecture of an optimization module in one embodiment;
fig. 13 is an internal structural view of a computer device in one embodiment.
Detailed Description
In order to make the objects, technical solutions and advantages of the present application more apparent, the present application will be further described in detail with reference to the accompanying drawings and examples. It should be understood that the specific embodiments described herein are for purposes of illustration only and are not intended to limit the present application.
The earth resistivity model modeling method provided by the application can be applied to an application environment shown in figure 1. Wherein the terminal 102 communicates with the server 104 via a network. The user may send the electric field and electromagnetic data capable of characterizing the earth resistivity, obtained based on a quadrupole method, the electric field and electromagnetic data capable of characterizing the earth resistivity, obtained based on a controllable source audio earth method, and the electric field and electromagnetic data capable of characterizing the earth resistivity, obtained based on an earth electromagnetic method, to the server 104 via the mobile terminal 102. The server 104 performs the steps of the earth resistivity modeling method based on the electric field and electromagnetic data described above to model the earth resistivity. The terminal 102 may be, but not limited to, various personal computers, notebook computers, smartphones, tablet computers, and portable wearable devices, and the server 104 may be implemented by a stand-alone server or a server cluster composed of a plurality of servers.
In one embodiment, as shown in fig. 2, there is provided a method for modeling earth resistivity, which is described by taking an example that the method is applied to the server in fig. 1, and includes the following steps:
and S20, measuring first earth resistivity data of the shallow layer by using a quadrupole method.
S40, measuring second earth resistivity data of the middle layer by adopting a controllable source audio magnetotelluric method.
S60, measuring third earth resistivity data of the deep layer by adopting an earth electromagnetic method.
S80, inverting the first earth resistivity data to obtain a first inversion objective function.
S100, inverting the second earth resistivity data to obtain a second inversion objective function.
S120, inverting the third earth resistivity data to obtain a third reflection objective function.
S140, acquiring a comprehensive inversion objective function according to the first inversion objective function, the second inversion objective function and the third inversion objective function.
And S160, carrying out soil parameter inversion on the comprehensive inversion objective function by adopting a differential evolution algorithm, and taking the comprehensive inversion objective 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 depth division.
The quadrupole method is a measurement method for deducing the underground resistivity distribution by utilizing the space potential difference formed by current backflow, has higher sensitivity to shallow earth resistivity, and is suitable for measuring shallow earth resistivity data.
Specifically, a measurement schematic diagram of the quadrupole method is shown in FIG. 3, P 1 And P 4 Is an electrode, P 2 And P 3 The measuring electrode is characterized in that the measuring electrode is a voltage electrode, a current electrode and a voltage electrode are respectively arranged on two sides of the measuring electrode by taking a measuring center point as a center, d is the distance between the electrodes, and b is the depth of the measuring electrode.
The specific measurement process for measuring the first earth resistivity data of the shallow layer by adopting the quadrupole method comprises the following steps:
the measurement locations are surveyed to determine a measurement center point.
Maintaining the measurement center pointThe arrangement is unchanged, and the wiring in the north-south direction (or any direction) is used for conducting the polar distance d i The pole distance is selected according to the principle of 1-2-3-5-7.
The same measurement is performed in an east-west direction (or perpendicular to the previous measurement direction) wiring mode and pole pitch arrangement, and the difference of apparent resistivity in the two directions is observed.
The measurement of m is completed according to the steps 1 After the pole pitch arrangement, according to a first inversion objective function F 1 Performing one-dimensional inversion of soil parameters:
wherein ρ is a (d i ) And ρ m (d i ) The i-th set of measurements, the forward and measured values of apparent resistivity at depth di, respectively. The method for calculating the forward value of the apparent resistivity by the quadrupole method is disclosed in other documents and will not be described herein.
In addition, when the pole distance exceeds 100m, the measured data is inaccurate due to the mutual inductance effect of the wires, and meanwhile, the quadrupole method is difficult to penetrate when encountering a high-resistance layer, so that the maximum measured pole distance of the quadrupole method is generally set to be about 100m, namely, the maximum measuring depth of the quadrupole method is 100m, and the earth resistivity cannot be measured when the depth exceeds 100 m.
The controllable source audio magnetotelluric method is developed to overcome the defect of low signal strength of magnetotelluric method and has the features of high signal-to-noise ratio and high interference resistance. But the controllable source audio magnetotelluric method is greatly influenced by frequency and is easy to interfere when measuring high-frequency shallow data, so that the controllable source audio magnetotelluric method is not as accurate as a quadrupole method aiming at the shallow data. In general, the sounding of the controllable source audio magnetotelluric method can reach 2.5km-3km, so in the embodiment of the invention, the four-pole method is adopted to measure the first earth resistivity data of the shallow layer, and the controllable source audio magnetotelluric method is adopted to measure the second earth resistivity data of the middle layer.
Specifically, fig. 4 illustrates a measurement scheme of a controllable source audio magnetotelluric method.
The measuring step of the second earth resistivity data of the middle layer by adopting the controllable source audio magnetotelluric method comprises the following steps:
the transmitter and the emission sources a and B are arranged and the electrodes and poles of the controllable source audio magnetotelluric method are arranged at the measurement point.
The lowest measurement frequency of the pole address is determined, thereby determining the maximum measurement depth.
Measuring apparent resistivity corresponding to a measurement point as a function of frequency f i Variable data, extracting electric field E in frequency range in receiving system x And magnetic field H y Calculate the frequency f i The corresponding apparent resistivity:
wherein Z is apparent resistivity, T is period, f i The frequency of the ith measurement.
Take out m 2 Characteristic points according to the second inversion objective function F 2 Performing one-dimensional inversion of soil parameters:
wherein ρ is a (f i ) And ρ m (f i ) The frequencies f of the i-th feature points respectively i The corresponding forward value and measured value of apparent resistivity are disclosed in other documents, and are not described 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 detection.
And moving the measuring points and repeating the steps until the measurement of all the measuring points is completed.
Magnetotelluric methods can theoretically reach a detection depth of 100 km, which depends on the lowest frequency that the measurement field device can receive. However, since the method needs to measure the data of the corresponding natural magnetic field and electric field at low frequency, the shallow and middle layer earth resistivity data measured by the method is easily submerged in the measuring site, and thus, in the embodiment of the invention, only the deep third earth resistivity data is measured by the magnetotelluric method.
Specifically, the arrangement of the magnetotelluric method is shown in fig. 5. The step of measuring deep layer earth resistivity data by using an magnetotelluric method comprises the following steps:
electrodes and magnetic poles are arranged at the measuring center point, non-polarized electrodes are respectively arranged in the four directions of southwest and northwest, the electrodes take a receiver as the center, and potential differences V in the horizontal direction x and the vertical direction y are measured x And V y And the magnetic poles are arranged in both x and y directions of the second quadrant and the fourth quadrant shown in the figure to measure the two-quadrant component H x And H y Respectively the components of the magnetic field. The electric fields in the x and y directions measured by the non-polarized electrodes are calculated as follows:
wherein E is x For the electric field in the x direction E y For the electric field in the y-direction, 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 a frequency band range, and calculating frequency f i The corresponding apparent resistivity:
wherein Z is 1 (f i ) Is the apparent resistivity in the xy direction, and Z 2 (f i ) Is the apparent resistivity in the yx direction.
Take out m 3 Characteristic points according to a third reflection objective functionF 3 Performing one-dimensional inversion of soil parameters:
wherein ρ is a (f i ) And ρ m (f i ) Respectively the ith characteristic point frequency f i The corresponding forward and measured values of apparent resistivity. The calculation method of the forward value of apparent resistivity of the magnetotelluric method has been disclosed in other documents and will not be described here.
And jointly determining the inversion result and the measurement data, and judging whether the requirements are met or not through the sounding. If not, repeating the above-mentioned magnetotelluric method measurement steps.
Adding the first inversion objective function, the second inversion objective function and the third inversion objective function to obtain a comprehensive inversion objective function, wherein the comprehensive inversion objective function minF can be expressed as:
minF＝F 1 +F 2 +F 3
in one embodiment, as shown in FIG. 6, the step of obtaining the synthetic inversion objective function from the first inversion objective function, the second inversion objective function, and the third inversion objective function includes:
s141, configuring a first weight parameter of the first inversion objective function, a second weight parameter of the second inversion objective function and a third weight parameter of the third inversion objective function according to the sounding.
The purpose of configuring the first, second and third weight parameters according to the sounding is to bias the earth resistivity model towards a more accurate method for the depth measurement for values of resistivity at different depths.
S142, carrying out weighted summation on the first inversion objective function, the second inversion objective function and the third inversion objective function to obtain the comprehensive inversion objective function.
Specifically, the configuration of the first weight parameter, the second weight parameter, and the third weight parameter may be expressed by the following formula to obtain the synthetic inversion objective function minF:
min F＝ω 1 F 1 +ω 2 F 2 +ω 3 F 3
wherein omega 1 、ω 2 And omega 3 Respectively 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 sounding. In particular, the first weight parameter is related to pole pitch, and the second weight parameter and the third weight parameter are related to frequency.
In one embodiment, the value of the first weight parameter decreases with increasing depth measurement, the value of the second weight parameter increases with increasing depth measurement and then decreases, and the value of the third weight parameter increases with increasing depth measurement and remains unchanged when the depth measurement reaches a maximum depth measurement.
In one embodiment, the earth resistivity modeling method further comprises:
s130, configuring an initial value of the first weight parameter as 1, configuring an initial value of the second weight parameter as 0 and configuring an 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 sounding includes:
when the sounding gradually increases to the first threshold, the first weight parameter is reduced, the second weight parameter is increased, and when the sounding increases to the first threshold, the first weight parameter is reduced to 0, and the second weight parameter is increased to 1;
And when the sounding is gradually increased within the range from the first threshold value to the second threshold value, decreasing the second weight parameter, and increasing the third weight parameter, and when the sounding is increased to the second threshold value, decreasing the second weight parameter to 0, and increasing the third weight parameter to 1.
Typical depth measurement ranges of the four-pole method are 0-100 m, typical depth measurement ranges of the controllable source audio magnetotelluric method are 10-2500 m, and typical depth measurement ranges of the magnetotelluric method are 2.5 km-100 km. Therefore, typical sounding ranges of the four-pole method, the controllable source audio magnetotelluric method and the magnetotelluric method are mutually overlapped to a certain extent, and the overlapped part is called a transition layer, for example, the transition layers of the shallow layer and the middle layer are 10 m-100 m. The shallow layer transition layer result in the earth resistivity model is more biased to a quadrupole method with higher measurement precision aiming at a shallow layer, the middle layer transition layer result is more biased to a controllable source audio earth electromagnetic method with higher measurement precision aiming at a middle layer, and the deep layer transition layer result is more biased to an earth electromagnetic method with higher measurement precision aiming at a deep layer, so that a more accurate wide area earth resistivity model is finally obtained.
Specifically, in shallow layer, the earth resistivity model mainly refers to the measurement result of the quadrupole method, the model output result is related to the first weight parameter and the second weight parameter, the first weight parameter initial value is set to 0, and the second weight parameter initial value is set to 1. As the sounding increases, the first weight parameter gradually decreases and the second weight parameter gradually increases. When the sounding 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, as the sounding continues to be increased, the second weight parameter is gradually reduced, the third weight parameter is increased from 0, the earth resistivity model mainly refers to the measurement result of the controllable source audio earth electromagnetic method, and the model output result is related to the second weight parameter and the third weight parameter. When the sounding is increased to a second threshold value for distinguishing the middle layer from the deep layer, the second weight parameter is reduced to 0, and the third weight parameter is increased to 1, at this time, along with the continuous increase of the sounding, the third weight parameter remains unchanged, 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 measured depth, and the first threshold value and the second threshold value are set according to the actual measured environment.
In one embodiment, as shown in fig. 7, the step S160 of inverting the soil parameters of the synthetic inversion objective function using the differential evolution algorithm includes:
s161, initializing a population by setting an inversion initial value of the comprehensive inversion objective function.
The population is formed by taking each soil parameter in the comprehensive inversion objective 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, depth, inversion precision 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 the algebra of evolution, the initial population g=0, and the individual s can be expressed as:
s162, iteratively executing the steps of performing mutation operation on individuals in the current population to generate variant individuals, performing cross operation on the current population and the variant individuals to generate experimental individuals, selecting excellent individuals between the experimental individuals and the individuals in the current population to form the next-generation population, and outputting soil parameters in the latest-generation population as soil parameters after inversion until termination conditions are met.
Specifically, the step of performing a mutation operation on the individuals in the current population to generate mutated individuals may be:
And randomly selecting three different individuals in the population, differentially weighting vectors of any two individuals, and overlapping the weighted vectors with vectors of the rest individuals to obtain the variant individuals.
Specifically, the process of generating variant individuals can be expressed as follows:
wherein r is 1 ,r 1 ,r 1 E (1, 2., N) and unlike i, k is a scale factor.
And then carrying out cross operation on the current population and the variant individuals to generate experimental individuals.
The above crossover operation can be expressed as follows:
wherein,for variant individuals, rand (j) is [0,1]Random numbers are uniformly distributed among the two; CR E [0,1 ]]Is the crossover probability; rnbr (i) is a random integer between {1,2,. }.
Fine individuals are selected between the test individuals and the individuals in the current population to form the next generation population.
And if the termination condition is judged to be met, outputting the soil parameters in the latest generation of population as the soil parameters after inversion.
If the algorithm reaches the preset inversion precision, the algorithm considers that the termination condition is met and outputs a feasible solution; if the algorithm exceeds the preset operation amount and does not reach the preset inversion precision, the method considers that the termination condition is met and outputs the current optimal solution. The preset inversion precision comprises inversion precision F1 of a first inversion objective function, inversion precision F2 of a second inversion objective function, inversion precision F3 of a third inversion objective function and inversion precision F of a comprehensive inversion objective function.
The preset operand includes the iteration times and running time of the algorithm.
In one embodiment, F 1 、F 2 And F 3 Are set to 6% and F is set to 18%. When F 1 、F 2 And F 3 And when the difference evolution algorithm is smaller than 6% and F is smaller than 18%, the difference evolution algorithm is considered to meet the termination condition, and a feasible solution is output.
If it is determined that the termination condition is not satisfied, S162 is iteratively performed.
Because the four-pole method is an electric method, the controllable source audio magnetotelluric method and the magnetotelluric method are electromagnetic methods, inversion theory of the electric method and the electromagnetic method cannot be applied across methods, and the normalization of different data to a unified size is a 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 of the comprehensive inversion objective function by adopting the differential evolution algorithm can enable the modeling process of the earth resistivity model to have higher efficiency, namely the algorithm can achieve convergence by fewer iteration times.
It should be understood that, although the steps in the flowcharts of fig. 2 or 6 or 7 are shown in order as indicated by the arrows, these steps are not necessarily performed in order as indicated by the arrows. The steps are not strictly limited to the order of execution unless explicitly recited herein, and the steps may be executed in other orders. Moreover, at least a portion of the steps in fig. 2 or 6 or 7 may include a plurality of steps or stages, which are not necessarily performed at the same time, but may be performed at different times, and the order of the execution of the steps or stages is not necessarily sequential, but may be performed in turn or alternately with at least a portion of the steps or stages in other steps or others.
In addition, the earth resistivity model modeling method provided by the embodiment of the application can be applied to extra-high voltage direct current earth electrode site earth resistivity modeling.
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 a grounding channel for rated current in monopolar earth operation or unbalanced current in bipolar asymmetric operation. When the direct current system runs on the monopole earth, strong rated current of the system flows into the earth through the grounding electrode to raise the surface potential, meanwhile, the current is dispersed in the earth to cause the surrounding soil to heat, and the stable running of the grounding electrode is seriously affected when the temperature is too high. The selection of the grounding electrode address during the DC engineering design stage is very important. Because of the strict requirements of the key performance index of the direct current grounding electrode, the grounding electrode is required to have a large enough current-dispersing area. The current land resources are tense, and the size of most conventional direct current grounding poles is large, so that the problems of difficult direct current grounding pole site selection and difficult land characterization are increasingly outstanding.
Therefore, it is a practical idea to extend the size of the grounding electrode in the vertical direction in practical engineering. When the earth resistivity of the electrode address deep layer is lower, current can be led to the deep layer of the underground, the method has certain advantages in the aspects of reducing the ground resistance and the step voltage, and meanwhile, the occupied area is reduced, so that the method is an effective ground electrode arrangement scheme. At present, the burial depth of the direct current grounding electrode can be in the range of hundreds to thousands of meters, and the electrical and thermal characteristics of the grounding electrode can be determined to a great extent by the earth resistivity of the electrode address within the range of 10 times of the burial depth.
However, at present, the conventional earth resistivity modeling method does not consider the distribution rule of earth resistivity at different depths, and cannot accurately perform uniform modeling on wide area earth resistivity covering tens of kilometers from the surface to the ground. Therefore, the embodiment of the application provides a combined measurement method and a geodetic resistivity modeling method which comprehensively use the geodetic method, the controllable source audio geodetic electromagnetic method and the geodetic electromagnetic method, so that the wide area geodetic resistivity covering tens of kilometers from the surface to the underground can be modeled more accurately in a unified way.
Specifically, the embodiments of the present invention are described in detail with respect to an extra-high voltage dc ground electrode address.
The earth electrode address requires the detection of earth resistivity at least 10km deep. According to different advantageous detection depths of the three methods, the earth resistivity test is divided into three different depth sections, and the three different methods are adopted for testing. The four-pole method is used for measuring the earth resistivity data of 0m-100m, the controllable source audio earth electromagnetic method is used for measuring the earth resistivity data of 0m-3km, and the earth electromagnetic method is used for measuring the earth resistivity data of 0-68 km. The pole address measurement is schematically shown in FIG. 8. Because of the limitation of field measurement conditions, the wiring of the four-pole method in the transverse and longitudinal directions adopts a three-longitudinal and one-transverse method, namely S1-S4 in the figure, and the maximum pole distance is 100m. The controllable source audio magnetotelluric method adopts a linear method arrangement form, namely L2-L6 in the figure. L2 and L3 are 70m apart, L3-L6 are 50m apart, each measuring line is 700m long, and the measuring point interval in each line is 12.5m. The local area of the grounding electrode address is basically covered. Depending on the field disturbance situation, it is difficult to arrange all magnetotelluric measurement points in a straight line. The magnetotelluric method therefore has four measurement points, M1-M4 in the figure.
The method in the embodiment of the application is based on the four-pole method, the controllable source audio magnetotelluric method and the magnetotelluric method, the shallow first magnetotelluric data, the middle second magnetotelluric data and the deep third magnetotelluric data are obtained, the comprehensive inversion objective function is obtained, and the differential evolution algorithm is adopted to carry out soil parameter inversion on the comprehensive inversion objective function. Setting the initial value of the number of the earth layers in the differential evolution algorithm to 8, if 'minF' is satisfied <18％，F 1 <6％，F 2 <6% and F 3 <Outputting a feasible solution if '6%' is not satisfied, if 'minF' is not satisfied<18％，F 1 <6％，F 2 <6% and F 3 <And 6% ", adding 1 to the number of soil layers to continue inversion until the preset algorithm stopping condition is met. According to the inversion thought, inverting the 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 the data of the four-pole method for the depth of 0-100m, the depth of 100-3000 m is more biased to the controllable source audio earth electromagnetic method, the depth of less than 3km is more biased to the earth electromagnetic method, 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 addresses
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 conducted on a direct current earth electrode. And carrying out the ground resistance test work of the grounding electrode by adopting a tripolar method. The measurement uses 220V alternating current power supply, the rectifier converts alternating current into direct current, and the maximum output is 7A. The voltage was measured at the end of the deep well earth electrode drainage cable, the measurement schematic diagram is shown in fig. 9, and the measurement results are shown in table 2.
Table 2 deep well ground electrode ground resistance measurements
By using the earth resistivity model obtained in the embodiment of the application, the earth resistance of the earth electrode is modeled and calculated in CDEGS simulation software, the calculation result is 0.138 omega, and compared with the 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 earth resistivity model modeling apparatus including: the system comprises a measurement module, a comprehensive inversion function establishment module and an optimization module, wherein:
a first earth resistivity acquisition module 20 for measuring shallow first earth resistivity data by a quadrupole method;
a second earth resistivity acquisition module 40 for measuring second earth resistivity data of the middle layer using a controllable source audio magnetotelluric method;
a third earth resistivity acquisition module 60 for measuring deep third earth resistivity data by an magnetotelluric method; the shallow layer, the middle layer and the deep layer are divided according to depth;
a first inversion function building module 80, configured to invert the first earth resistivity data to obtain a first inversion objective function;
a second inversion function building module 100, configured to invert the second earth resistivity data to obtain a second inversion objective function;
A third reflection function establishing module 120, configured to invert the third earth resistivity data to obtain a third reflection objective function;
the comprehensive inversion function building module 140 is configured to obtain a comprehensive inversion objective function according to the first inversion objective function, the second inversion objective function, and the third inversion objective function;
and the optimization module 160 is used for inverting the soil parameters of the comprehensive inversion objective function by adopting a differential evolution algorithm, and taking the comprehensive inversion objective function with the inverted soil parameters as a ground resistivity model.
In one embodiment, as shown in FIG. 11, the synthetic inversion function creation module includes:
the configuration unit 141 is configured to configure a first weight parameter of the first inversion objective function, a second weight parameter of the second inversion objective function, and a third weight parameter of the third inversion objective function according to the sounding.
And a superposition unit 142, configured to perform weighted summation on the first inversion objective function, the second inversion objective function, and the third inversion objective function, and obtain a comprehensive inversion objective function.
In one embodiment, the synthetic inversion function building block further comprises:
a weight initial value setting unit 130, configured to configure an initial value of the first weight parameter to be 1, configure an initial value of the second weight parameter to be 0, and configure an initial value of the third weight parameter to be 0.
In one embodiment, the configuration unit 141 includes:
the weight updating unit is used for reducing the first weight parameter and increasing the second weight parameter before the sounding gradually increases to the first threshold value, and when the sounding increases to the first threshold value, the first weight parameter is reduced to 0, and the second weight parameter is increased to 1;
and when the sounding is gradually increased within the range from the first threshold value to the second threshold value, decreasing the second weight parameter and increasing the third weight parameter, and when the sounding is increased to the second threshold value, decreasing the second weight parameter to 0 and increasing the third weight parameter 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 objective function.
The iteration execution unit 162 is configured to iteratively execute the steps of performing a mutation operation on individuals in the current population to generate mutated individuals, performing a crossover operation on the current population and the mutated individuals to generate experimental individuals, selecting good individuals between the experimental individuals and the individuals in the current population to form a next-generation population, and outputting soil parameters in the latest-generation population as soil parameters after inversion until a termination condition is satisfied.
In one embodiment, the variant generating unit includes:
and the differential weighting unit is used for randomly selecting three different individuals in the population, differentially weighting vectors of any two individuals, and then superposing the weighted vectors with the vectors of the rest individuals to obtain the variant individuals.
For specific limitations on the earth resistivity model modeling apparatus, reference may be made to the above limitations on the earth resistivity model modeling method, and no further description is given here. Each of the modules in the above-described earth resistivity model modeling apparatus may be implemented in whole or in part by software, hardware, and combinations thereof. The above modules may be embedded in hardware or may be independent of a processor in the computer device, or may be stored in software in a memory in the computer device, so that the processor may call and execute operations corresponding to the above modules.
In one embodiment, a computer device is provided, which may be a terminal, and the internal structure thereof 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 includes a non-volatile 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 the operating system and computer programs in the non-volatile storage media. The communication interface of the computer device is used for carrying out wired or wireless communication with an external terminal, and the wireless mode can be realized through WIFI, an operator network, NFC (near field communication) or other technologies. The computer program when executed by a processor implements a method for modeling a earth resistivity 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, can also be keys, a track ball or a touch pad arranged on the shell of the computer equipment, and can also be an external keyboard, a touch pad or a mouse and the like.
It will be appreciated by those skilled in the art that the structure shown in fig. 13 is merely a block diagram of a portion of the structure associated with the present application and is not limiting of the computer device to which the present application applies, and that a particular computer device may include more or fewer components than shown, or may combine some of the 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 stored therein a computer program, the processor when executing the computer program performing the steps of:
and S20, measuring first earth resistivity data of the shallow layer by using a quadrupole method.
S40, measuring second earth resistivity data of the middle layer by adopting a controllable source audio magnetotelluric method.
S60, measuring third earth resistivity data of the deep layer by adopting an earth electromagnetic method.
S80, inverting the first earth resistivity data to obtain a first inversion objective function.
S100, inverting the second earth resistivity data to obtain a second inversion objective function.
S120, inverting the third earth resistivity data to obtain a third reflection objective function.
S140, acquiring a comprehensive inversion objective function according to the first inversion objective function, the second inversion objective function and the third inversion objective function.
And S160, carrying out soil parameter inversion on the comprehensive inversion objective function by adopting a differential evolution algorithm, and taking the comprehensive inversion objective 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 depth division.
In one embodiment, the processor when executing the computer program further performs the steps of:
s130, configuring an initial value of the first weight parameter as 1, configuring an initial value of the second weight parameter as 0 and configuring an 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 objective function, a second weight parameter of the second inversion objective function and a third weight parameter of the third inversion objective function.
S142, carrying out weighted summation on the first inversion objective function, the second inversion objective function and the third inversion objective function to obtain the comprehensive inversion objective 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 objective function.
S162, iteratively executing the steps of performing mutation operation on individuals in the current population to generate variant individuals, performing cross operation on the current population and the variant individuals to generate experimental individuals, selecting excellent individuals between the experimental individuals and the individuals in the current population to form the next-generation population, and outputting soil parameters in the latest-generation population as soil parameters after inversion until termination conditions are met.
If it is determined that the termination condition is not satisfied, S162 is executed.
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 first earth resistivity data of the shallow layer by using a quadrupole method.
S40, measuring second earth resistivity data of the middle layer by adopting a controllable source audio magnetotelluric method.
S60, measuring third earth resistivity data of the deep layer by adopting an earth electromagnetic method.
S80, inverting the first earth resistivity data to obtain a first inversion objective function.
S100, inverting the second earth resistivity data to obtain a second inversion objective function.
S120, inverting the third earth resistivity data to obtain a third reflection objective function.
S140, acquiring a comprehensive inversion objective function according to the first inversion objective function, the second inversion objective function and the third inversion objective function.
And S160, carrying out soil parameter inversion on the comprehensive inversion objective function by adopting a differential evolution algorithm, and taking the comprehensive inversion objective 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 depth division.
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 objective function, a second weight parameter of the second inversion objective function and a third weight parameter of the third inversion objective function.
S142, carrying out weighted summation on the first inversion objective function, the second inversion objective function and the third inversion objective function to obtain the comprehensive inversion objective function.
In one embodiment, the computer program when executed by the processor further performs the steps of:
s130, configuring an initial value of the first weight parameter as 1, configuring an initial value of the second weight parameter as 0 and configuring an 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 objective function.
S162, iteratively executing the steps of performing mutation operation on individuals in the current population to generate variant individuals, performing cross operation on the current population and the variant individuals to generate experimental individuals, selecting excellent individuals between the experimental individuals and the individuals in the current population to form the next-generation population, and outputting soil parameters in the latest-generation population as soil parameters after inversion until termination conditions are met.
If it is determined that the termination condition is not satisfied, S162 is executed.
Those skilled in the art will appreciate that implementing all or part of the above described methods may be accomplished by way of a computer program stored on a non-transitory computer readable storage medium, which when executed, may comprise the steps of the embodiments of the methods described above. Any reference to memory, storage, database, or other medium used in embodiments provided herein may include at least one of non-volatile and volatile memory. The nonvolatile Memory may include Read-Only Memory (ROM), magnetic tape, floppy disk, flash Memory, optical Memory, or the like. Volatile memory can include random access memory (Random Access Memory, RAM) or external cache memory. By way of illustration, and not limitation, RAM can be in the form of a variety of forms, such as static random access memory (Static Random Access Memory, SRAM) or dynamic random access memory (Dynamic Random Access Memory, DRAM), and the like.
The technical features of the above embodiments may be arbitrarily combined, and all possible combinations of the technical features in the above embodiments are not described for brevity of description, however, as long as there is no contradiction between the combinations of the technical features, they should be considered as the scope of the description.
The above examples merely represent a few embodiments of the present application, which are described in more detail and are not to be construed as limiting the scope of the invention. It should be noted that it would be apparent to those skilled in the art that various modifications and improvements could be made without departing from the spirit of the present application, which would be within the scope of the present application. Accordingly, the scope of protection of the present application is to be determined by the claims appended hereto.
Claims (9)
1. A method of modeling a earth resistivity model, the method comprising:
measuring first earth resistivity data of the shallow layer by using a quadrupole method;
measuring second earth resistivity data of the middle layer by adopting a controllable source audio frequency earth electromagnetic method;
measuring third earth resistivity data of the 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 objective function;
inverting the second earth resistivity data to obtain a second inversion objective function;
inverting the third earth resistivity data to obtain a third reflection objective function;
configuring a first weight parameter of the first inversion objective function, a second weight parameter of the second inversion objective function and a third weight parameter of the third inversion objective function according to sounding; wherein the value of the first weight parameter decreases with increasing depth measurement, the value of the second weight parameter increases with increasing depth measurement and then decreases, and the value of the third weight parameter increases with increasing depth measurement, and remains unchanged when the depth measurement reaches a maximum depth measurement, so that the transition layer result between the shallow layer and the middle layer is more biased towards the controllable source audio magnetotelluric method measured for the middle layer, and the transition layer result between the middle layer and the deep layer is more biased towards the magnetotelluric method for the deep layer;
Carrying out weighted summation on the first inversion objective function, the second inversion objective function and the third inversion objective function to obtain a comprehensive inversion objective function;
and carrying out soil parameter inversion on the comprehensive inversion objective function by adopting a differential evolution algorithm, and taking the comprehensive inversion objective function with the inverted soil parameters as a ground resistivity model.
2. The method of claim 1, wherein the step of inverting the soil parameters of the synthetic inversion objective function using a differential evolution algorithm comprises:
initializing a population by setting an inversion initial value of the comprehensive inversion objective function, wherein the population is formed by taking each soil parameter in the comprehensive inversion objective function as an individual;
and iteratively executing the steps of carrying out mutation operation on individuals in the current population to generate variant individuals, carrying out cross operation on the current population and the variant individuals to generate experimental individuals, selecting excellent individuals between the experimental individuals and the individuals in the current population to form the next generation population, and outputting soil parameters in the latest generation population as the inverted soil parameters until termination conditions are met.
3. The method of claim 2, wherein the step of mutating individuals in the current population to generate mutated individuals comprises:
and randomly selecting three different individuals in the population, differentially weighting vectors of any two individuals, and overlapping the weighted vectors with vectors of the rest individuals to obtain the variant individuals.
4. The method as recited in claim 1, further comprising:
configuring an initial value of the first weight parameter to be 1; configuring an initial value of the second weight parameter to be 0; an initial value of the third weight parameter is configured to be 0.
5. The method of claim 4, wherein configuring the first weight parameter, the second weight parameter, and the third weight parameter according to sounding comprises:
decreasing the first weight parameter and increasing the second weight parameter before the sounding gradually increases to a first threshold, and decreasing the first weight parameter to 0 and increasing the second weight parameter to 1 when the sounding increases to the first threshold;
and when the sounding is gradually increased within the range from the first threshold value to the second threshold value, decreasing the second weight parameter and increasing the third weight parameter, and when the sounding is increased to the second threshold value, decreasing the second weight parameter to 0 and increasing the third weight parameter to 1.
6. A device for modeling a earth resistivity model, the device comprising:
the first earth resistivity acquisition module is used for measuring first earth resistivity data of the shallow layer by using 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 building module is used for inverting the first earth resistivity data to obtain a first inversion objective function;
the second inversion function building module is used for inverting the second earth resistivity data to obtain a second inversion objective function;
the third reflection function building module is used for inverting the third earth resistivity data to obtain a third reflection target function;
the comprehensive inversion function building module is used for obtaining comprehensive inversion target parameters according to the first inversion target function, the second inversion target function and the third inversion target function;
comprising the following steps:
the configuration unit is used for configuring a first weight parameter of the first inversion objective function, a second weight parameter of the second inversion objective function and a third weight parameter of the third inversion objective function according to sounding; wherein the value of the first weight parameter decreases with increasing depth measurement, the value of the second weight parameter increases with increasing depth measurement and then decreases, and the value of the third weight parameter increases with increasing depth measurement, and remains unchanged when the depth measurement reaches a maximum depth measurement, so that the transition layer result between the shallow layer and the middle layer is more biased towards the controllable source audio magnetotelluric method measured for the middle layer, and the transition layer result between the middle layer and the deep layer is more biased towards the magnetotelluric method for the deep layer;
The superposition unit is used for carrying out weighted summation on the first inversion objective function, the second inversion objective function and the third inversion objective function to obtain a comprehensive inversion objective function;
and the optimization module is used for inverting the soil parameters of the comprehensive inversion objective function by adopting a differential evolution algorithm, and taking the comprehensive inversion objective function with the inverted soil parameters as a ground resistivity model.
7. The apparatus of claim 6, wherein the synthetic inversion function building block further comprises:
a weight initial value setting unit configured to configure an initial value of the first weight parameter to 1, configure an initial value of the second weight parameter to 0, and configure an initial value of the third weight parameter to 0.
8. A computer device comprising a memory and a processor, the memory storing a computer program, characterized in that the processor implements the steps of the method of any one of claims 1 to 5 when the computer program is executed.
9. A computer readable storage medium, on which a computer program is stored, characterized in that the computer program, when being executed by a processor, implements the steps of the method of any of claims 1 to 5.
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 CN114047554A (en) | 2022-02-15 |
| CN114047554B true 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) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117270072B (en) * | 2023-09-19 | 2024-04-19 | 云南大学 | Gravity magnetic potential field imaging inversion method and system based on improved differential evolution algorithm |
Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| 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 |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8363509B2 (en) * | 2006-09-04 | 2013-01-29 | Daniele Colombo | Method for building velocity models for pre-stack depth migration via the simultaneous joint inversion of seismic, gravity and magnetotelluric data |
| US20090083006A1 (en) * | 2007-09-20 | 2009-03-26 | Randall Mackie | Methods and apparatus for three-dimensional inversion of electromagnetic data |
-
2021
- 2021-11-05 CN CN202111304750.4A patent/CN114047554B/en active Active
Patent Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| 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 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114047554A (en) | 2022-02-15 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Usui | 3-D inversion of magnetotelluric data using unstructured tetrahedral elements: applicability to data affected by topography | |
| CN107742015B (en) | DC induced polarization method three-dimensional numerical simulation method based on arbitrary dipole-dipole device | |
| CN112949134B (en) | Earth-well transient electromagnetic inversion method based on non-structural finite element method | |
| CN109100808A (en) | A kind of horizontal magnetic polarization field detection method of multi-thread source transient electromagnetic | |
| Ren et al. | 3-D direct current resistivity anisotropic modelling by goal-oriented adaptive finite element methods | |
| 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 | |
| CN107191181B (en) | Well periphery interface detection method based on electromagnetic scattering | |
| CN114047554B (en) | Earth resistivity model modeling method, apparatus, computer device and storage medium | |
| Zhang et al. | Adaptive finite element for 3D time-domain airborne electromagnetic modeling based on hybrid posterior error estimation | |
| CN115238550A (en) | Self-adaptive unstructured grid landslide rainfall geoelectric field numerical simulation calculation method | |
| CN113051779A (en) | Numerical simulation method of three-dimensional direct-current resistivity method | |
| Lin et al. | Three-dimensional inversion of CSAMT data in the presence of topography | |
| CN108873073A (en) | A kind of across hole resistivity tomography method of three-dimensional based on electrical method of network concurrency | |
| CN108388750B (en) | Novel radial grounding resistance calculation method | |
| CN112630842B (en) | Dipole differential device and method applied to direct-current resistivity method exploration | |
| Zhang et al. | 3D inversion of large-scale frequency-domain airborne electromagnetic data using unstructured local mesh | |
| CN112287545A (en) | Time-space division order conductivity modeling and simulating method for two-phase conductive medium | |
| CN113868919B (en) | Simplified method for electromagnetic wave logging while drilling 3D simulation | |
| CN113297526B (en) | Horizontal layered soil structure joint inversion method based on Wenner quadrupole and magnetotelluric data | |
| CN107942194B (en) | Regional area direct current fluctuation source positioning method based on sensitivity analysis | |
| CN113204905B (en) | Numerical simulation method for finite element by contact induced polarization method | |
| CN113238285B (en) | Resistivity calculation method, system and terminal for geophysical charging method exploration | |
| CN112882124B (en) | Method for delineating beneficial parts of three-dimensional uranium ore formation | |
| Lu et al. | 3D large-scale transient electromagnetic modeling using a Shift-and-Invert Krylov subspace method |
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 |