Search

CN-115300094-B - Pulse ablation area prediction device and method based on uncertain parameter quantification

CN115300094BCN 115300094 BCN115300094 BCN 115300094BCN-115300094-B

Abstract

The invention discloses a pulse ablation area prediction device and method based on uncertain parameter quantification, and belongs to the technical field of pulse ablation. The device comprises an acquisition module, an effective ablation boundary determination module and an ablation region determination module. The acquisition module is used for acquiring electrical parameters of a target object, wherein the electrical parameters reflect the electrical characteristics of tissue and blood of the target object. The effective ablation boundary determining module is used for determining an effective ablation boundary according to the electrical parameters and a pre-acquired ablation electric field distribution model. The ablation region determination module is used for determining a region enclosed by the effective ablation boundary and the tissue surface as an ablation region.

Inventors

  • WANG HUI
  • ZHAO FENG
  • GUO WENJUAN
  • ZHAO QIANCHENG
  • ZHANG WEI

Assignees

  • 上海商阳医疗科技有限公司

Dates

Publication Date
20260505
Application Date
20220815

Claims (15)

  1. 1. A pulse ablation zone prediction apparatus based on uncertain parameter quantification, the apparatus comprising: the acquisition module is used for acquiring electrical parameters of a target object, wherein the electrical parameters reflect the electrical characteristics of tissue and blood of the target object; The effective ablation boundary determining module is used for determining an effective ablation boundary according to the electrical parameters and a pre-acquired ablation electric field distribution model; an ablation region determining module, configured to determine a region enclosed by the effective ablation boundary and the tissue surface as an ablation region; The effective ablation boundary determination module includes: The first determining unit is used for determining the ablation depth according to the electrical parameters and a pre-acquired ablation electric field distribution model; a second determining unit, configured to determine the effective ablation boundary according to the ablation depth; The ablation electric field distribution model comprises the following formula: , wherein E is the field intensity of an electric field, the unit is V/cm or V/m, U is the ablation voltage, the unit is V, x is the coordinate value of the tissue in the depth direction, and the unit is cm or m; As a fitting coefficient, the unit is cm -1 or m -1 ; is a weight coefficient of the blood conductivity, As a function of the conductivity of the blood; to organize the weight coefficients of the initial conductivity, As a function of the initial conductivity of the tissue; as a weighting factor for the conductivity growth factor, As a function of the conductivity growth factor; for the weight coefficient of the electric field intensity corresponding to the center point of the transition region, As a function of the electric field strength corresponding to the center point of the transition region; the weighting coefficients for the coefficients are fitted for the dynamic conductivity, Fitting coefficients as a function of the dynamic conductivity; Are all the non-dimensional coefficients of the three-dimensional coefficients, Is in cm -1 or m -1 ; in the ablation electric field distribution model, , Wherein, the value of i is 1, 2, 3,4 and 5; When any one of the electrical parameters is a preset variation value and the other parameters are preset fixed values, according to the ablation depth variance obtained by the ablation electric field distribution model; Wherein, the Representing the corresponding variance when the blood conductivity is a preset change value, wherein the change range of the blood conductivity is 0.1-1S/m; representing the corresponding variance when the initial conductivity of the tissue is a preset change value, wherein the change range of the initial conductivity of the tissue is 0.05-0.9S/m; For the corresponding variance when the conductivity growth factor is a preset change value, the change range of the conductivity growth factor is 1-6; For the variance corresponding to the electric field intensity corresponding to the central point of the transition area when the electric field intensity corresponding to the central point of the transition area is a preset change value, the change range of the electric field intensity corresponding to the central point of the transition area is 200-1200V/cm; for the variance corresponding to the dynamic conductivity fitting coefficient when the dynamic conductivity fitting coefficient is a preset change value, the change range of the conductivity fitting coefficient is 0.00001-0.00003; And the sum of ablation depth variances obtained when any electrical parameter is a preset change value.
  2. 2. The apparatus of claim 1 wherein, in the ablative electric field distribution model, A kind of electronic device A kind of electronic device A kind of electronic device A kind of electronic device , Wherein, the The unit is S/m or S/cm, which is the blood conductivity; For the initial conductivity of the tissue, the unit is S/m or S/cm; is a conductivity growth factor, without units; The electric field intensity corresponding to the central point of the transition area is in units of V/cm or V/m; fitting coefficients for dynamic conductivity without units; ~ 、 ~ 、 ~ 、 ~ respectively, the fitting coefficients are predetermined in advance, ~ And ~ In units of cm -1 or m -1 , ~ And ~ Are each m/S or cm/S.
  3. 3. The device according to claim 1, characterized in that the second determination unit is specifically adapted to determine a corresponding field strength contour from the ablation depth, the corresponding field strength contour being taken as the effective ablation boundary.
  4. 4. The apparatus of claim 1, wherein the electrical parameter comprises blood conductivity, and wherein the acquisition module comprises: A first voltage applying unit for applying a first test voltage to a first electrode placed in a target area, the first electrode contacting blood of the target object and not contacting tissue of the target object; a first acquisition unit configured to acquire the first test voltage and a first current through the first electrode; and a third determining unit, configured to determine the blood conductivity according to the first voltage and the first current, and a first fitting function obtained in advance, where the first fitting function characterizes a relationship between the blood conductivity and the current and the voltage.
  5. 5. The apparatus of claim 4, wherein the electrical parameter further comprises an initial conductivity of the tissue, and wherein the acquisition module further comprises: A second voltage applying unit for applying a second test voltage to a second electrode placed in the target area, the second electrode contacting blood and tissue of the target object; a second acquisition unit configured to acquire the second test voltage and a second current through the second electrode; and a fourth determining unit, configured to determine the tissue initial conductivity according to the second voltage, the second current and the blood conductivity, and a second fitting function obtained in advance, where the second fitting function characterizes a relationship between tissue conductivity and voltage, current and blood conductivity.
  6. 6. The apparatus of claim 5, wherein the electrical parameters further comprise dynamic conductivity parameters that characterize a change in tissue conductivity with electric field strength during pulsed ablation, the acquisition module further comprising: The third acquisition unit is used for acquiring the mapping relation between the simulation current and the dynamic conductivity parameter; a third voltage applying unit for applying a pulse ablation voltage to a third electrode disposed at the target region, monitoring a current passing through the third electrode, and taking the stably output current as a third current; a fifth determining unit configured to determine a simulation current closest to the third current among a plurality of simulation currents stored in advance; And a sixth determining unit, configured to determine the dynamic conductivity parameter according to the closest simulation current and the mapping relationship.
  7. 7. The apparatus of claim 6, wherein the dynamic conductivity parameter comprises a conductivity growth factor and an electric field strength corresponding to a transition zone center point, the conductivity growth factor characterizing a magnitude of change in the tissue conductivity, the electric field strength corresponding to the transition zone center point being a mean value of the electric field strength corresponding to the tissue conductivity change start time and the electric field strength corresponding to the tissue conductivity change end time; The third acquisition unit includes: The first acquisition subunit is used for acquiring functions of the electric field intensity relations corresponding to the simulation current and the central point of the transition region under different conductivity growth factors; And the fitting subunit is used for fitting the function to obtain the mapping relation between the simulation current and the dynamic conductivity parameter.
  8. 8. The apparatus according to claim 7, wherein the fifth determining unit includes: the second acquisition subunit is used for acquiring the numerical value coincidence degree of the third current and a plurality of prestored simulation currents; and the first determination subunit is used for determining the simulation current corresponding to the minimum value in the numerical value matching degree as the closest simulation current.
  9. 9. The apparatus of claim 1, wherein the electrical parameters comprise dynamic conductivity fit coefficients, and wherein the acquisition module comprises: and the fourth acquisition unit is used for acquiring a first data set, wherein the first data set comprises a plurality of dynamic conductivity fitting coefficients, and the dynamic conductivity fitting coefficients are distributed in a Gaussian mode.
  10. 10. The apparatus of claim 1, wherein the electrical parameter further comprises an electric field strength threshold that characterizes an electric field strength of tissue cells when irreversible electroporation occurs, the acquisition module comprising: And a fifth acquisition unit, configured to acquire a second data set, where the second data set includes a plurality of electric field intensity thresholds, and the plurality of electric field intensity thresholds are in gaussian distribution.
  11. 11. The apparatus of claim 1, wherein the apparatus further comprises: the matching degree parameter acquisition module is used for acquiring a matching degree parameter according to the ablation depth and the target ablation depth, and the matching degree parameter characterizes the matching degree of the pulse parameter corresponding to the ablation depth and the pulse parameter for realizing the target ablation depth; and the judging module is used for judging whether to adjust the pulse ablation parameters according to the matching degree parameters.
  12. 12. The apparatus of claim 11, wherein the matching degree parameter comprises a mean value of the ablation depth, and wherein the determining module is specifically configured to: adjusting the pulse ablation parameters in response to the difference between the mean value and the target ablation depth exceeding a preset threshold, or The matching degree parameter comprises the probability of the tissue at the target ablation depth being ablated, and the judging module is specifically used for responding to the probability being less than 70 percent and adjusting the pulse ablation parameter.
  13. 13. A pulse ablation zone prediction method based on uncertain parameter quantification, the method comprising: Acquiring electrical parameters of a target object, wherein the electrical parameters reflect the electrical characteristics of tissue and blood of the target object; determining an effective ablation boundary according to the electrical parameters and an ablation electric field distribution model obtained in advance; Determining an area surrounded by the effective ablation boundary and the tissue surface as an ablation area; the method specifically comprises the following steps: determining the ablation depth according to the electrical parameters and a pre-acquired ablation electric field distribution model; determining the effective ablation boundary according to the ablation depth; The ablation electric field distribution model comprises the following formula: , wherein E is the field intensity of an electric field, the unit is V/cm or V/m, U is the ablation voltage, the unit is V, x is the coordinate value of the tissue in the depth direction, and the unit is cm or m; As a fitting coefficient, the unit is cm -1 or m -1 ; is a weight coefficient of the blood conductivity, As a function of the conductivity of the blood; to organize the weight coefficients of the initial conductivity, As a function of the initial conductivity of the tissue; as a weighting factor for the conductivity growth factor, As a function of the conductivity growth factor; for the weight coefficient of the electric field intensity corresponding to the center point of the transition region, As a function of the electric field strength corresponding to the center point of the transition region; the weighting coefficients for the coefficients are fitted for the dynamic conductivity, Fitting coefficients as a function of the dynamic conductivity; Are all the non-dimensional coefficients of the three-dimensional coefficients, Is in cm -1 or m -1 ; in the ablation electric field distribution model, , Wherein, the value of i is 1, 2, 3,4 and 5; When any one of the electrical parameters is a preset variation value and the other parameters are preset fixed values, according to the ablation depth variance obtained by the ablation electric field distribution model; Wherein, the Representing the corresponding variance when the blood conductivity is a preset change value, wherein the change range of the blood conductivity is 0.1-1S/m; representing the corresponding variance when the initial conductivity of the tissue is a preset change value, wherein the change range of the initial conductivity of the tissue is 0.05-0.9S/m; For the corresponding variance when the conductivity growth factor is a preset change value, the change range of the conductivity growth factor is 1-6; For the variance corresponding to the electric field intensity corresponding to the central point of the transition area when the electric field intensity corresponding to the central point of the transition area is a preset change value, the change range of the electric field intensity corresponding to the central point of the transition area is 200-1200V/cm; for the variance corresponding to the dynamic conductivity fitting coefficient when the dynamic conductivity fitting coefficient is a preset change value, the change range of the conductivity fitting coefficient is 0.00001-0.00003; And the sum of ablation depth variances obtained when any electrical parameter is a preset change value.
  14. 14. An electronic device comprising a memory, a processor, and a computer program stored on the memory and executable on the processor, wherein the processor implements the pulse ablation zone prediction method of claim 13 when the computer program is executed.
  15. 15. A computer readable storage medium having stored thereon a computer program, wherein the computer program when executed by a processor implements the pulse ablation zone prediction method of claim 13.

