Электрокардиология мәселесі - Forward problem of electrocardiology

Реалистік жүрек соғысын модельдеу.

ЭКГ-да көрсетілгендей, адамның жүрегі үшін қалыпты синус ырғағының схемасы.

The электрокардиологияның алға қойылған проблемасы зерттеу үшін есептеу-математикалық тәсіл болып табылады электрлік белсенділік туралы жүрек дене беті арқылы.^[1] Зерттеудің негізгі мақсаты ан электрокардиограмма (ЭКГ), оны анықтау үшін маңызды клиникалық маңызы бар жүрек патологиялары сияқты ишемия және инфаркт немесе тексеру үшін фармацевтикалық араласу. Олардың маңызды функционалдық мүмкіндіктерін және салыстырмалы түрде аз инвазивтілігін ескере отырып, электрокардиография әдістері клиникалық ретінде жиі қолданылады диагностикалық тесттер. Осылайша, ЭКГ-ны компьютерлік түрде көбейтуге көшу табиғи нәрсе, бұл организмдегі жүректің жүрісін математикалық модельдеуді білдіреді.^[1]

ЭКГ үшін алға бағытталған модельдің үш негізгі бөлігі:

жүректің электрлік белсенділігінің моделі;
экстракардиальды аймақты білдіретін электр қуаты ішіндегі электрлік потенциалдың диффузиясының моделі;
жүрек-денені біріктірудің кейбір ерекше шарттары.^[2]

Осылайша, ЭКГ алу үшін математикалық электрлік жүрек моделін қарастыру керек, диффузиялық модельмен бірге пассивті дирижер ішіндегі электрлік таралуын сипаттайтын торсық.^[1]

Біріктірілген модель әдетте а үш өлшемді модель терминдерімен көрсетілген дербес дифференциалдық теңдеулер. Мұндай модель әдетте көмегімен шешіледі ақырғы элемент әдісі шешімнің кеңістік эволюциясы үшін және жартылай жасырын сандық схемалар тарту ақырғы айырмашылықтар шешімнің уақыт эволюциясы үшін. Алайда, мұндай техниканың есептеу шығындары, әсіресе үш өлшемді модельдеу кезінде айтарлықтай жоғары. Осылайша, көбінесе жеңілдетілген модельдер қарастырылады, мысалы, жүректегі электрлік белсенділікті торсындағы проблемадан тәуелсіз. Нақты нәтижелер беру үшін жүректің және дененің үш өлшемді анатомиялық-шынайы модельдерін қолдану қажет.^[1]

Тағы бір ықтимал жеңілдету - а динамикалық модель үшеуінен жасалған қарапайым дифференциалдық теңдеулер.^[3]

Жүрек тіндерінің модельдері

Жүректің электрлік белсенділігі ағыннан туындайды иондар арқылы жасуша қабығы бойымен қозу толқынын анықтайтын жасушаішілік және жасушадан тыс кеңістіктер арасында жүрек бұлшықеті бұл жүректің қысылуын үйлестіретін және, демек, оны итеруге мүмкіндік беретін жүректің насостық әрекеті қан арқылы қанайналым жүйесі. Жүректің электрлік белсенділігін модельдеу а-да иондар ағынын модельдеуге байланысты микроскопиялық деңгей және қозу толқынының бұлшықет бойымен таралуы бойынша талшықтар үстінде макроскопиялық деңгей.^[1]^[4]

Макроскопиялық деңгейдегі математикалық модель арасында, Виллем Эйнтховен және Август Уоллер проекциясы тұрақты нүктенің айналасында айналатын дипольдің тұжырымдамалық моделі арқылы ЭКГ анықтады қорғасын ось жетекші жазбаларды анықтады. Содан кейін, а екі өлшемді көмегімен жүрек қызметін фронтальды жазықтықта қалпына келтіру мүмкін болды Эйнтховеннің аяқ-қолы алға шығады I, II және III теориялық негіз ретінде.^[5] Кейінірек айналмалы жүрек диполы жеткіліксіз болып саналды және оны алмастырды көпполярлы шектелген торс доменінің ішінде қозғалатын көздер. Осы дереккөздерді сандық бағалау әдістерінің негізгі жетіспеушілігі - олардың жүрек құбылыстарын шынайы имитациялау үшін өте маңызды бөлшектерінің болмауы.^[4]

Екінші жағынан, микроскопиялық модельдер бір жасушалардың мінез-құлқын бейнелеуге және оларды электрлік қасиеттерін ескере отырып қосуға тырысады.^[6]^[7]^[8] Бұл модельдер әр түрлі масштабтарға байланысты кейбір қиындықтарды тудырады, атап айтқанда, әсіресе қайта кіру немесе кең ауқымды құбылыстар үшін дене бетінің әлеуеті, ұяшықтардың ұжымдық мінез-құлқы әрбір ұяшыққа қарағанда импортты.^[4]

Жүректің электрлік белсенділігін модельдеудің үшінші нұсқасы - бұл модельге төменгі және жоғары деңгейдегі бөлшектерді қосатын «ортаңғы көзқарас» деп қарастыру. Бұл опция үздіксіз ұяшық деп аталатын ұяшықтар блогының әрекетін қарастырады, осылайша масштаб пен деталь проблемаларына жол бермейді. Алынған модель деп аталады бидомен үлгісі, оны көбінесе оны жеңілдетумен алмастырады монодомендік модель.^[4]

Бидомен үлгісі

Электрокардиографияның алға қойылған проблемасы үшін доменді және жазуды сипаттайтын адамның денесінің стильдендірілген көрінісі. Бөлінген екі аймақ және олардың шекарасы қарастырылады, олар жүрек пен оның айналасындағы адамның денесін бейнелейді.

