CN110334474A - A kind of quasi- true soil developing algorithm suitable for FDTD modeling - Google Patents

A kind of quasi- true soil developing algorithm suitable for FDTD modeling Download PDF

Info

Publication number
CN110334474A
CN110334474A CN201910645866.0A CN201910645866A CN110334474A CN 110334474 A CN110334474 A CN 110334474A CN 201910645866 A CN201910645866 A CN 201910645866A CN 110334474 A CN110334474 A CN 110334474A
Authority
CN
China
Prior art keywords
soil
true
quasi
fdtd
formula
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN201910645866.0A
Other languages
Chinese (zh)
Inventor
赵贤正
马春光
廖前进
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Dagang Oilfield Group Co Ltd
Original Assignee
Dagang Oilfield Group Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Dagang Oilfield Group Co Ltd filed Critical Dagang Oilfield Group Co Ltd
Priority to CN201910645866.0A priority Critical patent/CN110334474A/en
Publication of CN110334474A publication Critical patent/CN110334474A/en
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N27/00Investigating or analysing materials by the use of electric, electrochemical, or magnetic means
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V3/00Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
    • G01V3/12Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation operating with electromagnetic waves
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Theoretical Computer Science (AREA)
  • Remote Sensing (AREA)
  • Chemical & Material Sciences (AREA)
  • Electrochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • Geophysics (AREA)
  • Geology (AREA)
  • Chemical Kinetics & Catalysis (AREA)
  • Environmental & Geological Engineering (AREA)
  • Health & Medical Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • Electromagnetism (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

The invention discloses a kind of quasi- true soil developing algorithms suitable for FDTD modeling, firstly, constructing the electromagnetic parameter model of soil mixture based on the semiempirical formula for simulating true soil electromagnetic property;Secondly, generating soil random medium model according to the fractal characteristic that soil media is distributed.The beneficial effects of the present invention are: true soil modeling can be realized, the objective interaction process for reproducing electromagnetic wave and true soil, reduce the difference between FDTD Numerical Simulation Results and measured result, reliable training data is provided for the research of artificial intelligent Target identification technology, realizes the accurate analysis and explanation of electromagnetic wave lossless detection result.

Description

A kind of quasi- true soil developing algorithm suitable for FDTD modeling
Technical field
The present invention relates to a kind of quasi- true soil developing algorithm, specially a kind of quasi- true soil building suitable for FDTD modeling Algorithm belongs to the technical field of value simulation in terms of field of geophysical exploration electromagnetic wave lossless detection.
Background technique
When artificial intelligence technology is applied to electromagnetic wave lossless detection buried target body, it is desirable to provide a large amount of training with The processes such as Performance tuning, structural adjustment, parameter selection of the test data for neural network.Since Electromagnetic Simulation the data obtained has Have the characteristics that target position defines, known to model parameter, just has specific reference frame when carrying out neural net model establishing, it can To verify the forecasting accuracy of network.In addition, Electromagnetic Simulation, which is easy to batch, obtains the B total number evidence of different target parameter.Therefore, exist Before artificial intelligence technology is applied to measured data, first largely emulation data on the basis of carry out Primary Study be have it is important Meaning, electromagnetic wave lossless detection also is carried out to underground objective body, has established good preparation.
Currently, most common and most mature method is time-domain finite in the numerical simulation field of electromagnetic wave lossless detection Difference (FDTD) method.It generallys use simple geometry body and assumes the approximation method of uniform dielectric to construct simulation model, normal Satisfactory effect can be generally also obtained when advising numerical simulation.But carry out the purpose of electromagnetic wave lossless detection numerical simulation It is to provide training data for the research of artificial intelligent Target identification technology, needs to approach the emulation data of measured data, so Confidence level judge can be done to the actual performance of neural network.If simulation model is excessively ideal, the B total number generated is according in target area Domain can have apparent target signature, and clean image is presented in nontarget area thus with soil complicated during actual measurement Many disturbing factors brought by environment do not match.
Uniform dielectric model cannot simulate the random distribution characteristic of medium in soil.The soil particle of nature is rendered as not Uniform distribution characteristics, the uneven distribution of medium will lead to path and scattering signatures during Electromagnetic Wave Propagation in soil Change, image always is swept to electromagnetic wave B and causes to interfere.Convenient value only considers uniform dielectric when emulating, except in different medium circle Outside the position that medium parameter is mutated at face, it will not cause electromagnetic wave to generate any distortion and propagation path variations, thus generate Emulation data excessively idealize, it is difficult to meet the actual needs of artificial intelligence target identification training data.
Summary of the invention
The object of the invention is that providing a kind of quasi- true soil suitable for FDTD modeling to solve the above-mentioned problems Developing algorithm.
The present invention is through the following technical solutions to achieve the above objectives: a kind of quasi- true soil building suitable for FDTD modeling Algorithm, comprising the following steps:
Step 1, the semiempirical formula based on the true soil electromagnetic property of simulation, construct the electromagnetic parameter mould of soil mixture Type;
Step 2 sets the size of the net region FDTD as x, y, z, and generates the random three-dimensional data A of same size;
Step 3 carries out three dimensional fast Fourier transformation algorithm to matrix A, and result is translated, and is in zero frequency The center of three-dimensional matrice;
Step 4 sets the center point coordinate of matrix A as (ic, jc, kc), for each data point A (i, j, k), calculates it With the mould L of central point;I.e.
Step 5, according to approximate power law relation, matrix A is converted using formula;
Its formula are as follows:
Step 6 carries out inverse Fourier transform to matrix A again, and acquired results are normalized to and are quantified as 1 integer for arriving n, A series of corresponding above-mentioned soil complex dielectric permittivity numbers, and be that FDTD grid assigns corresponding medium parameter, structure in corresponding position Build quasi- true soil model.
Preferably, in the step 1, in the soil mixture model of building, soil by three kinds of sandy soil, silt and clay at To divide and constitutes, the proportion of these three ingredients is different can to indicate different soil types, and then also influences its electromagnetic property, meanwhile, Water content also plays significant impact to its electromagnetic property in soil.These comprehensive factors, the complex dielectric permittivity ε of soilmIt can by following formula :
Wherein, ε 'mWith ε 'm' be respectively complex dielectric permittivity real part and imaginary part;ρbWith ρsBe respectively soil gross density and The wherein density (unit: g/cm3) of sandy soil ingredient;β ' is relevant constant respectively with soil composition to β ";α is obtained by experience Constant, value 0.65, expression formula are as follows:
Wherein, S and C is the composition (0 < S < 1,0 < C < 1) of sandy soil and clay respectively.
ε ' in above formulafwWith ε "fwIt is the real part and imaginary part of Free water in soil, expression formula respectively are as follows:
Wherein, εw∞It is ε "fwThe limit in high frequency;τwFor the relaxation time of Free water;εw0It is opposite for the static state of water Dielectric constant, value 80.1;σeffFor effective conductivity, can be obtained by empirical equation:
σeff=0.0467+0.2204 ρb-0.4111S+0.6614C
After setting soil composition S and C, the electromagnetic parameter of the soil model is mainly by moisture content;Practical soil Water content is unevenly distributed in earth, thus can define the minimum value of water content and maximum value be respectively and;This is aqueous Range is equidistant takes n point for amount, and is substituted into the complex dielectric permittivity ε of soilmFormula can obtain a series of complex dielectric permittivities (i=1,2,3 ..., n).
Preferably, in the step 2, the model facetization of building is uniform square net, mesh spacing 1mm.
Preferably, in the step 6, the dual-mode antenna for constructing quasi- true soil model selects ideal dipole, center frequency Rate is 1.5GHz, and dual-mode antenna spacing is 20mm.
Preferably, in the step 2, the size of the constructed net region FDTD is 400 × 400 × 400.
The beneficial effects of the present invention are: should be suitable for FDTD modeling quasi- true soil developing algorithm design rationally, it can be achieved that True soil modeling, the objective interaction process for reproducing electromagnetic wave and true soil reduce FDTD Numerical Simulation Results and reality Survey the difference between result, provide reliable training data for the research of artificial intelligent Target identification technology, realization electromagnetic wave without Damage the accurate analysis and explanation of detection result.
Detailed description of the invention
Fig. 1 is three-dimensional modeling of the present invention and Numerical Simulation Results schematic diagram.
In figure: a, three-dimensional uniform dielectric soil model, b, three-dimensional random medium intend true soil model, c, three-dimensional uniform dielectric Soil model direct wave waveform, d, three-dimensional random medium intend true soil model direct wave waveform.
Specific embodiment
Following will be combined with the drawings in the embodiments of the present invention, and technical solution in the embodiment of the present invention carries out clear, complete Site preparation description, it is clear that described embodiments are only a part of the embodiments of the present invention, instead of all the embodiments.It is based on Embodiment in the present invention, it is obtained by those of ordinary skill in the art without making creative efforts every other Embodiment shall fall within the protection scope of the present invention.
Referring to Fig. 1, a kind of quasi- true soil developing algorithm suitable for FDTD modeling, comprising the following steps:
Step 1, the semiempirical formula based on the true soil electromagnetic property of simulation, construct the electromagnetic parameter mould of soil mixture Type;
In the soil mixture model of building, soil is made of three kinds of sandy soil, silt and clay ingredients, these three ingredients Proportion difference can indicate different soil types, and then also influence its electromagnetic property, meanwhile, water content is also to its electricity in soil Magnetic characteristic plays significant impact.These comprehensive factors, the complex dielectric permittivity ε of soilmIt can be obtained by following formula:
Wherein, ε 'mWith ε "mIt is the real part and imaginary part of complex dielectric permittivity respectively;ρbWith ρsBe respectively soil gross density and its The density (unit: g/cm3) of middle sandy soil ingredient;β ' is relevant constant respectively with soil composition to β ";α is resulting by experience Constant, value 0.65, expression formula are as follows:
Wherein, S and C is the composition (0 < S < 1,0 < C < 1) of sandy soil and clay respectively.
ε ' in above formulafwWith ε "fwIt is the real part and imaginary part of Free water in soil, expression formula respectively are as follows:
Wherein, εw∞It is ε "fwThe limit in high frequency;τwFor the relaxation time of Free water;εw0It is opposite for the static state of water Dielectric constant, value 80.1;σeffFor effective conductivity, can be obtained by empirical equation:
σeff=0.0467+0.2204 ρb-0.4111S+0.6614C
After setting soil composition S and C, the electromagnetic parameter of the soil model is mainly by moisture content;Practical soil Water content is unevenly distributed in earth, thus can define the minimum value of water content and maximum value be respectively and;This is aqueous Range is equidistant takes n point for amount, and is substituted into the complex dielectric permittivity ε of soilmFormula can obtain a series of complex dielectric permittivities (i=1,2,3 ..., n).
Programming mode are as follows:
import numpy as np
from scipy import fftpack
Step 2 sets the size of the net region FDTD as x, y, z, and the size of the constructed net region FDTD is 400 × 400 × 400, and the model facetization constructed is uniform square net, mesh spacing 1mm, and generate the random of same size Three-dimensional data A, the value range of numerical value are 0 to 1;
Programming mode are as follows:
Step 3 carries out three dimensional fast Fourier transformation algorithm to matrix A, and result is translated, and is in zero frequency The center of three-dimensional matrice;
Programming mode are as follows:
A=fftpack.fftn (A)
A=fftpack.fftshift (A)
Step 4 sets the center point coordinate of matrix A as (ic, jc, kc), for each data point A (i, j, k), calculates it With the mould L of central point;I.e.
Programming mode are as follows:
Step 5, according to approximate power law relation, matrix A is converted using formula;
Its formula are as follows:
Step 6 carries out inverse Fourier transform to matrix A again, needs to carry out translation to A again before inverse transformation, Acquired results are normalized to and are quantified as 1 integer for arriving n, correspond to a series of above-mentioned soil complex dielectric permittivity numbers, and in correspondence Position is that FDTD grid assigns corresponding medium parameter, constructs quasi- true soil model, and constructs the transmitting-receiving day of quasi- true soil model Line selection ideal dipole, centre frequency 1.5GHz, and dual-mode antenna spacing are 20mm.
Programming mode are as follows:
Working principle: when being suitable for the quasi- true soil developing algorithm of FDTD modeling using this, firstly, true based on simulation The semiempirical formula of soil electromagnetic property constructs the electromagnetic parameter model of soil mixture;Secondly, according to soil media distribution Fractal characteristic generates soil random medium model, it can be achieved that true soil modeling, the objective phase for reproducing electromagnetic wave and true soil Interaction process reduces the difference between FDTD Numerical Simulation Results and measured result, is artificial intelligent Target identification technology Research provides reliable training data, realizes the accurate analysis and explanation of electromagnetic wave lossless detection result.
It is obvious to a person skilled in the art that invention is not limited to the details of the above exemplary embodiments, Er Qie In the case where without departing substantially from spirit or essential attributes of the invention, the present invention can be realized in other specific forms.Therefore, no matter From the point of view of which point, the present embodiments are to be considered as illustrative and not restrictive, and the scope of the present invention is by appended power Benefit requires rather than above description limits, it is intended that all by what is fallen within the meaning and scope of the equivalent elements of the claims Variation is included within the present invention.Any reference signs in the claims should not be construed as limiting the involved claims.
In addition, it should be understood that although this specification is described in terms of embodiments, but not each embodiment is only wrapped Containing an independent technical solution, this description of the specification is merely for the sake of clarity, and those skilled in the art should It considers the specification as a whole, the technical solutions in the various embodiments may also be suitably combined, forms those skilled in the art The other embodiments being understood that.

Claims (5)

1. a kind of quasi- true soil developing algorithm suitable for FDTD modeling, it is characterised in that: the following steps are included:
Step 1, the semiempirical formula based on the true soil electromagnetic property of simulation, construct the electromagnetic parameter model of soil mixture;
Step 2 sets the size of the net region FDTD as x, y, z, and generates the random three-dimensional data A of same size;
Step 3 carries out three dimensional fast Fourier transformation algorithm to matrix A, and result is translated, and zero frequency is made to be in three-dimensional The center of matrix;
Step 4 sets the center point coordinate of matrix A as (ic, jc, kc), for each data point A (i, j, k), calculates it in The mould L of heart point;I.e.
Step 5, according to approximate power law relation, matrix A is converted using formula;
Its formula are as follows:
Step 6 carries out inverse Fourier transform to matrix A again, and acquired results are normalized to and are quantified as 1 integer for arriving n, corresponding A series of above-mentioned soil complex dielectric permittivity numbers, and be that FDTD grid assigns corresponding medium parameter in corresponding position, building is quasi- True soil model.
2. a kind of quasi- true soil developing algorithm suitable for FDTD modeling according to claim 1, it is characterised in that: described In step 1, in the soil mixture model of building, soil is made of three kinds of sandy soil, silt and clay ingredients, these three ingredients Proportion difference can indicate different soil types, and then also influence its electromagnetic property, meanwhile, water content is also to its electricity in soil Magnetic characteristic plays significant impact.These comprehensive factors, the complex dielectric permittivity ε of soilmIt can be obtained by following formula:
Wherein, ε 'mWith ε "mIt is the real part and imaginary part of complex dielectric permittivity respectively;ρbWith ρsIt is the gross density of soil respectively and wherein husky The density (unit: g/cm3) of native ingredient;β ' is relevant constant respectively with soil composition to β ";α is resulting often by experience Number, value 0.65, expression formula are as follows:
Wherein, S and C is the composition (0 < S < 1,0 < C < 1) of sandy soil and clay respectively.
ε ' in above formulafwWith ε "fwIt is the real part and imaginary part of Free water in soil, expression formula respectively are as follows:
Wherein, εw∞It is ε "fwThe limit in high frequency;τwFor the relaxation time of Free water;εw0For the static relative dielectric of water Constant, value 80.1;σeffFor effective conductivity, can be obtained by empirical equation:
σeff=0.0467+0.2204 ρb-0.4111S+0.6614C
After setting soil composition S and C, the electromagnetic parameter of the soil model is mainly by moisture content;In practical soil Water content is unevenly distributed, thus can define the minimum value of water content and maximum value be respectively and;By this water content model Enclose it is equidistant take n point, and substituted into the complex dielectric permittivity ε of soilmFormula can obtain a series of complex dielectric permittivity (i= 1,2,3,…,n)。
3. a kind of quasi- true soil developing algorithm suitable for FDTD modeling according to claim 1, it is characterised in that: described In step 2, the model facetization of building is uniform square net, mesh spacing 1mm.
4. a kind of quasi- true soil developing algorithm suitable for FDTD modeling according to claim 1, it is characterised in that: described In step 6, the dual-mode antenna for constructing quasi- true soil model selects ideal dipole, centre frequency 1.5GHz, and receives and dispatches day Line spacing is 20mm.
5. a kind of quasi- true soil developing algorithm suitable for FDTD modeling according to claim 1, it is characterised in that: described In step 2, the size of the constructed net region FDTD is 400 × 400 × 400.
CN201910645866.0A 2019-07-17 2019-07-17 A kind of quasi- true soil developing algorithm suitable for FDTD modeling Pending CN110334474A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910645866.0A CN110334474A (en) 2019-07-17 2019-07-17 A kind of quasi- true soil developing algorithm suitable for FDTD modeling

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910645866.0A CN110334474A (en) 2019-07-17 2019-07-17 A kind of quasi- true soil developing algorithm suitable for FDTD modeling

Publications (1)

Publication Number Publication Date
CN110334474A true CN110334474A (en) 2019-10-15

Family

ID=68145648

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910645866.0A Pending CN110334474A (en) 2019-07-17 2019-07-17 A kind of quasi- true soil developing algorithm suitable for FDTD modeling

Country Status (1)

Country Link
CN (1) CN110334474A (en)

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103969627A (en) * 2014-05-26 2014-08-06 苏州市数字城市工程研究中心有限公司 Ground penetrating radar large-scale three-dimensional forward modeling method based on FDTD
CN108761397A (en) * 2018-05-30 2018-11-06 中南大学 Polarization SAR model decomposition evaluation method based on electromagnetic scattering simulation

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103969627A (en) * 2014-05-26 2014-08-06 苏州市数字城市工程研究中心有限公司 Ground penetrating radar large-scale three-dimensional forward modeling method based on FDTD
CN108761397A (en) * 2018-05-30 2018-11-06 中南大学 Polarization SAR model decomposition evaluation method based on electromagnetic scattering simulation

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
NEIL R. PEPLINSKI等: ""Dielectric Properties of Soils in the 0.3-1.3-GHz Range"", 《IEEE TRANSACTIONS ON GEOSCIENCE AND REMOTE SENSING》 *
董娜: ""近地表土壤中无线通信特性研究"", 《中国优秀硕士学位论文全文数据库2013年第11期》 *
黄祥 等: ""各向异性表面生成算法及其应用"", 《华南理工大学学报(自然科学版)》 *

Similar Documents

Publication Publication Date Title
Fan et al. Multidomain pseudospectral time-domain simulations of scattering by objects buried in lossy media
CN103425816B (en) Matrix extraction method for rapidly acquiring electromagnetic scattering characteristics of metal bodies of revolution
McIntosh et al. SAR versus Sinc: What is the appropriate RF exposure metric in the range 1–10 GHz? Part II: Using complex human body models
Li et al. A three-dimensional model-based inversion algorithm using radial basis functions for microwave data
CN108959772A (en) Large-scale finite period array structure characteristic pattern analysis method
CN110852025A (en) Three-dimensional electromagnetic slow diffusion numerical simulation method based on hyperconvergence interpolation approximation
CN107271977A (en) High-precision SAR echo simulation methods based on mobile driving source FDTD algorithms
CN108761397A (en) Polarization SAR model decomposition evaluation method based on electromagnetic scattering simulation
CN110244351A (en) A kind of Uniform Construction inversion method of different constraint Geophysical Inverse Problems
CN111144013B (en) High-precision dielectric object scattering simulation method
Šušnjara et al. Stochastic-deterministic thermal dosimetry below 6 GHz for 5G mobile communication systems
Qureshi et al. Assessment of exposure to radio frequency electromagnetic fields from smart utility meters in GB; part II) numerical assessment of induced SAR within the human body
Gu et al. Hybrid method for full-wave simulations of forests at L-band
Roqueta et al. Microwave time-domain reflection imaging of steel fiber distribution on reinforced concrete
CN110334474A (en) A kind of quasi- true soil developing algorithm suitable for FDTD modeling
CN104062656B (en) Microwave imaging method based on vector surface element model
Stuchly et al. Exposure of man in the near-field of a resonant dipole: comparison between theory and measurements
Hart Mathematical modeling of electromagnetic interactions with biological systems
Psenakova Numerical modeling of electromagnetic field effects on the human body
Yang et al. Electromagnetic scattering of rough ground surface covered by multilayers vegetation
Han et al. Mid-and low-latitude ionospheric D region remote sensing by radio atmospherics—Part I: Forward modeling and field measurement validations
Hu et al. BFGS quasi-Newton method for solving electromagnetic inverse problems
CN107817492A (en) The imaging method and device of wide angle synthetic aperture radar
Oikawa et al. Analysis of lightning electromagnetic field on large‐scale terrain model using three‐dimensional MW‐FDTD parallel computation
Chen et al. SAR affected by shapes and electrical properties of the human head exposed to a cellular phone

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
RJ01 Rejection of invention patent application after publication

Application publication date: 20191015

RJ01 Rejection of invention patent application after publication