[go: up one dir, main page]

CN119818084A - SMV blood distribution calculation method and postoperative hepatic encephalopathy risk assessment method - Google Patents

SMV blood distribution calculation method and postoperative hepatic encephalopathy risk assessment method Download PDF

Info

Publication number
CN119818084A
CN119818084A CN202510310483.3A CN202510310483A CN119818084A CN 119818084 A CN119818084 A CN 119818084A CN 202510310483 A CN202510310483 A CN 202510310483A CN 119818084 A CN119818084 A CN 119818084A
Authority
CN
China
Prior art keywords
smv
blood
phase
stent
portal vein
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN202510310483.3A
Other languages
Chinese (zh)
Inventor
熊柱翔
熊霞
彭涵
郑庭辉
杨朝武
宋小燕
李雯
余春林
邱莫寒
张增荣
樊华
胡陈明
汪江先
杜龙环
陈家磊
夏波
杨礼
朱师良
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Sichuan Daheng Poultry Breeding Co ltd
Sichuan Muzhou Agriculture Co ltd
Sichuan Animal Science Academy
Original Assignee
Sichuan Daheng Poultry Breeding Co ltd
Sichuan Muzhou Agriculture Co ltd
Sichuan Animal Science Academy
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 Sichuan Daheng Poultry Breeding Co ltd, Sichuan Muzhou Agriculture Co ltd, Sichuan Animal Science Academy filed Critical Sichuan Daheng Poultry Breeding Co ltd
Priority to CN202510310483.3A priority Critical patent/CN119818084A/en
Publication of CN119818084A publication Critical patent/CN119818084A/en
Pending legal-status Critical Current

Links

Landscapes

  • Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)

Abstract

本发明公开了一种SMV血液分布计算方法及术后肝性脑病风险评估方法,涉及医学技术领域,包括步骤:S1、获取患者的CTA影像数据和超声数据;S2、对门静脉系统解剖和TIPS支架进行三维重建;S3、完成数值模拟;S4、定量分析SMV血液分布特征;S5、通过对TIPS手术前SMV血液分布计算,得到SMV血液分布特征,分析SMV血液流向门静脉左支和门静脉右支的比例,选择植入支架的位置,减少术后SMV血液进入TIPS支架的比例,评估TIPS手术后患者肝性脑病的风险。本发明可以分析SMV血液在门静脉系统的流动规律,并且能够获知任意位置处SMV血液的体积占比(即百分比);验证了“血氨假说”,提高了肝性脑病的预估准确度。

The present invention discloses a method for calculating SMV blood distribution and a method for assessing the risk of postoperative hepatic encephalopathy, which relates to the field of medical technology and includes the following steps: S1, obtaining CTA image data and ultrasound data of a patient; S2, performing three-dimensional reconstruction of the portal vein system anatomy and TIPS stent; S3, completing numerical simulation; S4, quantitatively analyzing the SMV blood distribution characteristics; S5, obtaining the SMV blood distribution characteristics by calculating the SMV blood distribution before TIPS surgery, analyzing the proportion of SMV blood flowing to the left branch of the portal vein and the right branch of the portal vein, selecting the position of implanting the stent, reducing the proportion of SMV blood entering the TIPS stent after surgery, and assessing the risk of hepatic encephalopathy in patients after TIPS surgery. The present invention can analyze the flow law of SMV blood in the portal vein system, and can obtain the volume proportion (i.e. percentage) of SMV blood at any position; the "blood ammonia hypothesis" is verified, and the accuracy of hepatic encephalopathy estimation is improved.

Description