Бидомендік модельдің негізгі болжамы - жүрек тінін екі омдық өткізгіш үздіксіз ортаға бөлуге болады, олар бір-бірімен байланысқан, бірақ жасуша қабықшасы арқылы бөлінген. Бұл ортаны жасушаішілік және жасушадан тыс аймақтар деп атайды, біріншісі жасушалық ұлпаларды, ал екіншісі жасушалар арасындағы кеңістікті білдіреді.^[2]^[1]

Иондық токтың динамикалық моделін қосқанда бидомендік модельдің стандартты формуласы келесідей^[2]

{ displaystyle { begin {case} nabla cdot ( mathbf { sigma} _ {i} nabla V_ {m}) + nabla cdot ( mathbf { sigma} _ {i} nabla u_ {e}) = A_ {m} солға (C_ {m} { frac { жартылай V_ {m}} { жартылай t}} + I_ {ион} (V_ {m}, ж) оңға) + I_ {app} & { text {in}} Omega _ {H} nabla cdot left (( mathbf { sigma} _ {i} + mathbf { sigma} _ {e}) nabla u_ {e} right) + nabla cdot ( mathbf { sigma} _ {i} nabla V_ {m}) = 0 & { text {in}} Omega _ {H} { frac { жарым-жартылай w} { жартылай {t}}} + g (V_ {m}, w) = 0 & { text {in}} Omega _ {H} end {case}}}

қайда ${ displaystyle V_ {m}}$ және ${ displaystyle u_ {e}}$ сәйкесінше трансмембраналық және жасушадан тыс потенциалдар болып табылады, ${ displaystyle I_ {ion}}$ иондық ток болып табылады, ол сондай-ақ қақпа айнымалысы деп аталады ${ displaystyle w}$ (жасушалық деңгейдегі иондық әрекеттерді есепке алу), және ${ displaystyle I_ {app}}$ бұл доменге қолданылатын сыртқы ток. Оның үстіне, ${ displaystyle mathbf { sigma} _ {i}}$ және ${ displaystyle mathbf { sigma} _ {e}}$ жасушаішілік және жасушадан тыс өткізгіштік тензорлары, ${ displaystyle A_ {m}}$ - бұл жасуша мембранасының және көлемнің қатынасы ${ displaystyle C_ {m}}$ бұл аудан бірлігіне мембрана сыйымдылығы. Мұнда домен ${ displaystyle Omega _ {H}}$ жүрек бұлшықетін білдіреді.^[2]

Бидомендік модельдің осы нұсқасының шекаралық шарттары жүректен тыс жасушаішілік потенциал ағыны жоқ деген болжам арқылы алынады, демек

{ displaystyle mathbf { sigma} _ {i} nabla V_ {m} cdot mathbf {n} + mathbf { sigma} _ {i} nabla u_ {e} cdot mathbf {n} = 0 quad quad { text {on}} Sigma}

қайда ${ displaystyle Sigma = жарым-жартылай Omega _ {H}}$ жүрек доменінің шекарасын және ${ displaystyle mathbf {n}}$ үшін сыртқы бірлік болып табылады ${ displaystyle Sigma}$ .^[2]

Монодомендік модель

Монодомендік модель - бұл кейбір физиологиялық болжамдарға қарамастан, ең болмағанда трансмембраналық потенциалға қатысты шынайы электрофизиологиялық құбылыстарды көрсете алатын бидомендік модельді жеңілдету. ${ displaystyle V_ {m}}$ .^[2]^[1]

Стандартты тұжырымдама - келесі белгісіз дифференциалдық теңдеу ${ displaystyle V_ {m}}$ трансмембраналық потенциал:

{ Displaystyle chi C_ {m} { frac { ішінара V_ {m}} { ішінара t}} - nabla cdot left ( mathbf { sigma} _ {i} { frac { lambda } {1+ lambda}} nabla V_ {m} right) + chi I_ {ion} = I_ {app} quad quad { text {in}} Omega _ {H}}

қайда ${ displaystyle lambda}$ - бұл жасушаішілік және жасушадан тыс өткізгіштік тензорларын байланыстыратын параметр.^[2]

Осы модель үшін пайдаланылатын шекаралық шарт болып табылады^[9]

{ displaystyle left ( mathbf { sigma} _ {i} nabla V_ {m} right) cdot mathbf {n} = 0 quad quad { text {on}} жарым-жартылай Omega _ {H}.}

Торс мата моделі

Алдыңғы электрокардиография мәселесінде торс пассивті өткізгіш ретінде қарастырылады және оның моделін Максвелл теңдеулері квазистатикалық болжам бойынша.^[1]^[2]

Стандартты тұжырымдама бір белгісіз скаляр өрісі, торс әлеуеті бар жартылай дифференциалдық теңдеуден тұрады ${ displaystyle u_ {T}}$ . Негізінен, торс модель келесі жалпыланған болып табылады Лаплас теңдеуі

{ displaystyle nabla cdot ( mathbf { sigma} _ {T} nabla u_ {T}) = 0 quad quad { text {in}}} Omega _ {T},}

қайда ${ displaystyle mathbf { sigma} _ {T}}$ - өткізгіштік тензоры және ${ displaystyle Omega _ {T}}$ бұл жүректі қоршап тұрған домен, яғни адамның денесі.^[2]

Шығу

Бидомендік модельге келетін болсақ, торс модель Максвелл теңдеулерінен және үздіксіздік теңдеуі кейбір болжамдардан кейін. Ең алдымен, дененің ішіндегі электрлік және магниттік белсенділік төмен деңгейде пайда болатындықтан, квазистатикалық болжамды қарастыруға болады. Осылайша, денені пассивті өткізгіш ретінде қарастыруға болады, яғни оның сыйымдылық, индуктивті және тарату әсерін елемеуге болады.^[1]

