System and device for generating needle distribution scheme before tumor ablation operation based on COMSOL

Abstract

The invention is suitable for the technical field of pulsed electric field treatment and computer simulation, and provides a COMSOL-based tumor ablation preoperative needle arrangement scheme generation system and device. In the embodiment of the invention, a three-dimensional simulation geometric model of a target tissue is constructed through COMSOL multi-physics field simulation software, and then an initial needle arrangement scheme is determined through the three-dimensional simulation geometric model. The method comprises the following steps: establishing a three-dimensional simulation geometric model, determining material parameters and boundary conditions of the three-dimensional simulation geometric model, constructing a multi-physical field tumor ablation simulation model, determining electric field energy distribution data during an ablation operation according to an initial needle arrangement scheme through the simulation model, and comparing the data with a target region to determine a target needle arrangement scheme. According to the method, the electric field and the thermal field are combined in construction of the multi-physics field tumor ablation simulation model, and the corresponding conductivity is set when the model material parameters are set, so that the problems of lack of multi-physics field coupling analysis and insufficient dynamic parameter optimization in the prior art are solved, and the accuracy of a needle arrangement scheme is improved.

Description

Technical Field

[0001] The present application relates to the field of pulsed electric field therapy and computer simulation technology, and in particular to a system and device for generating a needle placement plan before tumor ablation surgery based on COMSOL. Background Art

[0002] Cancer treatment remains a major global public health challenge, and new, safe and effective therapies are urgently needed to enable patients to achieve more predictable, long-term, cancer-free survival and a higher quality of life. With the increasing number of new cancer cases and deaths, the cancer prevention and control situation is becoming increasingly critical. Currently, surgery and chemotherapy remain the mainstays of treatment for most patients with advanced cancer. However, even after treatment, conventional therapies often fail to achieve optimal results due to drug resistance in tumor cells, resulting in poor prognosis, high recurrence rates, and even treatment failure. In recent years, pulsed electric field (PEF) tumor therapy, an emerging physical therapy, has demonstrated significant advantages in treating highly malignant tumors and improving prognosis by applying high-voltage pulsed electric fields to target tissue cells. This technology targets the cell membrane or the membranes of tumor organelles, including the nucleus and mitochondria, disrupting the intracellular physiological balance, inducing intracellular electrochemical effects and immunogenic cell death, initiating the patient's tumor-immune cycle, stimulating an anti-tumor immune response, and leading to cell death.

[0003] Tumor ablation (such as radiofrequency ablation, microwave ablation, and pulsed electric field ablation) is an important method for treating solid tumors, but its efficacy is highly dependent on the accuracy of the preoperative needle placement plan. However, existing technologies have the following problems: 1) Lack of multi-physics coupling analysis: Existing simulation tools mostly focus on a single physical field (such as electric field or thermal field), which makes it difficult to fully reflect the interaction between the electric field, thermal field, and biological tissue during the ablation process.

[0004] 2) Insufficient dynamic parameter optimization: Existing methods do not fully consider the dynamic changes of parameters such as tissue electrical conductivity and thermal conductivity (such as electroporation effect and thermal damage accumulation), resulting in a large deviation between simulation results and actual conditions. Summary of the Invention

[0005] In view of this, the present application provides a system and device for generating a preoperative needle placement plan for tumor ablation based on COMSOL to solve the problems of lack of multi-physics field coupling analysis and insufficient dynamic parameter optimization in the existing preoperative needle placement plan for tumor ablation.

[0006] In a first aspect, the present application provides a system for generating a needle placement plan before tumor ablation surgery based on COMSOL, the system comprising a three-dimensional model building module, a needle placement plan determination module, a multi-physics field coupling modeling module, and a needle placement plan adjustment module; The three-dimensional model construction module is used to determine the geometric information of the target area of the target tissue, and construct a three-dimensional simulation geometric model of the target tissue based on the geometric information using COMSOL multi-physics field simulation software, wherein the three-dimensional simulation geometric model includes the size parameters and geometric features of the target tissue; The needle arrangement plan determination module is used to determine the needle arrangement plan for the tumor ablation surgery of the target tissue according to the three-dimensional simulation geometric model, wherein the needle arrangement plan includes the specifications, quantity, spatial arrangement plan of the electrode needles, and the energy release size between the electrode needles; The multi-physics field coupling modeling module is used to determine the tissue material parameters and boundary conditions of the three-dimensional simulation geometric model, construct a multi-physics field tumor ablation simulation model based on the tissue material parameters and boundary conditions, and simulate and determine the electric field energy distribution data of the target tissue during tumor ablation surgery through the multi-physics field tumor ablation simulation model; The needle layout plan adjustment module is used to determine the comparison results of the electric field energy distribution data and the target area. If the comparison result meets expectations, the needle layout plan is determined as the target needle layout plan. If the comparison result does not meet expectations, the needle layout plan is adjusted, and the updated electric field energy distribution data is determined based on the adjusted needle layout plan. Until the comparison result between the updated electric field energy distribution data and the target area meets expectations, the adjusted needle layout plan is determined as the target needle layout plan.

