The present invention generally relates to medical systems, methods, and software. More specifically, the present invention provides adjustable open loop control systems, methods, and software for selectively heating tissues, particularly for the noninvasive treatment of urinary incontinence.
Urinary incontinence arises in both men and women with varying degrees of severity, and from different causes. In men, the condition frequently occurs as a result of prostatectomies which result in mechanical damage to the urinary sphincter. In women, the condition typically arises after pregnancy when musculoskeletal damage has occurred as a result of inelastic stretching of the structures supporting the genitourinary tract. Specifically, pregnancy can result in inelastic stretching of the pelvic floor, the external sphincter, and the tissue structures which support the bladder, urethra, and bladder neck region. In each of these cases, urinary leakage typically occurs when a patient's abdominal pressure increases as a result of stress, e.g., coughing, sneezing, laughing, exercise, or the like.
Treatment of urinary incontinence can take a variety of forms. Most simply, the patient can wear absorptive devices or clothing, which is often sufficient for minor leakage events. Alternatively or additionally, patients may undertake exercises intended to strengthen the muscles in the pelvic region, or may attempt a behavior modification intended to reduce the incidence of urinary leakage.
In cases where such non-interventional approaches are inadequate or unacceptable, the patient may undergo surgery to correct the problem. A wide variety of procedures have been developed to correct urinary incontinence in women. Several of these procedures are specifically intended to support the bladder neck region. For example, sutures, straps or other artificial structures are often looped around the bladder neck and affixed to the pelvis, the endopelvic fascia, the ligaments which support the bladder, or the like. Other procedures involve surgical injections of bulking agents, inflatable balloons, or other elements to mechanically support the bladder neck.
In work done related to the present invention, it has been proposed to treat urinary incontinence by selectively remodeling a portion of the pelvic support tissue, often so as to reposition the bladder and/or urogenital tract. U.S. Pat. No. 6,091,995 generally describes laparoscopic and other minimally invasive devices, methods, and systems for shrinking tissues, particularly for treatment of incontinence. U.S. Pat. Nos. 6,216,704; 6,558,381; and 6,546,934, describe noninvasive devices, methods, and systems for shrinking of tissues, often by cooling a surface of an intermediate tissue and directing energy through the cooled intermediate tissue to the target tissue so as to effect shrinkage. U.S. Pat. Nos. 6,156,060; 6,572,639; and 6,776,779, are directed to static devices and methods to shrink tissues for incontinence. Finally, U.S. Pat. No. 6,292,700 describes an endopelvic fascia treatment for incontinence in which a strength of a collagenous tissue increases, optionally without collagenous tissue contraction. U.S. patent application Ser. No. 10/759,732, filed Jan. 15, 2004, describes non-surgical incontinence treatment systems and methods. Each of these patents is assigned to the assignee of the present application, and their full disclosures are incorporated herein by reference.
While these recent proposals for treatment of incontinence represent significant advancements in the art, treatment of incontinence and other conditions related to insufficient collagenous tissue support could benefit from still further advances. For example, temperature sensing mechanisms such as tissue penetrating needles for feedback control may lead to burns on non-target healthy tissues. Temperature sensing needles may also not effect complete heating of target tissue due to a “tenting” effect caused by trapped air and fluid pockets which act to reduce thermal conductivity. For these reasons, it would be desirable to provide improved adjustable open loop control systems, methods, and software for selectively heating support tissues of the body. It would further be desirable if these improved systems and methods provide for truly noninvasive therapy for these support tissues, especially for the treatment of urinary incontinence in men and women. It would be still further desirable if these improved systems and methods provide a good ratio of both tissue treatment efficacy and safety while being less complex and costly to manufacture.
The present invention provides improved adjustable open loop power control systems, methods, and software for selectively heating fascia, tendons, and other support tissues of the body to a desired temperature range. In particular, the systems, methods, and software of the present invention control the delivery of a therapeutic energy that can heat and strengthen a collagenous structural support tissue within a pelvic support system. Advantageously, methods and systems of the present invention eliminate reliance on temperature sensors or tissue penetrating needles for control feedback, and as such provide a truly noninvasive therapy for support tissues, especially for the treatment of urinary incontinence in men and women. Such noninvasive systems are further simpler, more reliable and less costly to manufacture. It will further be appreciated that the present invention is not limited to incontinence therapy, but may also be applied to a variety of conditions such as bladder neck descent, hernias, cosmetic surgery, and the like. As discussed in more detail below, the present invention provides methods, systems, and computer implemented open loop power algorithms that yield enhanced efficacy through improved tissue treatment volumes while maintaining sufficient safety zones and minimizing complications, such as needle burns.
In one aspect of the present invention, a method for therapeutically heating a collagenous structural support tissue of a pelvic support system to a desired temperature range is provided. The method comprises delivering energy to the structural support tissue to heat the tissue to the desired temperature range by ramping up a power level for a first period of time. A first constant high power level is then maintained for a second period of time. The power level is then ramped down for a third period of time. A second constant lower power level is then maintained for a fourth period of time. This power application treatment yields favorable heat treatment temperatures maximizing predictability and efficacy while maintaining sufficient levels of safety.
A ramping up of the power level for the first period of time may comprise ramping up an initial starting power level of no greater than 22 watts, preferably no greater than 16 watts at a slope of no greater than 0.5 watts per second, preferably no greater than 0.25 watts per second. The first period of time may be in a range from 50 seconds to 220 seconds. The first constant high power dwell may be in a range from 34 watts to 40 watts, preferably no greater than 38 watts and the second period of time may be in a range from 60 seconds to 200 seconds. Ramping down of the power level for the third period of time may comprise ramping down the power level to a range from 29 watts to 33 watts at a slope in a range from 0.5 watts per second to 20 watts per second. The third transition period of time may be in a range from 1 second to 10 seconds, typically less than 3 seconds. The second constant low power dwell may be in a range from 29 watts to 33 watts, preferably 30 watts and the fourth period of time may be in a range from 15 seconds to 120 seconds.
Such open loop power methods result in heating the structural support tissue to the desired temperature range between 54° C. and 76° C. with improved predictability. The energy delivery patterns produce a mean minimum safety zone thickness in an intermediate tissue of at least 0.3 mm, preferably at least to 0.5 mm. The energy delivery patterns further produce a mean predominant safety zone thickness in an intermediate tissue of at least 0.5 mm, preferably at least 1.0 mm. The energy delivery patterns also provide enhanced efficacy by producing a tissue treatment volume in a range from 1 cubic centimeters to 5 cubic centimeters. An effective thermal capacity of the tissue treatment volume, denoted by capital letter Q herein, may be in a range from 40 joules/° C. to 87 joules/° C. A coefficient of thermal conductivity between a measured point in the tissue treatment volume and a non-treated tissue, denoted by the capital letter D herein, is in a range from 0.39 watts/° C. to 1.19 watts/° C. A coefficient of thermal conductivity between a measured point in the tissue treatment volume and an applicator body, denoted by the capital letter K herein, is in a range from 0.2 watts/° C. to 0.35 watts/° C.
The energy preferably comprises radio frequency energy, however other forms of heating energy may be adapted to the principles of the present invention, such as electro-resistive, sound, infra-red, radiation, and like energies which may be projected into a subsurface body of the tissue. In some embodiments, the structural support tissue may be cooled by conductive surface cooling. In such instances, a cooled electrode applicator may deliver at much higher power levels than a non-cooled electrode applicator since the tissue heating effect is the net of heating power less the heat removed by cooling. The energy may be delivered so as to effect shrinkage of the structural support tissue and/or to cause bulking and buttressing of the structural support tissue during healing. Tissue strengthening via shrinkage or tissue bulking/buttressing inhibit urinary incontinence or bladder neck descent, wherein the structural support tissue may comprise a collegenated tissue in an endopelvic fascia and the intermediate tissue may comprise vaginal mucosa. The structural support tissue may be accessed transvaginally or laparoscopically.
In another aspect of the present invention, a system for therapeutically heating a collagenous structural support tissue of a pelvic support system to a desired temperature range is provided. The system comprises an applicator body and a processor coupleable to the applicator body. The processor may be programmed to deliver energy to the structural support tissue with the applicator body by ramping up a power level for a first period of time, maintaining a high power dwell for a second period of time, ramping down the power level for a third period of time, and maintaining a low power dwell for a fourth period of time. The system may further comprise a power supply coupleable to the processor as well as a cooling source coupleable to the processor.
In yet another aspect of the present invention, a computer-readable storage medium having a computer-readable program embodied therein for directing operation of a computer system is provided. The computer system including a communications system, a processor, and a memory device. The computer-readable program includes instructions for therapeutically heating a collagenous structural support tissue of a pelvic support system to a desired temperature range in accordance with the any of the method steps described herein.
In still another aspect of the present invention, a method and device for heating living human tissue to a prescribed temperature range is provided. This is accomplished by application of heating energy in a particular pattern (e.g., power level versus time) such that the inherent ability of the specific tissue to absorb and dissipate heat interacts with the specific applied power pattern to yield the prescribed temperature range. As such, the need to invasively measure the tissue temperature and employ feedback control is thereby circumvented.
A further understanding of the nature and advantages of the present invention will become apparent by reference to the remaining portions of the specification and drawings.
The following drawings should be read with reference to the detailed description. Like numbers in different drawings refer to like elements. The drawings, which are not necessarily to scale, illustratively depict embodiments of the present invention and are not intended to limit the scope of the invention.
The present invention provides methods, systems, and software algorithms for controlling delivery of energy to a body's support tissue to enhance the structural support provided by the body's support tissues. The present invention may be directed to inducing controlled stiffening, contraction, or shrinkage of the structural support tissue of the body, typically being a collagenous tissue such as fascia, ligament, or the like.
For example, in one specific use, the present invention provides for treatment of urinary incontinence. The structural support tissue will be part of a pelvic support system that is responsible in some manner for control of urination, or for supporting such a tissue. The tissues of the pelvic support system generally maintain the position of the genitourinary tract, and particularly the position of urinary bladder, urethra, and the bladder neck coupling these structures. In general, endopelvic fascia may define a hammock-like structure which extends laterally between the left and right arcus tendineus fasciae pelvis (ATFP). These tendon structures may extend substantially between the anterior and posterior portions of the pelvis, so that the endopelvic fascia EF at least partially defines the pelvic floor.
The fascial tissue of the pelvic support system may comprise tissues referred to under different names by surgeons of different disciplines, and possibly even by different practitioners within a specialty. In fact, some surgeons may assign a collagenous support structure of the endopelvic fascia one name when viewed from a superior approach, and a different name when viewed from an inferior approach. Some of the endopelvic fascia may comprise two collagenous layers with a thin muscular layer therebetween, or may comprise a single collagenous layer. The hammock-like endopelvic fascia described herein may be damaged or missing, particularly after pregnancy, so that the support of the genitourinary tract is instead provided by a variety of fascial layers, muscular tissues, ligaments, and/or tendons within the pelvis. Hence, the treatment of the present invention may be directed at a variety of tissue structures defining the pelvic floor and/or diaphragm (including: anterior sacro-coccygeal ligament; arcus tendineus fasciae pelvis ATFP, the white line of the pelvis; fasciae of the obturator intemus muscle; the arcus tendineus levator ani or “picket fence” to the iliococcygeus portion of the levator ani muscle; bulbocavernosus muscle; ischiocavemosus muscle; urethrovaginal sphincter; m. compressor urethrae muscle; and m. sphincter urethrovaginal muscle which replaces deep perineal muscle); structures of the bladder and urethra (including: urethrovesical fascia; detrusor muscle; and the pubococcygeus muscle which relaxes to open the bladder neck, initiating micturation); structures of the vagina (including: vagino-uterine fascia, lamina propria—the dense connective tissue layer just under the epithelium; pubo-urethral or puboprostatic ligaments; pubo-vesicle ligament and posterior pubo-urethral or puboprostatic ligament; pubovesicle muscle, a smooth muscle that is integrated with the pubovesicle ligament; and pubocervical fascia which attaches to the ATFP); structures of the uterus (including: round ligament; sacrouterine ligament; and broad ligament); and structures of the bowel (including: rectal fascia and mackenrodt's ligament).
When the endopelvic fascia has excessive length or stretches excessively under a load, the fluid pressure within the bladder advances into the bladder neck and down the urethra more readily. Leakage may result in part because the endopelvic fascia allows the bladder, bladder neck, and/or urethra to drop below its desired position, at which fluid pressure within the bladder may actually help to seal the bladder neck. Stretching of the endopelvic fascia may also alter the timing of pressure pulse transmission to the urethra.
When a continent woman coughs, the pressure in the urethra will often increase more than one-tenth of a second prior to the increase in bladder pressure. In women with stress incontinence, the bladder pressure may rise first. For a continent woman having endopelvic fascia which stretches much less under the influence of a pressure pulse, the time delay between initiation of the pressure pulse and transferring sufficient force to urethra to effect closure may therefore be significantly less. By treating the endopelvic fascia to decrease its length and/or increase its stiffness, the descent time of the pelvic viscera during a cough will be shorter than an untreated, excessively long and/or excessively elastic tissue.
The support tissue may be treated non-surgically or it may be accessed for direct treatment in a variety of ways. When using a multi-electrode applicator, for example, the surface of the endopelvic fascia (or other tissue) may be accessed transvaginally by forming and displacing a flap from the vaginal wall with the assistance of a weighted speculum. Alternatively, the endopelvic fascia may be accessed laparoscopically. When using the noninvasive cooled electrode applicator, the tissue may be accessed directly by placing the applicator on the anterior vaginal wall.
Tissue contraction or stiffening results from controlled heating of the tissue by affecting the collagen molecules of the tissue. Contraction occurs as a result of heat-induced uncoiling and repositioning of the collagen β-pleated structure. By maintaining the times and temperatures set forth below, significant tissue contraction may be achieved. Tissue strengthening by controlled contraction or shrinkage is described in more detail in U.S. Pat. No. 6,836,688, assigned to the assignee of the present application and incorporated herein by reference. Stiffening results from the loss of elasticity of the tissue due to the formation of scar tissue and/or attachment of adjacent support tissues to each other as a result of controlled heating of the tissue. Tissue strengthening by bulking and buttressing of the structural support tissue due to the healing process and/or formation of scar tissue is described in more detail in U.S. Pat. No. 6,292,700, assigned to the assignee of the present application and incorporated herein by reference.
While the remaining description is generally directed to a system for treatment of urinary stress incontinence of a female patient, it will be appreciated that the present invention will find many other applications for selectively directing therapeutic heating energy into the tissues of a patient body. For example, treatment of other conditions may be effected by selective ablation, shrinking or stiffening of a wide variety of other tissues, including (but not limited to) the diaphragm, esophagus, the nasal concha, the abdominal wall, the breast supporting ligaments, the fascia and ligaments of the joints, the collagenous tissues of the skin, tumors, and the like.
While three active electrode segments are illustrated in
Control unit 20 includes a switch 36 that serves to activate and deactivate transmission of energy from a power source 38, such as a bipolar radiofrequency (RF) power source, to electrodes 12 on applicator 22. Switch 36 may be activated with activation and deactivation of the input device 49 on applicator 22, or the like. The power source 38 is coupled to the processor 32. In the three electrode configuration of
A cooling assembly 44 may optionally be coupled to processor 32 and applicator 22 and will be configured to pre-cool the tissue contacted by applicator 22 and/or cool the tissue during the delivery of the energy. A more complete description of some examples of cooling assembly 44 are described in commonly owned U.S. Pat. Nos. 6,091,995 and 6,480,746, the complete disclosures of which are incorporated herein by reference. As can be appreciated, cooling assembly 44 is optional and not all applicators of the present invention include cooling assembly 44.
Processor 32 may identify and display appropriate error messages pertaining to a variety of conditions, such as errors encountered during the diagnostic system tests, and the like. Some embodiments of processor 32 allow the user to set date and time, audio tone level, language selection for display on display device 28, power levels, desired treatment times, desired temperature goals, desired safety zone thickness, desired treatment volume, and the like. In some embodiments, such parameters may be preset. Processor 32 generates audio tones to prompt the user for actions and to indicate error and out of range conditions. A continuous or intermittent audio tone may be emitted by a speaker 26 associated with processor 32 at a steady rate when energy is applied. Processor 32 may generate a welcome screen showing a logo or other graphics desired by user of system 10. Processor 32 may display recoverable error condition messages and prompts the user to correct the cause. Unrecoverable error messages may be displayed on display device 28 and give appropriate error information.
Control unit 20 may be configured to complete a self-test each time the power source 38 is turned on. Control unit 20 allows processor 32 to complete its internal tests and display error messages accordingly. A fault in the power source output test can be diagnosed and displayed as an error condition. Processor 32 may be programmed to provide a clock signal for hardware detection of software operation. Processor 32 performs tests of internal subsystems, including but not limited to the analog and digital electronics. Control unit 20 provides a special test, diagnostics and service mode, which will allow the manufacturer or servicer of system 10 to bypass the normal diagnostic self-tests, be able to manually execute all functions and perform calibration and setup. This mode is generally not be accessible to the user.
At step 102, energy is delivered to the structural support tissue comprising collegenated tissue in an endopelvic fascia to heat the tissue to a desired temperature range by starting at a low power application level and slowly ramping up for a first period of time to a peak applied power (first constant power level). This allows the resulting “heat plume” beneath the tissue surface to develop into a characteristic form, rather than having the heat concentrate in a very small volume because it can not dissipate as fast as the energy is being applied. This allows the device to heat the maximum volume of tissue for any applied peak temperature as well as yields a more consistent temperature response curve. A ramping up of the power level for the first period of time may comprise ramping up an initial power level that is no greater than 22 watts, preferably no greater than 16 watts at a slope that is no greater than 0.5 watts per second, preferably no greater than to 0.25 watts per second. The first period of time may be in a range from 50 seconds to 220 seconds. Typically, the inherent conductive cooling rate applied by the device is kept up with the rate of energy application so that the surface safety zone is maintained at a maximum for any given peak temperature. This is because the rate of conductive cooling rises with temperature. As such, it is acceptable to go to higher power levels after temperature is built up.
At step 104, the peak applied power or first constant high power level is then maintained for a second period of time. The peak power level is low enough so that the heating rate does not outrun the inherent conductive cooling rate of the device, in order to maintain the safety zone at the tissue surface. The peak power level is high enough and sustained for long enough so that, together with the ramp period, it achieves the desired temperature range in an acceptable period of time. This peak power level will be above the equilibrium power rate required to sustain the target temperature range. The desired temperature range is determined by the need to heat as much volume as possible to tissue necrosis temperature (50° C.) and to achieve as high a volume of collagen shrinkage as possible, which is a time-temperature dependent effect, while keeping peak temperature below the maximum safe value for all patients. The actual practical peak temperature range may be determined by the variation in tissue thermal characteristics across the patient population, as these characteristics interact with a fixed power algorithm to achieve a range of outcomes. Typically, the first constant high power dwell may be in a range from 34 watts to 40 watts, preferably 35 watts and the second period of time may be in a range from 60 seconds to 200 seconds. Typically at the high power dwell level, an equilibrium temperature is above the desired peak temperature so that the desired temperature may be reached rapidly.
At step 106, once the peak power period (first constant high power level) is completed, the power level is then ramped down from the peak applied power for a third period of time to a level closer to the equilibrium power level required to maintain the target temperature range (second constant lower power level). Ramping down of the power level for the third period of time may comprise ramping down the power level to a range from 29 watts to 33 watts at a slope in a range from 0.5 watts per second to 20 watts per second. The third transition period of time may be in a range from 1 second to 10 seconds, generally less than 3 seconds. At step 108, the second constant lower power level is then maintained for a fourth period of time which is practical to achieve maximum dwell time near the desired peak temperature. The second constant low power dwell as noted above may be in a range from 29 watts to 33 watts, preferably 30 watts and the fourth period of time may be in a range from 15 seconds to 120 seconds. Due to the variability of tissue response, the range chosen for the second constant lower power level may cause some patient treatments to slightly rise in temperature while others may slightly fall in temperature and while others may remain constant. Typically at the low power dwell level, the equilibrium temperature is closer to the desired peak temperature.
The open loop control system of the present invention accounts for a variety of multi-variable components so as to achieve the desired heat treatment. For example, lower power levels, gentle power ramps, and lower maximum tissue temperatures have been found to increase safety zone thickness of intermediate tissue, such as the vaginal mucosa. On the other hand, higher power levels, increased time at a given power level, or higher maximum tissue temperatures result in rising tissue treatment volumes of the endopelvic fascia. Energy delivery will also depend in part on which tissue structure is being treated, how much tissue is disposed between the target tissue and the electrode, and the ability of the tissue to accept and store power. The power levels used in the present invention will also vary depending on the electrode size, electrode spacing, and whether or not cooling is used. For example, a cooled electrode applicator may deliver at much higher power levels than a non-cooled electrode applicator since the tissue heating effect is the net of heating power less the heat removed by cooling.
Hence, power levels, desired treatment times, desired temperature goals, desired safety zone thickness, desired treatment volume, a patient's anatomy, and applicator configurations are factors to be considered so as to improve efficacy while maintaining sufficient levels of safety. Generally, the energy delivery patterns produce a mean minimum safety zone thickness in an intermediate tissue of at least 0.3 mm, preferably at least 0.5 mm. The energy delivery patterns further produce a mean predominant safety zone thickness in an intermediate tissue of at least 0.5 mm, preferably at least 1.0 mm. The energy delivery patterns also provide enhanced efficacy by producing a tissue treatment volume in endopelvic fascia in a range from 1 cubic centimeters to 5 cubic centimeters.
The theory of open loop RF power control is now described with reference to
In these studies, heat losses due to the surroundings are assumed to be negligible. Therefore all heat flow is assumed to occur between the applicator head and the liver sample. Heat is applied through radio-frequency resistance heating. Cooling is by simple conduction due to the temperature difference between the tissue and the applicator head. For a given steady state RF power level and applicator head temperature, if the process is allowed to go on long enough, a characteristic equilibrium temperature will be reached. As the tissue is heated, cooling watts (LA) increase with tissue temperature, while heating rate (H) remains constant, until cooling rate equals heating rate and the equilibrium temperature is reached. This may be represented by the following equations:
LA=K(Ts−Th)
at equilibrium H=LA
H=LA=K(Tse−Th) or Tse=(H/K)+Th (Equation #1)
H and Th are set parameters, controlled by the applicator. K can be calculated from experimental data derived by running a heating cycle to equilibrium. In a typical test run, a steady state RF power application (H) of 20 watts with an applicator head temperature (Th) of 0° C., yields an equilibrium tissue temperature (Tse) of 77° C. in about 45 minutes. Calculating K from Equation #1 yields 0.26 watts/° C.
The value of K is determined by the nature of the applicator head and its materials as well as the power algorithm. Due to such factors, the value of K may be in a range from 0.2 watts/° C. to 0.35 watts/° C. in human tissues. It is a basic characteristic of the metal-to-tissue interface and is quite constant and similar between liver tissue samples and human tissue. Since the value of K can be arrived at through an equilibrium test (
Non-Equilibrium Factors
When a tissue sample is heated by the applicator 22, the equilibrium temperature will not be reached if one or both of the following conditions apply: (1) not enough time is allowed to elapse or (2) the equilibrium temperature exceeds the temperature at which the sample is altered or destroyed (in the treatment of incontinence, this corresponds to the occurrence of a “broad burn” which occurs around 80° C.). In the treatment of incontinence, the first condition prohibits reaching equilibrium as the lengthy treatment time is unpractical. The second condition should be avoided for safety reasons. Therefore the power algorithm selected after this characterization is complete, will be targeted to achieve a temperature range which tops out below 80° C.
In order to account for the time variation of temperature due to power input, we introduce another parameter Q which represents a thermal capacity of tissue volume being heated (joules/° C.). Q is a measure of the ability of the tissue to store heat. Conversely, Q together with the power input level, determines the rate of temperature rise. The basic geometry of the applicator head and the inherent characteristics of bipolar RF heating, cause the volume of tissue being heated to take on a very consistent characteristic geometry. This may be characterized by determining the 50° C. boundary. This is an elliptical cylindrical surface, where the tissue temperature is at 50° C. by the end of treatment. The 50° C. heat plume forms quickly during treatment and stabilizes, with further heating going into temperature rise within it. Outside the 50° C. heat plume, temperature drops off in a consistent steep gradient pattern. The 50° C. value is chosen because it is the approximate temperature at which tissue necrosis occurs.
Derivation of temperature versus time for constant RF power application will now be described. When operating the applicator on a liver tissue sample at a constant RF power level H, where heat losses are again assumed to be negligible, the instantaneous net heating which occurs is given by the following equations:
Hnet=H−LA or Hnet=H−K(Ts−Th)
The instantaneous rate of temperature rise is given by the following equations:
dTs/dt=Hnet/Q or dTs/dt=(H−K(Ts−Th))/Q
Solution of this differential equation is given in the form below:
Ts=Ce−(K/Q)t+((H+KTh)/K)(1−e−(K/Q)t)
By inspection of boundary conditions, C=To (the starting temperature of tissue (° C.)) and solving the equations for Ts, the tissue temperature as a function of time may be derived for constant power by the following equation:
Ts=T0e−(K/Q)t+((H+KTh)/K)(1−e−(K/Q)t) (Equation #2)
It will be appreciated that the for large values of time t, Equation #2 converges on Equation #1 as would be expected. If the Q value of a tissue sample is known, the curve of temperature versus time as the sample is heated by the applicator head at a constant RF power rate may be predicted by Equation #2.
To determine Q for the liver samples, a test heating at very low RF power may be performed, with the cooling system turned off. In this condition, the temperature of this sample rises in a straight sloped line, the slope of which is given by the following equations:
dTs/dt=H/Q or Q=H/slope
Knowing the H watts input, with a measurement of the slope of temperature versus time, Q may be in independently calculated. By the technique, it may be deduced that the effective Q for these liver samples is in a range from about 55 joules/° C. to 85 joules/° C. Determination of Q for human subjects varies over a wider range and will be accordingly deduced by other means as discussed below in the section on adapting the equations to human anatomy factors.
Derivation of temperature versus time for ramping RF power application will now be described. Since the actual open loop power algorithm used is segmented, with ramp-up and ramp-down periods as well as periods of level power, the tissue temperature as a function of time should be derived for a ramping power application. For RF power ramping at a slope P (watts/sec), at any given point in time, where Ho is the RF power at the beginning of the ramp, the instantaneous net heating which occurs is given by the following equation:
H=H0+Pt
The instantaneous rate of temperature rise is given by the following equations:
dTs/dt=H/Q or dTs/dt=(H−K(Ts−Th))/Q=((Ho+Pt)−KTs+KTh)/Q
Solution of this differential equation is given in the form below:
Ts=T0e(K/Q)t+((H0+KTh)/K)(1−e(K/Q)t)+(P/K)t−(PQ/K2)(1−e−(K/Q)t)
By inspection of boundary conditions, C=To (the starting temperature of tissue at the beginning of the RF power ramp (° C.)) and solving the equations for Ts, the tissue temperature as a function of time may be derived for a ramping RF power level by the following equation:
Ts=T0e−(K/Q)t+((H0+KTh)/K)(1−e−(K/Q)t)+(P/K)t−(PQ/K2)(1−e−(K/Q)t) (Equation #3)
Temperature as a function of time for any RF power algorithm which is in the form of a series of ramped and level RF power periods may now be expressed for the liver tissue samples. The preferred open loop RF power algorithm 110′ is depicted in
Derivation of temperature versus time for a segmented power curve will now be described. In particular, the temperature versus time curve which results from the application of algorithm 110′ in
TS=T0e−(K/Q)(t-t0)+((H+KTh)/K)(1−e−(K/Q)(t-t0))+(((H1−H0)/(t1−t0))/K)(t−t0)−((H1−H0)/(t1−t0))Q/K2)(1−e−(K/Q)(t-t0))
where H=H0+(H1−H0)(t−t0)/(t1−t0)
The temperature versus time curve which results from the application of algorithm 110′ in
TS=T1e−(K/Q)(t-t1)+((H1+KTh)/K)(1−e−(K/Q)(t-t1))
where H=H1
The temperature versus time curve which results from the application of algorithm 110′ in
Ts=T2e−(K/Q)(t-t2)+((H1+KTh)/K)(1−e−(K/Q)(t-t2))+(((H2−H1)/(t3−t2))/K)(t−t2)−((H2−H1)/(t3−t2))Q/K2)(1−e−(K/Q)(t-t2))
where H=H1+(H2−H1)(t−t2)/(t3−t2)
The temperature versus time curve which results from the application of algorithm 110′ in
Ts=T3e−(K/Q)(t-t3)+((H2+KTh)/K)(1−e−(K/Q)(t-t3))
where H=H2
The above equations are graphed for power levels typically used on the liver samples. The theoretically predicted form of the temperature response curve 112′ (depicted by the letter x) is superimposed on the driving RF open loop power algorithm 110′ (depicted by asterisks *) in
It will be appreciated that the values of H1 and H2 in
Human Anatomy Factors
Since the fat layer 118 varies in thickness and the intimacy of contact between the tissue structures across the space of retzius 120 is also variable, the losses of heat from the treatment zone 116 to the underlying tissues 122 is variable from patient to patient and even between sides of a particular patient. Considering the directions radially outward within the endopelvic fascia 116, the thermal gradients form a consistent pattern which is taken into account in the conception of the thermal capacity Q of the volume being treated. Therefore losses in those directions need not be considered. It is the variable losses in the increasing depth dimension that affect outcomes that should be considered. It will further be appreciated that blood circulation within the endopelvic fascia 116 is very minor, so heat transfer can be treated as conduction.
Transfer of heat from the tissue effects volume 124 to the rest of the body is given by the following equation:
LB=D(Ts−37)
LB represents the rate of loss of thermal energy from the treatment zone to the rest of the body. D represents a coefficient of thermal conductivity between the measured point (Ts) in the tissue effects volume 124 and the rest of the body (watts/° C.). The body temperature is naturally regulated at 37° C. The application of this equation is complicated by the fact that the applicator 22 may cool the treatment zone 116 below body temperature before beginning RF treatment to heat it up. Therefore, until Ts moves up and past body temperature, losses are negative as the body is in effect helping the RF in heating the treatment zone. After Ts passes above body temperature, do then losses become positive. At the instant the treatment temperature passes through body temperature, losses are momentarily zero.
Taking into account the rate of loss of thermal energy from the treatment zone to the body (LB), heating to equilibrium at constant RF power in vivo may now be expressed by the following equations:
LTOTAL=LA+LB=K(Ts−Th)+D(Ts−37)
at equilibrium H=LTOTAL
H=LTOTAL=K(Ts−Th)+D(Ts−37) or Tse=(H+KTh+37D)/(K+D) (Equation #4)
As discussed above, these equilibrium calculations are insufficient to describe power applications for incontinence treatment because the treatment procedure times for the equilibrium temperature to be reached are lengthy and unpractical. Therefore non-equilibrium effects in vivo need to be considered.
Taking into account the rate of loss of thermal energy from the treatment zone to the body (LB), derivation of temperature versus time for constant RF power application in vivo may now be expressed by the following equations:
Hnet=H−LA−LB or Hnet=H−K(Ts−Th)−D(Ts−37)
The instantaneous rate of temperature rise is given by the following equations:
dTs/dt=Hnet/Q or dT/ds=(H−(K+D)Ts+KTh+37D)/Q
Solution of this differential equation and application of the boundary condition to temperature=To yields:
Ts=T0e−((K+D)/Q)t+((H+KTh+37D)/(K+D))(1−e−((K+D)/Q)t) (Equation #5)
Taking into account the rate of loss of thermal energy from the treatment zone to the body (LB), derivation of temperature versus time for ramping RF power application in vivo may now be expressed by the following equation:
H=Ho+Pt
The instantaneous rate of temperature rise is given by the following equations:
dTs/dt=Hnet/Q or dTs/dt=(H0+Pt−(K+D)Ts+KTh+37D)/Q
Solution of this differential equation and application of the boundary condition to temperature=To yields:
Ts=T0e−((K+D)/Q)t+((H0+KTh+37D)/(K+D))(1−e−((K+D)/Q)t)+(P/(K+D))t−(PQ/(K+D)2)(1−e−((K+D)/Q)t) (Equation #6)
Taking into account the rate of loss of thermal energy from the treatment zone to the body (LB), derivation of temperature versus time for a segmented RF power application in vivo may now be expressed by the following heat transfer equation from time t0 to t1 (
Ts=T0e−((K+D)/Q)(t-t0)+((H0+KTh+37D)/(K+D))(1−e−((K+D)/Q)(t-t0))+(P/(K+D))(t−t0)−(PQ/(K+D)2)(1−e−((K+D)/Q)(t-t0))
where H=H0+(H1−H0)(t−t0)/(t1−t0)
The temperature versus time for a segmented RF power application in vivo may now be expressed by the following heat transfer equation from time t1 to t2 (
Ts=T1e−((K+D)/Q)(t-t1)+((H+KTh+37D)/(K+D))(1−e−((K+D)/Q)(t-t1))
where H=H1
The temperature versus time for a segmented RF power application in vivo may now be expressed by the following heat transfer equation from time t2 to t3 (
Ts=T2e−((K+D)/Q)(t-t2)+((H1+KTh+37D)/(K+D))(1−e−((K+D)/Q)(t-t2))+(P/(K+D))(t−t2)−(PQ/(K+D)2)(1−e−((K+D)/Q)(t-t2))
where H=H1+(H2−H1)(t−t2)/(t3−t2)
The temperature versus time curve for a segmented RF power application in vivo may now be expressed by the following heat transfer equation from time t3 to t4 (
Ts=T3e−((K+D)/Q)(t-t3)+((H+KTh+37D)/(K+D))(1−e−((K+D)/Q)(t-t3))
where H=H2
Although the effective thermal capacity Q and coefficient D for the temperature versus time equations for an individual treatment can not be predicted, actual response curves from multiple treatments can be used to develop a profile of the range of these figures. The coefficient K is the least variable, and remains very close to 0.26+/−0.01 watts/° C., as previously calculated. K is a function of the interface between the brass applicator head and the tissue, and varies little from treatment to treatment. Q and D however are highly dependent on anatomical variations among patients and even between sides on the same patient.
Thermal capacity Q would be expected to vary with the thickness of the endopelvic fascia 116, which typically varies from 6 to 9 mm as shown in
dTs/dt=(H−K(Ts−Th))/Q or Q=(H−37K+KTh)/(dTs/dt) (Equation #7)
H and Th are controlled variables. By observing the slope of the temperature line and the instantaneous power being applied as the temperature rises through body temperature, the required information to calculate effective Q for that particular treatment may be obtained.
Thermal transfer coefficient between the treatment zone and the rest of the body D would be expected to vary with the intimacy of contact between the endopelvic fascia 116 and the underlying muscular tissue 122. This varies with the thickness of the fat layer 118 behind the endopelvic fascia 116 and the size of the space of retzius 120. The peak temperature equations described above may be plotted using a plotting program, such that the theoretical curves resulting from the various values of the three coefficients, K, Q, and D may be readily observed. With the values of K set at nominal, and the value of Q determined from an actual observed temperature plot, iterative curve fitting may be utilized to estimate the value of D which would produce the observed in vivo temperature plot. By this means, the values for Q and D may be determined for the open loop power treatments of the present invention, for which temperature was monitored. The derived Q value was in a range from 40 joules/° C. to 87 joules/° C. and the derived D value was in a range from 0.39 watts/° C. to 1.19 watts/° C.
These extreme values have been applied to the heat transfer equations to calculate projected highest and lowest peak temperatures which will be seen with the open loop algorithm of the present invention. This process adds a very large measure of conservatism, because it is assumed that the three most extreme values coincide in the same treatment. Since the coinciding of these extreme values in a single treatment is itself an extremely unlikely occurrence, the calculated values for the temperature extremes are very rare. The inherent mathematics of the heat transfer equations leads to a non-normal distribution of peak temperatures. The distribution is skewed towards the lower temperatures. Therefore, although the average values for the coefficients are known, average peak temperature may not be reliably calculated. Fortunately, only the range is needed to select an algorithm. Since it has been observed with appropriate algorithm characteristics, such as gradual power ramp up and sufficient dwell time, efficacy performance is very forgiving of peak temperature variation. Hence, the objective of this calculation becomes mainly the assurance that the expected range of peak temperatures falls within the allowable range for safety, while allowing the algorithm to run long enough, and while favoring the high end of the temperature range for efficacy. In sum, a tight temperature range is not necessary.
As shown in
Experimental
The following description of experimental studies provides some specific, non-limiting examples that are encompassed by the present invention.
A series of experiments using a Sol2 head to treat 10 mm thick pieces of bovine liver were performed using the open loop power algorithms 110 of the present invention. At the conclusion of each treatment the safety zone thickness 126 (
Experiments were run using two comparison treatment algorithms. A Sol1 power step algorithm was run with its initial 24 second treatment at 20 watts followed by treatment at 41 watts. This value was chosen to achieve a mean treatment time on liver close to that seen in a Sol1 clinical study (109 seconds). The Sol1 power step algorithm treated until tissue temperature reached 70° C. or 150 seconds of treatment had occurred.
The Sol2 open loop power algorithm 110 was used with power levels chosen to achieve mean treatment temperatures of 59° C., 68° C., or 75° C. Each series involved 25 treatments. The Sol2 algorithm starts at 15 watts and rises to full treatment power after 140 seconds. It remains at that power level for 150 seconds and then drops by 5 watts for a 20 second dwell period. The 59° C. mean treatment temperature matches the mean temperature observed in open loop Mexico feasibility patients. These patients received full treatment power for only 60 seconds. The higher liver treatment target of 68° C. was chosen to be slightly larger in order to provide a worst case estimate of the expected safety zones in human patients. The 75° C. treatments were done to provide a worst case lower limit to the safety zone thickness.
As illustrated in
A separate series of experiments were run to measure the K value and the Q value for use in comparing the theoretical curves to the observed tissue temperature curves in liver. The liver samples in all experiments were covered with saran wrap on the untreated side. A melamine thermal insulation box was placed over the tissue to minimize air losses. As discussed in detail above, the following equation:
Tse=(H/K)+Th
is used to calculate K. The equilibrium tissue temperature is measured by performing 30 to 50 minute treatments at low power values.
As discussed in detail above, the following equation is used to calculate Q:
dTs/dt=H/Q or Q=H/slope
The tissue temperature slope is measured using short treatment times with a non-cooled applicator head. The heating power is in the range from 5 to 7 watts which matches the net RF heating rate (35 watt heating less 28 to 30 watts of cooling) during an actual treatment. The Q value was measured by averaging the temperature rise curve in six different runs.
The composite curves for the 25 treatment runs are shown in
Although certain exemplary embodiments and methods have been described in some detail, for clarity of understanding and by way of example, it will be apparent from the foregoing disclosure to those skilled in the art that variations, modifications, changes, and adaptations of such embodiments and methods may be made without departing from the true spirit and scope of the invention. Therefore, the above description should not be taken as limiting the scope of the invention which is defined by the appended claims.