SMV blood distribution calculation method and postoperative hepatic encephalopathy risk assessment method
Technical Field
The invention relates to the technical field of medicine, in particular to an SMV blood distribution calculation method and a postoperative hepatic encephalopathy risk assessment method.
Background
Portal hypertension, i.e., abnormally elevated portal system pressure, is a major consequence of cirrhosis. It can lead to serious complications such as esophageal gastric fundus varices rupture bleeding and refractory hydrothorax and ascites. Transjugular Intrahepatic Portosystemic Shunt (TIPS) is a minimally invasive procedure, conducted under angiographic guidance. The purpose of this procedure is to implant a stent graft at the left, right or trunk branch of the portal vein, creating a shunt path from the portal vein to the hepatic vein. In this way, a portion of portal blood may enter the hepatic vein directly and then flow into the inferior vena cava, thereby reducing portal pressure.
TIPS (transjugular intrahepatic portosystemic shunt) is used as a minimally invasive and accurate operation, can obviously relieve portal hypertension, and is the most widely applied treatment method in clinic at present. However, hepatic Encephalopathy (HE) occurs in up to 35% -55% of patients within one year after TIPS stent implantation. Hepatic encephalopathy, also known as hepatic coma, is a syndrome characterized by dysfunction of the central nervous system. Its main symptoms include a decline in consciousness level, abnormal behavior, and coma. If not treated effectively in time, the illness state may be worsened, so that a plurality of systems such as nerves, digestion, blood, endocrine of a patient are dysfunctional, and the life quality of the patient is seriously affected. Thus, for patients with cirrhosis with portal hypertension, prevention of hepatic encephalopathy after TIPS surgery is a key factor in maintaining their long-term quality of life.
Although the exact mechanism by which TIPS post-operative hepatic encephalopathy occurs is not yet fully understood, it is widely believed to be closely related to post-operative changes in the hemodynamic system of the portal vein. Portal venous blood flow is formed by the convergence of superior mesenteric veins, which are rich in hormones such as insulin and glucagon that promote hepatocyte proliferation, and splenic veins, which carry nutrients and harmful substances such as blood ammonia that are absorbed from the intestinal tract. Under normal physiological conditions, these blood flows into the liver through the left and right portal veins, and ammonia is bioconverted by synthesizing urea and glutamine. However, TIPS surgery results in partial portal venous blood flow bypassing the liver directly into the hepatic vein, and particularly superior mesenteric venous blood rich in blood ammonia flows into the systemic circulation without liver treatment, eventually reaching the brain, disrupting the blood brain barrier, which is believed to be the primary cause of hepatic encephalopathy. Therefore, deep research on hemodynamic changes of the portal vein system before and after TIPS operation, especially on the distribution characteristics of venous blood on mesentery, has a vital meaning for preventing hepatic encephalopathy after TIPS operation.
Currently, ultrasound and portal angiography techniques are mainly used clinically to analyze the blood flow pattern of the portal system and the blood distribution characteristics of the superior mesenteric vein. The method provides important support for optimizing the TIPS technical scheme, remarkably improves the success rate and safety of the TIPS operation, and is beneficial to effectively reducing the risk of hepatic encephalopathy of patients. However, ultrasound examination may be limited by factors such as patient size and intestinal gas interference, and portal angiography presents certain risks and potential complications as an invasive examination. More importantly, both ultrasound and portal angiography can only perform qualitative analysis on the flow distribution of SMV and SV in the portal vein, and cannot accurately quantify the specific distribution of SMV in the left and right portal veins. Therefore, in order to more accurately evaluate hemodynamic changes before and after TIPS surgery, particularly quantitative distribution characteristics of SMV in the portal venous system, a more accurate and safe evaluation means is urgently needed.
Disclosure of Invention
The invention aims to overcome the defects of the prior art and provides an SMV blood distribution calculation method and a postoperative hepatic encephalopathy risk assessment method.
The aim of the invention is realized by the following technical scheme:
In a first aspect, the invention discloses a method for calculating SMV blood distribution, comprising the following steps:
s1, acquiring CTA image data and ultrasonic data of a patient;
S2, performing three-dimensional reconstruction on the portal vein system anatomy and the TIPS bracket;
S3, completing numerical simulation;
S4, quantitatively analyzing SMV blood distribution characteristics.
Based on the first aspect, the step S2 specifically comprises reconstructing a three-dimensional anatomical model of a patient before and after a TIPS operation based on CTA data of the patient, wherein the three-dimensional anatomical model comprises a portal vein trunk, a portal vein left branch, a portal vein right branch, a superior mesenteric vein and a splenic vein, the three-dimensional anatomical model of the patient after the TIPS operation further comprises a TIPS bracket, all inlets and outlets on the model are prolonged by 10 times of the diameter of the three-dimensional anatomical model, mesh division is carried out through ANSYS ICEM software, unstructured mesh types are selected, the mesh size is set to be 1/15 of the diameter of the portal vein, the mesh near the wall of a blood vessel is encrypted, 6 layers of boundary layers are set, the total thickness is 1/10 of the size of an outlet of the model, and the growth rate is 1.2.
Based on the first aspect, step S3 specifically includes passing through the N-S equationAnd continuity equationNumerical simulation was performed in whichRepresenting the velocity vector of the fluid,Indicating the pressure of the fluid and,Indicating the density of the blood and,Is a vector differential operator, and is expressed by a formulaCalculating dynamic viscosity coefficient of bloodWhereinThe shear rate is indicated by the term "shear rate",,,,The inlet boundary condition of the numerical simulation is obtained by conversion according to clinical actual ultrasonic data, the outlet is set as a resistance boundary condition, and the numerical simulation is obtained by calculation of portal vein pressure and flow.
Based on the first aspect, step S3 further comprises setting convergence criteria for the numerical simulation, including residuals for both the mass equation and the momentum equation below 10 -5.
Based on the first aspect, step S3 further comprises dividing the results of the numerical simulation into pre-TIPS operation and post-TIPS operation;
Before TIPS operation, the SMV blood proportion of the left branch of the portal vein is R LPV=QLPV-SMV/QSMV,QLPV-SMV=QLPV*RLPV-SMV, wherein R LPV is the proportion of SMV blood in the left branch of the portal vein, Q LPV-SMV is the flow of SMV blood in the left branch of the portal vein, Q SMV is the total flow of SMV blood in the system of the portal vein, Q LPV is the total blood flow of the left branch of the portal vein, and R LPV-SMV is the proportion of SMV blood in the left branch of the portal vein to the total blood;
The SMV blood proportion of the right branch of the portal vein is R RPV=QRPV-SMV/QSMV,QRPV-SMV=QRPV*RRPV-SMV, wherein R RPV is the proportion of SMV blood in the right branch of the portal vein, Q RPV-SMV is the flow of SMV blood in the right branch of the portal vein, Q RPV is the total blood flow of the right branch of the portal vein, and R RPV-SMV is the proportion of SMV blood in the right branch of the portal vein to the total blood;
After TIPS operation, the SMV blood ratio in the stent is R Stent=QStent-SMV/QSMV,QStent-SMV=QStent*RStent-SMV, wherein R Stent is the SMV blood ratio in the TIPS stent, Q Stent-SMV is the SMV blood flow in the TIPS stent, Q Stent is the total blood flow of the TIPS stent, and R Stent-SMV is the SMV blood ratio in the TIPS stent.
Based on the first aspect, the step S4 specifically comprises incorporating a liquid-liquid two-phase flow module, analyzing quantitative distribution characteristics of SMV blood in portal vein, and combining the first stepVolume of phaseIs defined asWhereinIs the firstThe volume fraction of the phase is given by the equationFirst, theThe effective density of the phase isWhereinIs the firstPhysical Density of the phases (1)The continuity equation of the phase isNumber of integrated stereoCalculation of equation(s) for (i) and (ii) the next phase continuity equationVolume fraction of phases, whereinIs the firstVelocity of the phase (c)The momentum equation for the phase is: Wherein Is the acceleration of gravity and,Is the firstThe stress strain tensor of a phase is defined asWhereinIs the firstThe shear viscosity of the phase is determined by the viscosity of the phase,Is the firstThe bulk viscosity of the phase, T represents the transpose,Is an identity matrix of the unit cell,Is externally applied toThe volumetric force of the phase is such that,Is externally applied toThe internal lift of the phase is such that,Is applied to the firstThe virtual mass force of the phase(s),Is the firstDrag force of phase, through formulaCalculating a drag force, whereinRepresenting the discrete phase(s),The continuous phase is represented by the term "continuous phase",In continuous phaseIs used to control the drag force of the (c) in the (c),Is in a discrete phaseIs used to control the drag force of the (c) in the (c),Is in a discrete phaseIs used for the volume fraction of (a),In continuous phaseIs used for the volume fraction of (a),Is in a discrete phaseIs used for the physical density of the (c) a,Is in a discrete phaseIs used for the speed of the (c) in the (c),In continuous phaseIs used for the speed of the (c) in the (c),The calculation method is that the particle surface relaxation parameter time is as followsWhereinRepresenting discrete phasesThe diameter of the liquid droplet is such that,Representing discrete phasesShear viscosity, drag force function of (2)The calculation method of (1) is thatWhereinIn order for the coefficient of drag to be the same,Drag coefficient for relative Reynolds numberThe calculation method of (1) is thatRelative Reynolds numberThe calculation method of (1) is thatWhereinIn continuous phaseIs used for the physical density of the (c) a,In continuous phaseIs used for the shear viscosity of the rubber composition.
The invention also discloses a postoperative hepatic encephalopathy risk assessment method based on SMV blood distribution calculation, which comprises the following steps of S5, assessing the risk of hepatic encephalopathy of a patient after TIPS operation by calculating the SMV blood distribution before the TIPS operation, wherein the SMV blood distribution calculation specifically adopts the SMV blood distribution calculation method.
Based on the second aspect, step S5 comprises calculating SMV blood distribution characteristics of the SMV blood distribution before TIPS operation, analyzing the proportion of SMV blood flowing to the left branch of the portal vein and the right branch of the portal vein, thereby selecting the position of the implanted stent, reducing the proportion of SMV blood entering the TIPS stent after operation, and evaluating the risk of hepatic encephalopathy of the patient after TIPS operation.
The beneficial effects of the invention are as follows:
1) The invention can quantitatively analyze the distribution characteristics of SMV blood, and preliminarily verify the 'blood ammonia hypothesis', thereby improving the estimation accuracy of hepatic encephalopathy.
2) The invention can analyze the flow rule of SMV blood in the portal vein system, can acquire the volume ratio (i.e. percentage) of SMV blood at any position, and selects a proper position to implant a bracket by analyzing the ratio of SMV blood before TIPS operation to the left branch of the portal vein and the right branch of the portal vein, thereby reducing the ratio of SMV blood entering the TIPS bracket after operation and evaluating the risk of hepatic encephalopathy of patients after TIPS operation.
Drawings
FIG. 1 is a schematic flow chart of an embodiment of the present invention;
FIG. 2 is a graph showing the results of SMV blood ratio in TIPS postoperative stents in hepatic encephalopathy group and non-hepatic encephalopathy group according to the embodiment of the present invention;
fig. 3 is a schematic diagram of a three-dimensional model reconstruction process according to an embodiment of the present invention.
Detailed Description
The technical solutions of the present invention will be clearly and completely described below with reference to the embodiments, and it is apparent that the described embodiments are only some embodiments of the present invention, but not all embodiments. All other embodiments, which can be made by a person skilled in the art without any inventive effort, are intended to be within the scope of the present invention, based on the embodiments of the present invention.
The invention discloses an SMV blood distribution calculating method, which comprises the following steps:
s1, acquiring CTA image data and ultrasonic data of a patient;
S2, performing three-dimensional reconstruction on the portal vein system anatomy and the TIPS bracket;
S3, completing numerical simulation;
S4, quantitatively analyzing SMV blood distribution characteristics.
Specifically, step S2 comprises reconstructing a three-dimensional anatomical model of a patient before and after a TIPS operation based on CTA data of the patient, wherein the three-dimensional anatomical model comprises a portal vein trunk, a portal vein left branch, a portal vein right branch, an superior mesenteric vein and a splenic vein, the three-dimensional anatomical model after the TIPS operation further comprises a TIPS bracket, wherein for simulation, all the inlets and outlets are prolonged by 10 times of the diameter of the three-dimensional anatomical model to ensure that viscous fluid at the inlets can fully develop and reduce the influence of a backflow area at the outlets, meshing is performed by ANSYS ICEM software, unstructured mesh types are selected, the mesh size is set to be 1/15 of the diameter of the portal vein, the mesh near the wall of a blood vessel is subjected to encryption treatment in consideration of boundary layer effect, 6 layers of boundary layers are set, the total thickness is about 1/10 of the size of the outlet of the model, the growth rate is 1.2, a three-dimensional model reconstruction process is schematically shown in FIG. 3, wherein PV represents the portal vein, SMV represents the superior mesenteric vein, SV represents the left gastric vein, LPV represents the left portal vein, RPV represents the right portal vein, Q SMV represents the bracket, Q SMV represents the total blood flow resistance of the system LGR R represents the left portal vein, and the resistance value of the left portal vein is the total flow resistance of the system.
Specifically, step S3 includes providing that the pulsation characteristic of the blood flow is not significant in the venous system, and that the blood flow in the portal vein is relatively stable during the ultrasound examination and collection of the portal vein flow, assuming that the blood is a continuously homogeneous incompressible non-Newtonian fluid with a dynamic viscosity coefficientBy N-S equation, reference is made to the Carreau-Yasuda model in a non-Newtonian fluidAnd continuity equationNumerical simulation was performed in whichRepresenting the velocity vector of the fluid,Indicating the pressure of the fluid and,Represents the density of blood, is set to 1050 kg/m3,Is a vector differential operator, and is expressed by a formulaCalculating dynamic viscosity coefficient of bloodDescription of dynamic viscosity coefficients of blood Using Carreau modelWith shear rateΛ=3.313 s, n= 0.3568,,The inlet boundary condition of the numerical simulation is obtained by conversion according to clinical actual ultrasonic data, the outlet is set as a resistance boundary condition, and the numerical simulation is obtained by calculation of portal vein pressure and flow.
Specifically, step S3 also includes that the numerical simulation of the hemodynamics is performed in computational fluid dynamics (Computational Fluid Dynamics, CFD) software Ansys Fluent. All calculations use a software default split solver (SEGREGATED SOLVER), i.e., a Pressure-Based solver (Pressure-Based). The convergence criterion is set to be below 10-5 for both mass and momentum equation residuals.
Specifically, step S3 further includes dividing the results of the numerical simulation into a pre-TIPS operation and a post-TIPS operation, respectively, wherein before the TIPS operation, the proportion of SMV blood in the left branch of the portal vein is R LPV=QLPV-SMV/QSMV,QLPV-SMV=QLPV*RLPV-SMV, wherein R LPV is the proportion of SMV blood in the left branch of the portal vein, Q LPV-SMV is the flow of SMV blood in the left branch of the portal vein, Q SMV is the total flow of SMV blood in the portal vein system (i.e., the total flow of SMV blood in the right branch of the portal vein and the left branch of the portal vein), Q LPV is the total flow of SMV blood in the left branch of the portal vein, and R LPV-SMV is the proportion of SMV blood in the left branch of the portal vein to the total blood. The SMV blood ratio of the right portal branch is R RPV=QRPV-SMV/QSMV,QRPV-SMV=QRPV*RRPV-SMV, wherein R RPV is the SMV blood ratio of the right portal branch, Q RPV-SMV is the SMV blood flow of the right portal branch, Q RPV is the total SMV blood flow of the right portal branch, and R RPV-SMV is the SMV blood ratio of the right portal branch to the total SMV blood. After TIPS operation, the SMV blood ratio in the stent is R Stent=QStent-SMV/QSMV,QStent-SMV=QStent*RStent-SMV, wherein R Stent is the SMV blood ratio in the TIPS stent, Q Stent-SMV is the SMV blood flow in the TIPS stent, Q Stent is the total blood flow of the TIPS stent, and R Stent-SMV is the SMV blood ratio in the TIPS stent.
Specifically, step S4 includes incorporating a liquid-liquid two-phase flow module for the purpose of exploring the quantitative distribution characteristics of SMV blood in the portal vein. The model uses an Eulerian method to calculate each phase respectively, and solves simultaneous equations of mass conservation and momentum conservation, so that a high-precision simulation result can be provided. In order to more effectively apply the euler two-fluid model to the liquid-liquid two-phase flow problem of the portal vein system, the computational model is moderately simplified. In the process, the continuity and momentum conservation equation of each phase is solved independently, all phases share the same pressure field, mass and energy exchange does not exist between the phases, the diameter of the liquid drop is fixed, and the like. These simplifications help reduce model complexity while maintaining sufficient accuracy to capture critical hydrodynamic features within the portal system.
Specifically, the volume fraction equation and the phase volume fraction concept are one of the features of the euler dual-fluid model, and under the framework of the model, each phase needs to ensure that the mass and momentum conservation law is satisfied, so that the accurate description of the behaviors of each phase in multiphase flow is ensured. Will be the firstVolume of phaseIs defined asWhereinIs the firstThe volume fraction of the phase is given by the equationFirst, theThe effective density of the phase isWhereinIs the firstPhysical density of the phases.
Specifically, the continuity equation, thThe continuity equation of the phase isNumber of integrated stereoCalculation of equation(s) for (i) and (ii) the next phase continuity equationVolume fraction of phases, whereinIs the firstThe speed of the phases.
Specifically, the momentum equation, thThe momentum equation for the phase is: Wherein Is the acceleration of gravity and,Is the firstThe stress strain tensor of a phase is defined asWhereinIs the firstThe shear viscosity of the phase is determined by the viscosity of the phase,Is the firstThe bulk viscosity of the phase, T represents the transpose,Is an identity matrix of the unit cell,Is externally applied toThe volumetric force of the phase is such that,Is externally applied toThe internal lift of the phase is such that,Is applied to the firstThe virtual mass force of the phase(s),Is the firstThe drag force of the phase, since the present invention employs a liquid-liquid two-phase flow module, the flow in the vessel does not need to take into account external volumetric forces, and the lift force is a very small amount relative to the drag force and therefore can be ignored. In the present invention, it is assumed that the blood flowing from the SMV is a first blood, and the blood flowing from the splenic vein is a second blood, wherein the first blood and the second blood are both identical in material property settings, the first blood is a first phase, and the second blood is a second phase. The material properties (density and viscosity) of the two phases are set to be the same, so the virtual mass force is also negligible. Finally only the effect of drag force has to be considered.
Specifically, by the formulaCalculating a drag force, whereinRepresenting the discrete phase(s),The continuous phase is represented by the term "continuous phase",In continuous phaseIs used to control the drag force of the (c) in the (c),Is in a discrete phaseIs used to control the drag force of the (c) in the (c),Is in a discrete phaseIs used for the volume fraction of (a),In continuous phaseIs used for the volume fraction of (a),Is in a discrete phaseIs used for the physical density of the (c) a,Is in a discrete phaseIs used for the speed of the (c) in the (c),In continuous phaseIs used for the speed of the (c) in the (c),The calculation method is that the particle surface relaxation parameter time is as followsWhereinRepresenting discrete phasesThe diameter of the liquid droplet is such that,Representing discrete phasesShear viscosity, drag force function of (2)There is a certain difference in the different models, but they all contain a relative Reynolds numberAdopts a Schiller model and a drag functionThe calculation method of (1) is thatWhereinIn order for the coefficient of drag to be the same,Drag coefficient for relative Reynolds numberThe calculation method of (1) is thatRelative Reynolds numberThe calculation method of (1) is thatWhereinIn continuous phaseIs used for the physical density of the (c) a,In continuous phaseIs used for the shear viscosity of the rubber composition.
The invention also discloses a postoperative hepatic encephalopathy risk assessment method based on SMV blood distribution calculation, wherein the flow chart of the steps is shown in figure 1, and the method comprises the following steps of S5, assessing the risk of hepatic encephalopathy of a patient after a TIPS operation through calculating the SMV blood distribution before the TIPS operation, wherein the SMV blood distribution calculation specifically adopts the SMV blood distribution calculation method. The method specifically comprises the steps of calculating SMV blood distribution characteristics before TIPS operation, analyzing the proportion of SMV blood flowing to left portal vein branch and right portal vein branch, selecting the position of an implanted stent, reducing the proportion of SMV blood entering the TIPS stent after operation, and evaluating the risk of hepatic encephalopathy of patients after TIPS operation.
By way of example, the invention allows the quantitative analysis of the distribution characteristics of SMV blood and the preliminary verification of the "blood ammonia hypothesis", which indicates that, after a large amount of blood ammonia enters the brain circulation, it is able to cross the blood brain barrier and perform a synthetic treatment in the brain, which consumes a large amount of energy, thus affecting other normal brain functions and eventually causing hepatic encephalopathy. Multiple animal experiments prove that the blood ammonia concentration in SMV blood is obviously higher than that in other blood vessels, so that excessive SMV blood directly flows into the systemic circulation without being detoxified by the liver, and the SMV blood is probably an important cause of high incidence of hepatic encephalopathy. The results of SMV blood ratios in TIPS postoperative stents in hepatic encephalopathy group (HE, 33 cases) and non-hepatic encephalopathy group (N-HE, 12 cases) are schematically shown in FIG. 2. It can be seen that the proportion of SMV blood in the scaffolds of patients with hepatic encephalopathy group after TIPS surgery was significantly higher than that of patients with non-hepatic encephalopathy group, and the difference was statistically significant (P < 0.001). This result suggests that a significant amount of SMV blood rich in blood ammonia bypasses the detoxification process of the liver directly into the systemic circulation, probably being an important mechanism leading to the occurrence of postoperative hepatic encephalopathy. If the SMV blood proportion directly entering the bracket can be reduced, more SMV blood flows into the liver to detoxify, and blood ammonia is eliminated, the occurrence of hepatic encephalopathy after TIPS operation can be reduced.
Illustratively, the invention can select a proper position to implant the bracket by analyzing the proportion of SMV blood before TIPS operation to flow to the left branch of the portal vein and the right branch of the portal vein, thereby reducing the proportion of SMV blood after operation to enter the TIPS bracket and evaluating the occurrence probability of hepatic encephalopathy of patients after TIPS operation. In patient a and patient B, the portal vein pressure gradient (PCG) remained consistent before and after TIPS surgery. However, the ratios of pre-operative SMV blood flow to the left and right portal veins were 74.14% and 25.86% (patient a), and 36.08% and 63.92% (patient B), respectively, showing significant differences. After implantation of a TIPS stent of 8mm diameter into the right portal vein in both patients, the SMV blood flow to the stent was observed to be 45.83% (patient a) and 80.54% (patient B), respectively. Notably, patient a did not develop any postoperative complications during the up to 989 days of follow-up period, in contrast to patient B who developed hepatic encephalopathy only 27 days post-operation. This suggests that for the case like patient B, if the TIPS stent is chosen to be placed on the left portal branch where SMV blood is less, it may be helpful to reduce the risk of hepatic encephalopathy. According to the simulation result prediction of the invention, if a patient adopts the method, the proportion of SMV blood flowing into the stent after operation can be reduced to 63.56 percent.
The foregoing is merely a preferred embodiment of the invention, and it is to be understood that the invention is not limited to the form disclosed herein but is not to be construed as excluding other embodiments, but is capable of numerous other combinations, modifications and environments and is capable of modifications within the scope of the inventive concept, either as taught or as a matter of routine skill or knowledge in the relevant art. And that modifications and variations which do not depart from the spirit and scope of the invention are intended to be within the scope of the appended claims.

