阅读说明：本技术 生物组织个性化血流灌注率评估方法及系统 (Method and system for evaluating personalized blood perfusion rate of biological tissue ) 是由 王逸飞 张康伟 张爱丽 徐学敏 孙建奇 于 2020-06-10 设计创作，主要内容包括：本申请涉及生物医学工程领域,公开了一种生物组织个性化血流灌注率评估方法及系统,基于周边血流对目标生物组织冷冻后融化过程的影响,通过边界对流换热原理建立计算局部组织个性化血流灌注率的算法,所述方法包括：获得所述目标生物组织经冷冻而形成的冰球的表面积；获得所述冰球表面的温度梯度及其随时间的变化根据所述冰球表面的温度梯度及其随时间的变化以及所述冰球的表面积,确定所述目标生物组织周围的个性化血流灌注率；输出所述目标生物组织周围的个性化血流灌注率。(The application relates to the field of biomedical engineering, and discloses a method and a system for evaluating personalized blood perfusion rate of biological tissues, wherein an algorithm for calculating the personalized blood perfusion rate of local tissues is established through a boundary convection heat transfer principle based on the influence of peripheral blood flow on the freezing and thawing process of target biological tissues, and the method comprises the following steps: obtaining the surface area of an ice ball formed by freezing the target biological tissue; obtaining the temperature gradient of the ice ball surface and the change of the temperature gradient with time According to the temperature gradient of the ice hockey surface and the change of the ice hockey surface along with time And the surface area of the ice ball, determining personalized blood perfusion around the target biological tissueRate; outputting the personalized perfusion rate of blood flow around the target biological tissue.)

1. A biological tissue personalized blood perfusion rate assessment method is characterized in that an algorithm for calculating a local tissue personalized blood perfusion rate is established through a boundary convection heat transfer principle based on the influence of peripheral blood flow on a target biological tissue after freezing and thawing process, and the method comprises the following steps:

obtaining the surface area of an ice ball formed by freezing the target biological tissue;

obtaining the temperature gradient of the ice ball surface and the change of the temperature gradient with time

According to the temperature gradient of the ice hockey surface and the change of the ice hockey surface along with timeAnd the surface area of the ice hockey ball, determining the personalized blood perfusion rate around the target biological tissue;

outputting the personalized perfusion rate of blood flow around the target biological tissue.

2. The method of claim 1, wherein the temperature gradient is based on the surface of the puckDegree and its variation with timeAnd the surface area of the ice hockey puck, the step of determining the personalized perfusion rate of blood flow around the target biological tissue,

in the formula, ω_bIs the personalized perfusion rate of blood flow around the target biological tissue, k is the thermal conductivity of the target biological tissue,is the temperature gradient of the surface of the puck and its variation with time, T_bf(t)Is the blood flow temperature, T, of the ice ball surface over time_(t)Is the temperature, p, of the ice ball surface over time_bIs the blood density of the target biological tissue, c_bIs the blood specific heat capacity of the target biological tissue, A is the surface area of the ice hockey, and Σ is the surface of the ice hockey, wherein T is_(t)And T_bf(t)Is preset.

3. The method of claim 1, further comprising: according to the temperature gradient of the ice hockey surface and the change of the ice hockey surface along with timeAnd the surface area of the ice hockey, determining the personalized blood convection coefficient around the target biological tissue; and

outputting the personalized blood convection coefficient around the target biological tissue.

4. A method according to claim 3, wherein said temperature gradient is dependent on the surface of said puck and its variation over timeAnd the surface area of the ice hockey, and determining the personalized blood convection coefficient around the target biological tissue, wherein the personalized blood convection coefficient around the target biological tissue is calculated by the following method:

wherein h is_bIs the personalized blood convection coefficient around the target biological tissue, k is the thermal conductivity of the target biological tissue,is the temperature gradient of the surface of the puck and its variation with time, T_bf(t)Is the blood flow temperature, T, of the ice ball surface over time_(t)Is the temperature of the surface of the ice hockey as a function of time, A is the surface area of the ice hockey, Σ is the surface of the ice hockey wherein T is the temperature of the ice hockey surface_(t)And T_bf(t)Is preset.

5. A biological tissue personalized blood perfusion rate evaluation system is characterized in that the system establishes an algorithm for calculating a local tissue personalized blood perfusion rate through a boundary convection heat transfer principle based on the influence of peripheral blood flow on a target biological tissue after freezing and thawing process, and the system comprises:

an ice hockey size acquisition unit for acquiring the surface area of an ice hockey formed by freezing the target biological tissue;

an ice hockey temperature field distribution acquisition unit for acquiring the temperature gradient of the ice hockey surface and the change of the temperature gradient with time

An individualized blood perfusion rate obtaining unit for obtaining the blood perfusion rate according to the surface of the ice hockeyTemperature gradient and its variation with timeAnd the surface area of the ice hockey ball, determining the personalized blood perfusion rate around the target biological tissue;

an output unit for outputting an individualized blood perfusion rate around the target biological tissue.