[0007] Optionally, determining tissue material parameters of the three-dimensional simulation geometric model includes: The tissue types of the tumor area and the target area in the target tissue, as well as the materials of the conductive part and the insulating part of the electrode needle are set, and the conductivity of each tissue type and material is set.

[0008] Optionally, the boundary conditions include electrode needle boundaries, tissue boundaries, current boundaries, and thermal damage boundaries; The electrode needle boundary is used to determine the electrode potential on the surface of the electrode needle, the tissue boundary is used to determine the insulation state of the outer boundary of the target tissue, the current boundary is used to determine the potential inside the target tissue, and the thermal damage boundary is used to determine the thermal damage state of the target tissue.

[0009] Optionally, setting the conductivity of the tissue type of the tumor area includes: The Heaviside function is used to characterize the dynamic change of the conductivity of the tissue type in the tumor area with the electric field strength, and the conductivity is determined to be , where A is the growth factor of the conductivity after tissue electroporation, is the middle point of the transition area, is half of the transition area, is the temperature growth factor, is the initial ambient temperature.

[0010] Optionally, the method for determining the thermal damage boundary includes: The Arrhenius first-order kinetic model was used to evaluate the thermal damage to tissues caused by temperature increase during tissue electroporation ablation. Determine a dimensionless thermal damage function Ω(t), and then determine the thermal damage state of the target tissue using the thermal damage function Ω(t) and a preset threshold, wherein: is the function of tissue temperature and time, is the activation energy, is the gas constant.

[0011] A second aspect of the present application provides a device for generating a needle placement plan before tumor ablation surgery based on COMSOL, the device comprising: a three-dimensional model construction unit, configured to determine geometric information of a target area of a target tissue, and to construct a three-dimensional simulation geometric model of the target tissue based on the geometric information using COMSOL multi-physics simulation software, wherein the three-dimensional simulation geometric model includes dimensional parameters and geometric features of the target tissue; a needle arrangement plan determining unit, configured to determine a needle arrangement plan for tumor ablation surgery on the target tissue based on the three-dimensional simulation geometric model, wherein the needle arrangement plan includes specifications, quantity, spatial arrangement plan of electrode needles, and energy release between electrode needles; a multi-physics field coupling modeling unit, configured to determine tissue material parameters and boundary conditions of the three-dimensional simulation geometric model, construct a multi-physics field tumor ablation simulation model based on the tissue material parameters and boundary conditions, and simulate and determine electric field energy distribution data of the target tissue during tumor ablation surgery using the multi-physics field tumor ablation simulation model; The needle layout plan adjustment unit is used to determine the comparison result of the electric field energy distribution data and the target area. If the comparison result meets the expectation, the needle layout plan is determined as the target needle layout plan. If the comparison result does not meet the expectation, the needle layout plan is adjusted, and the updated electric field energy distribution data is determined according to the adjusted needle layout plan. Until the comparison result of the updated electric field energy distribution data and the target area meets the expectation, the adjusted needle layout plan is determined as the target needle layout plan.

[0012] Optionally, determining tissue material parameters of the three-dimensional simulation geometric model in the multi-physics field coupling modeling unit includes: The tissue types of the tumor area and the target area in the target tissue, as well as the materials of the conductive part and the insulating part of the electrode needle are set, and the conductivity of each tissue type and material is set.

[0013] Optionally, the boundary conditions in the multi-physics field coupling modeling unit include electrode needle boundary, tissue boundary, current boundary and thermal damage boundary; The electrode needle boundary is used to determine the electrode potential on the surface of the electrode needle, the tissue boundary is used to determine the insulation state of the outer boundary of the target tissue, the current boundary is used to determine the potential inside the target tissue, and the thermal damage boundary is used to determine the thermal damage state of the target tissue.