Claims (8)

1.一种SMV血液分布计算方法,其特征在于,包括以下步骤:1. A method for calculating SMV blood distribution, comprising the following steps: S1、获取患者的CTA影像数据和超声数据;S1. Obtain the patient's CTA imaging data and ultrasound data; S2、对门静脉系统解剖和TIPS支架进行三维重建;S2, 3D reconstruction of the portal venous system anatomy and TIPS stent; S3、完成数值模拟;S3, complete the numerical simulation; S4、定量分析SMV血液分布特征。S4. Quantitative analysis of SMV blood distribution characteristics. 2.根据权利要求1所述的一种SMV血液分布计算方法,其特征在于,步骤S2具体包括:基于患者CTA数据重建TIPS手术前、后的三维解剖模型,包括门静脉主干、门静脉左支、门静脉右支、肠系膜上静脉和脾静脉,TIPS手术后的三维解剖模型还包括TIPS支架;模型上的所有出入口均按照其直径的10倍延长,通过Ansys ICEM软件进行网格划分,选择非结构化网格类型,并设定网格大小为门静脉直径的1/15;在血管壁面附近的网格进行加密处理,设置6层边界层,总厚度为模型出口尺寸的1/10,增长率为1.2。2. A SMV blood distribution calculation method according to claim 1, characterized in that step S2 specifically comprises: reconstructing a three-dimensional anatomical model before and after TIPS surgery based on the patient's CTA data, including the portal vein trunk, the left branch of the portal vein, the right branch of the portal vein, the superior mesenteric vein and the splenic vein, and the three-dimensional anatomical model after TIPS surgery also includes a TIPS stent; all entrances and exits on the model are extended by 10 times their diameters, meshed by Ansys ICEM software, selecting an unstructured mesh type, and setting the mesh size to 1/15 of the portal vein diameter; the mesh near the blood vessel wall is encrypted, and 6 boundary layers are set, with a total thickness of 1/10 of the model outlet size and a growth rate of 1.2. 3.根据权利要求2所述的一种SMV血液分布计算方法,其特征在于,步骤S3具体包括:通过N-S方程和连续性方程进行数值模拟,其中表示流体的速度矢量,表示流体的压强,表示血液的密度,为矢量微分运算符;通过公式计算血液的动力粘性系数,其中表示剪切率,λ=3.313s,n=0.3568,,数值模拟的入口边界条件根据临床实际超声数据变换得到,出口设置为阻力边界条件,通过门静脉压强和流量计算得到。3. The method for calculating SMV blood distribution according to claim 2, characterized in that step S3 specifically comprises: using the NS equation and the continuity equation A numerical simulation is performed, in which represents the velocity vector of the fluid, represents the pressure of the fluid, Represents the density of blood, is a vector differential operator; by the formula Calculate the dynamic viscosity of blood ,in represents the shear rate, λ=3.313s, n=0.3568, , ,The inlet boundary conditions of the numerical simulation are obtained according to the ,transformation of actual clinical ultrasound data, and the outlet is set as the ,resistance boundary conditions, which are calculated by the portal vein pressure and ,flow. 4.根据权利要求3所述的一种SMV血液分布计算方法,其特征在于,步骤S3还包括:设定数值模拟的收敛标准,包括质量方程和动量方程的残差均低于10-54. The SMV blood distribution calculation method according to claim 3, characterized in that step S3 further comprises: setting a convergence standard for numerical simulation, including that the residuals of the mass equation and the momentum equation are both lower than 10 -5 . 5.根据权利要求4所述的一种SMV血液分布计算方法,其特征在于,步骤S3还包括:将数值模拟的结果分为TIPS手术前和TIPS手术后;5. The SMV blood distribution calculation method according to claim 4, characterized in that step S3 further comprises: dividing the results of the numerical simulation into before TIPS surgery and after TIPS surgery; 在TIPS手术前,门静脉左支的SMV血液比例为:RLPV=QLPV-SMV/QSMV,QLPV-SMV=QLPV*RLPV-SMV,其中,RLPV是门静脉左支中SMV血液的比例,QLPV-SMV是门静脉左支中SMV血液的流量,QSMV是门静脉系统SMV血液的总流量,QLPV是门静脉左支的总血流量,RLPV-SMV是门静脉左支中SMV血液占总血液的比例;Before TIPS surgery, the proportion of SMV blood in the left branch of the portal vein was: R LPV =Q LPV-SMV /Q SMV , Q LPV-SMV =Q LPV *R LPV-SMV , where R LPV is the proportion of SMV blood in the left branch of the portal vein, Q LPV-SMV is the flow rate of SMV blood in the left branch of the portal vein, Q SMV is the total flow rate of SMV blood in the portal vein system, Q LPV is the total blood flow rate of the left branch of the portal vein, and R LPV-SMV is the proportion of SMV blood in the left branch of the portal vein to the total blood; 门静脉右支的SMV血液比例为:RRPV=QRPV-SMV/QSMV,QRPV-SMV=QRPV*RRPV-SMV,其中,RRPV是门静脉右支中SMV血液的比例,QRPV-SMV是门静脉右支中SMV血液的流量,QRPV是门静脉右支的总血流量,RRPV-SMV是门静脉右支中SMV血液占总血液的比例;The proportion of SMV blood in the right branch of the portal vein is: R RPV =Q RPV-SMV /Q SMV , Q RPV-SMV =Q RPV *R RPV-SMV , where R RPV is the proportion of SMV blood in the right branch of the portal vein, Q RPV-SMV is the flow rate of SMV blood in the right branch of the portal vein, Q RPV is the total blood flow rate of the right branch of the portal vein, and R RPV-SMV is the proportion of SMV blood in the right branch of the portal vein to the total blood; TIPS手术后,支架内的SMV血液比例为:RStent=QStent-SMV/QSMV,QStent-SMV=QStent*RStent-SMV,其中,RStent是TIPS支架中SMV血液的比例,QStent-SMV是TIPS支架中SMV血液的流量,QStent是TIPS支架的总血流量,RStent-SMV是TIPS支架中SMV血液占总血液的比例。After TIPS surgery, the ratio of SMV blood in the stent is: R Stent =Q Stent-SMV /Q SMV , Q Stent-SMV =Q Stent *R Stent-SMV , where R Stent is the ratio of SMV blood in the TIPS stent, Q Stent-SMV is the flow rate of SMV blood in the TIPS stent, Q Stent is the total blood flow rate of the TIPS stent, and R Stent-SMV is the ratio of SMV blood to total blood in the TIPS stent. 6.根据权利要求5所述的一种SMV血液分布计算方法,其特征在于,步骤S4具体包括:纳入了液-液两相流模块,分析SMV血液在门静脉内的定量分布特征;将第相的体积定义为,其中为第相的体积分数,其方程为,第相的有效密度为,其中是第相的物理密度;第相的连续性方程为,联立体积分数的方程和下一相的连续性方程计算第相的体积分数,其中为第相的速度;第相的动量方程为: ;其中是重力加速度,是第相的应力应变张量,定义为,其中是第相的剪切粘度,是第相的体积粘度,T表示转置,是单位矩阵,是外部作用在第相的体积力,是外部作用在第相的内部升力,是施加给第相的虚拟质量力,是第相的曳力,通过公式计算曳力,其中表示离散相,表示连续相,为连续相的曳力,为离散相的曳力,为离散相的体积分数,为连续相的体积分数,为离散相的物理密度,为离散相的速度,为连续相的速度,为颗粒表面弛豫参数时间,其计算方法为,其中表示离散相液滴的直径,表示离散相的剪切粘度,曳力函数的计算方法为,其中为曳力系数,为相对雷诺数,曳力系数的计算方法为,相对雷诺数的计算方法为,其中为连续相的物理密度,为连续相的剪切粘度。6. A method for calculating SMV blood distribution according to claim 5, characterized in that step S4 specifically comprises: incorporating a liquid-liquid two-phase flow module to analyze the quantitative distribution characteristics of SMV blood in the portal vein; Phase volume Defined as ,in For the The volume fraction of the phase is given by , The effective density of the phase is ,in It is The physical density of the phase; The continuity equation of the phase is , combined volume fraction The equation of the next phase and the continuity equation of the next phase are used to calculate the The volume fraction of the phase, For the Phase speed; The momentum equation for the phase is: ;in is the acceleration due to gravity, It is The stress-strain tensor of the phase is defined as ,in It is The shear viscosity of the phase, It is The bulk viscosity of the phase, T represents the transposition, is the identity matrix, It is the external effect The body force of the phase, It is the external effect The internal lift of the phase, is applied to the The virtual mass force of the phase, It is The drag force of the phase is given by the formula Calculate the drag force, where represents the discrete phase, represents the continuous phase, Continuous phase The drag force, For discrete phase The drag force, For discrete phase The volume fraction of Continuous phase The volume fraction of For discrete phase The physical density of For discrete phase speed, Continuous phase speed, is the particle surface relaxation parameter time, which is calculated as follows: ,in Represents discrete phase The diameter of the droplet, Represents discrete phase Shear viscosity, drag function The calculation method is ,in is the drag coefficient, is the relative Reynolds number, the drag coefficient The calculation method is , relative Reynolds number The calculation method is ,in Continuous phase The physical density of Continuous phase Shear viscosity. 7.一种基于SMV血液分布计算的术后肝性脑病风险评估方法,其特征在于,该方法包括步骤:S5、通过对TIPS手术前SMV血液分布计算,评估TIPS手术后患者肝性脑病的风险;所述SMV血液分布计算,具体采用如权利要求1-6任一一项所述的SMV血液分布计算方法。7. A method for assessing the risk of postoperative hepatic encephalopathy based on SMV blood distribution calculation, characterized in that the method comprises the steps of: S5, assessing the risk of hepatic encephalopathy in patients after TIPS surgery by calculating the SMV blood distribution before TIPS surgery; the SMV blood distribution calculation specifically adopts the SMV blood distribution calculation method as described in any one of claims 1-6. 8.根据权利要求7所述的一种术后肝性脑病风险评估方法,其特征在于,步骤S5包括:对TIPS手术前SMV血液分布计算得到SMV血液分布特征,分析SMV血液流向门静脉左支和门静脉右支的比例,从而选择植入支架的位置,减少术后SMV血液进入TIPS支架的比例,用于评估TIPS手术后患者肝性脑病的风险。8. A method for assessing the risk of postoperative hepatic encephalopathy according to claim 7, characterized in that step S5 comprises: calculating the SMV blood distribution characteristics before TIPS surgery, analyzing the proportion of SMV blood flowing to the left branch of the portal vein and the right branch of the portal vein, thereby selecting the location for implanting the stent and reducing the proportion of SMV blood entering the TIPS stent after surgery, which is used to assess the risk of hepatic encephalopathy in patients after TIPS surgery.
CN202510310483.3A 2025-03-17 2025-03-17 SMV blood distribution calculation method and postoperative hepatic encephalopathy risk assessment method Pending CN119818084A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202510310483.3A CN119818084A (en) 2025-03-17 2025-03-17 SMV blood distribution calculation method and postoperative hepatic encephalopathy risk assessment method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202510310483.3A CN119818084A (en) 2025-03-17 2025-03-17 SMV blood distribution calculation method and postoperative hepatic encephalopathy risk assessment method