6. The system of claim 5, wherein the personalized perfusion rate obtaining unit is configured to obtain the personalized perfusion rate of blood flow by determining a temperature gradient of the surface of the ice hockey puck and a change of the temperature gradient with timeAnd the surface area of the ice hockey ball, determining the personalized blood perfusion rate around the target biological tissue:

in the formula, ω_bIs the personalized perfusion rate of blood flow around the target biological tissue, k is the thermal conductivity of the target biological tissue,is the temperature gradient of the surface of the puck and its variation with time, T_bf(t)Is the blood flow temperature, T, of the ice ball surface over time_(t)Is the temperature, p, of the ice ball surface over time_bIs the blood density of the target biological tissue, c_bIs the blood specific heat capacity of the target biological tissue, A is the surface area of the ice hockey, and Σ is the surface of the ice hockey, wherein T is_(t)And T_bf(t)Is preset.

7. The system of claim 5, further comprising:

a personalized blood convection coefficient obtaining unit for obtaining the temperature gradient of the ice hockey surface and the change of the temperature gradient with timeAnd the surface area of the ice hockey, determining the personalized blood convection coefficient around the target biological tissue; and

the output unit is also used for outputting the personalized blood convection coefficient around the target biological tissue.

8. The system of claim 7, wherein the personalized blood convection coefficient acquisition unit around the target biological tissue calculates the personalized blood convection coefficient around the target biological tissue by:

wherein h is_bIs the personalized blood convection coefficient around the target biological tissue, k is the thermal conductivity of the target biological tissue,is the temperature gradient of the surface of the puck and its variation with time, T_bf(t)Is the blood flow temperature, T, of the ice ball surface over time_(t)Is the temperature A of the ice ball surface over time is the surface area of the ice ball, Σ is the surface of the ice ball wherein said T_(t)And T_bf(t)Is preset.

9. The system of claim 5, wherein the output unit is a display screen.

10. A biological tissue personalized perfusion rate assessment device, comprising:

a memory for storing computer executable instructions; and the number of the first and second groups,

a processor for implementing the steps in the method of any one of claims 1 to 4 when executing the computer-executable instructions.

11. A computer-readable storage medium having stored thereon computer-executable instructions which, when executed by a processor, implement the steps in the method of any one of claims 1 to 4.

Technical Field

The application relates to the field of biomedical engineering, in particular to the technical field of data processing combining medical technology and computer technology.

Background

In recent years, physical ablation of target biological tissues such as tumors has been widely used in the treatment of benign and malignant tumors of the liver, lung, kidney, and the like as a minimally invasive ablation method. Common physical ablation methods include cryoablation and thermal ablation.

Specifically, the multi-modal ablation therapy combines cryoablation and heating ablation, thoroughly destroys tumor cells by using huge thermal stress generated by rapid temperature change in the tumor, releases tumor antigens, activates specific immune response of an organism to the tumor, and is suitable for various solid tumors. Based on various animal high-metastatic tumor models and a large number of experimental results, the effect is proved, and the survival rate is far higher than that of surgical excision. Clinical studies have also shown that it can prolong disease-free survival of patients.

In particular, in order to maximize the effectiveness of multi-modal ablation therapy, the boundary temperature of the tumor needs to be precisely controlled to form a transition zone of active tumor antigen exposure. However, it is difficult to control the boundary temperature during the actual procedure, especially during the heat ablation stage. This is mainly because the individualized blood flow condition of the patient can have an unpredictable effect on the temperature of the ablation region. The existing multi-modal ablation operation planning model does not consider personalized blood flow conditions, and has great difference with actual conditions.

Therefore, a simpler and more effective approach for personalized assessment of blood flow conditions in biological tissues is needed.

Disclosure of Invention

The objective of the present application is to provide a method and a system for evaluating a personalized blood perfusion rate of a biological tissue, which are simple and effective, and are helpful for precisely controlling the boundary temperature of a target biological tissue, such as a tumor tissue, during thermal physical ablation, so as to fully exert the therapeutic effect thereof, and also can be used for evaluating the physiological characteristics of the biological tissue.

The application discloses a biological tissue personalized blood perfusion rate evaluation method, which is based on the influence of peripheral blood flow on the process of thawing a target biological tissue after freezing and establishes an algorithm for calculating the local tissue personalized blood perfusion rate according to the boundary convection heat transfer principle, and the method comprises the following steps:

obtaining the surface area of an ice ball formed by freezing the target biological tissue;

obtaining the temperature gradient of the ice ball surface and the change of the temperature gradient with time

According to the temperature gradient of the ice hockey surface and the change of the ice hockey surface along with timeAnd the surface area of the ice hockey ball, determining the personalized blood perfusion rate around the target biological tissue;

outputting the personalized perfusion rate of blood flow around the target biological tissue.

Preferably, the temperature gradient according to the ice hockey surface and the change of the temperature gradient with timeAnd the surface area of the ice hockey puck, the step of determining the personalized perfusion rate of blood flow around the target biological tissue,