[0014] Optionally, setting the conductivity of the tissue type of the tumor area in the multi-physics field coupling modeling unit includes: The Heaviside function is used to characterize the dynamic change of the conductivity of the tissue type in the tumor area with the electric field strength, and the conductivity is determined to be , where A is the growth factor of the conductivity after tissue electroporation, is the middle point of the transition area, is half of the transition area, is the temperature growth factor, is the initial ambient temperature.

[0015] Optionally, the method for determining the thermal damage boundary in the multi-physics field coupling modeling unit includes: The Arrhenius first-order kinetic model was used to evaluate the thermal damage to tissues caused by temperature increase during tissue electroporation ablation. Determine a dimensionless thermal damage function Ω(t), and then determine the thermal damage state of the target tissue using the thermal damage function Ω(t) and a preset threshold, wherein: is the function of tissue temperature and time, is the activation energy, is the gas constant.

[0016] In the embodiment provided in the present application, a three-dimensional simulation geometric model of the target tissue is first constructed according to the geometric information of the target area of the target tissue using COMSOL multi-physics simulation software, and then the initial needle layout scheme is determined using the three-dimensional simulation geometric model. The material parameters and boundary conditions of the three-dimensional simulation geometric model are then determined, and a multi-physics tumor ablation simulation model is constructed. The electric field energy distribution data of the target tissue during tumor ablation surgery is then determined using the simulation model according to the initial needle layout scheme. Finally, the electric field energy distribution data is compared with the above-mentioned target area, and the needle layout scheme is adjusted according to the comparison results, thereby determining the target needle layout scheme with the highest degree of overlap between the target area and the electric field energy distribution data. The present application combines the electric field and the thermal field in constructing the multi-physics tumor ablation simulation model, and sets the corresponding conductivity when setting the model material parameters. This solves the problems of lack of multi-physics coupling analysis and insufficient dynamic parameter optimization in the prior art, and improves the accuracy of the needle layout scheme. BRIEF DESCRIPTION OF THE DRAWINGS

[0017] Figure 1 A system module diagram provided for an embodiment of the present application; Figure 2 The target tissue profile provided in the embodiment of the present application; Figure 3 A model diagram provided for an embodiment of the present application; Figure 4 A conductivity-electric field strength relationship curve provided in an embodiment of the present application; Figure 5 Provide ablation effect diagram for the embodiment of this application; Figure 6 Another ablation effect diagram provided in an embodiment of the present application; Figure 7 A diagram of the device structure provided in an embodiment of the present application; Figure 8 A schematic diagram of the internal structure of a computer device provided in an embodiment of the present application. DETAILED DESCRIPTION

[0018] Exemplary embodiments will be described in detail herein, with examples illustrated in the accompanying drawings. In the following description, when referring to the drawings, identical numerals in different figures represent identical or similar elements, unless otherwise indicated. The embodiments described in the following exemplary embodiments are not intended to represent all embodiments consistent with the present application. Rather, they are merely examples of apparatus and methods consistent with certain aspects of the present application, as detailed in the appended claims.

[0019] The terms used in this application are for the purpose of describing specific embodiments only and are not intended to limit this application. As used in this application and the appended claims, the singular forms "a," "an," "the," and "the" are intended to include the plural forms, unless the context clearly indicates otherwise. It should also be understood that the term "and / or" as used herein refers to and encompasses any and all possible combinations of one or more of the associated listed items.

[0020] It should be understood that although the terms first, second, third, etc. may be used in this application to describe various information, such information should not be limited to these terms. These terms are only used to distinguish information of the same type from each other. For example, without departing from the scope of this application, first information may also be referred to as second information, and similarly, second information may also be referred to as first information. Depending on the context, the word "if" as used herein may be interpreted as "at the time of" or "when" or "in response to determining".

[0021] The present application provides a system for generating a preoperative needle placement plan for tumor ablation based on COMSOL to solve the problems of the lack of multi-physics field coupling analysis and insufficient dynamic parameter optimization in the existing preoperative needle placement plan for tumor ablation.

[0022] The following specific embodiments are used to describe the technical solution of the present application in detail. The following specific embodiments can be combined with each other, and the same or similar concepts or processes may not be described in detail in some embodiments.

[0023] like Figure 1 As shown in FIG, a module diagram of a system for generating a needle placement plan for tumor ablation surgery based on COMSOL is provided in this application. The functions and effects of each module are described below.

[0024] 1. 3D model construction module. This module is used to determine the geometric information of the target area of the target tissue and construct a 3D simulation geometric model of the target tissue based on the geometric information using COMSOL multi-physics field simulation software. The 3D simulation geometric model includes the dimensional parameters and geometric features of the target tissue.