Publications (1)

Publication Number Publication Date
CN119818084A true CN119818084A (en) 2025-04-15

Family

ID=95306856

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202510310483.3A Pending CN119818084A (en) 2025-03-17 2025-03-17 SMV blood distribution calculation method and postoperative hepatic encephalopathy risk assessment method

Country Status (1)

Country Link
CN (1) CN119818084A (en)

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20170319984A1 (en) * 2016-05-03 2017-11-09 Saudi Arabian Oil Company Processes for analysis and optimization of multiphase separators, particular in regards to simulated gravity separation of immiscible liquid dispersions
CN117079811A (en) * 2023-08-11 2023-11-17 华西精创医疗科技(成都)有限公司 TIPS (tungsten inert gas) stent shunt proportion calculation method and hepatic encephalopathy risk assessment method

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20170319984A1 (en) * 2016-05-03 2017-11-09 Saudi Arabian Oil Company Processes for analysis and optimization of multiphase separators, particular in regards to simulated gravity separation of immiscible liquid dispersions
CN117079811A (en) * 2023-08-11 2023-11-17 华西精创医疗科技(成都)有限公司 TIPS (tungsten inert gas) stent shunt proportion calculation method and hepatic encephalopathy risk assessment method

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
WENZHI GAO等: "Hemodynamic investigations on the portal hypertension and treatment of transjugular intrahepatic portosystemic shunt (TIPS) based on CFD simulation", 《JOURNAL OF BIOMECHANICS》, vol. 181, 22 January 2025 (2025-01-22), pages 112516 - 1 *
XIONG ZHUXIANG等: "A new computational fluid dynamics based noninvasive assessment of portacaval pressure gradient", 《JOUNAL OF BIOMECHAINCS》, vol. 167, 30 April 2024 (2024-04-30), pages 112086 - 1, XP087512976, DOI: 10.1016/j.jbiomech.2024.112086 *

