[go: up one dir, main page]

WO2016006633A1 - Dispositif de diagnostic, système de diagnostic, procédé de simulation et programme - Google Patents

Dispositif de diagnostic, système de diagnostic, procédé de simulation et programme Download PDF

Info

Publication number
WO2016006633A1
WO2016006633A1 PCT/JP2015/069666 JP2015069666W WO2016006633A1 WO 2016006633 A1 WO2016006633 A1 WO 2016006633A1 JP 2015069666 W JP2015069666 W JP 2015069666W WO 2016006633 A1 WO2016006633 A1 WO 2016006633A1
Authority
WO
WIPO (PCT)
Prior art keywords
upper airway
model
fluid
dimensional shape
unit
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.)
Ceased
Application number
PCT/JP2015/069666
Other languages
English (en)
Japanese (ja)
Inventor
智憲 岩崎
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Kagoshima University NUC
Original Assignee
Kagoshima University NUC
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Kagoshima University NUC filed Critical Kagoshima University NUC
Priority to JP2016532957A priority Critical patent/JP6579472B2/ja
Publication of WO2016006633A1 publication Critical patent/WO2016006633A1/fr
Anticipated expiration legal-status Critical
Ceased legal-status Critical Current

Links

Images

Classifications

    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/08—Measuring devices for evaluating the respiratory organs