[0025] In this module, two-dimensional slice data of target tissue can be obtained through medical imaging technologies such as CT, MRI or ultrasound. Then, image processing software such as ITK-SNAP, 3D Slicer, Mimics, etc. can be used to segment the medical images and extract the contours of the target tissue. Figure 2As shown. The segmented two-dimensional contour data is then converted into three-dimensional geometric information. The three-dimensional geometric information is then imported into the COMSOL software, and the three-dimensional simulation geometric model is constructed using COMSOL's geometric modeling tools (such as stretching, rotation, sweeping, etc.). In COMSOL, the size of the geometric model can be controlled by defining parameters (such as length, width, height, radius, etc.). These parameters can be set in the "Global Definition" or "Geometry" section. The model is further refined in COMSOL based on the geometric characteristics of the target tissue (such as curvature, thickness, symmetry, etc.). Boolean operations, chamfering, stretching and other tools can be used to improve the geometric features. After setting the material parameters and boundary conditions of the model, it can be converted into a multi-physics tumor ablation simulation model that can simulate ablation surgery.

[0026] 2. Needle placement scheme determination module. This module is used to determine the needle placement scheme for tumor ablation surgery on the target tissue based on the three-dimensional simulation geometric model, wherein the needle placement scheme includes the specifications, quantity, spatial arrangement scheme of the electrode needles, and the energy release between the electrode needles.

[0027] In this module, electrode needle specifications include parameters such as electrode type, length, diameter, and performance. This is because the appropriate electrode needle type, such as radiofrequency or microwave, must be selected based on the tumor's size and location. The length and diameter of the electrode needle are selected based on the tumor's geometry to ensure adequate coverage. The performance parameters of the electrode needle determine the energy delivery range of a single needle, such as the thermal field diameter and ablation volume.

[0028] The number of electrode needles required can be estimated by the ratio of the total tumor volume to the ablation volume per needle. Spatial arrangement options include multi-needle parallel and multi-needle cross arrangements, which can be determined based on the shape and size of the tumor. For example, for medium-sized, regularly shaped tumors, a multi-needle parallel arrangement can be used. For large or irregularly shaped tumors, a cross arrangement can be used to improve energy field uniformity.

[0029] 3. Multi-physics coupling modeling module. This module is used to determine the tissue material parameters and boundary conditions of the three-dimensional simulation geometric model, construct a multi-physics tumor ablation simulation model based on the tissue material parameters and boundary conditions, and use the multi-physics tumor ablation simulation model to simulate and determine the electric field energy distribution data of the target tissue during tumor ablation surgery.

[0030] In this module, tissue material parameters include tissue parameters and electrode needle material parameters. Normal tissue and tumor tissue can be set according to the location of the tumor. For example, when the tumor is near the liver, the normal tissue is set to liver tissue. When it is near the pancreas, the normal tissue is set to liver tissue and the tumor tissue is set to the corresponding tumor tissue. The structure is as follows: Figure 3 As shown in the figure, each tissue has a traditional numerical model that can be used for modeling. Material parameters include the conductive and insulating parts of the electrode needle and the material. The corresponding materials are set according to the electrode needles determined in the needle layout plan.

[0031] Boundary conditions define the physical constraints and environmental interactions of the simulation model. For example, a constant temperature boundary is set to simulate the surface of tissue in contact with a constant temperature environment, such as setting the skin surface temperature to 37°C. A blood perfusion boundary is used to add blood perfusion terms to the tissue domain to describe the effect of blood flow on heat removal.

[0032] In another embodiment, determining the tissue material parameters of the three-dimensional simulation geometric model includes: The tissue types of the tumor area and the target area in the target tissue, as well as the materials of the conductive part and the insulating part of the electrode needle are set, and the conductivity of each tissue type and material is set.

[0033] From a tissue perspective, determining whether electroporation occurs in cells inside the tissue is mainly done by judging whether the electric field inside the target area is greater than the threshold for tissue electroporation, or by judging whether the tissue conductivity value has changed. In this embodiment, corresponding conductivity is set for each tissue and material, and the accuracy of the needle layout scheme is judged by monitoring the conductivity. Among them, the conductivity of the materials of normal tissue, the conductive part of the electrode needle, and the insulating part are all fixed values, which are set according to the tissue type of normal tissue and the material of the electrode needle. The conductivity setting method for the tissue type of the tumor area (hereinafter referred to as tumor tissue) is as follows: Since conductivity is affected not only by electric field changes but also by temperature, an expression of conductivity, temperature, and electric field is constructed and introduced into the material properties. This embodiment uses a smooth Heaviside function to characterize the dynamic change of tissue conductivity with electric field strength. The function of tissue conductivity with respect to electric field strength and temperature is expressed as , where A is the growth factor of the conductivity of the tissue in the tumor area after electroporation, is the middle value of the electric field strength interval corresponding to the conductivity change interval; It is half of the difference between the maximum and minimum values of the electric field intensity corresponding to the conductivity change range, such as Figure 4 As shown, assuming that the electric field strength range corresponding to the conductivity change area is 200-600, then Then it is (200+600) / 2=400, Then it is (600-200) / 2=200; is the temperature growth factor, is the initial temperature of the environment, A, , Both can be pre-set fixed values.

