CN106503302B - A kind of method and device for establishing blood electromagnetic simulation model - Google Patents
A kind of method and device for establishing blood electromagnetic simulation model Download PDFInfo
- Publication number
- CN106503302B CN106503302B CN201610862038.9A CN201610862038A CN106503302B CN 106503302 B CN106503302 B CN 106503302B CN 201610862038 A CN201610862038 A CN 201610862038A CN 106503302 B CN106503302 B CN 106503302B
- Authority
- CN
- China
- Prior art keywords
- blood
- simulation model
- electromagnetic simulation
- layer
- lipid
- 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
Classifications
-
- 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
Abstract
The present invention relates to lipids detection technical field, in particular to a kind of method and device for establishing blood electromagnetic simulation model.The method for establishing blood electromagnetic simulation model includes: step a: establishing blood electromagnetic simulation model；Step b: according to the volume of blood in the concentration of blood lipid and blood electromagnetic simulation model, the number of blood lipid particle in blood electromagnetic simulation model is calculated；Blood lipid particle: being randomly dispersed in blood electromagnetic simulation model by step c according to calculated result, and by the number of control blood lipid particle, establishes the variable blood electromagnetic simulation model of serum lipid concentrations.The present invention can be used for analyzing the interaction mechanism of blood lipid and electromagnetic wave, obtains the electromagnetic response characteristic under different serum lipid concentrations, provides important theory support for the development of noninvasive lipids detection technology.Meanwhile the present invention also has many advantages, such as that use cost is low, simulation accuracy is high.
Description
Technical field
The present invention relates to lipids detection technical field, in particular to a kind of method for establishing blood electromagnetic simulation model and dress
It sets.
Background technique
Blood lipid is neutral fat (triglycerides and cholesterol) and lipoid (phosphatide, glycolipid, sterol, steroids) in blood
General name, be widely present in human body.They are the basic metabolism necessary materials of life cells.Dyslipidemia person is often with more
Kind of cardiovascular risk factors result even in some diseases for seriously endangering human healths, as atherosclerosis, coronary heart disease,
Pancreatitis etc..In order to improve health of people living standard, detection is carried out to blood lipid and is very important.
Noninvasive lipids detection method is a kind of reflection using electromagnetic wave, transmission the method that obtains the serum lipid concentrations of patient,
There is no need to acquire the blood of patient, with noninvasive, easy, quick etc. advantage.Establishing lipids detection model is that research is noninvasive
The premise of lipids detection technology.Currently, in terms of blood lipid Modeling Research, researcher mainly use experiment statistics method to blood lipid into
Row modeling, i.e., by Animal Experimental Study under different feeding mode its Blood Lipid situation, thus establish hyperlipemia etc. doctor
Learn model.For example, Chinese patent CN103299950A discloses a kind of machin hyperlipemia and Atherosclerosis Model
Method for building up, this method complete high blood according to the changes of biochemical indexes of machin mainly by feeding high lipid food to machin
The foundation of rouge disease model.Another Chinese patent CN102907357A discloses a kind of building side of zebra fish hyperlipemia model
Method, this method make zebra fish histochemical stain or fluorescent staining, and obtain related figure by feeding zebra fish using yolk powder
Picture is analyzed and is counted to image/microwell plate, and the hyperlipemia model of zebra fish is finally established.
However, the medical models such as hyperlipemia that the above-mentioned method using experiment statistics is established are only applicable to research blood lipid
Changing rule, it is impossible to be used in the interaction mechanism of analysis blood lipid and electromagnetic wave, therefore can not be noninvasive lipids detection technology
Theory support is provided.
Summary of the invention
The present invention provides a kind of method and devices for establishing blood electromagnetic simulation model, it is intended at least to a certain extent
Solve one of above-mentioned technical problem in the prior art.
To solve the above-mentioned problems, the present invention provides the following technical scheme that
A method of establishing blood electromagnetic simulation model, comprising:
Step a: blood electromagnetic simulation model is established；
Step b: according to the volume of blood in the concentration of blood lipid and blood electromagnetic simulation model, blood Electromagnetic Simulation mould is calculated
The number of blood lipid particle in type；
Blood lipid particle: being randomly dispersed in blood electromagnetic simulation model by step c according to calculated result, and passes through control blood
The number of fat granule establishes the variable blood electromagnetic simulation model of serum lipid concentrations.
The technical solution that the embodiment of the present invention is taken further include: before the step a further include: determine blood Electromagnetic Simulation mould
The overall structure of type；The overall structure of the blood electromagnetic simulation model is cylindrical body, and the blood electromagnetic simulation model includes
Vascular lamina, blood layer and blood lipid particle layer, the distributing position of the vascular lamina, blood layer and blood lipid particle layer are respectively as follows: blood vessel
Layer is located at outermost layer, and the length of vascular lamina is h, and the internal diameter of vascular lamina is d, outer diameter D, vascular lamina with a thickness of r12=(D-d)/
2；Blood layer is located in vascular lamina, and the radius of blood layer is r=d/2；Blood lipid particle layer is distributed in blood layer.
The technical solution that the embodiment of the present invention is taken further include: in the step b, in the blood electromagnetic simulation model
Blood lipid population purpose calculation formula are as follows:
In above-mentioned formula, n is the number of blood lipid particle, and N is the number of neutral fat molecule and lipoid molecule, and ρ is blood lipid
Concentration, V are blood volume, M=ρ × V/1000, NA=6.02 × 1023, V=1000 × π (d/2)2×h。
The technical solution that the embodiment of the present invention is taken further include: described to establish what serum lipid concentrations can be changed in the step c
Blood electromagnetic simulation model specifically includes:
Step c1: n random numbers equally distributed between (- 1,1) are generated
Step c2: according to the central-limit theorem of probability theory, generating the random value of a Normal Distribution N (μ, σ),
Generate the formula of random value are as follows:
In above-mentioned formula,For uniform random number, μ is the desired value of normal distribution random number, and σ is normal distribution
The mean square deviation of random number, n are the numbers of uniform random number needed for generating normal distribution random number；
Step c3: blood lipid particle each in blood electromagnetic simulation model is determined according to the formula of above-mentioned generation random value
The variable blood electromagnetic simulation model of serum lipid concentrations is established in position.
The technical solution that the embodiment of the present invention is taken further include: after the step c further include: utilize triple debye-
Drude model is fitted the electromagnetic parameter of blood layer and vascular lamina in blood electromagnetic simulation model, and fitting result is led
Enter in blood electromagnetic simulation model, numerical value calculating is carried out to blood electromagnetic simulation model.
A kind of another technical solution that the embodiment of the present invention is taken are as follows: device for establishing blood electromagnetic simulation model, comprising:
First model building module: for establishing blood electromagnetic simulation model；
Blood lipid number of particles computing module: the body for blood in the concentration and blood electromagnetic simulation model according to blood lipid
Product calculates the number of blood lipid particle in blood electromagnetic simulation model；
Second model building module: for blood lipid particle to be randomly dispersed in blood electromagnetic simulation model according to calculated result
In, and by the number of control blood lipid particle, establish the variable blood electromagnetic simulation model of serum lipid concentrations.
The technical solution that the embodiment of the present invention is taken further includes structural calculation module, and the structural calculation module is for determining
The overall structure of blood electromagnetic simulation model；The overall structure of the blood electromagnetic simulation model is cylindrical body, the blood electricity
Magnetic simulation model includes vascular lamina, blood layer and blood lipid particle layer, the distribution position of the vascular lamina, blood layer and blood lipid particle layer
It sets and is respectively as follows: vascular lamina and is located at outermost layer, the length of vascular lamina is h, and the internal diameter of vascular lamina is d, outer diameter D, the thickness of vascular lamina
Degree is r12=(D-d)/2；Blood layer is located in vascular lamina, and the radius of blood layer is r=d/2；Blood lipid particle layer is distributed in blood
In layer.
The technical solution that the embodiment of the present invention is taken further include: the blood lipid number of particles computing module calculates blood electromagnetism
Blood lipid population purpose formula in simulation model are as follows:
In above-mentioned formula, n is the number of blood lipid particle, and N is the number of neutral fat molecule and lipoid molecule, and ρ is blood lipid
Concentration, V are blood volume, M=ρ × V/1000, NA=6.02 × 1023, V=1000 × π (d/2)2×h。
The technical solution that the embodiment of the present invention is taken further include: it is variable that second model building module establishes serum lipid concentrations
The mode of blood electromagnetic simulation model include: to generate n random numbers equally distributed between (- 1,1)；
According to the central-limit theorem of probability theory, the random value of a Normal Distribution N (μ, σ) is generated, the public affairs of random value are generated
Formula are as follows:
The position that blood lipid particle each in blood electromagnetic simulation model is determined according to the formula of above-mentioned generation random value, builds
The variable blood electromagnetic simulation model of vertical serum lipid concentrations；In above-mentioned formula,For uniform random number, μ is normal distribution
The desired value of random number, σ are the mean square deviations of normal distribution random number, and n is uniformly distributed needed for generating normal distribution random number
The number of random number.
The technical solution that the embodiment of the present invention is taken further includes parameter calculating module, and the parameter calculating module is for utilizing
Triple debye-drude models are fitted the electromagnetic parameter of blood layer and vascular lamina in blood electromagnetic simulation model, and will
Fitting result imports in blood electromagnetic simulation model, carries out numerical value calculating to blood electromagnetic simulation model.
Compared with the existing technology, the beneficial effect that the embodiment of the present invention generates is: the embodiment of the present invention establishes blood
The method and device of electromagnetic simulation model is in modeling process, according to blood body in the concentration of blood lipid and blood electromagnetic simulation model
Long-pending size calculates the number of blood lipid particle in blood electromagnetic simulation model, by controlling the number of blood lipid particle, to establish
The variable blood electromagnetic simulation model of serum lipid concentrations；And using triple debye-drude models in blood electromagnetic simulation model
The electromagnetic parameter of different tissues layer is fitted, to set its electromagnetic parameter.The blood Electromagnetic Simulation established through the invention
Model can be used for analyzing the interaction mechanism of blood lipid and electromagnetic wave, and the electromagnetic response obtained under different serum lipid concentrations is special
Property, important theory support is provided for the development of noninvasive lipids detection technology.Meanwhile the present invention also has use cost low, imitative
The advantages that very precision is high.
Detailed description of the invention
Fig. 1 is the flow chart of the method for establishing blood electromagnetic simulation model of the embodiment of the present invention；
Fig. 2 is the overall structure diagram of the blood electromagnetic simulation model of the embodiment of the present invention；
Fig. 3 is the blood electromagnetic simulation model for the different serum lipid concentrations established using blood lipid particle；
Fig. 4 is the structural schematic diagram of the device for establishing blood electromagnetic simulation model of the embodiment of the present invention.
Specific embodiment
In order to make the objectives, technical solutions, and advantages of the present invention clearer, with reference to the accompanying drawings and embodiments, right
The present invention is further elaborated.It should be appreciated that described herein, specific examples are only used to explain the present invention, not
For limiting the present invention.
Referring to Fig. 1, being the flow chart of the method for establishing blood electromagnetic simulation model of the embodiment of the present invention.The present invention is real
Apply the method for establishing blood electromagnetic simulation model of example the following steps are included:
Step 100: determining the overall structure of blood electromagnetic simulation model；
In step 100, the integrally-built method of determination of blood electromagnetic simulation model includes:
Step 110: according to human anatomy principle, studying blood vessel and blood in the distribution situation of human body, to blood vessel and blood
Liquid carries out abstract processing, determines the overall structure of blood electromagnetic simulation model；
In step 110, blood electromagnetic simulation model includes three vascular lamina, blood layer and blood lipid particle layer parts；This
In inventive embodiments, the overall structure of the blood electromagnetic simulation model is cylindrical body.Specifically as shown in Fig. 2, implementing for the present invention
The overall structure diagram of the blood electromagnetic simulation model of example.In other embodiments of the present invention, blood electromagnetic simulation model
Overall structure can also be other shapes.
Step 120: determining the distribution position of vascular lamina, blood layer and blood lipid particle layer in blood electromagnetic simulation model respectively
It sets；
In the step 120, the distributing position of vascular lamina, blood layer and blood lipid particle layer in blood electromagnetic simulation model point
Not are as follows: vascular lamina is located at outermost layer, and the length of vascular lamina is h, and the internal diameter of vascular lamina is d, outer diameter D, therefore vascular lamina
With a thickness of r12=(D-d)/2.Blood layer is located in vascular lamina, therefore the radius of blood layer is r=d/2.Blood lipid particle layer then divides
Cloth is in blood layer.
Step 200: blood vessel Electromagnetic Simulation mould is established according to the distributing position of vascular lamina, blood layer and blood lipid particle layer respectively
Type, blood electromagnetic simulation model and blood lipid particle electromagnetic simulation model；
In step 200, blood vessel electromagnetic simulation model, blood electromagnetic simulation model and blood lipid particle electromagnetic simulation model
Establish mode respectively include:
Step 210: establish blood vessel electromagnetic simulation model: using origin as the center of circle, establish length be h, radius be respectively d/2 and
The two cylindrical bodies are carried out additive operation by the cylindrical body of D/2, so that an inner hollow is obtained, with a thickness of r12=(D-d)/2
Blood vessel electromagnetic simulation model；
Step 220: establish blood electromagnetic simulation model: using origin as the center of circle, establishing length is h, and radius is the cylinder of d/2
Body, the cylindrical body are blood electromagnetic simulation model；
Step 230: establishing blood lipid particle electromagnetic simulation model；
In step 230, since mainly by neutral fat molecule, (triglycerides and gallbladder are solid for the blood lipid in blood of human body
Alcohol) and lipoid molecule (phosphatide, glycolipid, sterol, steroids) composition, and the mean radius one of neutral fat molecule and lipoid molecule
As be only 50nm.If establishing blood electromagnetic simulation model according to this size, and Electromagnetic Simulation is carried out with this size, then needed net
Lattice division is very fine, causes calculation amount very huge.In order to reduce calculation amount, accelerate simulation velocity, while not influencing to imitate again
True precision, the embodiment of the present invention in modeling process, by blood neutral fat molecule and lipoid molecule be equivalent to iipomicron
Son, Equivalent calculation method are as follows:
It being calculated to simplify, it is assumed that the shape of neutral fat molecule and lipoid molecule in blood of human body is sphere, and
Size is the same, radius rr, in electromagnetic simulation model, the shape for the blood lipid particle electromagnetic simulation model that will be established
It is set as sphere, radius Rr, then have
In formula (1), N is the number of neutral fat molecule and lipoid molecule.Due to neutral fat molecule and lipoid point
The radius of son is generally 50nm, i.e. rr=50nm is set as in blood electromagnetic simulation model when by the radius of blood lipid particle
0.5mm, i.e. RrN=1 × 10 can be calculated according to formula (1) in=0.5mm12, i.e., in modeling process, by every 1 × 1012It is a
Blood lipid molecule is equivalent to the blood lipid particle that a diameter is about 0.5mm, to greatly reduce simulation calculation amount.
Step 300: according to the size of blood volume in the concentration of blood lipid and blood electromagnetic simulation model, calculating blood electromagnetism
The number of blood lipid particle in simulation model；
In step 300, blood lipid population purpose calculation method specifically includes:
Step 310: calculating the volume of blood in blood electromagnetic simulation model: since blood electromagnetic simulation model is cylinder
Body, volume may be expressed as: V=1000 × π (d/2)2× h, wherein V is the volume of cylindrical body, and unit L, d/2 are cylindrical body
Radius, h is the length of cylindrical body, and the unit of d/2 and h are m.
Step 320: assuming that the mean concentration of blood lipid is ρ, unit mmol/L, when the volume of blood is V=1000 × π
(d/2)2The amount of × h, the substance of blood lipid molecule (i.e. neutral fat molecule and lipoid molecule) are then M=ρ × V/1000, wherein M
Unit be mol.
Step 330: according to Avgadro constant it is found that the substance of every 1mol contains 6.02 × 1023A molecule, i.e. NA=
6.02×1023, therefore, when serum lipid concentrations are ρ, and blood volume is V, the number of blood lipid particle is in blood electromagnetic simulation model
According to formula (2), the iipomicron in blood electromagnetic simulation model corresponding when different serum lipid concentrations can be calculated
The number of son.It therefore, can be by the number n of change blood lipid particle, to establish serum lipid concentrations in blood electromagnetic simulation model
Variable blood electromagnetic simulation model.
Step 400: blood lipid particle being pressed by normal state random distribution principle according to calculated result, it is imitative to be randomly dispersed in blood electromagnetism
In true mode, and by the number of control blood lipid particle, the variable blood electromagnetic simulation model of serum lipid concentrations is established；
In step 400, the concrete mode for establishing the variable blood electromagnetic simulation model of serum lipid concentrations includes:
Step 410: generating n random numbers equally distributed between (- 1,1)；
Step 420: according to the central-limit theorem in probability theory, generating a Normal Distribution N by formula (3)
The random value of (μ, σ):
In formula (3),For uniform random number, μ is the desired value of normal distribution random number, and σ is normal distribution
The mean square deviation of random number, n are the numbers of uniform random number needed for generating normal distribution random number.
Step 430: determining the position of blood lipid particle each in blood electromagnetic simulation model according to formula (3), finally build
The variable blood electromagnetic simulation model of vertical serum lipid concentrations, as shown in figure 3, being the different serum lipid concentrations established using blood lipid particle
Blood electromagnetic simulation model.
Step 500: using triple debye-drude models to different tissues layer (blood, blood in blood electromagnetic simulation model
Pipe) electromagnetic parameter be fitted, and by fitting result import blood electromagnetic simulation model in, to blood electromagnetic simulation model into
Line number value calculates；
In step 500, after establishing blood electromagnetic simulation model, also need in the model blood layer and vascular lamina setting
Electromagnetic parameter, electromagnetic parameter include dielectric constant and conductivity.Since blood and blood vessel belong to dispersive medium, dielectric constant
It can change with frequency with conductivity.During Electromagnetic Simulation, traditional Cole-Cole model may be only available for single
The emulation of frequency, it is helpless in broadband.Therefore, the embodiment of the present invention is using debye-drude model come to blood layer
It is fitted with the dielectric constant and conductivity of vascular lamina, relevant parameter is imported into corresponding organized layer again after fitting
In.Electromagnetic parameter fit approach specifically includes:
Step 510: obtaining the real part of the complex dielectric permittivity of blood layer and vascular lamina respectively by quadruple Cole-Cole model
Value and imaginary values；
In step 510, the blood quadruple Cole-Cole model that the embodiment of the present invention utilizes REMCOM company to provide, such as
Shown in formula (4), calculating find out in 10Hz-20GHz band limits the value of real part of the complex dielectric permittivity of blood layer and vascular lamina with
Imaginary values, it is dielectric constant that wherein value of real part is corresponding, and corresponding imaginary values are conductivity, and value of real part and imaginary values are distinguished
Origin software is imported as input value.
Step 520: triple debye-drude models are established, using triple debye-drude models to acquired blood
The value of real part of the complex dielectric permittivity of layer and vascular lamina is fitted respectively with imaginary values；Fitting formula is as follows:
Step 530: determining the parameter that blood layer and vascular lamina need to be fitted in triple debye-drude models；To triple
Debye-drude model carries out rational expression decomposition, calculates separately real part (ε '), the corresponding multinomial of imaginary part (ε "), calculation formula
It is as follows:
Determine that needing the parameter being fitted includes: ε∞=a, Δ ε1=b, τ1=c, Δ ε2=d, τ2=e, Δ ε3=g, τ3=h,
σi=k；
Step 540: corresponding nonlinear curve y1, y2 are established in Origin software,
Y1=ε ' (8)
Y2=ε " (9)
Wherein,
Y1=a+b/ (1+ (2 × 3.1415926 × f × c) ^2)+d/ (1+ (2*3.1415926*f*e) ^2)+
g/(1+(2×3.1415926×f×h)^2) (10)
Y2=(b × 2 × 3.1415926 × f × c)/(1+ (2 × 3.1415926 × f × c) ^2)+(d × 2 ×
3.1415926×f×e)/(1+(2×3.1415926×f×e)^2)+
(g×2×3.1415926×f×h)/(1+(2×3.1415926×f×h)^2)+k/(2×3.1415926×f
×8.85418×pow(10,-12)) (11)
Step 550: iteration being assigned to a, b, c, d, e, g, h, k parameter respectively using the initial value in Cole-Cole model and is intended
The initial value of conjunction carries out parameter fitting using Origin software；
Step 560: the parametric fitting results of blood layer and vascular lamina are substituted into the triple pole models of Debye-Drude respectively
In, then error analysis is carried out using parametric fitting results of the Origin software to blood layer and vascular lamina；
Step 570: the parametric fitting results of blood layer and vascular lamina being imported in blood electromagnetic simulation model, to blood electricity
Magnetic simulation model carries out numerical value calculating.
Referring to Fig. 4, being the structural schematic diagram of the device for establishing blood electromagnetic simulation model of the embodiment of the present invention.This hair
The device for establishing blood electromagnetic simulation model of bright embodiment includes structural calculation module, the first model building module, iipomicron
Subnumber mesh computing module, the second model building module and parameter calculating module.
Structural calculation module: for determining the overall structure of blood electromagnetic simulation model according to human anatomy principle；Its
In, the integrally-built method of determination of blood electromagnetic simulation model specifically: firstly, studying blood vessel according to human anatomy principle
With blood in the distribution situation of human body, abstract processing is carried out to blood vessel and blood, determines the entirety of blood electromagnetic simulation model
Structure, wherein blood electromagnetic simulation model includes three blood vessel, blood and blood lipid particle parts.Then, respectively determine blood vessel,
The distributing position of blood and blood lipid particle in blood electromagnetic simulation model；Wherein, blood vessel, blood and blood lipid particle are in blood electricity
Distributing position in magnetic simulation model is respectively as follows: blood vessel and is located at outermost layer, and the length of blood vessel is h, and the internal diameter of blood vessel is d, outside
Diameter is D, thus blood vessel with a thickness of r12=(D-d)/2.Blood is located in blood vessel, therefore the radius of blood is r=d/2.Blood lipid
Particle is then distributed in blood.
First model building module: for establishing blood vessel electricity respectively according to the distributing position of blood vessel, blood and blood lipid particle
Magnetic simulation model, blood electromagnetic simulation model and blood lipid particle electromagnetic simulation model；Specifically, the first model building module includes
Vascular pattern establishes unit, Blood Model establishes unit and blood lipid particle model establishes unit；
Vascular pattern establishes unit: for establishing blood vessel electromagnetic simulation model: using origin as the center of circle, establishing length is h, partly
Diameter is respectively the cylindrical body of d/2 and D/2, the two cylindrical bodies are carried out additive operation, to obtain an inner hollow, thickness
For r12The blood vessel electromagnetic simulation model of=(D-d)/2；
Blood Model establishes unit: for establishing blood electromagnetic simulation model: using origin as the center of circle, establishing length is h, partly
Diameter is the cylindrical body of d/2, which is blood electromagnetic simulation model；
Blood lipid particle model establishes unit: for establishing blood lipid particle electromagnetic simulation model: due to the blood in blood of human body
Rouge is mainly by neutral fat molecule and lipoid molecular composition, and the mean radius of neutral fat molecule and lipoid molecule is generally only
For 50nm.If establishing blood electromagnetic simulation model according to this size, and Electromagnetic Simulation is carried out with this size, then needs to draw grid
It is point very fine, cause calculation amount very huge.In order to reduce calculation amount, accelerate simulation velocity, while not influencing emulation essence again
Degree, the embodiment of the present invention in modeling process, by blood neutral fat molecule and lipoid molecule be equivalent to blood lipid particle, etc.
It is as follows to imitate calculation method:
It being calculated to simplify, it is assumed that the shape of neutral fat molecule and lipoid molecule in blood of human body is sphere, and
Size is the same, radius rr, in electromagnetic simulation model, the shape for the blood lipid particle electromagnetic simulation model that will be established
It is set as sphere, radius Rr, then have
In formula (1), N is the number of neutral fat molecule and lipoid molecule.Due to neutral fat molecule and lipoid point
The radius of son is generally 50nm, i.e. rr=50nm is set as in blood electromagnetic simulation model when by the radius of blood lipid particle
0.5mm, i.e. RrN=1 × 10 can be calculated according to formula (1) in=0.5mm12, i.e., in modeling process, by every 1 × 1012It is a
Blood lipid molecule is equivalent to the blood lipid particle that a diameter is about 0.5mm, to greatly reduce simulation calculation amount.
Blood lipid number of particles computing module: for blood volume in the concentration and blood electromagnetic simulation model according to blood lipid
Size calculates the number of blood lipid particle in blood electromagnetic simulation model；Blood lipid number of particles computing module specifically includes blood body
Product computing unit, serum lipid concentrations computing unit and blood lipid number of particles computing unit；
Blood volume computing unit: for calculating the volume of blood in blood electromagnetic simulation model: since blood electromagnetism is imitative
True mode is cylindrical body, and volume may be expressed as: V=1000 × π (d/2)2× h, wherein V is the volume of cylindrical body, and unit is
L, d/2 are the radius of cylindrical body, and h is the length of cylindrical body, and the unit of d/2 and h are m.
Serum lipid concentrations computing unit: for calculating serum lipid concentrations, it is assumed that the mean concentration of blood lipid is ρ, unit mmol/
L, when the volume of blood is V=1000 × π (d/2)2× h, the substance of blood lipid molecule (i.e. neutral fat molecule and lipoid molecule)
Amount be then M=ρ × V/1000, wherein the unit of M be mol.
Blood lipid number of particles computing unit: for being calculated in blood electromagnetic simulation model according to blood volume and serum lipid concentrations
The number of blood lipid particle；According to Avgadro constant it is found that the substance of every 1mol contains 6.02 × 1023A molecule, i.e. NA=
6.02×1023, therefore, when serum lipid concentrations are ρ, and blood volume is V, the number of blood lipid particle is in blood electromagnetic simulation model
According to formula (2), the iipomicron in blood electromagnetic simulation model corresponding when different serum lipid concentrations can be calculated
The number of son.It therefore, can be by the number n of change blood lipid particle, to establish serum lipid concentrations in blood electromagnetic simulation model
Variable blood electromagnetic simulation model.
Second model building module: random to divide for blood lipid particle to be pressed normal state random distribution principle according to calculated result
It is imitative to establish the variable blood electromagnetism of serum lipid concentrations in blood electromagnetic simulation model, and through the number of control blood lipid particle for cloth
True mode；Wherein, the second model building module establishes the concrete mode of the variable blood electromagnetic simulation model of serum lipid concentrations and includes:
1: generating n random numbers equally distributed between (- 1,1)；
2: according to the central-limit theorem in probability theory, generating a Normal Distribution N (μ, σ) by formula (3)
Random value:
In formula (3),For uniform random number, μ is the desired value of normal distribution random number, and σ is normal distribution
The mean square deviation of random number, n are the numbers of uniform random number needed for generating normal distribution random number.
3: determining the position of blood lipid particle each in blood electromagnetic simulation model according to formula (3), finally establish blood lipid
The blood electromagnetic simulation model of varying concentrations.
Parameter calculating module: for utilizing triple debye-drude models to different tissues in blood electromagnetic simulation model
The electromagnetic parameter of layer (blood, blood vessel) is fitted, and fitting result is imported in blood electromagnetic simulation model, to blood electromagnetism
Simulation model carries out numerical value calculating；Wherein, electromagnetic parameter includes dielectric constant and conductivity, since blood and blood vessel belong to color
Dispersion media, dielectric constant and conductivity can change with frequency.During Electromagnetic Simulation, traditional Cole-Cole
Model may be only available for the emulation of single-frequency, helpless in broadband.Therefore, the embodiment of the present invention uses debye-
Drude model come to blood layer and vascular lamina dielectric constant and conductivity be fitted, again by relevant parameter after fitting
It imported into corresponding organized layer.Electromagnetic parameter fit approach specifically includes:
1: obtaining the value of real part and void of the complex dielectric permittivity of blood layer and vascular lamina respectively by quadruple Cole-Cole model
Portion's value；The blood quadruple Cole-Cole model that the embodiment of the present invention utilizes REMCOM company to provide is calculated as shown in formula (4)
The value of real part and imaginary values of the complex dielectric permittivity of blood layer and vascular lamina in 10Hz-20GHz band limits are found out, wherein value of real part
Corresponding is dielectric constant, and corresponding imaginary values are conductivity, and value of real part and imaginary values are directed respectively into Origin software and made
For input value.
2: triple debye-drude models are established, using triple debye-drude models to acquired blood layer and blood
The value of real part of the complex dielectric permittivity of tube layer is fitted respectively with imaginary values；Fitting formula is as follows:
3: determining the parameter that blood layer and vascular lamina need to be fitted in triple debye-drude models；To triple debye-
Drude model carries out rational expression decomposition, calculates separately real part (ε '), the corresponding multinomial of imaginary part (ε "), calculation formula is as follows:
Determine that needing the parameter being fitted includes: ε∞=a, Δ ε1=b, τ1=c, Δ ε2=d, τ2=e, Δ ε3=g, τ3=h,
σi=k；
4: corresponding nonlinear curve y1, y2 are established in Origin software,
Y1=ε ' (8)
Y2=ε " (9)
Wherein,
Y1=a+b/ (1+ (2 × 3.1415926 × f × c) ^2)+d/ (1+ (2*3.1415926*f*e) ^2)+
g/(1+(2×3.1415926×f×h)^2) (10)
Y2=(b × 2 × 3.1415926 × f × c)/(1+ (2 × 3.1415926 × f × c) ^2)+(d × 2 ×
3.1415926×f×e)/(1+(2×3.1415926×f×e)^2)+
(g×2×3.1415926×f×h)/(1+(2×3.1415926×f×h)^2)+k/(2×3.1415926×f
×8.85418×pow(10,-12)) (11)
5: the first of iterative fitting is assigned to a, b, c, d, e, g, h, k parameter respectively using the initial value in Cole-Cole model
Initial value carries out parameter fitting using Origin software；
6: the parametric fitting results of blood layer and vascular lamina are substituted into respectively in the triple pole models of Debye-Drude, so
Error analysis is carried out using parametric fitting results of the Origin software to blood layer and vascular lamina afterwards；
7: the parametric fitting results of blood layer and vascular lamina being imported in blood electromagnetic simulation model, to blood Electromagnetic Simulation
Model carries out numerical value calculating.
The method and device for establishing blood electromagnetic simulation model of the embodiment of the present invention is in modeling process, according to blood lipid
The size of blood volume in concentration and blood electromagnetic simulation model calculates the number of blood lipid particle in blood electromagnetic simulation model,
By controlling the number of blood lipid particle, the variable blood electromagnetic simulation model of serum lipid concentrations is established, and utilize triple debye-
Drude model is fitted the electromagnetic parameter of different tissues layer in blood electromagnetic simulation model, to set its electromagnetic parameter.
The blood electromagnetic simulation model established through the invention can be used for analyzing the interaction mechanism of blood lipid and electromagnetic wave, obtain not
With the electromagnetic response characteristic under serum lipid concentrations, important theory support is provided for the development of noninvasive lipids detection technology.Meanwhile
The present invention also has many advantages, such as that use cost is low, simulation accuracy is high.
The foregoing description of the disclosed embodiments enables those skilled in the art to implement or use the present invention.
Various modifications to these embodiments will be readily apparent to those skilled in the art, as defined herein
General Principle can be realized in other embodiments without departing from the spirit or scope of the present invention.Therefore, of the invention
It is not intended to be limited to the embodiments shown herein, and is to fit to and the principles and novel features disclosed herein phase one
The widest scope of cause.
Claims (8)
1. a kind of method for establishing blood electromagnetic simulation model characterized by comprising
Step a: blood electromagnetic simulation model is established；
Step b: it according to the volume of blood in the concentration of blood lipid and blood electromagnetic simulation model, calculates in blood electromagnetic simulation model
The number of blood lipid particle；
Blood lipid particle: being randomly dispersed in blood electromagnetic simulation model by step c according to calculated result, and passes through control iipomicron
The number of son establishes the variable blood electromagnetic simulation model of serum lipid concentrations；
Wherein, in the step b, blood lipid population purpose calculation formula in the blood electromagnetic simulation model are as follows:
In above-mentioned formula, l is the number of blood lipid particle, and N is the number of neutral fat molecule and lipoid molecule, NAFor Avobenzene gal moral
Sieve constant, and NA=6.02 × 1023, M is the amount of the substance of neutral fat molecule and lipoid molecule, and M=ρ × V/1000, ρ are
Serum lipid concentrations, V are blood volume, and V=1000 × π (d/2)2× h, wherein d is the internal diameter of vascular lamina, and h is the length of vascular lamina
Degree.
2. the method according to claim 1 for establishing blood electromagnetic simulation model, which is characterized in that before the step a also
Comprise determining that the overall structure of blood electromagnetic simulation model；The overall structure of the blood electromagnetic simulation model is cylindrical body, institute
Stating blood electromagnetic simulation model includes vascular lamina, blood layer and blood lipid particle layer, the vascular lamina, blood layer and blood lipid particle layer
Distributing position be respectively as follows: vascular lamina and be located at outermost layer, the length of vascular lamina is h, and the internal diameter of vascular lamina is d, outer diameter D, blood
Tube layer with a thickness of r12=(D-d)/2；Blood layer is located in vascular lamina, and the radius of blood layer is r=d/2；Blood lipid particle layer point
Cloth is in blood layer.
3. the method according to claim 1 for establishing blood electromagnetic simulation model, which is characterized in that in the step c,
The variable blood electromagnetic simulation model of serum lipid concentrations of establishing specifically includes:
Step c1: n random numbers equally distributed between (- 1,1) are generated
Step c2: according to the central-limit theorem of probability theory, generating the random value of a Normal Distribution N (μ, σ), generates
The formula of random value are as follows:
In above-mentioned formula,For uniform random number, μ is the desired value of normal distribution random number, and σ is that normal distribution is random
Several mean square deviations, n are the numbers of uniform random number needed for generating normal distribution random number；
Step c3: the position of blood lipid particle each in blood electromagnetic simulation model is determined according to the formula of above-mentioned generation random value
It sets, establishes the variable blood electromagnetic simulation model of serum lipid concentrations.
4. the method according to claim 3 for establishing blood electromagnetic simulation model, which is characterized in that after the step c also
It include: to be carried out using triple debye-drude models to the electromagnetic parameter of blood layer and vascular lamina in blood electromagnetic simulation model
Fitting, and fitting result is imported in blood electromagnetic simulation model, numerical value calculating is carried out to blood electromagnetic simulation model.
5. a kind of device for establishing blood electromagnetic simulation model characterized by comprising
First model building module: for establishing blood electromagnetic simulation model；
Blood lipid number of particles computing module: for the volume of blood in the concentration and blood electromagnetic simulation model according to blood lipid, meter
Calculate the number of blood lipid particle in blood electromagnetic simulation model；
Second model building module: for blood lipid particle to be randomly dispersed in blood electromagnetic simulation model according to calculated result,
And by the number of control blood lipid particle, the variable blood electromagnetic simulation model of serum lipid concentrations is established；
Wherein, the blood lipid number of particles computing module calculates blood lipid population purpose formula in blood electromagnetic simulation model are as follows:
In above-mentioned formula, l is the number of blood lipid particle, and N is the number of neutral fat molecule and lipoid molecule, NAFor Avobenzene gal moral
Sieve constant, and NA=6.02 × 1023, M is the amount of the substance of neutral fat molecule and lipoid molecule, and M=ρ × V/1000, ρ are
Serum lipid concentrations, V are blood volume, and V=1000 × π (d/2)2× h, wherein d is the internal diameter of vascular lamina, and h is the length of vascular lamina
Degree.
6. the device according to claim 5 for establishing blood electromagnetic simulation model, which is characterized in that further include Structure Calculation
Module, the structural calculation module are used to determine the overall structure of blood electromagnetic simulation model；The blood electromagnetic simulation model
Overall structure be cylindrical body, the blood electromagnetic simulation model includes vascular lamina, blood layer and blood lipid particle layer, the blood vessel
The distributing position of layer, blood layer and blood lipid particle layer is respectively as follows: vascular lamina and is located at outermost layer, and the length of vascular lamina is h, vascular lamina
Internal diameter be d, outer diameter D, vascular lamina with a thickness of r12=(D-d)/2；Blood layer is located in vascular lamina, and the radius of blood layer is
R=d/2；Blood lipid particle layer is distributed in blood layer.
7. the device according to claim 5 for establishing blood electromagnetic simulation model, which is characterized in that second model is built
The mode that formwork erection block establishes the variable blood electromagnetic simulation model of serum lipid concentrations includes: to generate n to be uniformly distributed between (- 1,1)
Random numberAccording to the central-limit theorem of probability theory, a Normal Distribution N (μ, σ) is generated
Random value generates the formula of random value are as follows:
The position that blood lipid particle each in blood electromagnetic simulation model is determined according to the formula of above-mentioned generation random value, establishes blood
The variable blood electromagnetic simulation model of lipid concentration；In above-mentioned formula,For uniform random number, μ is that normal distribution is random
Several desired values, σ are the mean square deviations of normal distribution random number, and n is to be uniformly distributed needed for generating normal distribution random number at random
Several numbers.
8. the device according to claim 7 for establishing blood electromagnetic simulation model, which is characterized in that further include that parameter calculates
Module, the parameter calculating module be used for using triple debye-drude models to blood layer in blood electromagnetic simulation model and
The electromagnetic parameter of vascular lamina is fitted, and fitting result is imported in blood electromagnetic simulation model, to blood Electromagnetic Simulation mould
Type carries out numerical value calculating.
Priority Applications (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610862038.9A CN106503302B (en) | 2016-09-28 | 2016-09-28 | A kind of method and device for establishing blood electromagnetic simulation model |
| PCT/CN2016/103033 WO2018058720A1 (en) | 2016-09-28 | 2016-10-24 | Method and apparatus for establishing electromagnetic blood simulation model |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610862038.9A CN106503302B (en) | 2016-09-28 | 2016-09-28 | A kind of method and device for establishing blood electromagnetic simulation model |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106503302A CN106503302A (en) | 2017-03-15 |
| CN106503302B true CN106503302B (en) | 2019-08-23 |
Family
ID=58290050
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610862038.9A Active CN106503302B (en) | 2016-09-28 | 2016-09-28 | A kind of method and device for establishing blood electromagnetic simulation model |
Country Status (2)
| Country | Link |
|---|---|
| CN (1) | CN106503302B (en) |
| WO (1) | WO2018058720A1 (en) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108877445A (en) * | 2018-07-19 | 2018-11-23 | 天津大学 | A method of dividing blood distribution in shape reconstruct ear-lobe model based on DLA |
| CN110729050A (en) * | 2019-09-01 | 2020-01-24 | 天津大学 | Method for establishing three-dimensional earlobe model for blood sugar detection |
| CN113536626B (en) * | 2021-06-25 | 2023-05-26 | 北京航空航天大学 | DGTD electromagnetic transient simulation method based on Cole-Cole model |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102525443A (en) * | 2010-09-16 | 2012-07-04 | 西门子公司 | Method and system for non-invasive assessment of coronary artery disease |
| US8380541B1 (en) * | 2011-09-25 | 2013-02-19 | Theranos, Inc. | Systems and methods for collecting and transmitting assay results |
| CN104104450A (en) * | 2014-06-20 | 2014-10-15 | 中国科学院深圳先进技术研究院 | Human body communication channel modeling method based on non-uniform medium and system thereof |
| CN104720894A (en) * | 2015-02-11 | 2015-06-24 | 中山大学附属第一医院 | Rationality analyzing method for blood vessel operation mode |
| US9087147B1 (en) * | 2014-03-31 | 2015-07-21 | Heartflow, Inc. | Systems and methods for determining blood flow characteristics using flow ratio |
| CN105205211A (en) * | 2015-08-20 | 2015-12-30 | 电子科技大学 | Modeling method for three-dimensional electromagnetic simulation model of surface channel type mixing schottky diode |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP6260989B2 (en) * | 2013-12-05 | 2018-01-17 | 富士通株式会社 | Shape data generation apparatus, shape data generation method, and shape data generation program |
-
2016
- 2016-09-28 CN CN201610862038.9A patent/CN106503302B/en active Active
- 2016-10-24 WO PCT/CN2016/103033 patent/WO2018058720A1/en active Application Filing
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102525443A (en) * | 2010-09-16 | 2012-07-04 | 西门子公司 | Method and system for non-invasive assessment of coronary artery disease |
| US8380541B1 (en) * | 2011-09-25 | 2013-02-19 | Theranos, Inc. | Systems and methods for collecting and transmitting assay results |
| US9087147B1 (en) * | 2014-03-31 | 2015-07-21 | Heartflow, Inc. | Systems and methods for determining blood flow characteristics using flow ratio |
| CN104104450A (en) * | 2014-06-20 | 2014-10-15 | 中国科学院深圳先进技术研究院 | Human body communication channel modeling method based on non-uniform medium and system thereof |
| CN104720894A (en) * | 2015-02-11 | 2015-06-24 | 中山大学附属第一医院 | Rationality analyzing method for blood vessel operation mode |
| CN105205211A (en) * | 2015-08-20 | 2015-12-30 | 电子科技大学 | Modeling method for three-dimensional electromagnetic simulation model of surface channel type mixing schottky diode |
Non-Patent Citations (1)
| Title |
|---|
| 多电极电磁血液流速仪仿真建模研究;吴学礼等;《河北科技大学学报》;20141231;第35卷(第6期);第569-576页 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106503302A (en) | 2017-03-15 |
| WO2018058720A1 (en) | 2018-04-05 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Lai et al. | Ionically crosslinked complex gels loaded with oleic acid-containing vesicles for transdermal drug delivery | |
| CN106503302B (en) | A kind of method and device for establishing blood electromagnetic simulation model | |
| Ding et al. | Printability study of bioprinted tubular structures using liquid hydrogel precursors in a support bath | |
| Xie et al. | The effect and mechanism of transdermal penetration enhancement of Fu’s cupping therapy: new physical penetration technology for transdermal administration with traditional Chinese medicine (TCM) characteristics | |
| Taghibakhshi et al. | Three-dimensional modeling of avascular tumor growth in both static and dynamic culture platforms | |
| Li et al. | The design and evaluation of bionic porous bone scaffolds in fluid flow characteristics and mechanical properties | |
| Lu et al. | Three-dimensional microwave head imaging with GPU-based FDTD and the DBIM method | |
| Abedi et al. | A simulation-based methodology of developing 3D printed anthropomorphic phantoms for microwave imaging systems | |
| Pérez-Aliacar et al. | Predicting cell behaviour parameters from glioblastoma on a chip images. A deep learning approach | |
| Dustler et al. | Binary implementation of fractal Perlin noise to simulate fibroglandular breast tissue | |
| Chen et al. | Microwave-induced rapid shape change of 4D printed vegetable-based food | |
| Azhdari et al. | Diffusion, viscoelasticity and erosion: analytical study and medical applications | |
| Fishler et al. | A microfluidic model of biomimetically breathing pulmonary acinar airways | |
| Yu et al. | A fractional anomalous diffusion model and numerical simulation for sodium ion transport in the intestinal wall | |
| EP4146394A1 (en) | Microfluidic device | |
| Chen et al. | A heat and mass transfer model of peanut convective drying based on a two-component structure | |
| Wang et al. | Regularization solver guided FISTA for electrical impedance tomography | |
| Foroughimehr et al. | The Impact of Base Cell Size Setup on the Finite Difference Time Domain Computational Simulation of Human Cornea Exposed to Millimeter Wave Radiation at Frequencies above 30 GHz | |
| Karageorgos et al. | Prediction of viscoelastic properties of enzymatically crosslinkable tyramine–modified hyaluronic acid solutions using a dynamic monte carlo kinetic approach | |
| Huclova | Modeling of cell suspensions and biological tissue for computational electromagnetics | |
| Mai et al. | Simulation study of microwave ablation of porous lung tissue | |
| Annabi et al. | Voices of biotech research | |
| Ferrero et al. | Design and Characterization of an RF Applicator for In Vitro Tests of Electromagnetic Hyperthermia | |
| Cherubini et al. | Systems Biology Modeling of Nonlinear Cancer Dynamics | |
| Baravalle et al. | Procedural bread making |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | 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 |