Квазистатикалық болжам бойынша Максвелл теңдеулері болып табылады^[1]

{ displaystyle { begin {case} nabla cdot mathbf {E} & = { frac { rho} { epsilon _ {0}}} nabla times mathbf {E} & = 0 nabla cdot mathbf {B} & = 0 nabla times mathbf {B} & = mu _ {0} mathbf {J} end {case}}}

және үздіксіздік теңдеуі^[1]

{ displaystyle nabla cdot mathbf {J} = 0.}

Оның бұралуы нөлге тең болғандықтан, электр өрісін скалярлық потенциал өрісінің градиентімен, торс потенциалымен көрсетуге болады.

{ displaystyle mathbf {E} = - nabla u_ {T}}

(1)

мұндағы теріс белгі токтың әлеуетті аймақтардан жоғарыдан төменге қарай ағатынын білдіреді.^[1]

Сонымен, жалпы ток тығыздығын өткізгіштік токпен және басқа әр түрлі қолданылатын токтармен өрнектеуге болады, осылайша үздіксіздік теңдеуінен^[1]

{ displaystyle nabla cdot mathbf {J} = nabla cdot ( mathbf {J} _ {app} + mathbf { sigma} _ {T} mathbf {E}) = 0}

(2)

Содан кейін, ауыстыру (1) ішінде (2)

{ displaystyle nabla cdot ( mathbf { sigma} _ {T} nabla u_ {T}) = nabla cdot mathbf {J} _ {app} = I_ {v}}

онда ${ displaystyle I_ {v}}$ - бұл көлем бірлігіне келетін ток.^[1]

Ақырында, жүректен басқа дененің ішінде ток көзі жоқ болғандықтан, көлем бірлігіндегі ток нөлге теңестірілуі мүмкін, бұл жалпыланған Лаплас теңдеуін береді, бұл дененің ішіндегі диффузиялық есептің стандартты тұжырымын білдіреді^[1]

{ displaystyle nabla cdot ( mathbf { sigma} _ {T} nabla u_ {T}) = 0 quad quad { text {in}}} Omega _ {T}.}

Шектік шарт

Шекаралық шарттар денені қоршайтын орталардың, яғни дененің айналасындағы ауаның қасиеттерін ескереді. Әдетте, ауа нөлдік өткізгіштікке ие, яғни ток денеден тыс ағып кете алмайды. Бұл келесі теңдеуде аударылған^[1]

{ displaystyle mathbf { sigma} _ {T} nabla u_ {T} cdot mathbf {n} _ {T} = 0 quad quad { text {on}} Gamma _ {T}}

қайда ${ displaystyle mathbf {n} _ {T}}$ - бұл торсыққа және сыртқа қалыпты өлшем бірлігі ${ displaystyle Gamma _ {T}}$ дененің шекарасы, бұл торстың бетін білдіреді.^[1]^[2]

Торстың өткізгіштігі

Әдетте, торсот изотропты өткізгіштікке ие деп саналады, яғни ток барлық бағытта бірдей жүреді. Дегенмен, торс бос немесе біртекті конвертте емес, бірақ эксперименталды түрде алуға болатын әр түрлі өткізгіштік коэффициенттерімен сипатталатын әр түрлі мүшелерді қамтиды. Сүйектер мен өкпелерді қарастыратын денеде өткізгіштік параметрлерінің қарапайым мысалы келесі кестеде келтірілген.^[2]

Торстың өткізгіштік мәні.^[2]
${ displaystyle sigma _ {T} ^ { text {lungs}}}$	${ displaystyle sigma _ {T} ^ { text {bone}}}$	${ displaystyle sigma _ {T} ^ { text {қалған аймақтар}}}$
(S / см)	(S / см)	(S / см)
${ displaystyle 2.4 times 10 ^ {- 4}}$	${ displaystyle 4 times 10 ^ {- 5}}$	${ displaystyle 6 times 10 ^ {- 4}}$

Жүрек-торсық модельдері

Электрлік белсенділік моделі мен дененің моделі арасындағы байланыстыру эпикардта, яғни жүрек пен дененің арасындағы шекара бетінде қолайлы шекаралық жағдайлар арқылы жүзеге асырылады.^[1]^[2]

Егер жүрек электр моделі мен дененің моделі олардың арасындағы шектеулі немесе жетілмеген ақпарат алмасуымен бөлек шешілсе, екі домен арасындағы керемет электрлік беріліс қарастырылған болса немесе біріктірілуі мүмкін болса, жүрек-торс моделі толық байланыстырылуы мүмкін.^[2]

Толығымен біріктірілген жүрек-торс модельдері

Жүрек пен дененің толық байланысы жүрек пен торсық арасында электр өткізгіштігінің тамаша жағдайын тудырады. Бұл жасушадан тыс потенциал мен дененің әлеуеті арасындағы байланысты орнататын келесі екі теңдеуді ескере отырып жасалады^[2]

{ displaystyle { begin {aligned} u_ {e} & = u_ {T} quad quad & { text {on}} Sigma ( mathbf { sigma} _ {e} nabla u_ { e}) cdot mathbf {n} _ {e} & = - ( mathbf { sigma} _ {T} nabla u_ {T}) cdot mathbf {n} _ {T} quad quad & { text {on}} Sigma. end {aligned}}}

Бұл теңдеулер эпикардия бойындағы потенциалдың да, ағымның да үздіксіздігін қамтамасыз етеді.^[2]