[0034] Through the above method, the conductivity of each tissue material parameter used for modeling can be determined.

[0035] In another embodiment, the boundary conditions include an electrode needle boundary, a tissue boundary, a current boundary, and a thermal damage boundary; The electrode needle boundary is used to determine the electrode potential on the surface of the electrode needle, the tissue boundary is used to determine the insulation state of the outer boundary of the target tissue, the current boundary is used to determine the potential inside the target tissue, and the thermal damage boundary is used to determine the thermal damage state of the target tissue.

[0036] In this example, the electrode needles act as carriers for pulsed electric field application. Dirichlet boundary conditions are applied to the electrode needle boundaries, with the positive needle surface receiving the input pulse and the negative needle surface grounded. Neumann boundary conditions are applied to the outer boundary of the tumor tissue, simulating an insulating environment. This combination confines the electric field and thermal effects to the interior of the tumor tissue, providing conservative simulation results for safety assessments.

[0037] The current boundary adopts the Maxwell equations theory, and the potential inside the tumor tissue satisfies the external electric field. ,in, is the relative dielectric constant inside the tumor tissue; Vacuum dielectric constant; = is the electrical conductivity of the tumor tissue. Since the electrical conductivity inside the tumor tissue is affected not only by tissue electroporation but also by Joule heating inside the tumor tissue, the first pulse is used to evaluate the electrical conductivity of the tumor tissue as a function of the electric field. As multiple pulses are applied, the temperature rise inside the tumor tissue will also increase the electrical conductivity of the tumor tissue.

[0038] The thermal damage boundary is evaluated by the Arrhenius first-order kinetic model to evaluate the thermal damage caused by the temperature increase during tissue electroporation ablation. Determine a dimensionless thermal damage function Ω(t), and then determine the thermal damage state of the target tissue using the thermal damage function Ω(t) and a preset threshold, wherein: is the function of tissue temperature and time, is the activation energy, is the gas constant.

[0039] In this embodiment, it is necessary to first calculate the bioheat transfer process in the tumor tissue, which can be calculated using the Pennes bioheat transfer model. Calculate to determine the internal temperature T of the tissue. And determine the function of temperature and time change based on the relationship between temperature T and time ,in, is the tissue specific heat capacity, is the tissue thermal conductivity, is the blood density, is the blood perfusion rate, is the specific heat of blood, is the arterial blood temperature, is tissue metabolic heat, The Joule heat generated by the pulse can be detected and determined in real time on the target tissue.

[0040] Since the temperature inside human tissue exceeds 43-45°C, it will cause cell death, which is mainly attributed to the damage to cell structure and protein denaturation caused by the temperature increase. This example uses the Arrhenius first-order kinetic model to evaluate the tissue thermal damage caused by the temperature increase during tissue electroporation and ablation. The dimensionless thermal damage function Ω(t) is defined to characterize the degree of thermal damage to the tissue under the action of the pulsed electric field. ,in, is the activation energy; is the gas constant and is a fixed constant. When the thermal damage function Ω(t) ≤ 0.53, it indicates that the tissue has not been thermally damaged; when 0.53 ≤ Ω(t) ≤ 1, it indicates that the tissue has been reversibly thermally damaged; and when Ω(t) > 1, it indicates that the tissue has been irreversibly thermally damaged.

[0041] The electrode needle boundary, tissue boundary, current boundary and thermal damage boundary of the three-dimensional simulation geometric model can be determined by the above method.

[0042] After determining the tissue material parameters and boundary conditions of the 3D simulation geometry model, the model was input into COMSOL software to generate a multi-physics tumor ablation simulation model capable of performing ablation simulation. The initial needle placement plan was input into the multi-physics tumor ablation simulation model for ablation simulation, generating corresponding electric field energy distribution data.