Description

Pulse ablation area prediction device and method based on uncertain parameter quantification Technical Field The invention relates to the technical field of pulse ablation, in particular to a pulse ablation area prediction device and method based on uncertain parameter quantification. Background Pulse ablation is a non-thermal ablation modality and is an effective method for treating persistent atrial fibrillation and ventricular arrhythmias. By applying an ultra-fast electric field to the target tissue, irreversible nanoscale pores are formed in the tissue cells, allowing the cell contents to leak out, destroying the stability of the cells, and thus leading to cell death. The therapeutic effect of pulse ablation mainly depends on pulse parameters, and excessive ablation can be caused by too high pulse dose, and incomplete treatment and easy recurrence can be caused by too low pulse dose. To ensure proper pulse delivery, it is necessary to determine the ablation zone and set pulse parameters based on the determined ablation zone. Therefore, accurately determining the ablation region is a key step in ensuring the pulse ablation effect. In the related art, an ablation area is determined according to the position of the tissue to be ablated, and the electrical characteristics of the physiological tissue of the patient are not combined, so that the determined ablation area is inaccurate, and the pulse ablation effect is affected. Disclosure of Invention The invention aims to overcome the defect that the pulse ablation area determination is inaccurate for different patients in the prior art, and provides a pulse ablation area prediction device and method based on uncertain parameter quantification. The invention solves the technical problems by the following technical scheme: In a first aspect, the present invention provides a pulse ablation zone prediction apparatus based on uncertain parameter quantification, the apparatus comprising: the acquisition module is used for acquiring electrical parameters of a target object, wherein the electrical parameters reflect the electrical characteristics of tissue and blood of the target object; The effective ablation boundary determining module is used for determining an effective ablation boundary according to the electrical parameters and a pre-acquired ablation electric field distribution model; and the ablation area determining module is used for determining an area enclosed by the effective ablation boundary and the tissue surface as an ablation area. Optionally, the effective ablation boundary determination module includes: The first determining unit is used for determining the ablation depth according to the electrical parameters and a pre-acquired ablation electric field distribution model; and the second determining unit is used for determining the effective ablation boundary according to the ablation depth. Optionally, the ablation electric field distribution model includes the following formula: wherein E is the field intensity, U is the ablation voltage, x is the coordinate value of the depth direction of the tissue, In order to fit the coefficients of the coefficients, Is a weight coefficient of the blood conductivity,As a function of the conductivity of the blood; to organize the weight coefficients of the initial conductivity, As a function of the initial conductivity of the tissue; as a weighting factor for the conductivity growth factor, As a function of the conductivity growth factor; for the weight coefficient of the electric field intensity corresponding to the center point of the transition region, As a function of the electric field strength corresponding to the center point of the transition region; the weighting coefficients for the coefficients are fitted for the dynamic conductivity, As a function of the fitting coefficients for the dynamic conductivity. Optionally, in the ablative electric field distribution model, A kind of electronic device A kind of electronic device A kind of electronic device A kind of electronic device , Wherein, the For the electrical conductivity of the blood,For the initial conductivity of the tissue,As a factor of the increase in the electrical conductivity, For the electric field strength corresponding to the center point of the transition region,For the dynamic conductivity fit coefficients,~、~、~、~Respectively, predetermined fitting coefficients. Optionally, in the ablative electric field distribution model, Wherein, the value of i is 1, 2, 3,4 and 5; When any one of the electrical parameters is a preset variation value and the other parameters are preset fixed values, according to the ablation depth variance obtained by the ablation electric field distribution model; Wherein, the Representing the corresponding variance when the blood conductivity is a preset change value, wherein the change range of the blood conductivity is 0.1-1S/m; representing the corresponding variance when the initial conductivity of the tiss