Осы шекаралық шарттарды қолдана отырып, жүректің электрлік белсенділігі үшін бидоменді немесе монодомендік модельді ескере отырып, екі түрлі толық байланысқан жүрек-торсық модельдерін алуға болады. Сандық тұрғыдан алғанда, екі модель есептеу үшін өте қымбат және есептеу шығындары ұқсас.^[2]

Альтернативті шекаралық шарттар

Жүрек пен дененің арасындағы тамаша электр байланысын білдіретін шекаралық жағдайлар ең көп қолданылатын және классикалық болып табылады. Алайда, жүрек пен дененің арасында перикардия, электр берілісіне ерекше әсер ететін серозды сұйықтықты қамтитын екі қабатты қабырға. Ескере отырып сыйымдылық ${ displaystyle C_ {p}}$ және қарсылық ${ displaystyle R_ {p}}$ перикардтың әсерін, осы әсерді ескеретін альтернативті шекаралық шарттарды келесідей тұжырымдауға болады^[10]

{ displaystyle { begin {aligned} R_ {p} mathbf { sigma} _ {e} nabla u_ {e} cdot mathbf {n} & = R_ {p} C_ {p} { frac { жартылай (u_ {T} -u_ {e})} { жартылай t}} + (u_ {T} -u_ {е}) quad quad & { text {on}} Sigma mathbf { sigma} _ {e} nabla u_ {e} cdot mathbf {n} & = mathbf { sigma} _ {T} nabla u_ {T} cdot mathbf {n} quad quad & { text {on}} Sigma end {aligned}}}

Бидомен үлгісімен тұжырымдау

Толығымен байланыстырылған жүрек-торс моделі бидомен үлгісі жүректің электрлік қызметі үшін оның толық түрінде болады^[2]

{ displaystyle { begin {case} nabla cdot ( mathbf { sigma} _ {i} nabla V_ {m}) + nabla cdot ( mathbf { sigma} _ {i} nabla u_ {e}) = A_ {m} солға (C_ {m} { frac { ішінара V_ {m}} { бөлшек t}} + I_ {ион} (V_ {m}, ж) оңға) + I_ {app} & { text {in}} Omega _ {H} nabla cdot left (( mathbf { sigma} _ {i} + mathbf { sigma} _ {e}) nabla u_ {e} right) + nabla cdot ( mathbf { sigma} _ {i} nabla V_ {m}) = 0 & { text {in}} Omega _ {H} { frac { жарым-жартылай w} { жартылай {t}}} + g (V_ {m}, w) = 0 & { text {in}} Omega _ {H} nabla cdot ( mathbf { sigma} _ {T} nabla u_ {T}) = 0 & { text {in}} Omega _ {T} ( mathbf { sigma} _ {i} nabla V_ {m}) cdot mathbf {n} + ( mathbf { sigma} _ {i} nabla u_ {e}) cdot mathbf {n} = 0 & { text {on}} Sigma u_ {e} = u_ {T} & { text {on}} Sigma ( mathbf { sigma} _ {e} nabla u_ {e}) cdot mathbf {n} = ( mathbf { sigma} _ {T} nabla u_ {T}) cdot mathbf {n} & { text {on}} Sigma ( mathbf { sigma} _ {T} nabla u_ {T}) cdot mathbf {n} = 0 & { text {on}} Gamma _ {T} end {case}}}

мұндағы алғашқы төрт теңдеу дербес дифференциалдық теңдеулер бидомен моделін, иондық модель мен дененің моделін, ал қалғандары бидомейн мен торс модельдерінің шекаралық шарттарын және олардың арасындағы байланыс шарттарын білдіреді.^[2]

Монодомендік модельмен тұжырымдау

Толығымен байланыстырылған жүрек-торс моделі монодомендік модель өйткені жүректің электрлік белсенділігі бидомен мәселесінен де күрделі. Шынында да, түйісу шарттары торс әлеуетін жасушадан тыс потенциалмен байланыстырады, оны монодомендік модель есептемейді. Сонымен, теңдеуін де қолдану керек бидомен үлгісі (монодомендік модель алынған сол жорамалдар бойынша):^[2]

{ displaystyle nabla cdot left (( mathbf { sigma} _ {i} + mathbf { sigma} _ {e}) nabla u_ {e} right) + nabla cdot ( mathbf { sigma} _ {i} nabla V_ {m}) = 0 quad quad { text {in}} Omega _ {H}.}

Осылайша, байланыстыру шарттарын өзгертудің қажеті жоқ, ал жүрек-торстың толық моделі екі түрлі блоктан тұрады:^[2]

Алдымен монодомендік модель әдеттегі шекаралық шартпен шешілуі керек:

{ displaystyle { begin {case} chi C_ {m} { frac { ішінара V_ {m}} { ішінара t}} - nabla cdot left ( mathbf { sigma} _ {i} { frac { lambda} {1+ lambda}} nabla V_ {m} right) + chi I_ {ion} = I_ {app} & { text {in}} Omega _ {H} солға ( mathbf { sigma} _ {i} nabla V_ {m} оңға) cdot mathbf {n} = 0 & { text {on}}} ішінара Omega _ {H}. end {істер}}}

Одан кейін, жасушадан тыс потенциалды, торс модельді және байланыс шарттарын есептеуді қосатын модель шешілуі керек:

{ displaystyle { begin {case} nabla cdot left (( mathbf { sigma} _ {i} + mathbf { sigma} _ {e}) nabla u_ {e} right) + nabla cdot ( mathbf { sigma} _ {i} nabla V_ {m}) = 0 & { text {in}} Omega _ {H} nabla cdot ( mathbf { sigma} _ {T} nabla u_ {T}) = 0 & { text {in}} Omega _ {T} mathbf { sigma} _ {T} nabla u_ {T} cdot mathbf {n} = 0 & { text {on}} Gamma _ {T} u_ {e} = u_ {T} & { text {on}} Sigma ( mathbf { sigma} _ {e} nabla u_ {e}) cdot mathbf {n} = ( mathbf { sigma} _ {T} nabla u_ {T}) cdot mathbf {n} & { text {on}} Sigma соңы {істер}}}