[0043] 4. Needle placement plan adjustment module. This module is used to determine the comparison result between the electric field energy distribution data and the target area. If the comparison result meets expectations, the needle placement plan is determined as the target needle placement plan. If the comparison result does not meet expectations, the needle placement plan is adjusted and updated electric field energy distribution data is determined based on the adjusted needle placement plan. When the comparison result between the updated electric field energy distribution data and the target area meets expectations, the adjusted needle placement plan is determined as the target needle placement plan.

[0044] By comparing the electric field energy distribution data with the target area and performing visualization processing, the ablation range can be obtained, such as Figure 5 As shown, it can be seen that the overlap between the two is not high, indicating that the needle layout scheme needs to be adjusted. Based on the comparison results, the specifications, quantity, spatial arrangement scheme of the electrode needles and the energy release between the electrode needles can be adjusted. For example, when the left side of the tumor is not covered, you can choose to adjust the position of the needle to the left, or when the position cannot be adjusted, you can choose to increase the energy to expand the range. The adjusted needle layout scheme is then input into the multi-physics field tumor ablation simulation model for simulated ablation. When it is not covered, it is adjusted again until the entire tumor is completely covered by the electric field energy. The adjusted ablation range is shown as follows. Figure 6 As shown in FIG, it can be seen that the tumor area has been completely covered by the ablation range. At this time, the adjusted needle arrangement plan is determined as the target needle arrangement plan for the tumor.

[0045] So far, completed Figure 1 Description of each module shown.

[0046] In an embodiment of the present application, a three-dimensional simulation geometric model of the target tissue is first constructed according to the geometric information of the target area of the target tissue using COMSOL multi-physics simulation software, and then the initial needle layout scheme is determined using the three-dimensional simulation geometric model. The material parameters and boundary conditions of the three-dimensional simulation geometric model are then determined to construct a multi-physics tumor ablation simulation model, and then the electric field energy distribution data of the target tissue during tumor ablation surgery is determined using the simulation model according to the initial needle layout scheme. Finally, the electric field energy distribution data is compared with the above-mentioned target area, and the needle layout scheme is adjusted according to the comparison results, thereby determining the target needle layout scheme with the highest degree of overlap between the target area and the electric field energy distribution data. The present application combines the electric field and thermal field in constructing the multi-physics tumor ablation simulation model, and sets the corresponding conductivity when setting the model material parameters. This solves the problems of lack of multi-physics coupling analysis and insufficient dynamic parameter optimization in the prior art, and improves the accuracy of the needle layout scheme.

[0047] like Figure 7 As shown, the present application also provides a device for generating a needle placement plan before tumor ablation surgery based on COMSOL, the device comprising: A three-dimensional model construction unit 701 is used to determine geometric information of a target area of a target tissue, and to construct a three-dimensional simulation geometric model of the target tissue based on the geometric information using COMSOL multi-physics field simulation software, wherein the three-dimensional simulation geometric model includes dimensional parameters and geometric features of the target tissue; A needle placement plan determining unit 702 is configured to determine a needle placement plan for tumor ablation surgery on the target tissue based on the three-dimensional simulated geometric model, wherein the needle placement plan includes specifications, quantity, spatial arrangement plan of electrode needles, and energy release between electrode needles; a multi-physics coupling modeling unit 703, configured to determine tissue material parameters and boundary conditions of the three-dimensional simulation geometric model, construct a multi-physics tumor ablation simulation model based on the tissue material parameters and boundary conditions, and simulate and determine electric field energy distribution data of the target tissue during tumor ablation surgery using the multi-physics tumor ablation simulation model; The needle layout plan adjustment unit 704 is used to determine the comparison result of the electric field energy distribution data and the target area. If the comparison result meets the expectation, the needle layout plan is determined as the target needle layout plan. If the comparison result does not meet the expectation, the needle layout plan is adjusted, and the updated electric field energy distribution data is determined according to the adjusted needle layout plan. Until the comparison result of the updated electric field energy distribution data and the target area meets the expectation, the adjusted needle layout plan is determined as the target needle layout plan.

[0048] In another embodiment, determining the tissue material parameters of the three-dimensional simulation geometric model in the multi-physics field coupling modeling unit includes: The tissue types of the tumor area and the target area in the target tissue, as well as the materials of the conductive part and the insulating part of the electrode needle are set, and the conductivity of each tissue type and material is set.

