CN118380155B - Method, device, equipment, medium and product for calibrating passive efficiency of human body irradiation - Google Patents
Method, device, equipment, medium and product for calibrating passive efficiency of human body irradiation Download PDFInfo
- Publication number
- CN118380155B CN118380155B CN202410563575.8A CN202410563575A CN118380155B CN 118380155 B CN118380155 B CN 118380155B CN 202410563575 A CN202410563575 A CN 202410563575A CN 118380155 B CN118380155 B CN 118380155B
- Authority
- CN
- China
- Prior art keywords
- model
- respiratory tract
- dimensional
- sub
- fluid
- 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 48
- 239000002245 particle Substances 0.000 claims abstract description 149
- 210000002345 respiratory system Anatomy 0.000 claims abstract description 139
- 238000004088 simulation Methods 0.000 claims abstract description 97
- 239000012530 fluid Substances 0.000 claims abstract description 87
- 230000008021 deposition Effects 0.000 claims abstract description 86
- 238000009826 distribution Methods 0.000 claims abstract description 56
- 238000005259 measurement Methods 0.000 claims abstract description 28
- 230000000241 respiratory effect Effects 0.000 claims abstract description 27
- 230000000704 physical effect Effects 0.000 claims abstract description 24
- 239000013078 crystal Substances 0.000 claims abstract description 13
- 229910052732 germanium Inorganic materials 0.000 claims abstract description 13
- GNPVGFCGXDBREM-UHFFFAOYSA-N germanium atom Chemical compound [Ge] GNPVGFCGXDBREM-UHFFFAOYSA-N 0.000 claims abstract description 12
- 230000008676 import Effects 0.000 claims abstract 3
- 239000000443 aerosol Substances 0.000 claims description 37
- 210000000621 bronchi Anatomy 0.000 claims description 35
- 230000005855 radiation Effects 0.000 claims description 30
- 210000001519 tissue Anatomy 0.000 claims description 28
- 230000002285 radioactive effect Effects 0.000 claims description 27
- 238000004364 calculation method Methods 0.000 claims description 26
- 210000004872 soft tissue Anatomy 0.000 claims description 16
- 238000004590 computer program Methods 0.000 claims description 14
- 210000000214 mouth Anatomy 0.000 claims description 13
- 230000029058 respiratory gaseous exchange Effects 0.000 claims description 11
- 210000003437 trachea Anatomy 0.000 claims description 11
- 210000001331 nose Anatomy 0.000 claims description 10
- 238000005295 random walk Methods 0.000 claims description 10
- 239000013598 vector Substances 0.000 claims description 10
- 230000007613 environmental effect Effects 0.000 claims description 8
- 230000008569 process Effects 0.000 claims description 8
- 238000003860 storage Methods 0.000 claims description 8
- 210000003928 nasal cavity Anatomy 0.000 claims description 7
- 230000005484 gravity Effects 0.000 claims description 6
- 238000012937 correction Methods 0.000 claims description 4
- 230000006870 function Effects 0.000 claims description 4
- 238000006073 displacement reaction Methods 0.000 claims description 3
- 210000004072 lung Anatomy 0.000 claims description 3
- 230000010352 nasal breathing Effects 0.000 claims description 3
- 208000001705 Mouth breathing Diseases 0.000 claims 2
- 238000004980 dosimetry Methods 0.000 claims 2
- 230000000877 morphologic effect Effects 0.000 claims 2
- 238000001514 detection method Methods 0.000 abstract description 5
- 238000002474 experimental method Methods 0.000 abstract description 3
- 230000001678 irradiating effect Effects 0.000 description 17
- 238000004891 communication Methods 0.000 description 11
- 238000004458 analytical method Methods 0.000 description 9
- 238000013461 design Methods 0.000 description 9
- 238000010586 diagram Methods 0.000 description 8
- 210000003800 pharynx Anatomy 0.000 description 7
- 230000008878 coupling Effects 0.000 description 6
- 238000010168 coupling process Methods 0.000 description 6
- 238000005859 coupling reaction Methods 0.000 description 6
- 238000000342 Monte Carlo simulation Methods 0.000 description 4
- 230000008859 change Effects 0.000 description 4
- 230000000694 effects Effects 0.000 description 4
- 239000007789 gas Substances 0.000 description 4
- 210000003123 bronchiole Anatomy 0.000 description 3
- 238000002591 computed tomography Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 239000000463 material Substances 0.000 description 3
- 239000000203 mixture Substances 0.000 description 3
- XKRFYHLGVUSROY-UHFFFAOYSA-N Argon Chemical compound [Ar] XKRFYHLGVUSROY-UHFFFAOYSA-N 0.000 description 2
- CURLTUGMZLYLDI-UHFFFAOYSA-N Carbon dioxide Chemical compound O=C=O CURLTUGMZLYLDI-UHFFFAOYSA-N 0.000 description 2
- 229910000831 Steel Inorganic materials 0.000 description 2
- 230000002146 bilateral effect Effects 0.000 description 2
- 239000003814 drug Substances 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 238000000338 in vitro Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000004224 protection Effects 0.000 description 2
- HCWPIIXVSYCSAN-UHFFFAOYSA-N radium atom Chemical compound [Ra] HCWPIIXVSYCSAN-UHFFFAOYSA-N 0.000 description 2
- 239000010959 steel Substances 0.000 description 2
- JFALSRSLKYAFGM-UHFFFAOYSA-N uranium(0) Chemical compound [U] JFALSRSLKYAFGM-UHFFFAOYSA-N 0.000 description 2
- ZCYVEMRRCGMTRW-UHFFFAOYSA-N 7553-56-2 Chemical compound [I] ZCYVEMRRCGMTRW-UHFFFAOYSA-N 0.000 description 1
- 235000020847 Body for Life Nutrition 0.000 description 1
- 208000006545 Chronic Obstructive Pulmonary Disease Diseases 0.000 description 1
- 206010058467 Lung neoplasm malignant Diseases 0.000 description 1
- 208000028571 Occupational disease Diseases 0.000 description 1
- 206010041067 Small cell lung cancer Diseases 0.000 description 1
- 229910052770 Uranium Inorganic materials 0.000 description 1
- 230000001133 acceleration Effects 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 230000002411 adverse Effects 0.000 description 1
- 238000003915 air pollution Methods 0.000 description 1
- 229910052782 aluminium Inorganic materials 0.000 description 1
- 229910052786 argon Inorganic materials 0.000 description 1
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 229910052790 beryllium Inorganic materials 0.000 description 1
- 239000012472 biological sample Substances 0.000 description 1
- 229910052792 caesium Inorganic materials 0.000 description 1
- TVFDJXOCXUVLDH-UHFFFAOYSA-N caesium atom Chemical compound [Cs] TVFDJXOCXUVLDH-UHFFFAOYSA-N 0.000 description 1
- 229910002092 carbon dioxide Inorganic materials 0.000 description 1
- 239000001569 carbon dioxide Substances 0.000 description 1
- 229910017052 cobalt Inorganic materials 0.000 description 1
- 239000010941 cobalt Substances 0.000 description 1
- GUTLYIVDDKVIGB-UHFFFAOYSA-N cobalt atom Chemical compound [Co] GUTLYIVDDKVIGB-UHFFFAOYSA-N 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000005520 cutting process Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 239000012634 fragment Substances 0.000 description 1
- 238000005286 illumination Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000001727 in vivo Methods 0.000 description 1
- 239000011630 iodine Substances 0.000 description 1
- 229910052740 iodine Inorganic materials 0.000 description 1
- 210000000867 larynx Anatomy 0.000 description 1
- 201000005202 lung cancer Diseases 0.000 description 1
- 208000020816 lung neoplasm Diseases 0.000 description 1
- 201000005243 lung squamous cell carcinoma Diseases 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000005025 nuclear technology Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 210000000056 organ Anatomy 0.000 description 1
- 239000001301 oxygen Substances 0.000 description 1
- 229910052760 oxygen Inorganic materials 0.000 description 1
- 230000035515 penetration Effects 0.000 description 1
- 210000004879 pulmonary tissue Anatomy 0.000 description 1
- 239000000941 radioactive substance Substances 0.000 description 1
- 230000001950 radioprotection Effects 0.000 description 1
- 229910052705 radium Inorganic materials 0.000 description 1
- 229910052704 radon Inorganic materials 0.000 description 1
- SYUHGPGVQRZVTB-UHFFFAOYSA-N radon atom Chemical compound [Rn] SYUHGPGVQRZVTB-UHFFFAOYSA-N 0.000 description 1
- 230000036387 respiratory rate Effects 0.000 description 1
- 208000023504 respiratory system disease Diseases 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 208000000587 small cell lung carcinoma Diseases 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 229910052712 strontium Inorganic materials 0.000 description 1
- CIOAGBVUUVVLOB-UHFFFAOYSA-N strontium atom Chemical compound [Sr] CIOAGBVUUVVLOB-UHFFFAOYSA-N 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 229910052724 xenon Inorganic materials 0.000 description 1
- FHNFHKCVQCLJFQ-UHFFFAOYSA-N xenon atom Chemical compound [Xe] FHNFHKCVQCLJFQ-UHFFFAOYSA-N 0.000 description 1
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
- G16H50/50—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for simulation or modelling of medical disorders
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/50—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/42—Arrangements for detecting radiation specially adapted for radiation diagnosis
- A61B6/4208—Arrangements for detecting radiation specially adapted for radiation diagnosis characterised by using a particular type of detector
- A61B6/4258—Arrangements for detecting radiation specially adapted for radiation diagnosis characterised by using a particular type of detector for detecting non x-ray radiation, e.g. gamma radiation
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
- A61B6/5205—Devices using data or image processing specially adapted for radiation diagnosis involving processing of raw data to produce diagnostic data
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/58—Testing, adjusting or calibrating thereof
- A61B6/582—Calibration
- A61B6/585—Calibration of detector units
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C20/00—Chemoinformatics, i.e. ICT specially adapted for the handling of physicochemical or structural data of chemical particles, elements, compounds or mixtures
- G16C20/30—Prediction of properties of chemical compounds, compositions or mixtures
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C20/00—Chemoinformatics, i.e. ICT specially adapted for the handling of physicochemical or structural data of chemical particles, elements, compounds or mixtures
- G16C20/70—Machine learning, data mining or chemometrics
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H10/00—ICT specially adapted for the handling or processing of patient-related medical or healthcare data
- G16H10/60—ICT specially adapted for the handling or processing of patient-related medical or healthcare data for patient-specific data, e.g. for electronic patient records
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
- G16H50/20—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for computer-aided diagnosis, e.g. based on medical expert systems
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
- G16H50/30—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for calculating health indices; for individual health risk assessment
Landscapes
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Public Health (AREA)
- General Health & Medical Sciences (AREA)
- Biomedical Technology (AREA)
- Pathology (AREA)
- Veterinary Medicine (AREA)
- Data Mining & Analysis (AREA)
- Radiology & Medical Imaging (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Heart & Thoracic Surgery (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- High Energy & Nuclear Physics (AREA)
- Biophysics (AREA)
- Physics & Mathematics (AREA)
- Primary Health Care (AREA)
- Epidemiology (AREA)
- Databases & Information Systems (AREA)
- Optics & Photonics (AREA)
- Chemical & Material Sciences (AREA)
- Computing Systems (AREA)
- Theoretical Computer Science (AREA)
- Bioinformatics & Computational Biology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Crystallography & Structural Chemistry (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Artificial Intelligence (AREA)
- Software Systems (AREA)
- Evolutionary Computation (AREA)
- Dentistry (AREA)
- Oral & Maxillofacial Surgery (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本申请公开了一种人体内照射无源效率刻度方法、装置、设备、介质及产品，涉及人体内照射检测技术领域。所述方法是先重建得到目标人体的呼吸道三维模型，然后根据呼吸道内流体的物理属性，使用有限元分析软件计算得到在目标环境下且在人体呼吸时，流体在呼吸道三维模型内部的流场分布情况，以及计算仿真颗粒在流场分布情况下的运动情况，得到仿真颗粒在呼吸道内壁上的沉积信息，再然后将沉积信息导入到蒙特卡洛模拟软件中进行射线发射仿真，得到射线在目标探测器的蒙特卡洛模型中的锗晶体处的测量数据，最后根据测量数据完成对探测器的无源效率刻度任务，如此不但可有效提高仿真实验精确度，还可针对不同人群及探测器进行综合实验，有效提高实用性。
The present application discloses a method, device, equipment, medium and product for calibrating passive efficiency of human body irradiation, and relates to the technical field of human body irradiation detection. The method is to first reconstruct a three-dimensional respiratory model of a target human body, and then use finite element analysis software to calculate the flow field distribution of the fluid inside the three-dimensional respiratory model under the target environment and when the human body breathes according to the physical properties of the fluid in the respiratory tract, and calculate the movement of the simulated particles under the flow field distribution, and obtain the deposition information of the simulated particles on the inner wall of the respiratory tract, and then import the deposition information into the Monte Carlo simulation software to simulate the ray emission, and obtain the measurement data of the ray at the germanium crystal in the Monte Carlo model of the target detector, and finally complete the passive efficiency calibration task of the detector according to the measurement data, so that not only the accuracy of the simulation experiment can be effectively improved, but also comprehensive experiments can be carried out for different populations and detectors, effectively improving practicality.
Description
Technical Field
The invention belongs to the technical field of human body internal irradiation detection, and particularly relates to a method, a device, equipment, a medium and a product for passive efficiency calibration of human body internal irradiation.
Background
Air pollution in a radioactive workplace refers to tiny colloidal particles or gases of a radioactive nature suspended in the air, which may be released to a large extent to the surrounding environment, ultimately adversely affecting the human body. After uptake of radionuclides into the body, they will deposit and lodge in the tissues or organs to which they are attached, creating a dose of radiation to them, causing damage. Wherein the radioactive aerosol is critical for assessing respiratory tract exposure dose for professionals. The traditional calculation method assumes that the radionuclides in the aerosol enter the respiratory tract of the staff completely through respiration and adopts the default value of the international radioprotection committee (International Commission on Radiological Protection, abbreviated as ICRP) as the intra-nuclide irradiation conversion factor. However, in practice, the deposition position and deposition amount of the aerosol in the respiratory tract depend on the diameter size and the respiratory rate of the human body. Thus, it is simply assumed that the aerosol is deposited completely on the respiratory tract or a specific region, which may be greatly different from the actual situation.
For the radioactive aerosols discharged to the atmosphere in the related activities such as the current global atmospheric nuclear test, nuclear facility accident and nuclear technology application, the radioactive aerosols mainly commonly comprise argon (Ar-41), cobalt (Co-60), strontium (Sr-90), iodine (I-131), xenon (Xe-133), cesium (Cs-137), radon (Rn-222), radium (Ra-226) and/or uranium (U-235) and the like, and when the aerosols containing the radioactive substances are inhaled by a human body, the aerosols can be irradiated on corresponding deposition positions of the human body for a period of time, and certain radioactive nuclides can exist in the human body for life due to the factors such as long decay period. For the evaluation of the internal irradiation dose, the common detection methods are selected from the group consisting of in vitro direct measurement, biological sample analysis, and air sampling analysis. In practical situations, one or more measurement methods are selected according to specific situations to measure.
In general, in vitro direct measurement is the most effective method of making an internal radiation dose estimate by bringing a radiation detector close to the body to detect radiation emitted from within the body and which can pass through the body tissue. Can be used for detecting radionuclides capable of releasing characteristic X-rays, gamma rays, positrons and high-energy beta particles. The apparatus for direct measurement consists of one or more high efficiency detectors mounted in a low background environment.
Before the measurement device is used for measuring the human body, the detection efficiency of the related measurement equipment is firstly required to be calibrated, and the conventional efficiency calibration method comprises an active efficiency curve method and a passive efficiency calibration method. The active efficiency curve method has the problems of complex operation, time consumption, difficulty in error evaluation and the like, and is gradually replaced by the passive efficiency scale technology at present.
In the existing passive efficiency calibration technology, the efficiency calibration of the irradiation dose in the human body is mainly related to the human body model recommended based on ICRP, but in reality, the individual difference of the human body is obvious, the conditions of radiation pollution in different human bodies are different under the same environment, and the method for calibrating the efficiency of the detection instrument by only using the human body model is rough.
Disclosure of Invention
The invention aims to provide a method, a device, computer equipment, a computer readable storage medium and a computer program product for irradiating passive efficiency scales in a human body, which are used for solving the problem that the accuracy of the efficiency scales cannot be improved due to personalized differences of the human body and lack of reference to more physical factors in the conventional technology for irradiating passive efficiency scales in the human body.
In order to achieve the above purpose, the present invention adopts the following technical scheme:
In a first aspect, a method for illuminating passive efficiency scales in a human body is provided, comprising:
Reconstructing to obtain a respiratory tract three-dimensional model of a target human body, wherein the respiratory tract three-dimensional model comprises a front nasal passage sub-model, a rear nasal passage sub-model, a nasal opening sub-model, a throat sub-model and a main bronchus bifurcation structure sub-model;
According to the physical properties of fluid in the respiratory tract, calculating by using finite element analysis software to obtain the flow field distribution condition of the fluid in the respiratory tract in a target environment and when a human body breathes, wherein the fluid in the respiratory tract is a turbulent constant-temperature incompressible Newtonian fluid, and the environmental parameters of the target environment comprise a gravity field parameter and a temperature field parameter;
Adding simulation particles with physical properties of aerosol into the three-dimensional respiratory tract model, and calculating the motion condition of the simulation particles under the flow field distribution condition by using the finite element analysis software to obtain the deposition information of the simulation particles on the inner wall of the three-dimensional respiratory tract model, wherein the deposition information comprises deposition positions or deposition positions and deposition quantity;
The respiratory tract three-dimensional model is imported into Monte Carlo simulation software, a human tissue three-dimensional model for simulating radiation scattering and energy attenuation in a soft tissue structure is added around the respiratory tract three-dimensional model, and a Monte Carlo model of a target detector is further placed around the human tissue three-dimensional model;
The deposition information is imported into Monte Carlo simulation software to carry out ray emission simulation, and measurement data of rays at germanium crystals in a Monte Carlo model of the target detector are obtained;
and according to the measurement data, completing the passive efficiency calibration task of the target detector.
Based on the above-mentioned invention, a new scheme of coupling multiple physical fields finite element analysis to improve and finish the irradiation passive efficiency scale in human body is provided on the basis of traditional Monte Carlo simulation, namely, firstly, rebuild to obtain the respiratory tract three-dimensional model of the target human body, then, according to the physical property of fluid in respiratory tract, calculate to obtain the flow field distribution condition of fluid in the respiratory tract three-dimensional model under the target environment and when the human body breathes by using finite element analysis software, and calculate the motion condition of simulation particles under the flow field distribution condition, obtain the deposition information of simulation particles on the inner wall of respiratory tract, then, the deposition information is imported into Monte Carlo simulation software to perform ray emission simulation, obtain the measurement data of rays at germanium crystal in Monte Carlo model of the target detector, finally, finish the passive efficiency scale task of the detector according to the measurement data, thus, compared with the traditional scheme, the accuracy of simulation experiment can be effectively improved by comprehensively calculating specific model structure, gravitational field, temperature field, fluid flow field and the like, and the method can also be used for different people, the comprehensive experiment can be effectively improved, the practical application is facilitated.
In one possible design, reconstructing the three-dimensional model of the respiratory tract of the target human body includes:
acquiring CT data of the respiratory tract of a target human body;
Establishing an upper respiratory tract three-dimensional model of the target human body according to the respiratory tract CT data, wherein the upper respiratory tract three-dimensional model comprises a front nasal passage sub-model, a rear nasal passage sub-model, a nasal opening sub-model, a throat sub-model and a main air pipe sub-model;
Establishing a multi-stage bronchus sub-model which is communicated with the main bronchus sub-model and has a bilateral symmetry structure based on the main bronchus sub-model, obtaining a main bronchus bifurcation structure sub-model comprising the main bronchus sub-model and the multi-stage bronchus sub-model, and obtaining a respiratory tract three-dimensional model of the target human body and comprising the anterior nasal passage sub-model, the posterior nasal passage sub-model, the nasal opening sub-model, the laryngeal sub-model and the main bronchus bifurcation structure sub-model, wherein the included angle of each bronchus bifurcation position of the multi-stage bronchus sub-model is subjected to fillet length as follows Round corners of the steel plate are processed smoothly,Represents the cross-sectional diameter of the associated bronchi at the bifurcation of the bronchi.
In one possible design, adding simulation particles with physical properties of aerosol into the respiratory tract three-dimensional model, and calculating the motion condition of the simulation particles under the condition of flow field distribution by using the finite element analysis software to obtain the deposition information of the simulation particles on the inner wall of the respiratory tract three-dimensional model, wherein the method comprises the following steps:
Obtaining initial configuration information of simulation particles, wherein the simulation particles have physical properties of aerosol, the initial configuration information comprises initial coordinates, initial speed, particle density, particle diameter and particle mass flow, the initial coordinates are at a respiratory inlet cross section of the respiratory tract three-dimensional model, the initial speed is consistent with the initial fluid speed of fluid in the respiratory tract, and the particle mass flow refers to total mass of particles passing through the respiratory inlet cross section of the respiratory tract three-dimensional model in one second;
According to the initial configuration information of the simulation particles, the simulation particles are added into the three-dimensional model of the respiratory tract, and the motion condition of the simulation particles under the distribution condition of the flow field is calculated by using the finite element analysis software, so that the deposition information of the simulation particles on the inner wall of the three-dimensional model of the respiratory tract is obtained, wherein the deposition information comprises deposition positions or deposition positions and deposition quantity.
In one possible design, the calculating the motion of the simulated particles in the flow field distribution using the finite element analysis software includes:
simulating random movements of the simulated particles in the interior space of the three-dimensional model of the respiratory tract using a random walk model;
Or for simulated particles having a size dimension greater than the fluid molecules in the respiratory tract, using a random walk model and a Stokes-Cunningham drag model to simulate random movement of the corresponding particles in the interior space of the three-dimensional model of the respiratory tract.
In one possible design, importing the deposition information into the monte carlo simulation software for radiation emission simulation includes:
Importing the deposition position in the deposition information into Monte Carlo simulation software to serve as a position of a radioactive aerosol which becomes a radioactive source and generates radioactive rays after being deposited on the inner wall of the respiratory tract three-dimensional model;
the emission direction of the radioactive rays generated by the radioactive source is calculated by simulation using a random trigonometric function equation.
In one possible design, the physical properties of the fluid within the respiratory tract include fluid density, fluid viscosity, fluid initial velocity, fluid turbulence kinetic energy, and fluid turbulence dissipation ratio;
and/or the flow field distribution condition comprises a fluid flow velocity distribution condition and a fluid pressure distribution condition.
The second aspect provides a passive efficiency scale device for human body internal irradiation, which comprises a three-dimensional model reconstruction unit, a flow field distribution calculation unit, a particle motion calculation unit, a simulation model arrangement unit, a ray emission simulation unit and a scale task completion unit;
The three-dimensional model reconstruction unit is used for reconstructing and obtaining a respiratory tract three-dimensional model of a target human body, wherein the respiratory tract three-dimensional model comprises a front nose channel sub-model, a rear nose channel sub-model, a nose mouth sub-model, a throat sub-model and a main bronchus bifurcation structure sub-model;
The flow field distribution calculation unit is in communication connection with the three-dimensional model reconstruction unit and is used for calculating the flow field distribution condition of the fluid in the respiratory tract three-dimensional model under the target environment and when a human body breathes by using finite element analysis software, wherein the fluid in the respiratory tract is turbulent constant-temperature incompressible Newtonian fluid, and the environmental parameters of the target environment comprise gravity field parameters and temperature field parameters;
The particle motion calculation unit is in communication connection with the flow field distribution calculation unit and is used for adding simulation particles with physical properties of aerosol into the three-dimensional model of the respiratory tract, calculating the motion condition of the simulation particles under the flow field distribution condition by using the finite element analysis software, and obtaining the deposition information of the simulation particles on the inner wall of the three-dimensional model of the respiratory tract, wherein the deposition information comprises deposition positions or deposition positions and deposition quantity;
the simulation model arrangement unit is in communication connection with the three-dimensional model reconstruction unit and is used for importing the respiratory tract three-dimensional model into Monte Carlo simulation software, adding a human tissue three-dimensional model for simulating radiation scattering and energy attenuation in a soft tissue structure around the respiratory tract three-dimensional model, and placing a Monte Carlo model of a target detector around the human tissue three-dimensional model;
The radiation emission simulation unit is respectively in communication connection with the particle motion calculation unit and the simulation model arrangement unit and is used for importing the deposition information into the Monte Carlo simulation software to perform radiation emission simulation so as to obtain measurement data of radiation at a germanium crystal position in a Monte Carlo model of the target detector;
and the calibration task completion unit is in communication connection with the ray emission simulation unit and is used for completing the passive efficiency calibration task of the target detector according to the measurement data.
In a third aspect, the present invention provides a computer device comprising a memory, a processor and a transceiver in communication connection in sequence, wherein the memory is configured to store a computer program, the transceiver is configured to send and receive messages, and the processor is configured to read the computer program and perform the method of irradiating passive efficiency scales in a human body according to the first aspect or any of the possible designs of the first aspect.
In a fourth aspect, the present invention provides a computer readable storage medium having instructions stored thereon which, when executed on a computer, perform the method of irradiating passive efficiency scales in the human body as described in the first aspect or any of the possible designs of the first aspect.
In a fifth aspect, the present invention provides a computer program product comprising a computer program or instructions which, when executed by a computer, implement the method of irradiating passive efficiency scales in the human body as described in the first aspect or any of the possible designs of the first aspect.
The beneficial effect of above-mentioned scheme:
(1) The invention creatively provides a new scheme for improving and completing irradiation passive efficiency scale in human body by coupling finite element analysis of multiple physical fields on the basis of traditional Monte Carlo simulation, namely, firstly reconstructing to obtain a respiratory tract three-dimensional model of a target human body, then calculating by using finite element analysis software according to physical properties of fluid in the respiratory tract to obtain flow field distribution conditions of the fluid in the respiratory tract three-dimensional model in the target environment and when the human body breathes, calculating movement conditions of simulation particles under the flow field distribution conditions to obtain deposition information of the simulation particles on the inner wall of the respiratory tract, then introducing the deposition information into Monte Carlo simulation software to perform ray emission simulation to obtain measurement data of rays at germanium crystal in the Monte Carlo model of the target detector, and finally completing passive efficiency scale task of the detector according to the measurement data;
(2) The method can be used for detecting the efficiency scale of the internal irradiation detector during human respiration more efficiently and accurately, and the theory of the method can be used in other fields, such as medicine flow rate and atomized particle size control during atomized medicine inhalation treatment;
(3) The accuracy of the whole experiment is improved by coupling finite element analysis of multiple physical fields on the basis of traditional Monte Carlo simulation, and the optimal efficiency scale can be realized;
(4) According to the actual condition of the patient, when the patient is simulated to have chronic obstructive pulmonary disease, the bronchus fragments of part of respiratory tract can not effectively carry out smooth circulation of oxygen or carbon dioxide;
(5) The distribution of main deposition positions of substances such as aerosol particles in the respiratory tract under the action of human respiration can be accurately simulated, and reference basis is provided for the treatment of related radioactive occupational diseases in the later period;
(6) Compared with the traditional internal irradiation passive efficiency scale, the passive efficiency scale of the internal irradiation detector can be higher, namely the traditional passive efficiency scale is formed by assuming that the radioactive source is uniformly distributed in the respiratory tract or is singly and quantitatively arranged according to an empirical formula, so that the radiation source is inconsistent with the actual current situation in Monte Carlo analysis, a certain error is generated in the result, the passive efficiency scale development of the internal irradiation detector under the future dose is not facilitated, and the scheme is free of concern.
Drawings
In order to more clearly illustrate the embodiments of the invention or the technical solutions in the prior art, the drawings that are required in the embodiments or the description of the prior art will be briefly described, it being obvious that the drawings in the following description are only some embodiments of the invention, and that other drawings may be obtained according to these drawings without inventive effort for a person skilled in the art.
Fig. 1 is a flow chart of a method for illuminating passive efficiency scales in a human body according to an embodiment of the present application.
Fig. 2 is an exemplary diagram of a three-dimensional model of the respiratory tract provided in an embodiment of the present application.
Fig. 3 is a diagram showing a combination example of a three-dimensional model of the respiratory tract and a three-dimensional model of human tissue for irradiating passive efficiency scales in a human body according to an embodiment of the present application.
Fig. 4 is an exemplary diagram of a monte carlo model of a target detector for illuminating a passive efficiency scale in a human body provided in an embodiment of the present application.
Fig. 5 is a diagram illustrating a positional relationship between a three-dimensional model of the respiratory tract and a three-dimensional model of human tissue and a monte carlo model of a target detector according to an embodiment of the present application.
Fig. 6 is a schematic diagram of flow velocity distribution of an intra-respiratory tract flow field after finite element analysis according to an embodiment of the present application, wherein (a) in fig. 6 shows a front view and (b) in fig. 6 shows a side view.
Fig. 7 is a schematic diagram of the distribution of particles deposited on the inner wall of the respiratory tract after finite element analysis according to an embodiment of the present application.
Fig. 8 is a graph showing the distribution trend of particles with different particle diameters in the respiratory tract according to the embodiment of the present application.
Fig. 9 is a graph showing an example of radiation energy of the detector detecting different aerosol particles with different diameters according to an embodiment of the present application, wherein (a) in fig. 9 shows an example of radiation energy of the detector detecting Co-60-based aerosol particles with different diameters, and (b) in fig. 9 shows an example of radiation energy of the detector detecting Cs-137-based aerosol particles with different diameters.
Fig. 10 is a schematic structural diagram of a passive efficiency scale device for human body internal irradiation according to an embodiment of the present application.
Fig. 11 is a schematic structural diagram of a computer device according to an embodiment of the present application.
In the drawings, the composition comprises 11-nasal cavity, 12-oral cavity, 13-main trachea, 14-bronchus, 20-lung tissue, 30-soft tissue, 101-beryllium window, 102-aluminum shell, 103-germanium crystal and 104-vacuum environment cavity.
Detailed Description
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the present invention will be briefly described below with reference to the drawings and the description of the embodiments or the prior art, and it is apparent that the following descriptions about the structures of the drawings are only some embodiments of the present invention, and other embodiments may be obtained according to these embodiments without inventive effort to those of ordinary skill in the art. It should be noted that the description of these examples is for aiding in understanding the present invention, but is not intended to limit the present invention.
It should be understood that although the terms first and second, etc. may be used herein to describe various objects, these objects should not be limited by these terms. These terms are only used to distinguish one object from another. For example, a first object may be referred to as a second object, and similarly a second object may be referred to as a first object, without departing from the scope of example embodiments of the invention.
It should be understood that for the term "and/or" that may appear herein, it is merely one kind of association relation describing the associated object, it may be indicated that there may be three kinds of relations, for example, a and/or B, it may be indicated that there are three kinds of cases, for example, a alone, B alone or both a and B, and for example, A, B and/or C may be indicated that there is any one of A, B and C or any combination thereof, for the term "/and" that may appear herein, it is another kind of association object relation, it may be indicated that there may be two kinds of relations, for example, a/and B, it may be indicated that there are a alone or both a and B together, and in addition, for the character "/", that may appear herein, it is generally indicated that the associated object is an "or" relation.
Examples
As shown in fig. 1, the method for illuminating passive efficiency scales in a human body according to the first aspect of the present embodiment may be, but is not limited to, executed by a computer device having a certain computing resource, for example, executed by a computer device such as a cloud server. As shown in FIG. 1, the method for irradiating passive efficiency scales in a human body can include, but is not limited to, the following steps S1 to S6.
S1, reconstructing to obtain a three-dimensional respiratory tract model of a target human body, wherein the three-dimensional respiratory tract model comprises, but is not limited to, a front nasal passage sub-model, a rear nasal passage sub-model, a nasal opening sub-model, a throat sub-model, a main bronchus bifurcation structure sub-model and the like.
In the step S1, the target human body may be a standard human body of a certain group of people, or may be a personalized human body of a certain person. Specifically, the three-dimensional model of the respiratory tract of the target human body is obtained through reconstruction, including but not limited to the following steps S11-S13.
S11, acquiring CT (Computed Tomography electronic computer tomography) data of the respiratory tract of a target human body.
S12, establishing an upper respiratory tract three-dimensional model of the target human body according to the respiratory tract CT data, wherein the upper respiratory tract three-dimensional model comprises, but is not limited to, a front nasal passage sub-model, a rear nasal passage sub-model, a nasal opening sub-model, a throat sub-model, a main air pipe sub-model and the like.
In the step S12, the reconstruction of the three-dimensional model of the upper respiratory tract may be specifically performed by related medical software (e.g. MIMICS MEDICAL). In the construction process of the upper respiratory tract three-dimensional model, the transverse cutting treatment can be carried out on the nasal cavity opening at the front end according to the morphology characteristics of the dose and the human respiratory tract, so that the end face is in a geometric coordinate plane, and the method can be used for the kinematic finite element simulation in the later period of nasal respiration, meanwhile, the oral cavity is constructed into a lip tight state and an oral cavity micro-opening state according to the normal respiration state of the human body, the lip is divided into an upper lip and a lower lip, in the later period of simulation, the opening angle of the lip can be controlled, and the kinematic finite element simulation of the oral respiration of the human body in different states is simulated.
S13, establishing a multi-stage bronchus sub-model which is communicated with the main bronchus sub-model and has a bilateral symmetry structure based on the main bronchus sub-model, obtaining a main bronchus bifurcation structure sub-model comprising the main bronchus sub-model and the multi-stage bronchus sub-model, and obtaining a respiratory tract three-dimensional model of the target human body and comprising the anterior nasal passage sub-model, the posterior nasal passage sub-model, the nasal opening sub-model, the laryngeal sub-model and the main bronchus bifurcation structure sub-model, wherein the fillet length of an included angle of each bronchus bifurcation position of the multi-stage bronchus sub-model is as followsRound corners of the steel plate are processed smoothly,Represents the cross-sectional diameter of the associated bronchi at the bifurcation of the bronchi.
In the step S13, it is considered that the lung cancer such as small cell lung cancer or squamous cell carcinoma after radiation is mainly focused on the main bronchiole region (i.e., the main bronchus to bronchiole may be divided into 16-stage structures, the main bronchiole region is 0-8-stage structures, wherein the 0-stage is the main bronchus, and since the airway is long and curved in real circumstances, the CT scan data is used for the establishment), and thus the multistage bronchi model may be specifically 1-8-stage structures, i.e., 316 solid models may be established in total, so as to be used for simulating the respiratory model of a patient suffering from a certain respiratory disease in a specific environment. For example, the three-dimensional model of the respiratory tract may be as shown in fig. 2. In addition, it is considered that aerosol is easy to accumulate at the bifurcation of the airway along with respiration, and in a real human body tracheal model, each bifurcation of the bronchus needs to present a certain radian, and because ICRP lacks related data, in an actual three-dimensional model, rounding and rounding of the length of the included angle at the bifurcation of the bronchus is also required, and the data of the related structure are shown in the following table 1.
TABLE 1 parameter Table of multistage bronchial model
S2, calculating the flow field distribution condition of the fluid in the respiratory tract three-dimensional model by using finite element analysis software according to the physical properties of the fluid in the respiratory tract, wherein the fluid in the respiratory tract is a turbulent constant-temperature incompressible Newtonian fluid, and the environmental parameters of the target environment comprise gravitational field parameters and temperature field parameters.
In the step S2, the intra-airway fluid is air fluid, and considering that the air flow in the airway is affected by multiple factors such as the geometry of the airway, the viscosity and inertia of the air flow, the speed, direction and turbulence of the air flow entering the airway, and the final simulation result will directly affect the deposition of aerosol in the airway, so that the physical properties of the intra-airway fluid include, but are not limited to, the fluid density, the fluid viscosity, the fluid initial velocity, the fluid turbulence kinetic energy, the fluid turbulence dissipation rate, and the like, wherein the fluid density and the fluid viscosity are specifically set individually according to the simulation environment, for example, the fluid density is 1.225kg/m 3 and the fluid viscosity is 0.000017894 kg/(m.s), respectively, and the initial inlet velocity of the fluid initial velocity at the nasal cavity or oral inlet is set(Wherein,Indicating the tidal volume of the breath,Indicating the duration of an incoming call during a breath,Representing the cross-sectional area at the nasal cavity orifice), fluid turbulence kinetic energyThe calculation formula of (2) is(Wherein,Representing fluid density; representing time; Representing a fluid velocity component; Representing the spatial coordinates; Representing the jth component in the spatial coordinates; Represents the dynamic viscosity of the fluid; indicating turbulent viscosity by correction and Calculating relevant model parameters; Typically a value of 1.0 or 1.3; the term representing the generation of turbulent energy, typically calculated by modeling), fluid turbulence dissipation ratio The calculation formula of (2) is(Wherein,Representing the proportional relationship between the generation of a determined turbulent dissipation rate and the consumption of turbulent energy, generally inSet to 0.9 in the model; Representing the proportional relationship between the degree of dissipation determining the turbulence dissipation ratio and the magnitude of the turbulence energy, generally in Set to 1.44 in the model). It is also considered that when aerosol deposition simulation is performed, the simulation of the air flow and pressure in the respiratory tract during human breathing can provide basic information, and accurate boundary conditions and input parameters are provided for subsequent aerosol deposition simulation, so that a more reliable simulation result is obtained, and therefore, the step needs to be performed before the subsequent step S3. Further, the environmental parameters of the target environment include, but are not limited to, gravitational field parameters (e.g., set gravitational acceleration at 9.81m/s 2, direction parallel to the level 0 main gas pipe and vertically downward) and temperature field parameters (e.g., set environmental temperature at 298.15K), and the flow field distribution conditions include, but are not limited to, fluid flow velocity distribution conditions, fluid pressure distribution conditions, and the like, as shown in FIG. 6.
S3, adding simulation particles with physical properties of aerosol into the three-dimensional model of the respiratory tract, and calculating the motion condition of the simulation particles under the flow field distribution condition by using finite element analysis software to obtain the deposition information of the simulation particles on the inner wall of the three-dimensional model of the respiratory tract, wherein the deposition information comprises deposition positions or deposition positions and deposition quantity.
In the step S3, after the finite element analysis software is used to calculate the convergence of the continuous phase (i.e. air) flow to the flow field, particles are added to perform coupling solution. In particular, after particles such as aerosol enter breath, the particles can be deposited or adhered on the inner wall of the respiratory tract along with the change of the movement direction when the airflow flows through complex mechanisms such as an oral cavity, a pharynx, a larynx, a trachea, a bronchus and the like, and part of the particles can be directly captured by the inner wall of the respiratory tract under the influence of factors such as gravity, air resistance, vortex and the like. Different physical parameters are required to be set aiming at the influence of the key factors, and corresponding modification is carried out on particles with different particle sizes, so that an optimal simulation result is achieved. Specifically, a calculation mode can be adopted to calculate the deposition positions and the shares of the particles in the respiratory tract under different particle diameters in a simulation mode, so that a relatively accurate simulation result is obtained, a guarantee is provided for the subsequent passive efficiency scale, namely, preferably, the simulation particles with the physical properties of aerosol are added into the respiratory tract three-dimensional model, and the motion condition of the simulation particles under the condition of flow field distribution is calculated by using the finite element analysis software, so that the deposition information of the simulation particles on the inner wall of the respiratory tract three-dimensional model is obtained, wherein the deposition information comprises the following steps S31-S32.
S31, obtaining initial configuration information of simulation particles, wherein the simulation particles have physical properties of aerosol, the initial configuration information comprises, but is not limited to, initial coordinates, initial speed, particle density, particle diameter, particle mass flow and the like, the initial coordinates are at a respiratory inlet cross section of the respiratory tract three-dimensional model, the initial speed is consistent with the initial fluid speed of fluid in the respiratory tract, and the particle mass flow refers to the total mass of particles passing through the respiratory inlet cross section of the respiratory tract three-dimensional model in one second.
In the step S31, aerosol deposition simulation may be performed by respectively creating the simulation particles of 1um, 3um, 5um, 7um and 10um according to design requirements, wherein the particle density is 3000kg/m 3, and 15000 particles are injected into 0-1.5S according to the respiratory cycle of 3S of the human body to simulate aerosol inhaled by the human body.
S32, adding the simulation particles into the three-dimensional model of the respiratory tract according to the initial configuration information of the simulation particles, and calculating the motion condition of the simulation particles under the distribution condition of the flow field by using the finite element analysis software to obtain the deposition information of the simulation particles on the inner wall of the three-dimensional model of the respiratory tract, wherein the deposition information comprises but is not limited to deposition positions or deposition positions and deposition quantity.
In the step S32, since the influence of the fluid turbulence on the particle motion track is specifically considered when the coupling calculation is performed on the particles and the fluid, that is, the motion situation of the simulated particles under the flow field distribution situation is preferably calculated by using the finite element analysis software, including using a Random Walk model (Random Walk) to simulate the Random motion of the simulated particles in the internal space of the respiratory tract three-dimensional model. The random walk model can be expressed in three dimensions asWherein, the method comprises the steps of, wherein,Is shown in the first(Representing a positive integer) the position vector of the object at the step,Is shown in the firstThe position vector of the object in the step,Is a random vector which is independently and uniformly distributed and represents the firstRandom displacement vectors of steps.
In the step S32, the air resistance of the tiny particles moving in the gas is considered, especially under the condition of low reynolds number, when the particle size is large relative to the molecular size of the gas, a Stokes-Cunningham drag model (the calculation formula is thatWherein, the method comprises the steps of, wherein,Represents the average free path of air molecules, takes a value of 69.1 nm at normal temperature and normal pressure,The particle diameter is indicated as the diameter of the particles,A bottom representing a natural logarithm), i.e., preferably, calculating the motion of the simulated particles in the flow field distribution using the finite element analysis software, comprising simulating random motion of corresponding particles in the interior space of the three-dimensional model of the respiratory tract using a random walk model and a Stokes-Cunningham drag model for simulated particles having a dimension greater than the dimension of the size of the simulated particles of fluid molecules in the respiratory tract. In addition, specifically, the correction factors of the Stokes-Cunningham drag model may be set to 1.1737, 1.0579, 1.0347, 1.0248 and 1.0174 according to the particle diameters of 1um, 3um, 5um, 7um and 10um, respectively.
In said step S32, as shown in fig. 7 and 8, for the model wall boundary conditions, it is specifically defined that the particles are supposed to be deposited at this position when they come into contact with the wall boundary, that the pressure outlet interface is defined at the 8 th-stage bronchial end section, and that the particles passing through the 8 th-stage bronchial end section are supposed to be particles to enter the bronchiolar region, and that the simulation of passive efficiency scale is not performed by taking into account the late monte carlo. In addition, the deposition position is the space coordinates (X, Y, Z) of the position where the particles are deposited, and the deposition information can also comprise the deposition time.
S4, importing the three-dimensional model of the respiratory tract into Monte Carlo simulation software, adding a human tissue three-dimensional model for simulating the scattering and energy attenuation of rays in a soft tissue structure around the three-dimensional model of the respiratory tract, and placing a Monte Carlo model of a target detector around the three-dimensional model of the human tissue.
In the step S4, the three-dimensional model of human tissue specifically includes, but is not limited to, a three-dimensional model of pulmonary tissue, a three-dimensional model of adipose layer soft tissue, a three-dimensional model of skin layer soft tissue, and the like, as shown in fig. 3. The three-dimensional model of human tissue is composed of tissue equivalent materials, and detailed components of the related tissue equivalent materials are shown in the following table 2.
TABLE 2 composition of tissue equivalent material elements
The whole model is considered to be placed in a vacuum environment close to reality and used for simulating ionization of rays in air, and meanwhile, soft tissues should be attached to lung tissues as much as possible. Wherein the composition of the air is carbon-C0.000124, nitrogen-N0.755268, oxygen-O0.231781 and argon-Ar 0.012827.
In the step S4, the target detector is an in-vivo illumination measurement tool with passive efficiency calibration to be performed, and the monte carlo model thereof may be as shown in fig. 4. Specifically, the Monte Carlo model may be placed in the middle of the respiratory tract approximately 1-2 mm from the outer layer of soft tissue, as shown in FIG. 5.
S5, importing the deposition information into the Monte Carlo simulation software to perform ray emission simulation, and obtaining measurement data of rays at germanium crystals in a Monte Carlo model of the target detector.
In the step S5, the deposition information is specifically imported into Monte Carlo simulation software to perform radiation emission simulation, including but not limited to importing the deposition position in the deposition information into Monte Carlo simulation software to serve as a position of a radioactive aerosol which becomes a radioactive source and generates radioactive rays after being deposited on the inner wall of the three-dimensional model of the respiratory tract, and simulating and calculating the emission direction of the radioactive rays generated by the radioactive source by adopting a random trigonometric function equation. In order to ensure the accuracy of the simulation experiment, the simulation particles emit rays 5000000 times, and in order to prevent random errors, the simulation experiment uses the energy of characteristic X-rays of radioactive aerosol as the energy of a radioactive source, wherein Co-60:1.33MeV and Cs-137:0.661MeV are used. In addition, considering that in practical situations, the efficiency of the detector for detecting the radiation irradiated in the human body is generally low, and the detector generally needs to detect the radiation for a quite long time at the same position, for the simulation of Monte Carlo, the number of particles needs to be transmitted for at least millions of iterations according to the input coordinates so that the detector can detect effective information about the energy of the particles. According to the ray penetration law, gamma rays > X rays > alpha rays or beta rays are generally adopted, and in a simulation experiment, in order to improve the operation efficiency, gamma rays are selected as rays generated during the decay of the related radioactive aerosol to carry out the simulation experiment.
In the step S5, specific processes of obtaining measurement data of rays at germanium crystals in a Monte Carlo model of the target detector include, but are not limited to, counting deposition energy of rays at germanium crystals in the Monte Carlo model of the target detector, and plotting a change trend comparison chart according to particle size change or particle source energy change, as shown in FIG. 9.
S6, completing the passive efficiency calibration task of the target detector according to the measurement data.
In step S6, the specific process of completing the passive efficiency calibration task according to the measurement data is the prior art, for example, an energy-efficiency curve is obtained by energy spectrum contrast.
The method comprises the steps of providing a new scheme of improving and completing the irradiation passive efficiency scale in the human body by coupling the finite element analysis of multiple physical fields on the basis of the traditional Monte Carlo simulation based on the irradiation passive efficiency scale method in the human body described in the steps S1-S6, namely, firstly reconstructing to obtain a three-dimensional respiratory tract model of a target human body, then calculating by using finite element analysis software according to the physical properties of fluid in the respiratory tract to obtain the flow field distribution condition of the fluid in the three-dimensional respiratory tract model in the target environment and when the human body breathes, calculating the movement condition of simulation particles in the flow field distribution condition to obtain the deposition information of the simulation particles on the inner wall of the respiratory tract, then introducing the deposition information into Monte Carlo simulation software to perform radiation emission simulation to obtain the measurement data of the radiation at the germanium crystal in the Monte Carlo model of the target detector, and finally completing the passive efficiency scale task of the detector according to the measurement data.
As shown in fig. 10, a second aspect of the present embodiment provides a virtual device for implementing the method for irradiating passive efficiency scales in a human body according to the first aspect, which includes a three-dimensional model reconstruction unit, a flow field distribution calculation unit, a particle motion calculation unit, a simulation model arrangement unit, a radiation emission simulation unit, and a scale task completion unit;
The three-dimensional model reconstruction unit is used for reconstructing and obtaining a respiratory tract three-dimensional model of a target human body, wherein the respiratory tract three-dimensional model comprises a front nose channel sub-model, a rear nose channel sub-model, a nose mouth sub-model, a throat sub-model and a main bronchus bifurcation structure sub-model;
The flow field distribution calculation unit is in communication connection with the three-dimensional model reconstruction unit and is used for calculating the flow field distribution condition of the fluid in the respiratory tract three-dimensional model under the target environment and when a human body breathes by using finite element analysis software, wherein the fluid in the respiratory tract is turbulent constant-temperature incompressible Newtonian fluid, and the environmental parameters of the target environment comprise gravity field parameters and temperature field parameters;
The particle motion calculation unit is in communication connection with the flow field distribution calculation unit and is used for adding simulation particles with physical properties of aerosol into the three-dimensional model of the respiratory tract, calculating the motion condition of the simulation particles under the flow field distribution condition by using the finite element analysis software, and obtaining the deposition information of the simulation particles on the inner wall of the three-dimensional model of the respiratory tract, wherein the deposition information comprises deposition positions or deposition positions and deposition quantity;
the simulation model arrangement unit is in communication connection with the three-dimensional model reconstruction unit and is used for importing the respiratory tract three-dimensional model into Monte Carlo simulation software, adding a human tissue three-dimensional model for simulating radiation scattering and energy attenuation in a soft tissue structure around the respiratory tract three-dimensional model, and placing a Monte Carlo model of a target detector around the human tissue three-dimensional model;
The radiation emission simulation unit is respectively in communication connection with the particle motion calculation unit and the simulation model arrangement unit and is used for importing the deposition information into the Monte Carlo simulation software to perform radiation emission simulation so as to obtain measurement data of radiation at a germanium crystal position in a Monte Carlo model of the target detector;
and the calibration task completion unit is in communication connection with the ray emission simulation unit and is used for completing the passive efficiency calibration task of the target detector according to the measurement data.
The working process, working details and technical effects of the foregoing device provided in the second aspect of the present embodiment may refer to the method for irradiating passive efficiency scales in a human body described in the first aspect, which are not described herein again.
As shown in fig. 11, a third aspect of the present embodiment provides a computer device for performing the method for irradiating passive efficiency calibration in a human body according to the first aspect, which includes a memory, a processor, and a transceiver that are sequentially communicatively connected, where the memory is configured to store a computer program, the transceiver is configured to send and receive a message, and the processor is configured to read the computer program, and perform the method for irradiating passive efficiency calibration in a human body according to the first aspect. By way of specific example, the Memory may include, but is not limited to, random-Access Memory (RAM), read-Only Memory (ROM), flash Memory (Flash Memory), first-in first-out Memory (First Input First Output, FIFO), and/or first-in last-out Memory (First Input Last Output, FILO), etc., and the processor may be, but is not limited to, a microprocessor of the STM32F105 family. In addition, the computer device may include, but is not limited to, a power module, a display screen, and other necessary components.
The working process, working details and technical effects of the foregoing computer device provided in the third aspect of the present embodiment may refer to the method for irradiating passive efficiency scales in a human body described in the first aspect, which are not described herein again.
A fourth aspect of the present embodiment provides a computer readable storage medium storing instructions comprising the method of irradiating passive efficiency scales in a human body according to the first aspect, i.e. instructions stored on the computer readable storage medium, which when run on a computer, perform the method of irradiating passive efficiency scales in a human body according to the first aspect. The computer readable storage medium refers to a carrier for storing data, and may include, but is not limited to, a floppy disk, an optical disk, a hard disk, a flash Memory, and/or a Memory Stick (Memory Stick), where the computer may be a general purpose computer, a special purpose computer, a computer network, or other programmable devices.
The working process, working details and technical effects of the foregoing computer readable storage medium provided in the fourth aspect of the present embodiment may refer to the method for irradiating passive efficiency scales in a human body as described in the first aspect, which are not described herein.
A fifth aspect of the present embodiment provides a computer program product comprising a computer program or instructions which, when executed by a computer, implement the method of irradiating passive efficiency scales in a human body as described in the first aspect. Wherein the computer may be a general purpose computer, a special purpose computer, a computer network, or other programmable apparatus.
Finally, it should be noted that the above description is only of the preferred embodiments of the present invention and is not intended to limit the scope of the present invention. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present invention should be included in the protection scope of the present invention.
Claims (5)
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202410563575.8A CN118380155B (en) | 2024-05-08 | 2024-05-08 | Method, device, equipment, medium and product for calibrating passive efficiency of human body irradiation |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202410563575.8A CN118380155B (en) | 2024-05-08 | 2024-05-08 | Method, device, equipment, medium and product for calibrating passive efficiency of human body irradiation |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN118380155A CN118380155A (en) | 2024-07-23 |
| CN118380155B true CN118380155B (en) | 2025-02-14 |
Family
ID=91901564
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202410563575.8A Active CN118380155B (en) | 2024-05-08 | 2024-05-08 | Method, device, equipment, medium and product for calibrating passive efficiency of human body irradiation |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN118380155B (en) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN119400347B (en) * | 2024-12-31 | 2025-04-18 | 三亚中心医院(海南省第三人民医院、三亚中心医院医疗集团总院) | System and method for monitoring use of atomized medicine |
| CN119920479A (en) * | 2025-04-02 | 2025-05-02 | 中国人民解放军军事科学院军事医学研究院 | A full-scale human respiratory rhythmic breathing three-dimensional geometric model and its construction method |
| CN119920478A (en) * | 2025-04-02 | 2025-05-02 | 中国人民解放军军事科学院军事医学研究院 | A numerical model for predicting the deposition distribution of heterogeneous particles in the respiratory tract of animals and its construction method |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114707416A (en) * | 2022-04-18 | 2022-07-05 | 成都理工大学 | A kind of human body radiation dose detection method, device, system and computer equipment |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS5999276A (en) * | 1982-11-27 | 1984-06-07 | Hitachi Ltd | Radioactivity measurement method and device |
| WO2023220604A2 (en) * | 2022-05-11 | 2023-11-16 | University Of Iowa Research Foundation | Individualized whole-lung deposition model |
| CN117970414A (en) * | 2023-12-21 | 2024-05-03 | 中国辐射防护研究院 | Nondestructive gamma measurement and hot spot reconstruction method adopting multiple NaI detectors |
-
2024
- 2024-05-08 CN CN202410563575.8A patent/CN118380155B/en active Active
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114707416A (en) * | 2022-04-18 | 2022-07-05 | 成都理工大学 | A kind of human body radiation dose detection method, device, system and computer equipment |
Non-Patent Citations (1)
| Title |
|---|
| 人体上呼吸道气流结构与气溶胶沉积数值模拟;赵秀国;徐新喜;孙栋;谭树林;牛福;刘亚军;任旭东;;系统仿真学报;20120808(第08期);第1582-1586页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN118380155A (en) | 2024-07-23 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN118380155B (en) | Method, device, equipment, medium and product for calibrating passive efficiency of human body irradiation | |
| Schulz et al. | Particle deposition in the respiratory tract | |
| Xi et al. | Parametric study on mouth–throat geometrical factors on deposition of orally inhaled aerosols | |
| Bailey | The new ICRP model for the respiratory tract | |
| Harrison et al. | Radiation doses and risks from internal emitters | |
| CN114707416B (en) | Method, device and system for detecting irradiation dose in human body and computer equipment | |
| Dubsky et al. | Imaging regional lung function: A critical tool for developing inhaled antimicrobial therapies | |
| Yu et al. | Assessment of environmental radon hazard using human respiratory tract models | |
| Sakoda et al. | Lung dosimetry of inhaled radon progeny in mice | |
| Fleming et al. | A technique for simulating radionuclide images from the aerosol deposition pattern in the airway tree | |
| Schroeter et al. | A computer model of lung morphology to analyze SPECT images | |
| Balásházy et al. | Non-linear relationship of cell hit and transformation probabilities in a low dose of inhaledradon progenies | |
| Verbanck et al. | Inspiratory and expiratory aerosol deposition in the upper airway | |
| CN112255663B (en) | Method and system for measuring alpha potential exposure of living radon daughter | |
| Zhang et al. | An analysis of dependency of counting efficiency on worker anatomy for in vivo measurements: whole-body counting | |
| Parker | The dilemma of lung dosimetry | |
| Dong et al. | A Hybrid Sourceless Efficiency Calibration Method for Internal Irradiation Dose Calculation | |
| Karami et al. | A method for assessing subject-specific counting efficiency of whole-body monitoring systems for radioiodine measurements | |
| Ishikawa et al. | Deposition and clearance for radon progeny in the human respiratory tract | |
| Rabi et al. | CFD modeling of radon progeny transport and deposition in the human respiratory tract | |
| Takahashi et al. | Construction of a Computational Tracheobronchial Model to Simulate Radon Progeny Deposition with Computational Fluid Dynamics | |
| DONG et al. | New Method for Sourceless Efficiency Calibration for Internal Irradiation Dose Calculation Based on Coupled Cfpd and Mc Simulations | |
| Delache et al. | Positron emission tomography-based comparison of methods for exposing macaques to respiratory pathogens | |
| CN119960002A (en) | Method and device for determining virtual scale of lung counter and storage medium | |
| Fernández-Tena et al. | Performance of different materials and their coating for use in in-vitro tests to simulate particle deposition in 3D printed lung models |
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 |