Біріктірілген жүрек-торс модельдері

Толық байланыстырылған жүрек-торсық модельдері өте егжей-тегжейлі модельдер, бірақ оларды есептеу өте қымбат.^[2] Мүмкін болатын жеңілдету деп аталады біріктірілмеген болжам онда жүрек жүректен толығымен электрлік оқшауланған болып саналады.^[2] Математикалық тұрғыдан, бұл ток эпикардия арқылы жүректен торсыққа қарай өте алмайтындығына әсер етеді, атап айтқанда^[2]

{ displaystyle mathbf { sigma} _ {e} nabla u_ {e} cdot mathbf {n} = 0 quad quad { text {on}} Sigma.}

Бұл теңдеуді толық байланыстырылған модельдердің шекаралық шарттарына қолдана отырып, электрлік модельдерді есептеу шығындарын төмендететін торс модельдерінен бөлек шешуге болатын екі біріктірілген жүрек-торс модельдерін алуға болады.^[2]

Бидомендік модельмен біріктірілген жүрек-торсық моделі

Жүректің электр белсенділігін бейнелеу үшін бидоменді қолданатын толық байланысқан жүрек-торсық моделінің біріктірілмеген нұсқасы екі бөлінген бөліктен тұрады:^[2]

The бидомен үлгісі оның оқшауланған түрінде

{ displaystyle { begin {case} nabla cdot ( mathbf { sigma} _ {i} nabla V_ {m}) + nabla cdot ( mathbf { sigma} _ {i} nabla u_ {e}) = A_ {m} солға (C_ {m} { frac { ішінара V_ {m}} { бөлшек t}} + I_ {ион} (V_ {m}, ж) оңға) + I_ {app} & { text {in}} Omega _ {H} nabla cdot left (( mathbf { sigma} _ {i} + mathbf { sigma} _ {e}) nabla u_ {e} right) + nabla cdot ( mathbf { sigma} _ {i} nabla V_ {m}) = 0 & { text {in}} Omega _ {H} { frac { жарым-жартылай w} { жартылай {t}}} + g (V_ {m}, w) = 0 & { text {in}} Omega _ {H} ( mathbf { sigma} _ {i} nabla V_ {m}) cdot mathbf {n} + ( mathbf { sigma} _ {i} nabla u_ {e}) cdot mathbf {n} = 0 & { text {on }} Sigma mathbf { sigma} _ {e} nabla u_ {e} cdot mathbf {n} = 0 & { text {on}} Sigma. End {case}}}

Дененің диффузиялық моделі, оның стандартты тұжырымдамасында, потенциалды үздіксіздік шартымен

{ displaystyle { begin {case} nabla cdot ( mathbf { sigma} _ {T} nabla u_ {T}) = 0 & { text {in}} Omega _ {T} mathbf { sigma} _ {T} nabla u_ {T} cdot mathbf {n} = 0 & { text {on}} Gamma _ {T} u_ {T} = u_ {e} & { мәтін {on}} Sigma end {case}}}

Модомендік модельмен біріктірілген жүрек-торсық моделі

Монодомендік модельді қолданатын толық байланысқан жүрек-торс моделі жағдайындағыдай, сонымен қатар сәйкес келмейтін модельде жасушадан тыс потенциалды есептеу керек. Бұл жағдайда үш түрлі және тәуелсіз мәселелерді шешу керек:^[2]

Монодомендік модель әдеттегі шекаралық шартпен:

{ displaystyle { begin {case} chi C_ {m} { frac { ішінара V_ {m}} { ішінара t}} - nabla cdot left ( mathbf { sigma} _ {i} { frac { lambda} {1+ lambda}} nabla V_ {m} right) + chi I_ {ion} = I_ {app} & { text {in}} Omega _ {H} солға ( mathbf { sigma} _ {i} nabla V_ {m} оңға) cdot mathbf {n} = 0 & { text {on}}} ішінара Omega _ {H}. end {істер}}}

Эпикардтағы шекаралық шарты бар жасушадан тыс потенциалды есептеу проблемасы, жасуша ішілік ток ағынын белгілемейді:

{ displaystyle { begin {case} nabla cdot left (( mathbf { sigma} _ {i} + mathbf { sigma} _ {e}) nabla u_ {e} right) + nabla cdot ( mathbf { sigma} _ {i} nabla V_ {m}) = 0 & { text {in}} Omega _ {H} ( mathbf { sigma} _ {i} + mathbf { sigma} _ {e}) nabla u_ {e} cdot mathbf {n} = - ( mathbf { sigma} _ {i} nabla V_ {m}) cdot mathbf {n } & { text {on}} Sigma end {case}}}

Эпикардтағы потенциалды үздіксіздік шекарасының шартымен диффузиялық модель:

{ displaystyle { begin {case} nabla cdot ( mathbf { sigma} _ {T} nabla u_ {T}) = 0 & { text {in}} Omega _ {T} mathbf { sigma} _ {T} nabla u_ {T} cdot mathbf {n} = 0 & { text {on}} Gamma _ {T} u_ {T} = u_ {e} & { мәтін {on}} Sigma end {case}}}

Электрокардиограмма бойынша есептеу

ЭКГ-де прекордиальды жетекші