in the formula, ω_bIs the personalized perfusion rate of blood flow around the target biological tissue, k is the thermal conductivity of the target biological tissue,is the temperature gradient of the surface of the puck and its variation with time, T_bf(t)Is the blood flow temperature, T, of the ice ball surface over time_(t)Is the temperature, p, of the ice ball surface over time_bIs the blood density of the target biological tissue, c_bIs the blood specific heat capacity of the target biological tissue, A is the surface area of the ice hockey, and Σ is the surface of the ice hockey, wherein T is_(t)And T_bf(t)Is preset.

Preferably, the method further comprises: according to the temperature gradient of the ice hockey surface and the change of the ice hockey surface along with timeAnd the surface area of the ice hockey, determining the personalized blood convection coefficient around the target biological tissue; and

outputting the personalized blood convection coefficient around the target biological tissue.

Preferably, the temperature gradient according to the ice hockey surface and the change of the temperature gradient with timeAnd the surface area of the ice hockey, and determining the personalized blood convection coefficient around the target biological tissue, wherein the personalized blood convection coefficient around the target biological tissue is calculated by the following method:

wherein h is_bIs the personalized blood convection coefficient around the target biological tissue, k is the thermal conductivity of the target biological tissue,is the temperature gradient of the surface of the puck and its variation with time, T_bf(t)Is the blood flow temperature, T, of the ice ball surface over time_(t)Is that the surface of the ice hockey is always in contact withA is the surface area of the ice hockey, Σ is the surface of the ice hockey wherein T is the temperature of the ice hockey_(t)And T_bf(t)Is preset.

The application also discloses individualized blood flow perfusion rate evaluation system of biological tissue, based on the influence of peripheral blood flow to the process of thawing target biological tissue after freezing, establish the algorithm of calculating individualized blood flow perfusion rate of local tissue through boundary convection heat transfer principle, the system includes:

an ice hockey size acquisition unit for acquiring the surface area of an ice hockey formed by freezing the target biological tissue;

an ice hockey temperature field distribution acquisition unit for acquiring the temperature gradient of the ice hockey surface and the change of the temperature gradient with time

A personalized blood perfusion rate acquisition unit for acquiring the blood perfusion rate according to the temperature gradient of the ice ball surface and the change of the ice ball surface along with the timeAnd the surface area of the ice hockey ball, determining the personalized blood perfusion rate around the target biological tissue;

an output unit for outputting an individualized blood perfusion rate around the target biological tissue.

Preferably, the personalized blood perfusion rate obtaining unit is used for obtaining the personalized blood perfusion rate according to the temperature gradient of the ice hockey surface and the change of the ice hockey surface along with the timeAnd the surface area of the ice hockey ball, determining the personalized blood perfusion rate around the target biological tissue:

in the formula, ω_bIs the one around the target biological tissueA sexually transfused blood flow rate, k being the thermal conductivity of the target biological tissue,is the temperature gradient of the surface of the puck and its variation with time, T_bf(t)Is the blood flow temperature, T, of the ice ball surface over time_(t)Is the temperature, p, of the ice ball surface over time_bIs the blood density of the target biological tissue, c_bIs the blood specific heat capacity of the target biological tissue, A is the surface area of the ice hockey, and Σ is the surface of the ice hockey, wherein T is_(t)And T_bf(t)Is preset.

Preferably, the system further comprises:

a personalized blood convection coefficient obtaining unit for obtaining the temperature gradient of the ice hockey surface and the change of the temperature gradient with timeAnd the surface area of the ice hockey, determining the personalized blood convection coefficient around the target biological tissue; and

the output unit is also used for outputting the personalized blood convection coefficient around the target biological tissue.

Preferably, the personalized blood convection coefficient acquisition unit around the target biological tissue calculates the personalized blood convection coefficient around the target biological tissue by:

wherein h is_bIs the personalized blood convection coefficient around the target biological tissue, k is the thermal conductivity of the target biological tissue,is the temperature gradient of the surface of the puck and its variation with time, T_bf(t)Is that the surface of said puck changes with timeTemperature of blood flow, T_(t)Is the temperature A of the ice ball surface over time is the surface area of the ice ball, Σ is the surface of the ice ball wherein said T_(t)And T_bf(t)Is preset.

Preferably, the output unit is a display screen.

The application also discloses individualized blood flow perfusion rate assessment equipment of biological tissue includes:

a memory for storing computer executable instructions; and the number of the first and second groups,

a processor for implementing the steps in the method as described hereinbefore when executing the computer-executable instructions.

The present application also discloses a computer-readable storage medium having stored therein computer-executable instructions which, when executed by a processor, implement the steps in the method as described above.

The system and the method for evaluating the personalized blood perfusion rate of the biological tissue, provided by the embodiment of the specification, can evaluate the personalized blood perfusion rate through the size or the temperature state of an ice ball formed after cryoablation without adding additional temperature measurement damage and equipment, are simple and effective, are beneficial to accurately controlling the boundary temperature of the target biological tissue, and further give full play to the curative effect of multi-modal ablation.