Definitions

  • the present invention relates to a diagnostic apparatus, a diagnostic system, a simulation method, and a program.
  • Obstructive Sleep Apnea Syndrome is a respiratory disease that causes various respiratory effects by obstructing the upper respiratory tract during sleep and causing respiratory problems.
  • Obstructive sleep apnea syndrome (hereinafter simply referred to as sleep apnea syndrome) is a cause of serious traffic accidents because it does not provide sufficient sleep at night and causes excessive sleepiness during the day. Attention has been paid.
  • Non-Patent Document 1 discloses a method for diagnosing sleep apnea syndrome by constructing a fluid model of an upper airway on a computer based on a three-dimensional CT image and performing fluid analysis using the fluid model.
  • This method based on the pressure distribution in the upper airway and the air velocity distribution, it is possible to narrow down the ventilation trouble site to some extent.
  • Non-Patent Document 1 is a few effective methods for narrowing down the cause site of sleep apnea syndrome, but due to the following disadvantages, effective treatment results are not always obtained. .
  • the simulation model of the upper airway is modeled as a rigid body that does not deform. For this reason, fluid analysis in consideration of elastic deformation of the upper airway is difficult. It is also difficult to reproduce the upper airway obstruction.
  • the main cause of sleep apnea syndrome is narrowing of the upper respiratory tract due to tongue depression due to gravity in the supine position.
  • the tongue which is a tissue around the upper airway, is not modeled. For this reason, the narrowing of the upper airway due to the drop of the tongue due to gravity cannot be reproduced by simulation. Therefore, with this technique, it becomes difficult to specify that the cause of sleep apnea syndrome is narrowing of the upper respiratory tract due to tongue depression due to gravity in the supine position.
  • the present invention has been made to solve the above-described problems, and is a diagnostic device, a diagnostic system, a simulation method, and a program that can obtain better treatment results with respect to respiratory diseases more reliably.
  • the purpose is to provide.
  • a diagnostic apparatus provides: Using a numerical analysis technique, based on the three-dimensional image data inside the subject's maxillofacial portion, a three-dimensional shape model of the upper airway and tissues around the upper airway, and air fluid in the upper airway A model generation unit for generating a model; A physical property value specific to the three-dimensional shape model and the fluid model generated by the model generation unit is given, and the upper airway, the surrounding tissue, and the air fluid in the upper airway accompanying the subject's breathing A simulation execution unit that calculates information on the ventilation state of the upper airway that fluctuates according to respiration by performing a simulation of structural coupling analysis; Is provided.
  • a model display unit that displays a three-dimensional shape model of the upper airway and the surrounding tissue generated by the model generation unit;
  • the simulation execution unit By performing a simulation of the fluid-structure interaction analysis using the three-dimensional model of the upper airway and the surrounding tissue whose three-dimensional shape has been changed by the changing unit, information on the ventilation state of the upper airway is obtained. calculate, It is good as well.
  • the three-dimensional model of the nasal cavity of the upper airway is a rigid body, and performs fluid analysis over one respiratory period to calculate information on the ventilation state of the nasal cavity, and at the boundary between the nasal cavity and the pharynx
  • a fluid analysis unit for calculating a cross-sectional average pressure The simulation execution unit Using the cross-sectional average pressure at the boundary between the nasal cavity and the pharynx over one respiratory cycle calculated by the fluid analysis unit as an initial condition, a model of the three-dimensional shape of the upper airway and the surrounding tissue excluding the nasal cavity is used. , By calculating a simulation of the fluid structure coupling analysis, to calculate information on the ventilation state of the upper airway, It is good as well.
  • the simulation execution unit In a state where the three-dimensional shape model of the surrounding tissue is deformed by gravity, the fluid-structure interaction analysis is performed, and information on the ventilation state of the upper airway is calculated. It is good as well.
  • the simulation execution unit As information on the ventilation state of the upper airway, Calculating information on the flow of air in the upper airway or information on deformation of the upper airway and a region around the upper airway; It is good as well.
  • simulation execution unit As information on the air flow in the upper airway, Calculating pressure distribution or flow velocity distribution in the upper airway, It is good as well.
  • the diagnostic system according to the second aspect of the present invention is: An imaging device for imaging three-dimensional image data inside the maxillofacial portion of the subject; A diagnostic apparatus according to the present invention using the three-dimensional image data imaged by the imaging apparatus; Is provided.
  • a simulation method provides: A computer uses a numerical analysis method to calculate a model of a three-dimensional shape of the upper airway and a tissue around the upper airway based on the three-dimensional image data inside the maxillofacial portion of the subject, A model generation process for generating a fluid model of the air; A computer gives specific physical property values to the three-dimensional shape model and the fluid model generated in the model generation step, and the upper airway, the surrounding tissue, and the upper airway associated with the subject's breathing A simulation execution step of calculating information on the ventilation state of the upper airway, which fluctuates according to breathing, by performing a simulation of fluid-structure interaction analysis of air; Is provided.
  • a program is: Computer Using a numerical analysis method, based on the three-dimensional image data inside the subject's maxillofacial region, a model of the three-dimensional shape of the upper airway and the tissue around the upper airway, and the air in the upper airway A model generation unit for generating a fluid model; A physical property value specific to the three-dimensional shape model and the fluid model generated by the model generation unit is given, and the upper airway, the surrounding tissue, and the air fluid in the upper airway accompanying the subject's breathing A simulation execution unit that calculates information on the ventilation state of the upper airway that fluctuates according to respiration by performing a simulation of structural coupling analysis; To function as.
  • simulation of fluid-structure interaction analysis is performed using a model of the three-dimensional shape of not only the upper airway but also the surrounding tissue.
  • this simulation it is possible to calculate not only the flow of air in the upper airway accompanying breathing but also the deformation of the upper airway and surrounding tissues. Since the state of the upper airway and the surrounding tissue during actual sleep can be reproduced by simulation, the cause of the respiratory system disease can be more accurately identified. As a result, good treatment results can be obtained for respiratory diseases more reliably.
  • FIG. 9A is a diagram showing an example of a cross section of the upper airway when not lying on the back
  • FIG. 9B is a diagram showing an example of a cross section of the upper airway when lying on the back.
  • An example of pressure distribution in the upper pharyngeal airway is shown by fluid-structure interaction analysis accompanying respiration.
  • An example of the displacement distribution of the upper pharyngeal airway by the fluid-structure interaction analysis accompanying respiration is shown.
  • It is a flowchart of a diagnostic process. It is a flowchart of a model production
  • Embodiment 1 FIG. First, a first embodiment of the present invention will be described.
  • the diagnostic system 100 is used to reproduce on the computer the ventilation state of the upper respiratory tract of the human body that is the subject in order to identify the cause of sleep apnea syndrome. It is done.
  • the diagnostic system 100 is based on the three-dimensional image data inside the subject's maxillofacial region acquired by an X-ray CT (Computer Tomography) device or the like in order to reproduce the ventilation state of the subject's upper airway.
  • a model of a three-dimensional shape of the upper airway in the maxillofacial portion and a tissue around the upper airway and a fluid model of air in the upper airway are generated.
  • the upper respiratory tract means from the nose to the nasal cavity, nasopharynx, pharynx, and larynx.
  • the diagnosis system 100 performs a simulation of fluid-structure interaction analysis using a three-dimensional shape model of the upper airway and the surrounding tissue of the upper airway and a fluid model of the air in the upper airway.
  • the upper airway, the three-dimensional shape model of the surrounding tissue, and the fluid model of the air in the upper airway are given unique physical properties, and the simulation of fluid-structure interaction analysis makes the upper airway fluctuating according to respiration, It is possible to clarify the interaction between the surrounding tissue deformation and the air flow due to breathing in the upper airway. In addition, it is possible to accurately reproduce both the air flow (pressure distribution, flow velocity distribution) in the upper airway and the deformation of the upper airway due to the deformation of surrounding tissue during breathing. By reproducing these, the cause site of sleep apnea syndrome can be identified more accurately.
  • the diagnostic system 100 includes an imaging device 1 and a computer 2.
  • the imaging device 1 and the computer 2 are connected via a communication network. Through this communication network, data can be transmitted and received between the imaging apparatus 1 and the computer 2.
  • the imaging apparatus 1 is an X-ray CT apparatus.
  • the imaging apparatus 1 captures a three-dimensional X-ray CT image of the subject's maxillofacial portion.
  • the imaging device 1 includes an X-ray tube and a detector. When the X-ray CT scan is executed, X-rays irradiated from the X-ray tube and transmitted through the subject are detected by the detector. The detection result of the detector is stored in the imaging device 1 as raw data.
  • the imaging apparatus 1 generates cross-sectional data (slice image data) of the maxillofacial portion of the subject by performing image reconstruction processing based on the stored raw data. Furthermore, the imaging device 1 generates three-dimensional image data inside the maxillofacial portion of the subject based on the slice image data. As described above, three-dimensional image data inside the maxillofacial portion of the subject is obtained. The obtained three-dimensional image data is transmitted to the computer 2 through the communication network as DICOM (Digital Imaging and Communications Communications in Medicine) data.
  • DICOM Digital Imaging and Communications Communications in Medicine
  • DICOM data refers to data generated based on the DICOM standard.
  • the DICOM standard is mainly used as a format for medical image data.
  • the DICOM data is composed of the three-dimensional image data and incidental information according to the DICOM standard.
  • the incidental information is attribute information of image data such as patient information, imaging condition information, image information, and display information, and is embedded as tag information in DICOM data.
  • the computer 2 generates a model of a three-dimensional shape of the upper airway in the maxillofacial portion of the subject and the tissue around the upper airway based on the three-dimensional X-ray CT image data included in the received DICOM data.
  • This model is generated based on the finite element method.
  • This model is generated by dividing an object into a finite number of elements bounded by nodes. That is, the generated upper airway and the tissue model around the upper airway are so-called mesh models configured by joining a plurality of elements in a mesh shape at nodes.
  • the element may be divided on the basis of voxel data constituting the medical image data.
  • the computer 2 gives unique physical property values to the mesh model of the three-dimensional shape of the tissue around the upper airway and its maxillofacial surface, and performs a fluid-structure interaction analysis simulation accompanying the breathing of the subject, Get information about the ventilation status of the upper airway.
  • FIG. 2 shows the hardware configuration of the computer 2 in FIG.
  • the computer 2 includes a control unit 31, a main storage unit 32, an external storage unit 33, an operation unit 34, a display unit 35, and a communication unit 36.
  • the main storage unit 32, the external storage unit 33, the operation unit 34, the display unit 35, and the communication unit 36 are all connected to the control unit 31 via the internal bus 30.
  • the control unit 31 includes a CPU (Central Processing Unit) and the like.
  • the CPU executes the program 39 stored in the external storage unit 33, thereby realizing each component of the computer 2 shown in FIG.
  • the main storage unit 32 is composed of RAM (Random-Access Memory) or the like.
  • the main storage unit 32 is loaded with a program 39 stored in the external storage unit 33.
  • the main storage unit 32 is used as a work area (temporary data storage area) of the control unit 31.
  • the external storage unit 33 includes a non-volatile memory such as a flash memory, a hard disk, a DVD-RAM (Digital Versatile Disc Random-Access Memory), a DVD-RW (Digital Versatile Disc ReWritable).
  • a program 39 to be executed by the control unit 31 is stored in advance. Further, the external storage unit 33 supplies data used when executing the program 39 to the control unit 31 in accordance with an instruction from the control unit 31, and stores the data supplied from the control unit 31.
  • the operation unit 34 includes a pointing device such as a keyboard and a mouse, and an interface device that connects the keyboard and the pointing device to the internal bus 30. Information regarding the content operated by the operator is input to the control unit 31 via the operation unit 34.
  • the display unit 35 includes a CRT (Cathode Ray Tube), an LCD (Liquid Crystal Display), an organic EL (ElectroLuminescence), or the like.
  • CTR Cathode Ray Tube
  • LCD Liquid Crystal Display
  • organic EL ElectroLuminescence
  • the communication unit 36 includes a serial interface or a parallel interface.
  • the communication unit 36 is connected to the imaging device 1 via the communication network and receives the 3D X-ray CT image data transmitted from the imaging device 1.
  • control unit 31 includes a control unit 31, a main storage unit 32, an external storage unit 33, an operation unit 34, a display unit 35, a communication unit 36, etc. as hardware resources. It demonstrates its function by being used as a.
  • the computer 2 having a hardware configuration as shown in FIG. 2 has a storage unit 10, a data acquisition unit 11, a model generation unit 12, a fluid analysis unit 13, A simulation execution unit 14 and an output unit 15 are provided.
  • the storage unit 10 corresponds to the external storage unit 33 in FIG. 2 in the hardware configuration shown in FIG.
  • the storage unit 10 stores various data.
  • One of the data stored by the storage unit 10 is DICOM data 21.
  • the data acquisition unit 11 corresponds to the control unit 31 and the communication unit 36 in the hardware configuration shown in FIG.
  • the data acquisition unit 11 receives the three-dimensional X-ray CT image data (DICOM data) of the maxillofacial portion transmitted from the imaging device 1.
  • the data acquisition unit 11 stores the received DICOM data 21 in the storage unit 10.
  • the model generation unit 12 corresponds to the control unit 31 in the hardware configuration illustrated in FIG.
  • the model generation unit 12 inputs DICOM data 21 stored in the storage unit 10, that is, three-dimensional image data inside the subject's maxillofacial portion. Based on the DICOM data 21, the model generation unit 12 generates a three-dimensional shape model of the upper airway and the tissue around the upper airway and a fluid model of air in the upper airway using the finite element method.
  • the three-dimensional shape model of the upper airway and the tissues around the upper airway is constructed based on the governing equations indicating the structure and shape deformation of each tissue.
  • Tissues around the upper respiratory tract include soft tissues such as bones, tongue, soft palate, and maxillofacial muscles.
  • the fluid model is constructed based on a fluid governing equation (for example, Bernoulli's equation, Naviestoke's equation, etc.) using information on the air flow such as the pressure and flow velocity of the fluid as variables.
  • a fluid governing equation for example, Bernoulli's equation, Naviestoke's equation, etc.
  • FIG. 3 shows an example of a three-dimensional image of the upper airway, tongue, and soft palate.
  • A1 is the upper airway
  • A2 is the soft palate
  • A3 is the tongue.
  • the tongue A3 and the soft palate A2 are tissues around the upper airway adjacent to the upper airway A1.
  • the tongue A3 is known as the main cause of sleep apnea syndrome. This is because the tongue A3 falls down and squeezes the upper airway when sleeping in the supine position.
  • FIG. 4 shows a three-dimensional image of bone in addition to the upper airway, soft palate, and tongue.
  • A4 is a bone.
  • the mandible and cervical vertebrae exist around the upper airway, tongue, and soft palate.
  • their three-dimensional shape models are generated as tissues around the upper airway.
  • FIG. 5 shows an example of a three-dimensional image of soft tissues such as muscles and fat around the jaw.
  • soft tissues such as muscles and fat around the jaw.
  • a three-dimensional shape model is generated as tissue around the upper airway. This is because enlargement of the soft tissue around the jaw greatly affects respiration.
  • FIG. 6 shows an example of a three-dimensional mesh model of soft tissue portions around the upper airway, tongue, soft palate, bone, and jaw.
  • nodes connecting elements are common on the boundary surface of each tissue.
  • the upper airway and its surrounding area are formed so as to be regarded as one three-dimensional shape model. In this way, it is possible to reproduce, by the three-dimensional shape model, how the adjacent tissues are deformed due to the deformation of each tissue.
  • the model generation unit 12 generates a three-dimensional model by dividing the upper airway into a nasal cavity portion and a portion below the pharynx.
  • the nasal cavity portion is generated as STL (Standard (Triangulated Language) data.
  • STL data is data representing three points of a triangular patch used to approximate a solid, and the object surface is represented as a polygon by triangular patches (facets).
  • FIG. 7 shows an example of a three-dimensional image of the upper airway.
  • the part above the line L is the nasal cavity part.
  • the portion above the line L is stored in the storage unit 10 as the nasal cavity portion STL data 22.
  • the model generation unit 12 generates a three-dimensional mesh model for the upper respiratory tract below the pharynx.
  • the generated three-dimensional mesh model of the upper airway below the pharynx and the three-dimensional mesh models of other tissues are stored in the storage unit 10 as the three-dimensional shape model data 23.
  • the fluid analysis unit 13 corresponds to the control unit 31 in the hardware configuration illustrated in FIG.
  • the fluid analysis unit 13 uses the nasal cavity portion STL data 22 to perform a fluid analysis over one respiratory cycle with the nasal cavity as a rigid body, and calculates information regarding the ventilation state of the nasal cavity. This is because the internal structure of the nasal cavity is complicated, there is almost no change in the shape of the airway, and fluid-structure interaction analysis is difficult. As information about the calculated ventilation state of the nasal cavity, for example, there is a pressure distribution in the nasal cavity.
  • the fluid analysis unit 13 calculates the cross-sectional average pressure P (t) at the boundary between the nasal cavity and the pharynx based on the pressure distribution in the nasal cavity.
  • FIG. 8 shows an example of the cross-sectional average pressure P (t) at the boundary between the nasal cavity and the pharynx over one respiratory cycle. As shown in FIG. 8, the cross-sectional average pressure P (t) repeats exhalation and inspiration with a period T. In expiration, air in the nasal cavity is positive with respect to the external pressure, and in inspiration, air in the nasal cavity is in negative pressure with respect to the external pressure.
  • the simulation execution unit 14 corresponds to the control unit 31 in the hardware configuration illustrated in FIG.
  • the simulation execution unit 14 gives specific physical property values to the three-dimensional shape model and the fluid model generated by the model generation unit 12, and the upper airway, surrounding tissues, and air fluid in the upper airway accompanying the subject's breathing
  • information on the ventilation state of the upper airway that varies with breathing is calculated.
  • the fluid structure coupled analysis simulation is a numerical analysis simulation for analyzing a phenomenon in which a solid structure is deformed by a force exerted by a fluid flow.
  • Sleep apnea syndrome is often the root cause of tissue surrounding the upper respiratory tract. For this reason, if the simulation of fluid-structure interaction analysis including the tissues around the upper airway is performed to reproduce the deformation of the tissues around the upper airway during sleep, the cause of sleep apnea syndrome Can be identified more accurately.
  • Specific physical property values of the three-dimensional shape model include Young's modulus, linear expansion coefficient, Poisson's ratio, shear elastic modulus, etc. of each tissue.
  • the viscosity coefficient of air, density, bulk modulus, Reynolds number, etc. is there.
  • the physical property values those stored in the storage unit 10 as the physical property value data 25 are used.
  • the physical property values those known as the physical property values of the tissue may be used, but those measured for each subject may be used.
  • the simulation execution unit 14 uses a calculated cross-sectional average pressure P (t) at the boundary between the nasal cavity and the pharynx over one respiratory cycle as an initial condition, and generates a three-dimensional mesh model of the upper airway excluding the nasal cavity and surrounding tissues. It is used to simulate the fluid-structure interaction analysis of the upper respiratory tract and surrounding tissues associated with the subject's breathing.
  • the simulation execution unit 14 performs a simulation of fluid-structure interaction analysis of the upper airway and the surrounding tissue that accompanies the subject's breathing in a state where the three-dimensional model of the surrounding tissue is deformed by gravity.
  • FIG. 9A shows an example of a cross section of the upper airway when not lying on the back
  • FIG. 9B shows an example of a cross section of the upper airway when lying on the back. ing.
  • the upper airway is narrower when lying on its back. This is because the tongue part is sunk by gravity and compresses the upper airway.
  • the simulation execution unit 14 performs the simulation of the fluid structure coupled analysis in consideration of the gravity. Therefore, the ventilation of the upper airway is accurately reproduced while the upper airway deforming according to the posture of the subject is accurately reproduced. The state can be calculated.
  • the simulation execution unit 14 calculates information regarding the air flow in the upper airway and the deformation state of the upper airway and the three-dimensional shape model of the tissue around the upper airway.
  • the simulation execution unit 14 calculates the pressure distribution in the upper airway, the flow velocity distribution, or the displacement distribution of the three-dimensional shape model of the tissue around the upper airway and the upper airway as information on the air flow in the upper airway. These pieces of information are stored in the storage unit 10 as simulation result data 26.
  • the output unit 15 corresponds to the control unit 31 and the display unit 35 in the hardware configuration of FIG.
  • the output unit 15 displays the pressure distribution in the upper airway, the flow velocity distribution, or the deformation state of the upper airway.
  • FIG. 10 shows an example of the pressure distribution in the pharyngeal airway part by the fluid structure coupled analysis accompanying respiration.
  • FIG. 11 shows an example of the displacement distribution of the upper pharyngeal airway portion by the fluid structure coupled analysis accompanying respiration.
  • the higher the pressure is the darker the color is displayed.
  • the larger the displacement is the darker the color is displayed.
  • the pressure on the upper airway is higher.
  • the output unit 15 can also display and output the upper airway deformed by the fluid-structure interaction analysis and the three-dimensional shape of the tissue around the upper airway.
  • step S1 an imaging process is performed in the imaging apparatus 1 (step S1).
  • step S1 three-dimensional image data inside the subject's maxillofacial portion is obtained.
  • the three-dimensional image data is sent to the computer 2 and stored in the storage unit 10 of the computer 2 as DICOM data 21 by the data acquisition unit 11.
  • the computer 2 performs a model generation process (step S2).
  • the model generation unit 12 performs the three-dimensional analysis of the upper airway and surrounding tissues based on the three-dimensional image data (DICOM data 21) inside the subject's maxillofacial portion.
  • a shape model and a fluid model in the upper airway are generated, and the data of the generated fluid model is stored in the storage unit 10 as the three-dimensional shape model data 23 (step S21).
  • the model generation unit 12 generates the three-dimensional model of the nasal cavity part and the three-dimensional model of the part other than the nasal cavity for the upper airway separately.
  • the model generation unit 12 generates STL data of the nasal cavity portion and stores it as the nasal cavity portion STL data 22 in the storage unit 10 (step S22).
  • the fluid analysis unit 13 performs a fluid analysis of the nasal cavity part based on the nasal cavity part STL data 22, calculates the pressure distribution in the nasal cavity, and stores it in the storage unit 10 as the intranasal pressure distribution data 24 (Ste S23).
  • the fluid analysis is performed on the assumption that the nasal cavity portion is rigid and does not deform.
  • the fluid analysis unit 13 calculates an average cross-sectional pressure P (t) between the nasal cavity and the pharynx based on the pressure distribution obtained by the fluid analysis, and also includes it in the intranasal pressure distribution data 24 to store it. 10 (step S24).
  • the computer 2 performs a simulation execution process (step S3).
  • the simulation execution unit 14 refers to the physical property value data 25 stored in the storage unit 10, and sets the physical property values of the upper airway and surrounding tissues to the respective three-dimensional shape models. (Step S31).
  • the physical property value of the tongue is set for the three-dimensional shape model of the tongue
  • the physical property value of the soft palate is set for the three-dimensional shape model of the soft palate
  • the physical property value of the bone is set for the three-dimensional shape model of the bone.
  • the physical property value of the soft tissue is set for the three-dimensional shape model of the soft tissue around the jaw.
  • the simulation execution unit 14 sets the direction of gravity (step S32). This setting is performed by an operation input from the operation unit 34. For example, when diagnosing the ventilation state of the upper airway in the supine position, the occipital side of the subject is set as the lower side.
  • the simulation executing unit 14 performs a fluid structure coupled analysis simulation (step S33).
  • the simulation is performed in a state where the three-dimensional shape model of the upper airway and surrounding tissue and the fluid model of the upper airway are integrated.
  • the average cross-sectional pressure P (t) between the nasal cavity and pharynx in the upper airway is used as an initial condition.
  • the fluid analysis that analyzes the air flow in the upper airway and the structural analysis that analyzes the deformation of the upper airway and the surrounding tissue take into account the mutual effects of the three-dimensional shape model and the fluid model.
  • the pressure distribution, flow velocity distribution, and upper airway displacement distribution in the upper airway in one breath are calculated.
  • FIG. 15 shows the flow of the simulation process of the fluid structure interaction analysis performed in step S33.
  • the simulation execution unit 14 performs a gravity calculation so that the three-dimensional shape model of the surrounding tissue is deformed (step S41).
  • the three-dimensional shape model including the upper airway and surrounding tissues is in the state shown in FIG.
  • FIG. 16 schematically shows the upper airway in the supine position and the tissue around the upper airway.
  • the upper airway A1 is vertically sandwiched between the tongue A3 and the lower bone A4. Furthermore, due to the drop of the tongue A3 due to the gravity G, the upper airway A1 is deformed and becomes narrower than the awake state.
  • the simulation execution unit 14 first performs a fluid analysis of the fluid model A5 in the upper airway A1 to obtain the pressure P in the upper airway (step S42).
  • the simulation execution part 14 calculates
  • the simulation execution unit 14 performs structural analysis of the upper airway A1 and surrounding tissue using the load condition as a boundary condition, and the displacement of the inner wall surface of the upper airway in a state where the gravity G and the pressure P are balanced (for example, , [Delta] x in FIG. 16 and the like are obtained (step S43).
  • the simulation execution unit 14 updates the three-dimensional shape of the model of the upper airway and the tissue around the upper airway based on the displacement of the inner wall surface obtained by the structural analysis.
  • the simulation execution unit 14 determines whether or not the termination condition is satisfied (step S44).
  • the end condition may be that the load condition and the displacement of the wall surface converge within an allowable range.
  • step S44 the simulation executing unit 14 performs a fluid analysis based on the inner wall surface displaced by the structural analysis (step S42). In this way, steps S41 ⁇ S42 are repeated. By repeating this, the load condition (pressure P) and the displacement ( ⁇ x) of the inner wall surface converge to constant values.
  • a negative pressure P is generated in the upper airway during breathing and in the inspiratory phase. This negative pressure P further narrows the upper airway and displaces the boundary surface.
  • the final load that converged within the allowable range of each boundary element was repeatedly performed by fluid analysis using the fluid model in the upper airway and structural analysis using the model of the upper airway and surrounding tissue.
  • a set of conditions (pressure P) and displacement ( ⁇ x) of the inner wall surface is obtained as pressure distribution and displacement distribution in the upper airway.
  • the simulation result data 26 stored in the storage unit 10 includes the pressure distribution in the upper airway, the displacement distribution of the upper airway, and the deformed three-dimensional shape of the tissue around the upper airway when the termination condition is satisfied. It is data.
  • the output unit 15 displays and outputs the simulation result (step S34).
  • the display unit 35 displays a pressure distribution and a flow velocity distribution in the upper airway, a deformed three-dimensional model of the tissue around the upper airway, and the like.
  • the simulation result displayed on the output unit 15 it is possible to know the pressure distribution during breathing in the upper airway of the subject, the deformation of the upper airway, the deformation of surrounding tissues, and the like. Based on these pieces of information, it is possible to detect which part of the upper airway is constricted. If the upper airway is narrowed, it is easy to identify the cause of sleep apnea syndrome.
  • the output unit 15 can display the deformation of the three-dimensional shape of the tissue model around the upper airway along with the pressure distribution in the upper airway or the displacement distribution of the upper airway. In this way, it becomes possible to more accurately identify which part of the surrounding tissue is causing the stenosis in the upper airway.
  • the nasal cavity is considered as a cause of sleep apnea syndrome.
  • sleep apnea syndrome Such as nasal congestion and nasal catarrh.
  • the site to be treated is the nose.
  • adenoid is considered as a cause of sleep apnea syndrome.
  • excision of adenoids is an appropriate treatment method.
  • the tonsils of the palate are considered as the cause of sleep apnea syndrome.
  • the site to be treated is the soft palate.
  • the jaw When the pressure is extremely high at B4, the jaw is considered to be a cause of sleep apnea syndrome.
  • the site to be treated is the lower jaw. Correction or weight loss of the upper teeth and lower teeth is the main treatment method.
  • the cause of sleep apnea syndrome is not always one.
  • a plurality of parts may be the cause part. If the simulation of the fluid-structure interaction analysis according to this embodiment is performed, it is easy to find out that a plurality of parts cause sleep apnea syndrome.
  • simulation of fluid-structure interaction analysis is performed using a model of a three-dimensional shape of not only the upper airway but also the surrounding tissue.
  • Embodiment 2 FIG. Next, a second embodiment of the present invention will be described.
  • the diagnostic system according to Embodiment 2 of the present invention is used for planning a treatment plan for sleep apnea syndrome, which is a respiratory disease.
  • the output unit 15 displays and outputs the three-dimensional shape model of the upper airway generated by the model generation unit 12 and the tissue around the upper airway.
  • the computer 2 is different from the first embodiment in that the changing unit 16 is provided.
  • the output unit 15 displays and outputs the upper airway and the tissue around the copper based on the nasal cavity portion STL data 22 and the three-dimensional shape model data 23 stored in the storage unit 10 in response to the operation input.
  • the user operates the operation unit 34 (for example, by operating the mouse) while looking at the 3D model of the upper airway displayed on the output unit 15 and the tissue around the upper airway (for example, by operating the mouse). Specify an area.
  • the operation input input from the operation unit 34 is input to the control unit 31, and the control unit 31 performs processing according to the operation input. This processing corresponds to the changing unit 16 in FIG.
  • an area D shown in FIG. 19 is designated as a specific area by an operation input from the operation unit 34.
  • the change unit 16 changes the three-dimensional shape mesh model of the upper airway and the tissue around the upper airway so as to cut out the designated region D according to the operation input of the operation unit 34.
  • the output unit 15 displays the three-dimensional shape mesh model of the changed upper airway and the tissue around the upper airway.
  • the changing unit 16 changes the three-dimensional mesh model of the upper airway and the tissue around the upper airway according to the operation input.
  • FIG. 20 shows the upper airway after excision of the designated region D and the tissue around the upper airway.
  • the changed three-dimensional model is stored in the storage unit 10 as the three-dimensional shape model data 23.
  • the STL data of the changed nasal cavity part is stored in the storage unit 10 as the nasal cavity part STL data 22.
  • the simulation executing unit 14 uses the three-dimensional shape mesh model of the upper airway and the tissue around the upper airway edited by the changing unit 16 and stored in the storage unit 10.
  • the fluid structure coupled analysis of the upper airway associated with breathing and the tissue surrounding the upper airway is simulated to calculate information on the air flow in the upper airway and information on the deformation of the upper airway.
  • the calculated simulation result is stored in the storage unit 10 as simulation result data 26 and is displayed and output by the output unit 15. With reference to this information, it is possible to check the ventilation state of the upper airway after the designated area D is deleted.
  • the fluid-structure interaction analysis is simulated using the model of the three-dimensional shape of the upper airway changed by the changing unit 16 and the surrounding tissue.
  • this simulation it is possible to predict information about the ventilation state of the upper airway after treatment. If it is possible to predict the upper respiratory tract and the surrounding tissue during sleep after surgery by simulation, the optimal amount of tissue resection can be obtained before surgery, and an appropriate treatment plan can be established. Can stand. As a result, good treatment results for respiratory diseases can be obtained more reliably.
  • the diagnosis system 100 is used to change the three-dimensional shape of each tissue suspected as a causal site and perform simulation to analyze how the ventilation state of the upper airway changes. By doing so, it is also possible to accurately specify a complex cause site.
  • the fluid-structure interaction analysis was performed using a weakly coupled method (time difference method).
  • the strongly coupled method integrated solution solves the governing equations of the fluid and the structure exactly at the same time.
  • fluid-structure interaction analysis may be performed.
  • the simulation method of the fluid-structure interaction analysis is not limited to the above-described method, and various methods can be applied.
  • the finite element method is used to generate the three-dimensional shape model of the upper airway and the tissue around the upper airway and the fluid model of the air in the upper airway. Absent.
  • a numerical analysis method such as a finite difference method, a boundary element method, a finite volume method, etc.
  • An air fluid model may be generated.
  • the imaging apparatus 1 and the X-ray CT apparatus are used, but the present invention is not limited to this.
  • An MRI (Magnetic Resonance Imaging) apparatus or an ultrasonic diagnostic apparatus may be used as the imaging apparatus 1.
  • one 3D image data is generated from a plurality of 3D image data obtained from an X-ray CT apparatus, an MRI apparatus, and an ultrasonic diagnostic apparatus, and a model of a 3D shape of each tissue is generated from the generated image data. You may make it produce
  • the present invention is not limited to this.
  • An evaluation system can be used for diagnosis and treatment as long as it is a respiratory disease and is related to the shape of the upper respiratory tract.
  • the present invention can also be used to identify the cause of symptoms such as hypertension.
  • the central part that performs processing of the computer 2 composed of the control unit 31, the main storage unit 32, the external storage unit 33, the operation unit 34, the display unit 35, the communication unit 36, the internal bus 30, and the like is a dedicated system Regardless, it can be realized using a normal computer system.
  • a computer program for executing the above operation is stored and distributed in a computer-readable recording medium (flexible disk, CD-ROM, DVD-ROM, etc.), and the computer program is installed in the computer.
  • the computer 2 that executes the above-described processing may be configured.
  • the computer 2 may be configured by storing the computer program in a storage device included in a server device on a communication network such as the Internet and downloading it by a normal computer system.
  • the functions of the computer 2 are realized by sharing an OS (operating system) and an application program or by cooperation between the OS and the application program, only the application program portion may be stored in a recording medium or a storage device. .
  • the computer program may be posted on a bulletin board (BBS, “Bulletin” Board System) on the communication network, and the computer program may be distributed via the network.
  • BBS bulletin board
  • the computer program may be started and executed in the same manner as other application programs under the control of the OS, so that the above-described processing may be executed.
  • 1 imaging device 2 computer, 10 storage unit, 11 data acquisition unit, 12 model generation unit, 13 fluid analysis unit, 14 simulation execution unit, 15 output unit, 16 change unit, 21 DICOM data, 22 nasal cavity part STL data, 23 3D shape model data, 24 nasal pressure distribution data, 25 physical property value data, 26 simulation result data, 30 internal bus, 31 control unit, 32 main storage unit, 33 external storage unit, 34 operation unit, 35 display unit, 36 Communication department, 39 program, 100 diagnostic system.

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Biomedical Technology (AREA)
  • Molecular Biology (AREA)
  • Pulmonology (AREA)
  • Biophysics (AREA)
  • Pathology (AREA)
  • Engineering & Computer Science (AREA)
  • Physiology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Medical Informatics (AREA)
  • Physics & Mathematics (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • Apparatus For Radiation Diagnosis (AREA)
  • Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)

Abstract

Selon la présente invention, une unité de génération de modèle (12) utilise un procédé d'analyse numérique pour générer, sur la base de données d'image en trois dimensions de l'intérieur d'une partie maxillo-faciale d'un cobaye, un modèle de forme tridimensionnelle des voies respiratoires supérieures et des tissus entourant les voies respiratoires supérieures, et un modèle de fluide de l'air dans les voies respiratoires supérieures. Une unité d'exécution de simulation (14) affecte des valeurs de propriété physique unique au modèle de forme tridimensionnelle et au modèle de fluide générés par l'unité de génération de modèle (12), et calcule des informations relatives à l'état du passage d'air dans les voies respiratoires supérieures, qui varie en fonction de la respiration, par réalisation d'une simulation d'une analyse couplée de structure de fluide des voies respiratoires supérieures, des tissus environnants, et de l'air à l'intérieur des voies respiratoires supérieures accompagnant la respiration du cobaye.
PCT/JP2015/069666 2014-07-10 2015-07-08 Dispositif de diagnostic, système de diagnostic, procédé de simulation et programme Ceased WO2016006633A1 (fr)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2016532957A JP6579472B2 (ja) 2014-07-10 2015-07-08 顎顔面部の流体構造連成解析を用いた気道通気状態解析システム

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2014142076 2014-07-10
JP2014-142076 2014-07-10

Publications (1)

Publication Number Publication Date
WO2016006633A1 true WO2016006633A1 (fr) 2016-01-14

Family

ID=55064262

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2015/069666 Ceased WO2016006633A1 (fr) 2014-07-10 2015-07-08 Dispositif de diagnostic, système de diagnostic, procédé de simulation et programme

Country Status (2)

Country Link
JP (1) JP6579472B2 (fr)
WO (1) WO2016006633A1 (fr)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2018159759A1 (fr) * 2017-03-01 2018-09-07 国立大学法人 鹿児島大学 Système d'étalonnage d'état de ventilation des voies respiratoires et système pour prédire la déformation des voies respiratoires pendant le sommeil
WO2018193955A1 (fr) * 2017-04-18 2018-10-25 国立大学法人 鹿児島大学 Système de test de fonction de déglutition utilisant une caméra 3d
JP6456461B1 (ja) * 2017-11-17 2019-01-23 株式会社アルム 鼻弁狭窄診断装置、および鼻弁狭窄診断システム
JP2019526291A (ja) * 2016-08-30 2019-09-19 太暘科技國際股▲ふん▼有限公司 睡眠時気道の距離調整方法、ソフトウェア、システム及び記録媒体
CN110916667A (zh) * 2019-11-20 2020-03-27 四川大学华西医院 一种偏瘫患者仰卧位姿势智能评估床
WO2020232502A1 (fr) * 2019-05-20 2020-11-26 Commonwealth Scientific And Industrial Research Organisation Prédiction et intervention en lien avec une apnée obstructive du sommeil

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2005537068A (ja) * 2002-08-30 2005-12-08 ユニバーシティー オブ フロリダ 呼吸仕事量を予測するための方法および装置
JP2014104361A (ja) * 2012-11-27 2014-06-09 Ge Medical Systems Global Technology Co Llc スカウト画像内でローカライザを自動的に決定するための方法およびシステム

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2005537068A (ja) * 2002-08-30 2005-12-08 ユニバーシティー オブ フロリダ 呼吸仕事量を予測するための方法および装置
JP2014104361A (ja) * 2012-11-27 2014-06-09 Ge Medical Systems Global Technology Co Llc スカウト画像内でローカライザを自動的に決定するための方法およびシステム

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
TOMONORI IWASAKI ET AL.: "Clinical Application of Upper Airway Ventilation Simulation Using Fluid-Structure-Coupled Analysis in Children with Obstructive Sleep Apnea Syndrome", JOURNAL OF THE JAPANESE ASSOCIATION FOR DENTAL SCIENCE, vol. 33, March 2014 (2014-03-01), pages 39 - 43 *

Cited By (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2019526291A (ja) * 2016-08-30 2019-09-19 太暘科技國際股▲ふん▼有限公司 睡眠時気道の距離調整方法、ソフトウェア、システム及び記録媒体
JPWO2018159759A1 (ja) * 2017-03-01 2019-12-26 国立大学法人 鹿児島大学 気道通気状態キャリブレーションシステム及び睡眠時の気道変形予測システム
WO2018159759A1 (fr) * 2017-03-01 2018-09-07 国立大学法人 鹿児島大学 Système d'étalonnage d'état de ventilation des voies respiratoires et système pour prédire la déformation des voies respiratoires pendant le sommeil
JP7075131B2 (ja) 2017-03-01 2022-05-25 国立大学法人 鹿児島大学 気道通気状態キャリブレーションシステム及び睡眠時の気道変形予測システム
US11432742B2 (en) 2017-03-01 2022-09-06 Kagoshima University Airway ventilation state calibration system and system for predicting airway deformation during sleep
WO2018193955A1 (fr) * 2017-04-18 2018-10-25 国立大学法人 鹿児島大学 Système de test de fonction de déglutition utilisant une caméra 3d
JPWO2018193955A1 (ja) * 2017-04-18 2020-02-27 国立大学法人 鹿児島大学 ３ｄカメラを用いた摂食嚥下機能検査システム
JP6456461B1 (ja) * 2017-11-17 2019-01-23 株式会社アルム 鼻弁狭窄診断装置、および鼻弁狭窄診断システム
US12211209B2 (en) 2019-05-20 2025-01-28 Commonwealth Scientific And Industrial Reseach Organisation Prediction and intervention of obstructive sleep apnoea
WO2020232502A1 (fr) * 2019-05-20 2020-11-26 Commonwealth Scientific And Industrial Research Organisation Prédiction et intervention en lien avec une apnée obstructive du sommeil
AU2020277895B2 (en) * 2019-05-20 2025-03-27 Commonwealth Scientific And Industrial Research Organisation Prediction and intervention of obstructive sleep apnoea
CN110916667A (zh) * 2019-11-20 2020-03-27 四川大学华西医院 一种偏瘫患者仰卧位姿势智能评估床
CN110916667B (zh) * 2019-11-20 2022-08-19 四川大学华西医院 一种偏瘫患者仰卧位姿势智能评估床

Also Published As

Publication number Publication date
JP6579472B2 (ja) 2019-09-25
JPWO2016006633A1 (ja) 2017-05-25

Similar Documents

Publication Publication Date Title
Bates et al. Assessing the relationship between movement and airflow in the upper airway using computational fluid dynamics with motion determined from magnetic resonance imaging
JP6579472B2 (ja) 顎顔面部の流体構造連成解析を用いた気道通気状態解析システム
Zhao et al. Computational fluid dynamics for the assessment of upper airway response to oral appliance treatment in obstructive sleep apnea
De Backer et al. Functional imaging using computational fluid dynamics to predict treatment success of mandibular advancement devices in sleep-disordered breathing
Haskell et al. Effects of mandibular advancement device (MAD) on airway dimensions assessed with cone-beam computed tomography
Chiang et al. Three-dimensional airway evaluation in 387 subjects from one university orthodontic clinic using cone beam computed tomography
JP2016526466A (ja) 呼吸パラメータの決定
Zheng et al. Computational fluid dynamics simulation of the upper airway response to large incisor retraction in adult class I bimaxillary protrusion patients
De Backer et al. Novel imaging techniques using computer methods for the evaluation of the upper airway in patients with sleep-disordered breathing: a comprehensive review
TWI413511B (zh) The method of forming the visualization map of the pressure distribution of the flow field of the sleep apnea and the computer
Chang et al. Fluid structure interaction simulations of the upper airway in obstructive sleep apnea patients before and after maxillomandibular advancement surgery
Van Holsbeke et al. Anatomical and functional changes in the upper airways of sleep apnea patients due to mandibular repositioning: a large scale study
JP7075131B2 (ja) 気道通気状態キャリブレーションシステム及び睡眠時の気道変形予測システム
Stratemann et al. Three-dimensional analysis of the airway with cone-beam computed tomography
Chen et al. Analyses of aerodynamic characteristics of the oropharynx applying CBCT: obstructive sleep apnea patients versus control subjects
US20210196217A1 (en) Dynamic 3-d anatomical mapping and visualization
Moreddu et al. Computational fluid dynamics in the assessment of nasal obstruction in children
Van Gaver et al. Functional imaging improves patient selection for mandibular advancement device treatment outcome in sleep-disordered breathing: a prospective study
Lucey et al. Measurement, reconstruction, and flow-field computation of the human pharynx with application to sleep apnea
Dultra et al. Pharyngeal airspace of asthmatic individuals and those suffering from obstructive sleep apnea syndrome: Study by CBCT
Chen et al. Comparison of anatomic and aerodynamic characteristics of the upper airway among edentulous mild, moderate, and severe obstructive sleep apnea in older adults
Fletcher et al. The effect of geniglossal advancement on airway flow using a computational flow dynamics model
Patel et al. Modeling congenital nasal pyriform aperture stenosis using computational fluid dynamics
Farazi et al. A 3D dynamic biomechanical swallowing model for training and diagnosis of dysphagia
Barbour et al. Analysis of upper airway flow dynamics in Robin sequence infants using 4-D computed tomography and computational fluid dynamics

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 15818462

Country of ref document: EP

Kind code of ref document: A1

ENP Entry into the national phase

Ref document number: 2016532957

Country of ref document: JP

Kind code of ref document: A

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 15818462

Country of ref document: EP

Kind code of ref document: A1