Толығымен байланыстырылған немесе біріктірілмеген жүрек-торсық модельдерін шешу адамның денесінің әр нүктесінде, атап айтқанда, дененің бүкіл бетінде жүректен пайда болатын электрлік потенциалды алуға мүмкіндік береді. Торсадағы электродтардың орналасуын анықтай отырып, осындай нүктелердегі потенциалдың уақыттық эволюциясын табуға болады. Содан кейін электрокардиограммалар есептеуге болады, мысалы, келесі формулаларды ескере отырып, 12 стандартты нұсқауларға сәйкес^[2]

{ displaystyle { begin {case} I & = u_ {T} (L) -u_ {T} (R) II & = u_ {T} (F) -u_ {T} (R) III & = u_ {T} (F) -u_ {T} (L) aVR & = { frac {3} {2}} (u_ {T} (R) -u_ {w}) aVL & = { frac { 3} {2}} (u_ {T} (L) -u_ {w}) aVF & = { frac {3} {2}} (u_ {T} (F) -u_ {w}) V_ {i} & = u_ {T} (V_ {i}) - u_ {w} quad i = 1, нүкте, 6 соңы {жағдай}}}

қайда ${ displaystyle u_ {W} = (u_ {T} (L) + u_ {T} (R) + u_ {T} (F)) / 3}$ және ${ displaystyle L, R, F, {V_ {i} } _ {i = 1, нүктелер, 6}}$ бұл электродтардың стандартты орналасуы.^[2]

Сандық әдістер

Жүрек-торсық модельдері терминдермен көрсетілген дербес дифференциалдық теңдеулер оның белгісіздері кеңістіктің де, уақыттың да қызметі. Олар өз кезегінде иондық модельмен біріктіріледі, ол әдетте жүйесінде көрсетілген қарапайым дифференциалдық теңдеулер. Сол мәселелерді шешу үшін сандық сұлбаларды қолдануға болады. Әдетте ақырғы элемент әдісі кеңістікті дискретизациялау үшін қолданылады, ал уақыт бойынша дискреттеу үшін жартылай айқын емес ақырлы айырмашылық схемалары қолданылады.^[1]^[2]

Байланыстырылмаған жүрек-торсық моделі сандық тәсілмен қарапайым, өйткені жүректің электрлік моделін торсықтан бөлек шешуге болады, осылайша олардың әрқайсысын шешудің классикалық сандық әдістерін қолдануға болады. Бұл бидомен және монодомен модельдерін мысалы а-ға дейін шешуге болатындығын білдіреді кері дифференциалдау формуласы уақытты дискретизациялау үшін, ал жасушадан тыс потенциал мен торс потенциалын есептеуге арналған есептерді тек ақырғы элементтер әдісін қолдану арқылы оңай шешуге болады, өйткені олар уақытқа тәуелді емес.^[1]^[2]

Толығымен біріктірілген жүрек-торс модельдері, керісінше, күрделі және сандық модельдерді қажет етеді. Мысалы, жүректің жүріс-тұрысын электрлік модельдеу үшін бидомен үлгісін қолданатын толық жүрек-торсық моделін шешуге болады. доменнің ыдырауы Дирихлет-Нейман доменінің ыдырауы сияқты әдістер.^[2]^[11]

Торстың геометриялық моделі

Көптеген органдарды қоса алғанда, үш өлшемді торс модель.^[12]

Модельдеу үшін және электрокардиограмма толық байланыстырылған немесе біріктірілмеген модельдерді қолдана отырып, адам денесін үш өлшемді қайта құру қажет. Бүгінгі күні диагностикалық бейнелеу әдістері МРТ және КТ адамның анатомиялық бөліктерін егжей-тегжейлі қалпына келтіруге және сәйкесінше дененің геометриясын алуға мүмкіндік беретін жеткілікті дәл кескіндер бере алады.^[13] ішкі органдармен, сүйек құрылымымен және бұлшық еттерімен егжей-тегжейлі сипатталған, дененің үш өлшемді моделін құру үшін пайдалы деректер жиынтығы.^[1]

Электрокардиограмманың динамикалық моделі

Нәтижелер егжей-тегжейлі сипатталған болса да, үш өлшемді модельді шешу әдетте өте қымбатқа түседі. Мүмкін болатын жеңілдету - үш қарапайым дифференциалдық теңдеуге негізделген динамикалық модель.^[3]

Жүрек ырғағының квази кезеңділігі үш өлшемді траектория арқылы көбейеді, тартымды шекті цикл ${ displaystyle (x, y)}$ ұшақ. P, Q, R, S және T болатын ЭКГ-нің негізгі шыңдары бекітілген бұрыштарда сипатталған ${ displaystyle theta _ {P}, theta _ {Q}, theta _ {R}, theta _ {S} { text {and}} theta _ {T}}$ , олар келесі үш ODE береді^[3]