The present specification describes a number of technical features distributed throughout the various technical aspects, and if all possible combinations of technical features (i.e. technical aspects) of the present specification are listed, the description is made excessively long. In order to avoid this problem, the respective technical features disclosed in the above summary of the invention of the present application, the respective technical features disclosed in the following embodiments and examples, and the respective technical features disclosed in the drawings may be freely combined with each other to constitute various new technical solutions (which are considered to have been described in the present specification) unless such a combination of the technical features is technically infeasible. For example, in one example, the feature a + B + C is disclosed, in another example, the feature a + B + D + E is disclosed, and the features C and D are equivalent technical means for the same purpose, and technically only one feature is used, but not simultaneously employed, and the feature E can be technically combined with the feature C, then the solution of a + B + C + D should not be considered as being described because the technology is not feasible, and the solution of a + B + C + E should be considered as being described.

Drawings

FIG. 1 is a schematic flow chart of a method for personalized perfusion rate assessment of biological tissue according to a second embodiment of the present application;

FIG. 2 is a schematic diagram of the structure of a system for personalized perfusion rate assessment of biological tissue according to a first embodiment of the present application;

FIG. 3 is a schematic diagram of a probe of a system for personalized perfusion rate assessment of biological tissue according to a first embodiment of the present application;

FIG. 4 is a schematic diagram of experimental results of a system for personalized perfusion rate assessment of biological tissue according to a first embodiment of the present application;

fig. 5 is a schematic diagram of experimental results of a system for evaluating a personalized blood perfusion rate of a biological tissue according to a first embodiment of the present application.

Detailed Description

In the following description, numerous technical details are set forth in order to provide a better understanding of the present application. However, it will be understood by those skilled in the art that the technical solutions claimed in the present application may be implemented without these technical details and with various changes and modifications based on the following embodiments.

Description of partial concepts:

the blood flow condition refers to the abundance degree of blood vessels at the ablation site, the blood flow velocity, and the like, and the blood flow condition may include: blood perfusion rate, blood convection coefficient, etc.

Multi-modality ablation, comprising a cryoablation stage followed by a thermal ablation stage, wherein cryoablation is performed by means including, but not limited to, liquid nitrogen ablation, argon-helium ablation, and cryogenic alcohol ablation, and thermal ablation is performed by means including, but not limited to, radio frequency ablation, microwave ablation, and focused ultrasound ablation.

Ice ball: refers to the frozen tissue formed by freezing the target biological tissue, for example, tumor tissue.

The rewarming temperature, which is the temperature of the iceball after a period of rewarming, can be measured by the tip of the treatment probe, see fig. 3.

Personalized perfusion rate, refers to the amount of blood that a target biological tissue, such as a tumor, flows into the tissue per unit time.

Personalized blood convection coefficient refers to the heat exchange capability between blood and the surface of a target biological tissue such as a tumor.

To make the objects, technical solutions and advantages of the present application more clear, embodiments of the present application will be described in further detail below with reference to the accompanying drawings. The first embodiment of the present application relates to a method for evaluating an individualized blood perfusion rate of a biological tissue, which establishes an algorithm for calculating an individualized blood perfusion rate of a local tissue based on the influence of peripheral blood flow on a target biological tissue after freezing and thawing process by using a boundary convective heat transfer principle, where the individualized blood perfusion rate includes a target biomass, for example, an individualized blood perfusion rate of a tumor tissue, and the flow of the method is as shown in fig. 1, and includes:

step 110: the surface area of the puck is obtained.

Note that this step is preceded by a tissue freezing unit that forms the target biological tissue into an ice ball. As described above, the target biological tissue may be tumor tissue.

It is noted that the surface area of the puck can be obtained in different ways.

Preferably, the surface area of the ice hockey puck can be obtained by any one of an X-ray machine, CT machine, MRI machine or ultrasound machine.

Preferably, the surface area of the ice hockey can also be converted by a relational expression between the surface area and the rewarming temperature of the ice hockey. Specifically, the rewarming temperature of the ice ball can be obtained, the radial and axial dimensions corresponding to the ice ball are determined according to the rewarming temperature, and the surface area corresponding to the ice ball is determined according to the radial and axial dimensions.

Preferably, the rewarming temperature of the ice hockey after a period of time can be obtained by a thermocouple inside the treatment probe.

The advantage of this is that not only the accuracy of rewarming temperature can be guaranteed, but also additional injuries to the patient are avoided.

Preferably, the length of the rewarming process is selectable when obtaining the rewarming temperature of the ice hockey. For example, the rewarming temperature may be selected to be 5 minutes after rewarming. The rewarming period can also be longer or shorter, provided that the temperature at this point is above freezing and above 1 ℃ below normal body temperature.

Preferably, the ice hockey has the following correspondence between radial and axial dimensions and corresponding rewarming temperatures:

l_radial＝-0.1408T₅+43.807, and/or

l_axial＝-0.2816T₅+27.113。