[0049] In another embodiment, the boundary conditions in the multi-physics field coupling modeling unit include an electrode needle boundary, a tissue boundary, a current boundary, and a thermal damage boundary; The electrode needle boundary is used to determine the electrode potential on the surface of the electrode needle, the tissue boundary is used to determine the insulation state of the outer boundary of the target tissue, the current boundary is used to determine the potential inside the target tissue, and the thermal damage boundary is used to determine the thermal damage state of the target tissue.

[0050] In another embodiment, setting the conductivity of the tissue type of the tumor region in the multi-physics field coupling modeling unit includes: The Heaviside function is used to characterize the dynamic change of the conductivity of the tissue type in the tumor area with the electric field strength, and the conductivity is determined to be , where A is the growth factor of the conductivity after tissue electroporation, is the middle point of the transition area, is half of the transition area, is the temperature growth factor, is the initial ambient temperature.

[0051] In another embodiment, the method for determining the thermal damage boundary in the multi-physics field coupling modeling unit includes: The Arrhenius first-order kinetic model was used to evaluate the thermal damage to tissues caused by temperature increase during tissue electroporation ablation. Determine a dimensionless thermal damage function Ω(t), and then determine the thermal damage state of the target tissue using the thermal damage function Ω(t) and a preset threshold, wherein: is the function of tissue temperature and time, is the activation energy, is the gas constant.

[0052] The above-mentioned embodiment of the present invention provides a system for generating a needle placement plan before tumor ablation surgery based on COMSOL, and based on this system, provides a device for generating a needle placement plan before tumor ablation surgery based on COMSOL. Through the above-mentioned system and device, the problems of lack of multi-physical field coupling analysis and insufficient dynamic parameter optimization in the existing preoperative needle placement plan for tumor ablation surgery are solved, and the accuracy of the needle placement plan is improved.

[0053] This embodiment also discloses a computer device, such as Figure 8 As shown, the computer device includes a processor and a memory, wherein the memory stores at least one instruction, and the at least one instruction is loaded and executed by the processor to implement any of the above-mentioned methods on the COMSOL-based tumor ablation pre-needle plan generation system.

[0054] In addition, in the above-mentioned example of the COMSOL-based device for generating a needle layout plan for tumor ablation surgery, the logical division of each program module is only an example. In actual applications, the above-mentioned functions can be assigned to different program modules as needed, for example, for the convenience of corresponding hardware configuration requirements or software implementation. That is, the internal structure of the COMSOL-based device for generating a needle layout plan for tumor ablation surgery is divided into different program modules to complete all or part of the functions described above.

[0055] The above description is only a preferred embodiment of the present application and is not intended to limit the present application. Any modifications, equivalent replacements, improvements, etc. made within the spirit and principles of the present application shall be included in the scope of protection of the present application.

Claims

1. A system for generating needle placement plans before tumor ablation surgery based on COMSOL, characterized by: The system includes a three-dimensional model building module, a needle arrangement plan determination module, a multi-physics field coupling modeling module and a needle arrangement plan adjustment module; The three-dimensional model construction module is used to determine the geometric information of the target area of the target tissue, and construct a three-dimensional simulation geometric model of the target tissue based on the geometric information using COMSOL multi-physics field simulation software, wherein the three-dimensional simulation geometric model includes the size parameters and geometric features of the target tissue; The needle arrangement plan determination module is used to determine the needle arrangement plan for the tumor ablation surgery of the target tissue according to the three-dimensional simulation geometric model, wherein the needle arrangement plan includes the specifications, quantity, spatial arrangement plan of the electrode needles, and the energy release size between the electrode needles; The multi-physics field coupling modeling module is used to determine the tissue material parameters and boundary conditions of the three-dimensional simulation geometric model, construct a multi-physics field tumor ablation simulation model based on the tissue material parameters and boundary conditions, and simulate and determine the electric field energy distribution data of the target tissue during tumor ablation surgery through the multi-physics field tumor ablation simulation model; The needle layout plan adjustment module is used to determine the comparison results of the electric field energy distribution data and the target area. If the comparison result meets expectations, the needle layout plan is determined as the target needle layout plan. If the comparison result does not meet expectations, the needle layout plan is adjusted, and the updated electric field energy distribution data is determined based on the adjusted needle layout plan. Until the comparison result between the updated electric field energy distribution data and the target area meets expectations, the adjusted needle layout plan is determined as the target needle layout plan.

2. The system according to claim 1, wherein: The determining of tissue material parameters of the three-dimensional simulation geometric model further comprises: The tissue types of the tumor area and the target area in the target tissue, as well as the materials of the conductive part and the insulating part of the electrode needle are set, and the conductivity of each tissue type and material is set.