{ displaystyle { begin {case} x '& = alfa z- omega y y' & = alpha y + omega x z '& = - sum _ {i in {P, Q, R, S, T }} a_ {i} Delta theta _ {i} { text {exp}} left (- { frac { Delta theta _ {i} ^ {2}} {2b_ {i} ^ {2}}} оң жақ) - (z-z_ {0}) end {жағдайлар}}}

бірге ${ displaystyle alpha = 1 - { sqrt {x ^ {2} + y ^ {2}}}}$ , ${ displaystyle Delta theta _ {i} = ( theta - theta _ {i}) { text {mod}} (2 pi)}$ , ${ displaystyle theta = { text {atan}} 2 (y, x)}$

Сияқты теңдеулерді классикалық сандық алгоритмдер арқылы оңай шешуге болады Рунге-кутта әдістері ODE үшін.^[3]

Сондай-ақ қараңыз

Әдебиеттер тізімі

^ ^а ^б ^c ^г. ^e ^f ^ж ^сағ ^мен ^j ^к ^л ^м ⁿ ^o ^б ^q ^р ^с ^т ^сен Пуллан, Эндрю Дж.; Бист, Мартин Л .; Ченг, Лео К. (2005). Жүректің электрлік белсенділігін математикалық модельдеу: жасушадан дене бетіне және қайтадан. Әлемдік ғылыми. ISBN 978-9812563736.
^ ^а ^б ^c ^г. ^e ^f ^ж ^сағ ^мен ^j ^к ^л ^м ⁿ ^o ^б ^q ^р ^с ^т ^сен ^v ^w ^х ^ж ^з ^аа ^аб ^ак ^{жарнама} ^ае ^аф Булакия, Муриэль; Cazeau, Serge; Фернандес, Мигель А .; Жербо, Жан-Фредерик; Земземи, Неджиб (24 желтоқсан 2009). «Электрокардиограмманы математикалық модельдеу: сандық зерттеу» (PDF). Биомедициналық инженерия шежіресі. 38 (3): 1071–1097. дои:10.1007 / s10439-009-9873-0. PMID 20033779. S2CID 10114284.
^ ^а ^б ^c ^г. МакШарри, П.Е .; Клиффорд, Г.Д .; Тарасенко, Л .; Смит, Л.А. (наурыз 2003). «Синтетикалық электрокардиограмма сигналдарын генерациялаудың динамикалық моделі». Биомедициналық инженерия бойынша IEEE транзакциялары. 50 (3): 289–294. дои:10.1109 / TBME.2003.808805. PMID 12669985. S2CID 544816.
^ ^а ^б ^c ^г. Сызықтар, Г.Т .; Буист, М.Л .; Гротум, П .; Пуллан, А.Ж .; Сунднес, Дж .; Tveito, A. (1 шілде 2002). «Жүрек электрофизиологиясындағы алға қойылған есептердің математикалық модельдері және сандық әдістері». Ғылымдағы есептеу және көрнекілік. 5 (4): 215–239. дои:10.1007 / s00791-003-0101-4. S2CID 123211416.
^ Эйнтховен, Уоллер (1903). «Die galvanometrische Registrierung des menschlichen Elektrokardiogramms, zugleich eine Beurteilung der Anwendung des Capillarelektrometers in der Physiologie». Pfluegers Arch. 99 (9–10): 472–480. дои:10.1007 / BF01811855. S2CID 10400303.
^ Хенрикес, СС .; Плони, Р. (1987). «Резистивтік үзілістердің бір жүрек талшығындағы толқын пішіні мен жылдамдығына әсері». Мед. Биол. Eng. Есептеу. 25 (4): 428–438. дои:10.1007 / BF02443364. PMID 3450994. S2CID 3038844.
^ Мюллер-Борер, Б.Дж .; Эрдман, Д.Дж .; Бьюкенен, Дж. (1994). «Анатомиялық модельденген қарыншалар тіндеріндегі электр байланысы және импульсті көбейту». IEEE Транс. Биомед. Eng. 41 (5): 445–454. дои:10.1109/10.293219. PMID 8070804. S2CID 14407776.
^ Хрен, Р .; Нененен, Дж .; Хоракек, Б.М. (1998). «Акциздеу кезінде имитацияланған эпикардтық потенциалдық карталар миокардтың талшықты құрылымын көрсетеді». Энн. Биомед. Eng. 26 (6): 1022–1035. дои:10.1114/1.73. PMID 9846940. S2CID 9978399.
^ Кинер, Джеймс; Sneyd, James (2009). Математикалық физиология 2009: жүйелік физиология II (2-ші редакциялау). Спрингер. ISBN 978-1-4939-3709-7.
^ Булакия, Муриэль; Фернандес, Мигель А .; Жербо, Жан-Фредерик; Земземи, Неджиб (2007). «Электрокардиограмманы сандық модельдеуге қарай». Жүректі функционалды бейнелеу және модельдеу. Информатика пәнінен дәрістер. Спрингер. 4466: 240–249. дои:10.1007/978-3-540-72907-5_25. ISBN 978-3-540-72906-8.
^ Фернандес, Мигель А .; Земземи, Неджиб (1 шілде 2010). «Жүректің есептеу электрофизиологиясындағы және ЭКГ сандық модельдеуіндегі ажыратылған уақыттық жүру схемалары» (PDF). Математикалық биология. 226 (1): 58–75. дои:10.1016 / j.mbs.2010.04.003. ISSN 0025-5564. PMID 20416327.
^ Феррер, Ана; Себастьян, Рафаэль; Санчес-Кинтана, Дамиан; Родригес, Хосе Ф .; Годой, Эдуардо Дж .; Мартинес, Лаура; Saiz, Javier (2 қараша 2015). «Жүрекшені белсендіруді модельдеуге арналған адамның атриясы мен торсының егжей-тегжейлі анатомиялық-электрофизиологиялық модельдері». PLOS ONE. 10 (11): e0141573. дои:10.1371 / journal.pone.0141573. ISSN 1932-6203. PMC 4629897. PMID 26523732.
^ Шпитцер, V .; Аккерман, Дж .; Шерзингер, А.Л .; Уитлок, Д. (1996 ж. 1 наурыз). «Көрінетін ер адам: техникалық есеп». Американдық медициналық информатика қауымдастығының журналы. 3 (2): 118–130. дои:10.1136 / jamia.1996.96236280. PMC 116294. PMID 8653448.

[Pullan-2005-1] а ^б ^c ^г. ^e ^f ^ж ^сағ ^мен ^j ^к ^л ^м ⁿ ^o ^б ^q ^р ^с ^т ^сен Пуллан, Эндрю Дж.; Бист, Мартин Л .; Ченг, Лео К. (2005). Жүректің электрлік белсенділігін математикалық модельдеу: жасушадан дене бетіне және қайтадан. Әлемдік ғылыми. ISBN 978-9812563736.

[Boulakia-2009-2] а ^б ^c ^г. ^e ^f ^ж ^сағ ^мен ^j ^к ^л ^м ⁿ ^o ^б ^q ^р ^с ^т ^сен ^v ^w ^х ^ж ^з ^аа ^аб ^ак ^{жарнама} ^ае ^аф Булакия, Муриэль; Cazeau, Serge; Фернандес, Мигель А .; Жербо, Жан-Фредерик; Земземи, Неджиб (24 желтоқсан 2009). «Электрокардиограмманы математикалық модельдеу: сандық зерттеу» (PDF). Биомедициналық инженерия шежіресі. 38 (3): 1071–1097. дои:10.1007 / s10439-009-9873-0. PMID 20033779. S2CID 10114284.

[McSharry-2003-3] а ^б ^c ^г. МакШарри, П.Е .; Клиффорд, Г.Д .; Тарасенко, Л .; Смит, Л.А. (наурыз 2003). «Синтетикалық электрокардиограмма сигналдарын генерациялаудың динамикалық моделі». Биомедициналық инженерия бойынша IEEE транзакциялары. 50 (3): 289–294. дои:10.1109 / TBME.2003.808805. PMID 12669985. S2CID 544816.

[Lines_-2002-4] а ^б ^c ^г. Сызықтар, Г.Т .; Буист, М.Л .; Гротум, П .; Пуллан, А.Ж .; Сунднес, Дж .; Tveito, A. (1 шілде 2002). «Жүрек электрофизиологиясындағы алға қойылған есептердің математикалық модельдері және сандық әдістері». Ғылымдағы есептеу және көрнекілік. 5 (4): 215–239. дои:10.1007 / s00791-003-0101-4. S2CID 123211416.

[5] Эйнтховен, Уоллер (1903). «Die galvanometrische Registrierung des menschlichen Elektrokardiogramms, zugleich eine Beurteilung der Anwendung des Capillarelektrometers in der Physiologie». Pfluegers Arch. 99 (9–10): 472–480. дои:10.1007 / BF01811855. S2CID 10400303.

[Henriquez-1987-6] Хенрикес, СС .; Плони, Р. (1987). «Резистивтік үзілістердің бір жүрек талшығындағы толқын пішіні мен жылдамдығына әсері». Мед. Биол. Eng. Есептеу. 25 (4): 428–438. дои:10.1007 / BF02443364. PMID 3450994. S2CID 3038844.

[Muller-Borer-7] Мюллер-Борер, Б.Дж .; Эрдман, Д.Дж .; Бьюкенен, Дж. (1994). «Анатомиялық модельденген қарыншалар тіндеріндегі электр байланысы және импульсті көбейту». IEEE Транс. Биомед. Eng. 41 (5): 445–454. дои:10.1109/10.293219. PMID 8070804. S2CID 14407776.

[Hren-1998-8] Хрен, Р .; Нененен, Дж .; Хоракек, Б.М. (1998). «Акциздеу кезінде имитацияланған эпикардтық потенциалдық карталар миокардтың талшықты құрылымын көрсетеді». Энн. Биомед. Eng. 26 (6): 1022–1035. дои:10.1114/1.73. PMID 9846940. S2CID 9978399.

[9] Кинер, Джеймс; Sneyd, James (2009). Математикалық физиология 2009: жүйелік физиология II (2-ші редакциялау). Спрингер. ISBN 978-1-4939-3709-7.

[Boulakia-2007-10] Булакия, Муриэль; Фернандес, Мигель А .; Жербо, Жан-Фредерик; Земземи, Неджиб (2007). «Электрокардиограмманы сандық модельдеуге қарай». Жүректі функционалды бейнелеу және модельдеу. Информатика пәнінен дәрістер. Спрингер. 4466: 240–249. дои:10.1007/978-3-540-72907-5_25. ISBN 978-3-540-72906-8.

[Fernández-11] Фернандес, Мигель А .; Земземи, Неджиб (1 шілде 2010). «Жүректің есептеу электрофизиологиясындағы және ЭКГ сандық модельдеуіндегі ажыратылған уақыттық жүру схемалары» (PDF). Математикалық биология. 226 (1): 58–75. дои:10.1016 / j.mbs.2010.04.003. ISSN 0025-5564. PMID 20416327.

[12] Феррер, Ана; Себастьян, Рафаэль; Санчес-Кинтана, Дамиан; Родригес, Хосе Ф .; Годой, Эдуардо Дж .; Мартинес, Лаура; Saiz, Javier (2 қараша 2015). «Жүрекшені белсендіруді модельдеуге арналған адамның атриясы мен торсының егжей-тегжейлі анатомиялық-электрофизиологиялық модельдері». PLOS ONE. 10 (11): e0141573. дои:10.1371 / journal.pone.0141573. ISSN 1932-6203. PMC 4629897. PMID 26523732.

[Spitzer-13] Шпитцер, V .; Аккерман, Дж .; Шерзингер, А.Л .; Уитлок, Д. (1996 ж. 1 наурыз). «Көрінетін ер адам: техникалық есеп». Американдық медициналық информатика қауымдастығының журналы. 3 (2): 118–130. дои:10.1136 / jamia.1996.96236280. PMC 116294. PMID 8653448.

[1]

[2]

[3]

[4]

[5]

[6]

[7]

[8]

[9]

[10]

[11]

[12]

[13]