Step 120: obtaining the temperature gradient of the ice ball surface and the change of the ice ball surface along with the time

Preferably, the temperature gradient of the surface of the ice hockey puck and its change over time after the end of the cryotherapy phase can be calculated according to the solid heat transfer equation

Specifically, the temperature gradient of the surface of the ice hockey puck and its change with time were calculated in the following manner

Where T is the target biological tissue temperature, k is the thermal conductivity of the target biological tissue, T is time, ρ is the density of the target biological tissue, and c is the specific heat capacity of the target biological tissue. Therefore, the temperature field distribution is calculated, and the temperature gradient of the ice hockey surface and the change of the ice hockey surface along with time can be obtained.

Step 130, according to the temperature gradient of the ice hockey surface and the change of the ice hockey surface along with the timeAnd the surface area of the ice ball, determining the personalized blood perfusion rate around the target biological tissue.

Preferably, the perfusion rate of blood flow around the target biological tissue may be determined according to the following:

in the formula, ω_bIs the personalized perfusion rate of blood flow around the target biological tissue, k is the thermal conductivity of the target biological tissue,is the temperature gradient of the surface of the puck and its variation with time, T_bf(t)Is the blood flow temperature, T, of the ice ball surface over time_(t)Is the temperature, p, of the ice ball surface over time_bIs the blood density of the target biological tissue, c_bIs the blood specific heat capacity of the target biological tissue, A is the surface area of the ice hockey, and Σ is the surface of the ice hockey, wherein T is_(t)And T_bf(t)Is preset.

Step 140: according to the temperature gradient of the ice hockey surface and the change of the ice hockey surface along with timeAnd the surface area of the ice hockey, determining the personalized blood convection coefficient around the target biological tissue, wherein,

wherein h is_bIs personality around target biological tissueThe blood convection coefficient, k is the thermal conductivity of the target biological tissue,is the temperature gradient of the surface of the puck and its variation with time, T_bf(t)The temperature of the blood flow on the surface of the ice ball can be preset, for example, the temperature of the blood flow on the surface of the ice ball at the end of freezing is preferably 8 ℃, and T_(t)Is the temperature of the surface of the puck over time, may be preset, for example, preferably, may be set to a temperature of 0 c of the puck surface at the end of freezing, a is the puck surface area, Σ is the puck surface area.

It is noted that in the embodiments of the present specification, since the frozen target biological tissue, i.e., the region covered by the ice hockey is slightly larger than the tissue, it is preferable to calculate the ice hockey boundary as the tissue boundary.

Preferably, for a cryoablation steady-state temperature field distribution, assuming that all blood vessels within the puck are destroyed by cryoablation and no blood flows inside the puck, in this case a thermal conduction model can be employed within the puck boundary as a biological heat transfer model, calculating the cryoablation steady-state temperature field distribution, which can be described by the following equation:

where T is the target biological tissue temperature, k is the thermal conductivity of the target biological tissue, T is time, ρ is the density of the target biological tissue, and c is the specific heat capacity of the target biological tissue. Preferably, the thermal conductivity (e.g. of the liver) k of the target biological tissue can be described by the following equation:

wherein T represents the target biological tissue temperature.

Preferably, the specific heat capacity of the target biological tissue (e.g. the specific heat capacity of the liver) c can be described by the following equation:

wherein T represents the temperature of the target biological tissue.

Further, the blood convection condition around the ice hockey puck can be expressed as

Wherein the content of the first and second substances,is the temperature gradient of the surface of the puck and its variation with time, T_bf(t)Is the blood flow temperature of the ice ball surface over time, k is the thermal conductivity of the target biological tissue, T_(t)Is the temperature of the surface of the puck over time, h_bIs the personalized blood convection coefficient.

Preferably, the personalized blood convection coefficient calculated in the above manner may be used as a constant for the re-freezing and rf ablation processes, provided that the conditions outside the pre-frozen target biological tissue are less affected by the core temperature of the pre-frozen target biological tissue.

Preferably, the personalized blood perfusion rate around the target biological tissue calculated in the above manner may be used as a constant for the annealing and rf ablation processes, provided that the conditions outside the pre-frozen target biological tissue are less affected by the core temperature of the pre-frozen target biological tissue.

Step 150: and outputting the personalized blood perfusion rate and the personalized blood convection coefficient around the target biological tissue.

Preferably, the personalized blood convection coefficient and the personalized blood perfusion rate around the target biological tissue are output through a display screen.

It is noted that, as mentioned above, the blood flow condition may comprise an individualized blood perfusion rate, and may further comprise an individualized blood convection coefficient. Therefore, if the blood flow condition only includes the personalized blood perfusion rate, the personalized blood perfusion rate may be outputted in this step.

It should be noted that, in the embodiments of the present specification, the order of obtaining the personalized blood perfusion rate and the personalized blood convection coefficient is not limited.

It is noted that, in the embodiments of the present specification, the target biological tissue may be a tumor tissue, but is not limited thereto, and may be other target biological tissues, such as liver, bladder, prostate, and the like.