3. The system according to claim 1, wherein: The boundary conditions include electrode needle boundary, tissue boundary, current boundary and thermal damage boundary; The electrode needle boundary is used to determine the electrode potential on the surface of the electrode needle, the tissue boundary is used to determine the insulation state of the outer boundary of the target tissue, the current boundary is used to determine the potential inside the target tissue, and the thermal damage boundary is used to determine the thermal damage state of the target tissue.

4. The system according to claim 2, wherein: Setting the conductivity of the tissue type of the tumor region includes: The Heaviside function is used to characterize the dynamic change of the conductivity of the tissue type in the tumor area with the electric field strength, and the conductivity is determined to be , where A is the growth factor of the conductivity after tissue electroporation, is the middle point of the transition area, is half of the transition area, is the temperature growth factor, is the initial ambient temperature.

5. The system according to claim 3, wherein: The method for determining the thermal damage boundary includes: The Arrhenius first-order kinetic model was used to evaluate the thermal damage to tissues caused by temperature increase during tissue electroporation ablation. Determine a dimensionless thermal damage function Ω(t), and then determine the thermal damage state of the target tissue using the thermal damage function Ω(t) and a preset threshold, wherein: is the function of tissue temperature and time, is the activation energy, is the gas constant.

6. A device for generating a needle placement plan before tumor ablation surgery based on COMSOL, characterized in that: The device comprises: a three-dimensional model construction unit, configured to determine geometric information of a target area of a target tissue, and to construct a three-dimensional simulation geometric model of the target tissue based on the geometric information using COMSOL multi-physics simulation software, wherein the three-dimensional simulation geometric model includes dimensional parameters and geometric features of the target tissue; a needle arrangement plan determining unit, configured to determine a needle arrangement plan for tumor ablation surgery on the target tissue based on the three-dimensional simulation geometric model, wherein the needle arrangement plan includes specifications, quantity, spatial arrangement plan of electrode needles, and energy release between electrode needles; a multi-physics field coupling modeling unit, configured to determine tissue material parameters and boundary conditions of the three-dimensional simulation geometric model, construct a multi-physics field tumor ablation simulation model based on the tissue material parameters and boundary conditions, and simulate and determine electric field energy distribution data of the target tissue during tumor ablation surgery using the multi-physics field tumor ablation simulation model; The needle layout plan adjustment unit is used to determine the comparison result of the electric field energy distribution data and the target area. If the comparison result meets the expectation, the needle layout plan is determined as the target needle layout plan. If the comparison result does not meet the expectation, the needle layout plan is adjusted, and the updated electric field energy distribution data is determined according to the adjusted needle layout plan. Until the comparison result of the updated electric field energy distribution data and the target area meets the expectation, the adjusted needle layout plan is determined as the target needle layout plan.

7. The device according to claim 6, characterized in that Determining the tissue material parameters of the three-dimensional simulation geometric model in the multi-physics field coupling modeling unit includes: The tissue types of the tumor area and the target area in the target tissue, as well as the materials of the conductive part and the insulating part of the electrode needle are set, and the conductivity of each tissue type and material is set.

8. The device according to claim 6, characterized in that The boundary conditions in the multi-physics field coupling modeling unit include electrode needle boundary, tissue boundary, current boundary and thermal damage boundary; The electrode needle boundary is used to determine the electrode potential on the surface of the electrode needle, the tissue boundary is used to determine the insulation state of the outer boundary of the target tissue, the current boundary is used to determine the potential inside the target tissue, and the thermal damage boundary is used to determine the thermal damage state of the target tissue.

9. The device according to claim 7, characterized in that Setting the conductivity of the tissue type of the tumor area in the multi-physics field coupling modeling unit includes: The Heaviside function is used to characterize the dynamic change of the conductivity of the tissue type in the tumor area with the electric field strength, and the conductivity is determined to be , where A is the growth factor of the conductivity after tissue electroporation, is the middle point of the transition area, is half of the transition area, is the temperature growth factor, is the initial ambient temperature.

10. The device according to claim 8, characterized in that The method for determining the thermal damage boundary in the multi-physics field coupling modeling unit includes: The Arrhenius first-order kinetic model was used to evaluate the thermal damage to tissues caused by temperature increase during tissue electroporation ablation. Determine a dimensionless thermal damage function Ω(t), and then determine the thermal damage state of the target tissue using the thermal damage function Ω(t) and a preset threshold, wherein: is the function of tissue temperature and time, is the activation energy, is the gas constant.

Images