Similar Documents

Publication Publication Date Title
Chen et al. A patient-specific study of type-B aortic dissection: evaluation of true-false lumen blood exchange
Sahbaee et al. The effect of contrast material on radiation dose at CT: part I. Incorporation of contrast material dynamics in anthropomorphic phantoms
WO2016182508A1 (en) Medical image processing methods and systems
US20180242917A1 (en) Systems and methods for determining image quality and models based on contrast injection protocols
US20190189286A1 (en) Systems and methods for performing computer-simulated evaluation of treatments on a target population
Modi et al. Optimal Application of Fractional Flow Reserve to Assess Serial Coronary Artery Disease: A 3D‐Printed Experimental Study With Clinical Validation
Wang et al. Dynamic targeting in cancer treatment
Mut et al. Effects of changing physiologic conditions on the in vivo quantification of hemodynamic variables in cerebral aneurysms treated with flow diverting devices
Hunter et al. Simulations of congenital septal defect closure and reactivity testing in patient-specific models of the pediatric pulmonary vasculature: a 3D numerical study with fluid-structure interaction
CN109758228A (en) Optimal stent placement method for transjugular portosystemic shunt surgery based on three-dimensional reconstruction
CN103823956B (en) Determine that contrast agent uses method, device and the imaging device of parameter
CN117637099A (en) Forecasting methods, devices, equipment and storage media
Becker et al. Non-invasive physiological assessment of coronary artery obstruction on coronary computed tomography angiography
CN119818084A (en) SMV blood distribution calculation method and postoperative hepatic encephalopathy risk assessment method
Gutierrez et al. Assessment of coronary artery aneurysms caused by Kawasaki disease using transluminal attenuation gradient analysis of computerized tomography angiograms
Sarkhosh et al. Prediction of time averaged wall shear stress distribution in coronary arteries’ bifurcation varying in morphological features via deep learning
Gognieva et al. One-dimensional mathematical model-based automated assessment of fractional flow reserve in a patient with silent myocardial ischemia
Javadi et al. Predicting hemodynamic performance of Fontan operation for Glenn physiology using computational fluid dynamics: ten patient-specific cases
Malek et al. Machine learning-based prediction of hemodynamic parameters in left coronary artery bifurcation: A CFD approach
Chivukula et al. Leveraging patient-specific simulated angiograms to characterize cerebral aneurysm hemodynamics using computational fluid dynamics
Aramburu et al. Patient-specific closed-loop model of the fontan circulation: Calibration and validation
Sun et al. Does lower limb exercise worsen renal artery hemodynamics in patients with abdominal aortic aneurysm?
Segalova et al. Computational modeling of shear-based hemolysis caused by renal obstruction
Qiu et al. Modelling the liver region as porous media to noninvasively measure portal vein pressure gradient (PPG) with numerical methods
Bartolo et al. Computational framework for the generation of one‐dimensional vascular models accounting for uncertainty in networks extracted from medical images

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