The second embodiment of the application relates to a biological tissue personalized blood perfusion rate evaluation system, which is used for establishing an algorithm for calculating a local tissue personalized blood perfusion rate through a boundary convection heat transfer principle based on the influence of peripheral blood flow on a target biological tissue after freezing and melting process. Preferably, taking the example of multi-mode ablation using liquid nitrogen ablation and then radiofrequency ablation as an example, the flow chart is shown in fig. 2, and includes:

an ice hockey size acquisition unit for acquiring the surface area of an ice hockey formed by freezing the target biological tissue;

an ice hockey temperature field distribution acquisition unit for calculating the temperature gradient of the ice hockey surface and the change thereof along with the time

An individualized blood perfusion rate obtaining unit for obtaining the blood perfusion rate according to the temperature gradient of the ice hockey surface and the change of the ice hockey surface along with the timeAnd the surface area of the ice ball, determining the personalized blood perfusion rate around the target biological tissue:

in the formula, ω_bIs the personalized perfusion rate of blood flow around the target biological tissue, k is the thermal conductivity of the target biological tissue,is the temperature gradient of the surface of the puck and its variation with time, T_bf(t)Is the blood flow temperature, T, of the ice ball surface over time_(t)Is the temperature, p, of the ice ball surface over time_bIs the blood density of the target biological tissue, c_bIs the blood specific heat capacity of the target biological tissue, A is the surface area of the ice hockey, and Σ is the surface of the ice hockey, wherein T is_(t)And T_bf(t)Is preset

An output unit for outputting an individualized blood perfusion rate around the target biological tissue.

Further, the system for evaluating a personalized blood perfusion rate of a biological tissue according to the present disclosure may further include:

a personalized blood convection coefficient obtaining unit for obtaining the temperature gradient of the ice hockey surface and the change of the temperature gradient with timeAnd the surface area of the ice hockey, determining the personalized blood convection coefficient around the target biological tissue; and

the output unit is also used for outputting the personalized blood convection coefficient around the target biological tissue.

The system for evaluating the personalized blood perfusion rate of biological tissues will be explained in detail below.

Specifically, the ice hockey size acquiring unit may determine the surface area of the ice hockey by different devices, for example, preferably, the ice hockey size acquiring unit may be one of the following: x-ray machine, CT machine, MRI machine or ultrasound machine.

Preferably, the surface area of the ice hockey can also be converted by a relational expression between the surface area and the rewarming temperature of the ice hockey. Specifically, the ice hockey size obtaining unit may obtain the rewarming temperature of the ice hockey, determine the radial and axial sizes corresponding to the ice hockey according to the rewarming temperature, and determine the surface area corresponding to the ice hockey according to the radial and axial sizes.

Preferably, as shown in fig. 3, the rewarming temperature of the ice hockey after a certain time can be obtained by a thermocouple inside the treatment probe.

The advantage of this is that not only the accuracy of rewarming temperature can be guaranteed, but also additional injuries to the patient are avoided.

Preferably, the length of the rewarming process is selectable when obtaining the rewarming temperature of the ice hockey. For example, the rewarming temperature may be selected to be 5 minutes after rewarming. The rewarming period can also be longer or shorter, provided that the temperature at this point is above freezing and above 1 ℃ below normal body temperature.

For example, it is preferable that the ice hockey has the following correspondence between radial and axial dimensions and corresponding rewarming temperatures:

l_radial＝-0.1408T₅+43.807, and/or

l_axial＝-0.2816T₅+27.113。

Preferably, the ice hockey temperature field distribution acquisition unit calculates the temperature gradient of the ice hockey surface according to the following mode

Where T is the target biological tissue temperature, k is the thermal conductivity of the target biological tissue, T is time, ρ is the density of the target biological tissue, and c is the specific heat capacity of the target biological tissue. Therefore, the temperature field distribution is calculated, and the temperature gradient of the ice hockey surface and the change of the ice hockey surface along with time can be obtained.

Preferably, the personalized blood convection coefficient acquisition unit determines the personalized blood convection coefficient around the target biological tissue according to the temperature gradient ^ T of the surface of the ice hockey and the surface area of the ice hockey by:

wherein h is_bIs the personalized blood convection coefficient around the target biological tissue, k is the thermal conductivity of the target biological tissue,is the temperature gradient of the surface of the puck and its variation with time, T_bf(t)The temperature of the blood flow on the surface of the ice ball can be preset, for example, the temperature of the blood flow on the surface of the ice ball at the end of freezing is preferably 8 ℃, and T_(t)Is the temperature of the surface of the puck over time, may be preset, for example, preferably, may be set to a temperature of 0 c of the puck surface at the end of freezing, a is the puck surface area, Σ is the puck surface area.

It is noted that in the embodiments of the present specification, since the frozen target biological tissue, i.e., the region covered by the ice hockey is slightly larger than the tissue, it is preferable to calculate the ice hockey boundary as the tissue boundary.

It is further noted that, in the examples of the present specification:

preferably, for a cryoablation steady-state temperature field distribution, assuming that all blood vessels within the puck are destroyed by cryoablation and no blood flows inside the puck, in this case a thermal conduction model can be employed within the puck boundary as a biological heat transfer model, calculating the cryoablation steady-state temperature field distribution, which can be described by the following equation:

where T is the target biological tissue temperature, k is the thermal conductivity of the target biological tissue, T is time, ρ is the density of the target biological tissue, and c is the specific heat capacity of the target biological tissue. Preferably, the thermal conductivity (e.g. of the liver) k of the target biological tissue can be described by the following equation:

wherein T represents the target biological tissue temperature.

Preferably, the specific heat capacity of the target biological tissue (e.g. the specific heat capacity of the liver) c can be described by the following equation:

wherein T represents the temperature of the target biological tissue.

Further, the blood convection condition around the ice hockey puck can be expressed as

Wherein the content of the first and second substances,is the temperature gradient of the surface of the puck and its variation with time, T_bf(t)Is the blood flow temperature of the ice ball surface over time, k is the thermal conductivity of the target biological tissue, T_(t)Is the temperature of the surface of the puck over time, h_bIs the personalized blood convection coefficient.

Preferably, the personalized blood convection coefficient obtaining unit may use the personalized blood convection coefficient calculated in the above manner as a constant of the annealing and rf ablation processes, assuming that a condition outside the pre-frozen target biological tissue is less affected by the central temperature of the pre-frozen target biological tissue.

Preferably, the personalized blood convection coefficient obtaining unit may use the personalized blood perfusion rate around the target biological tissue calculated in the above manner as a constant of the complex temperature and the rf ablation process, assuming that the condition of the outside of the pre-frozen target biological tissue is less affected by the central temperature of the pre-frozen target biological tissue.

Preferably, the output unit may be a display screen for outputting the personalized blood perfusion rate and the personalized blood convection coefficient around the target biological tissue. The output unit may also be a data transmission component for transmitting the personalized blood perfusion rate and the personalized blood convection coefficient to the designated target component.

The system for evaluating the personalized blood perfusion rate of the biological tissue provided by the embodiment of the specification can evaluate the personalized blood perfusion rate through the size or the temperature state of an ice ball formed after cryoablation without adding additional temperature measurement damage and equipment, is simple and effective, and is beneficial to accurately controlling the boundary temperature of non-tissues, such as tumor tissues, so that the curative effect of multi-modal ablation is more fully exerted.

The experimental results are as follows:

to verify the technical effect of the examples of the present specification, the applicant carried out in vivo experiments on two pigs, each weighing 55-60 kg. The relative position of the probe and a thermocouple outside the probe is determined through a CT machine, and the relative position is used for determining the size of the ice hockey.

The personalized blood convection coefficients of 3 positions of the pork liver calculated according to the mode are 1514, 1850 and 2077W/m respectively²K. They were used for the simulation of rewarming process with the results shown in table 1. The difference between the temperature of the probe tip of the rewarming probe and the measured temperature is small, the maximum error is about 5%, and the average error is about 3.8%.

This shows that the personalized blood convection coefficient calculated according to the above method is very close to the actual situation, and the simplified model can well simulate the freezing and rewarming processes.

Further, the calculated personalized perfusion rates were 0.048, 0.059 and 0.0651/s, respectively, within the reasonable range recognized in the art.

TABLE 1

Position of	Length of rewarming time(s)	Measured temperature (. degree. C.)	Simulation temperature (. degree. C.)
				1	486	25.3	26.53
2	396	26.81	27.794
				3	335	26.34	26.143

Through a plurality of groups of experimental data, the relationship between the personalized blood convection coefficient and the needle tip temperature of 5 minutes of rewarming is obtained by fitting:

wherein, T₅The tip temperature is shown after 5 minutes of rewarming.

Meanwhile, the size of the ice ball and the rewarming temperature can be converted by the following relational expression:

l_radial＝-0.1408T₅+43.807

l_axial＝-0.2816T₅+27.113

thus, the input to the relational expression calculation module 1022 may be the size of the ice ball formed after cryoablation or the temperature after a period of rewarming.

Further, fig. 4 shows a comparison graph of the personalized blood perfusion rate around the target biological tissue, the personalized blood perfusion rate commonly used in the prior art, and the actually measured thermocouple temperature outside the treatment needle in the simulation of the subsequent heating ablation. Wherein, the curve 41 represents the measured temperature curve, the curve 42 represents the model simulation temperature curve of the non-personalized blood flow condition, and the curve 43 represents the model simulation temperature curve of the personalized blood flow condition obtained by the embodiment of the present specification.

In this experimental position, the outer thermocouple was 9.6mm higher than the probe tip and 8.6mm from the central axis. In this embodiment, the temperature of the probe tip can be controlled at 95 ℃ for 15 minutes by controlling the RF voltage. The model models the entire treatment area including the iceball, and in the whole a Pennes biological heat transfer model was used, the personalized perfusion rate around the target biological tissue was 0.0591/s, whereas the personalized perfusion rate around the non-target biological tissue was set to 0.00641/s, which is commonly used in the liver.

After 15 minutes of treatment, the model simulation temperature curve 43 of the personalized blood flow condition obtained by the embodiment of the specification is closer to the measured temperature curve 41 than the model simulation temperature curve 42 of the non-personalized blood flow condition, which indicates that the personalized blood flow perfusion condition obtained by evaluation is close to the real condition.

Fig. 5 is a graph of three sets of individualized blood perfusion rates around the target biological tissue used in a simulation of subsequent thermal ablation versus actual measured thermocouple temperatures outside the treatment needle. Where curves 51, 53 and 55 represent simulated temperature curves and curves 52, 54 and 56 represent measured temperature curves.

In these 3 experimental positions, the external thermocouples of the examples of this specification were 3.3mm, 9.6mm and 10.4mm higher than the probe tips, respectively, and 5.5mm, 8.6mm and 5.1mm from the central axis, respectively. The treatment probe tips in the first and second positions were controlled at 95 deg.C and the treatment probe tips in the third position were controlled at 65 deg.C.

After 15 minutes of treatment, the maximum temperature difference between the simulated temperature curves 51, 53 and 55 and the measured temperature curves 52, 54 and 56 of the embodiment is within 0.5 ℃, which indicates that the requirement of temperature control accuracy is met.

The system and the method for evaluating the personalized blood perfusion rate of the biological tissue, which are provided by the embodiment of the specification, can evaluate the personalized blood perfusion rate through the size or the temperature state of an ice hockey formed after cryoablation under the condition of not adding additional temperature measurement damage and equipment, are simple and effective, are beneficial to accurately controlling the temperature of a tumor boundary, and thus the curative effect of multi-modal ablation is more fully exerted.

It should be noted that, as will be understood by those skilled in the art, the implementation functions of the modules shown in the embodiment of the system for evaluating a personalized blood perfusion rate of a biological tissue may be understood by referring to the related description of the method for evaluating a personalized blood perfusion rate of a biological tissue. The functions of the modules shown in the embodiment of the system for evaluating a personalized blood perfusion rate of biological tissue described above can be realized by a program (executable instructions) running on a processor, and can also be realized by specific logic circuits. The system for evaluating the personalized blood perfusion rate of the biological tissue according to the embodiment of the present invention may also be stored in a computer-readable storage medium if the system is implemented in the form of a software functional module and sold or used as a stand-alone product. Based on such understanding, the technical solutions of the embodiments of the present application may be essentially implemented or portions thereof contributing to the prior art may be embodied in the form of a software product stored in a storage medium, and including several instructions for causing a computer device (which may be a personal computer, a server, or a network device) to execute all or part of the methods described in the embodiments of the present application. And the aforementioned storage medium includes: various media capable of storing program codes, such as a usb disk, a removable hard disk, a Read Only Memory (ROM), a magnetic disk, or an optical disk. Thus, embodiments of the present application are not limited to any specific combination of hardware and software.

Accordingly, the present application also provides a computer storage medium, in which computer executable instructions are stored, and when executed by a processor, the computer executable instructions implement the method embodiments of the present application.

In addition, the present application also provides a biological tissue personalized blood perfusion rate assessment apparatus, which includes a memory for storing computer executable instructions, and a processor; the processor is configured to implement the steps of the method implementations described above when executing the computer-executable instructions in the memory. The Processor may be a Central Processing Unit (CPU), other general-purpose Processor, a Digital Signal Processor (DSP), an Application Specific Integrated Circuit (ASIC), or the like. The aforementioned memory may be a read-only memory (ROM), a Random Access Memory (RAM), a Flash memory (Flash), a hard disk, or a solid state disk. The steps of the method disclosed in the embodiments of the present invention may be directly implemented by a hardware processor, or implemented by a combination of hardware and software modules in the processor.

It is noted that, in the present patent application, relational terms such as first and second, and the like are used solely to distinguish one entity or action from another entity or action without necessarily requiring or implying any actual such relationship or order between such entities or actions. Also, the terms "comprises," "comprising," or any other variation thereof, are intended to cover a non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements does not include only those elements but may include other elements not expressly listed or inherent to such process, method, article, or apparatus. Without further limitation, the use of the verb "comprise a" to define an element does not exclude the presence of another, same element in a process, method, article, or apparatus that comprises the element. In the present patent application, if it is mentioned that a certain action is executed according to a certain element, it means that the action is executed according to at least the element, and two cases are included: performing the action based only on the element, and performing the action based on the element and other elements. The expression of a plurality of, a plurality of and the like includes 2, 2 and more than 2, more than 2 and more than 2.

All documents mentioned in this application are to be considered as being incorporated in their entirety into the disclosure of this application so as to be subject to modification as necessary. Further, it is understood that various changes or modifications may be made to the present application by those skilled in the art after reading the above disclosure of the present application, and such equivalents are also within the scope of the present application as claimed.