ISSN 1004-4140
CN 11-3017/P

双源CT噪声优化虚拟单能量成像技术在门静脉高压患者肝门静脉成像中的应用价值

柳维娜, 孙小丽, 段永利, 霍萌, 张明霞, 王仁贵

柳维娜,孙小丽,段永利,等. 双源CT噪声优化虚拟单能量成像技术在门静脉高压患者肝门静脉成像中的应用价值[J]. CT理论与应用研究(中英文),xxxx,x(x): 1-8. DOI: 10.15953/j.ctta.2024.084.
引用本文: 柳维娜,孙小丽,段永利,等. 双源CT噪声优化虚拟单能量成像技术在门静脉高压患者肝门静脉成像中的应用价值[J]. CT理论与应用研究(中英文),xxxx,x(x): 1-8. DOI: 10.15953/j.ctta.2024.084.
LIU W N, SUN X L, DUAN Y L, et al. Application Value of Dual-Source CT Noise-Optimized Virtual Monoenergetic Imaging Technology in Portal Vein Imaging of Patients with Portal Hypertension[J]. CT Theory and Applications, xxxx, x(x): 1-8. DOI: 10.15953/j.ctta.2024.084. (in Chinese).
Citation: LIU W N, SUN X L, DUAN Y L, et al. Application Value of Dual-Source CT Noise-Optimized Virtual Monoenergetic Imaging Technology in Portal Vein Imaging of Patients with Portal Hypertension[J]. CT Theory and Applications, xxxx, x(x): 1-8. DOI: 10.15953/j.ctta.2024.084. (in Chinese).

双源CT噪声优化虚拟单能量成像技术在门静脉高压患者肝门静脉成像中的应用价值

详细信息
    作者简介:

    柳维娜,女,首都医科大学附属北京世纪坛医院医学影像科主管技师,研究方向为门静脉高压患者CT门静脉成像研究,E-mail:lwn8810@bjsjth.cn

    通讯作者:

    王仁贵✉,男,首都医科大学附属北京世纪坛医院医学影像科主任医师,研究方向为弥漫性肺疾病及淋巴系统疾病的影像学诊断,E-mail:wangrg@bjsjth.cn

Application Value of Dual-Source CT Noise-Optimized Virtual Monoenergetic Imaging Technology in Portal Vein Imaging of Patients with Portal Hypertension

  • 摘要:

    目的:对比分析双源CT噪声优化虚拟单能量成像(VMI+)不同 keV值对门静脉高压患者门静脉成像质量的影响,探讨VMI+最佳 keV值。方法:对入组的37例门静脉高压患者的双源CT门静脉期图像进行后处理,采用VMI+技术重建40、45、50、55、60、65,70,75和80 keV共9组单能量图像,与自动生成的线性融合系数为0.6(M-0.6)的标准混合图像对比,分别测量10组图像的门静脉主干(MPV)和肝实质的CT值、SD值,分别计算出MPV的信噪比(SNR)和与肝实质的对比噪声比(CNR);在门静脉期最大密度投影(MIP)图像上分别测量MPV、脾静脉(SPV)、门静脉左支(IHLPV)、门静脉右支(IHRPV)的管腔直径;分别记录10组图像中病变血管的检出例数。客观评价采用方差分析,主观评价由两名副主任影像诊断医师采用五分制评价10组图像上MPV与背景肝实质的对比效果以及整体图像质量。结果:门静脉期40 ~ 65 keV 6组图像MPV的CT值、SNR、CNR以及主观评分均比M-0.6组高,差异有统计学意义;70 keV图像上MPV的CT值、SNR、CNR以及主观评分均与M-0.6组相近,差异无统计学意义;75 keV和80 keV两组图像上MPV的CT值,SNR、CNR以及主观评分均比M-0.6组低,差异有统计学意义;各组间MPV、IHLPV、IHRPV和SPV的管腔直径,以及病变血管的检出例数一致,差异无统计学意义;其中40 keV组上MPV的CT值、SNR和CNR 均是10组中最高的。结论:应用VMI+技术(40 ~ 65 keV)可以提高门脉高压患者门静脉期图像质量,有利于门静脉血管的显示,其中40 keV是最佳单能量能级。

    Abstract:

    Objective: This study aimed to compare and analyze the effects of different keV values of dual-source computed tomography (CT) noise-optimized virtual monoenergetic imaging (VMI+) on the quality of portal vein imaging in patients with portal hypertension and to determine the optimal keV value for VMI+. Methods: Postprocessing was performed on portal venous phase images of 37 patients with portal hypertension who underwent dual-source CT. The VMI+ technique was used to reconstruct nine groups of monoenergetic images at 40, 45, 50, 55, 60, 65, 70, 75, and 80 keV. Compared with the standard mixed images automatically generated with a linear blending coefficient of 0.6 (M-0.6), the CT values and standard deviation of the main portal vein (MPV) and liver parenchyma were measured for each group of images. The signal-to-noise ratio (SNR) of the MPV and the contrast-to-noise ratio (CNR) with the liver parenchyma were calculated. On the maximum intensity projection images of the portal venous phase, the luminal diameters of the MPV, splenic vein (SPV), left branch of the portal vein (IHLPV), and right branch of the portal vein (IHRPV) were measured, and the number of detected cases of diseased vessels in each image group was recorded. Objective evaluation was conducted using analysis of variance. Subjective evaluation was performed by two deputy chief radiologists who used a five-point scale to assess the contrast effect of the MPV against the background liver parenchyma and the overall image quality in the 10 groups of images. Results: The CT value, SNR, CNR of the MPV, and subjective image scores measured for the 40 ~ 65 keV portal vein groups were higher than those for the M-0.6 group, and the difference was statistically significant. The CT value, SNR, CNR of the MPV, and subjective image scores measured for 70 keV images were similar to those of the M-0.6 group, with no statistical difference. The CT value, SNR, CNR of the MPV, and subjective image scores measured for the 75 and 80 keV groups were lower than those for the M-0.6 group, and the difference was statistically significant. The diameters of the MPV, IHLPV, IHRPV, and SPV were the same among all the groups, and the same number of pathological vessels were detected, with no statistical significance. Among them, the 40 keV group had the highest CT values, SNR, and CNR of the MPV among all 10 groups. Conclusion: The application of VMI+ technology (40 ~ 65 keV) in dual-source CT significantly improves the portal-vein image quality of patients with portal hypertension, which is beneficial for displaying portal vein vessels. The optimal monoenergetic energy level is 40 keV.

  • 随着勘探开发逐步深入,勘探开发对象多转向薄层或薄互层类型油气藏[1-2],而薄互层调谐效应影响了AVO变化规律[3],导致薄互层油藏AVO应用困难。

    前人针对薄互层储层在地震响应特征、频率特征及其厚度关系方面都开展了大量研究[4-7]。Widess[8]首先基于零相位地震子波研究了薄层厚度与调谐振幅之间关系,当薄层厚度为1/4波长时,反射波波形为入射波波形的一阶导数,突破了求取薄层厚度的极限,同时将厚度小于1/4波长的地层定义为薄层。

    此后围绕薄互层厚度与地震响应之间关系,当薄互层存在时,由于调谐效应、层间多次波及衰减等因素的影响,会导致其反射特征和单界面反射有很大的区别,尤其对AVO响应有很大差异。Simmons和Backus[9]指出转换波对薄层地震响应有一定的影响,仅用Zoeppritz方程无法准确得到真实薄层的AVO响应特征,同时也分析了转换波对薄层AVO响应特征的影响;Liu和Schmitt [10]推导了薄储层振幅随偏移距变化的AVO响应精确解,分析了储层厚度、泊松比等因素对薄层AVO规律的影响;Pan等[11]指出传统的Zoeppritz方程不适合分析薄层的反射和透射响应,利用Breshkovskikh方程分析了弹性介质中薄层的AVO特征,指出声波和弹性波介质中薄层响应存在差异,并且衰减严重的地层,其AVO曲线接近于Zoeppritz方程的结果;赵伟等[12]分析了薄互层调谐效应对AVO的影响,当薄层厚度小于1/16波长时,薄层调谐效应则会导致AVO的变化规律发生反转;周丽等[13]设计了不同含油气性薄层砂体模型来进行波动方程正演模拟,获得了薄层地震正演AVO曲线,分析了塔河油田碎屑岩储层含不同流体时的AVO特征;陈珊等[14]对薄层进行Brekhovskikh正演,发现不同入射角下振幅与厚度关系不同,与叠后全叠加振幅相比,截距属性和薄层厚度之间的线性关系更强,说明截距属性比全叠加振幅更有利于薄层厚度预测;常锁亮等[15]基于薄层反射系数谱理论中Brekhovskikh方程,计算并总结围岩变化对煤层AVO曲线、属性及道集的影响;郭智奇等[16]分析了薄层反射AVO响应与薄层厚度有关,薄层厚度增加时,对其含油和含气饱和度的变化,梯度属性比截距属性更敏感;王云专等[17]对不同薄层类型进行地震反射特征分析研究,得出薄层厚度与峰值频率的定量关系式,并对其进行高阶近似,得到适用于不同类型薄层的峰值频率随厚度变化的公式,并探讨了厚度和子波峰值频率对薄层峰值频率的升频作用;陈小宏等[18]将分频技术与AVO研究密切结合是有效解决薄互层气藏地震AVO识别,讨论了地震AVO与地震频率的关系,指出进行分频AVO分析的依据。

    本文主要从理论推导和正演模拟出发,以薄层地震反射特征为研究对象,系统开展了薄层AVO变化规律研究。结合薄层耦合AVO反射表达式,构建不同厚度的薄层概念模型,利用反射率法进行叠前AVO正演,探究不同频带宽度对叠前AVO稳定性的影响。在此基础上研究了薄层分频AVO 特征,通过叠前道集优势信息重构可有效增强烃类检测的稳定性,实际资料也验证了该方法的有效性。

    设平面波从顶层入射到下面各层中,则顶层中位移−应力关系[19-20]为:

    $$ \left( {\begin{array}{*{20}{c}} {{u^{\left( 1 \right)}}} \\ {{w^{\left( 1 \right)}}} \\ {{\sigma _{zz}}^{\left( 1 \right)}} \\ {{\sigma _{zx}}^{\left( 1 \right)}} \end{array}} \right) = {\boldsymbol{A}}_{ij}^{\left( 2 \right)}\left( {\begin{array}{*{20}{c}} {{u^{\left( 3 \right)}}} \\ {{w^{\left( 3 \right)}}} \\ {{\sigma _{zz}}^{\left( 3 \right)}} \\ {{\sigma _{zx}}^{\left( 3 \right)}} \end{array}} \right) \text{,} $$ (1)

    式中,上标1、2和3为第1层、第2层和第3层介质;$ u $$ w $为位移值,$ {\sigma _{zz}} $为垂直应力;$ {\sigma _{zx}} $为水平方向应力。$ {{\boldsymbol{A}}_{ij}} $为传播矩阵,上标2为传播矩阵$ {{\boldsymbol{A}}_{ij}} $仅和介质2有关。

    反射系数和透射系数计算公式为:

    $$ \begin{aligned}&{R_{pp}} = \frac{{A_2^1}}{{A_1^1}} \text{,} {R_{ps}} = \frac{{B_2^1}}{{A_1^1}}\frac{{{\alpha _1}}}{{{\beta _1}}} \text{,}\\ &{T_{pp}} = \frac{{A_1^3}}{{A_1^1}}\frac{{{\alpha _1}}}{{{\alpha _3}}} \text{,} {T_{ps}} = \frac{{B_1^3}}{{A_1^1}}\frac{{{\alpha _1}}}{{{\beta _3}}} , \end{aligned}$$ (2)

    式中,$ {R_{pp}} $$ {R_{ps}} $$ {T_{pp}} $$ {T_{ps}} $分别为纵波反射系数、横波反射系数、纵波透射系数和横波透射系数;$ \alpha $$ \beta $分别为横波和纵波速度。

    当薄层厚度等于1/4地震波长时,地震波场出现相长干涉,使振幅出现极大值,该现象称为薄层的调谐效应。为了进一步分析薄层厚度、地震主频与其叠前AVO特征关系,这里设计单砂体楔状地质模型,楔状模型弹性参数见表1,分别采用5~40 Hz的窄频带子波和5~80 Hz的宽频带子波(图1)对单砂体楔状模型进行叠前AVO正演模拟,正演过程中,首先采用薄层反射系数公式计算不同厚度条件下,随角度变化的反射系数,然后利用不同带宽子波进行褶积运算,叠加后形成不同带宽的正演模拟结果(图2)。

    表  1  楔状模型弹性参数信息表
    Table  1.  Wedge model elastic parameter information
    AVO
    类型
    泥岩围岩 砂岩
    纵波速度/(m/s) 横波速度/(m/s) 密度 纵波速度 横波速度 密度
    II类 4150 2220 2.63 4300 2610 2.40
    下载: 导出CSV 
    | 显示表格
    图  1  不同频段地震子波
    Figure  1.  Seismic subgroups in different frequency bandwidths
    图  2  楔状模型正演模拟结果
    Figure  2.  Wedge model forward simulation results

    从不同带宽地震正演结果上看,随着砂岩储层厚度的增加,地震振幅能量逐渐增加,当储层厚度达到调谐厚度(窄频子波调谐厚度为47.5 m,宽频子波调谐厚度为25.5 m)时,地震振幅出现极大值;随着砂岩储层厚度进步增加,地震振幅能量逐渐减弱直至能量趋于稳定,其中宽频带正演对薄层识别能力上要优于窄频带正演。为了进一步精细分析不同带宽地震叠前响应随厚度变化特征,这里分别以20 m间隔抽取正演叠前AVO道集并拾取振幅随角度变化曲线(图3图4)。

    图  3  窄频带叠前道集振幅随角度变化曲线
    Figure  3.  Narrow bandwidth pre-stack angle gather amplitude curve
    图  4  宽频带叠前道集振幅随角度变化曲线
    Figure  4.  Wide bandwidth pre-stack angle gather amplitude curve

    图3为窄频带叠前道集振幅随角度变化曲线结果,当地层厚度大于调谐厚度时,其砂岩顶底干涉程度弱,砂岩顶面AVO特征基本稳定且差异不大;当地层厚度小于调谐厚度时,由于顶底界面的干涉作用,薄层叠前响应较为复杂,AVO在远角度伴随能量反转,影响AVO属性的准确提取。

    图4为宽频带叠前道集振幅随角度变化曲线,这里宽频地震叠前道集对薄层识别能力得到显著提升,调谐厚度也由47.5 m变为25.5 m,当地层厚度大于调谐厚度时,其模拟结果与窄频带相似AVO曲线相对比较稳定,远角度未出现能量反转现象,都表现为稳定的II类AVO特征,当地层厚度小于调谐厚度时,AVO道集远角度仍有能量反转现象。

    对比窄频带和宽频带模拟结果,窄频带由于带宽窄顶底干涉现象严重,在调谐厚度附近AVO曲线“V”字型特征更显著,而对于宽频带子波模拟结果,由于地震高频成分增加,AVO曲线“V”字型现象得到压制。因此频带宽度决定着地震薄层的分辨能力和AVO曲线稳定性,为了压制薄层AVO的调谐效应,应该采用宽频地震资料进行AVO分析。

    由于地震频带宽度直接影响AVO的稳定性,尤其不利于薄层的AVO分析,因此需要借助高分辨率时频分析方法精细解析地震数据的频带特征。目前常用的时频分析方法包括短时傅里叶变换、小波变换、匹配追踪等。本文研究采用Wigner-Ville分布进行时频分析研究[21]。这里假设任一信号$f\left( t \right) \in {L^2}\left( R \right)$可以分解为:

    $$ f\left( t \right) = \sum\limits_{n = 0}^{ + \infty } {\left\langle {{R^n}f,{g_{{\gamma _n}}}} \right\rangle } {g_{{\gamma _n}}}\text{,} $$ (3)

    其中:${g_{{\gamma _n}}}$为第n次匹配因子;${R^n}f$为第n次匹配余项;${\gamma _n} = \left( {{s_n},{u_n},{\xi _n}} \right)$是尺度参数${s_n}$、位移参数${u_n}$、频率相位参数${\xi _n}$的函数,且满足${g_{{\gamma _n}}} = \displaystyle\frac{1}{{\sqrt {{s_n}} }}g\left( {\frac{{t - {u_n}}}{{{s_n}}}} \right){{{\mathrm{exp}}}({i{\xi _n}t}})$为最佳匹配余项。

    依据Wigner分布的扩张和平移特性,对于任一信号可以表示为:

    $$\begin{gathered} W\left( {f,f} \right)\left( {t,\omega } \right) = \\ \dfrac{1}{{2{\text{π}} }}\int_{ - \infty }^{ + \infty } {f\left( {t + \dfrac{\tau }{2}} \right)} {f^ * }\left( {t - \dfrac{\tau }{2}} \right){{\mathrm{exp}}({ - i\omega \tau }})d\tau \text{,} \\ \end{gathered}$$ (4)

    对于Wigner分布可以表示为:

    $$ {W_{{g_\gamma }}}\left( {t,\omega } \right) = {W_g}\left( {\frac{{t - u}}{s},s\left( {\omega - \xi } \right)} \right)\text{。} $$ (5)

    因此去掉交叉项后信号的时频平面分布可表示为:

    $$ \begin{gathered} Ef\left( {t,\omega } \right) = \\ \sum\limits_{n = 0}^{ + \infty } {{{\left| {\left\langle {{R^n}f,{g_{{\gamma _n}}}} \right\rangle } \right|}^2}} {W_g}\left[ {\frac{{t - {u_n}}}{{{s_n}}},{s_n}\left( {\omega - {\xi _n}} \right)} \right]。 \\ \end{gathered} $$ (6)

    Wigner-Ville分布具有较高时频分辨率,其时间和频率分辨率之积能够达到测不准原理的极限分辨率,Wigner-Ville分布能量聚集性好,且没有交叉项的干扰。

    通过宽频地震楔状模型正演发现,在调谐厚度内AVO曲线在远偏移距存在能量反转现象不利于烃类检测分析,为此需要开展薄层AVO稳定性研究。分频AVO能量扫描可以较全面地解析地震数据频带特征和分布规律。

    本文利用Wigner-Ville分布时频分析方法开展叠前道集分频能量扫描分析。图5为宽频带叠前道集分频扫描结果(20 m砂厚),通过拾取砂岩顶面的AVO曲线进行分析发现,当截止频率低于60 Hz时,砂岩顶面AVO随着角度增加而增加,AVO曲线基本稳定;当截至频率超过60 Hz时,砂岩顶面AVO曲线开始出现不稳定现象,在远角度产生能量反转,因此可以选择AVO曲线稳定变化优势频带信息,即10~60 Hz的频率成分进行AVO曲线重构。

    图  5  宽频带叠前道集分频扫描结果
    Figure  5.  Wide bandwidth pre-stack angle gather crossover frequency scanning results

    图6是10~60 Hz和10~65 Hz叠前信息分别进行AVO曲线重构后对比效果。当重构截至频率在60 Hz以内时,叠前道集砂岩顶面远道未出现能量反转现象,AVO道集比较稳定;当截至频率超过60 Hz时,叠前道集砂岩顶面远道开始出现能量反转,不利于叠前道集AVO属性信息的提取及烃类检测研究,为此需要通过叠前分频扫描确定优势频带信息以确保AVO曲线变化的稳定性,以降低储层烃类检测多解性。

    图  6  宽频带叠前道集重构结果对比
    Figure  6.  Comparison of reconstruction results of wide bandwidth pre-stack angle gather

    Y气田位于东海陆架盆地西湖凹陷西次凹内,构造总体上表现为形态比较完整的NE-SW向背斜构造。Y气田自下而上发育新生界古近系、新近系及第四系;渐新统花港组是主要含气层段,其中主产层为花港组上段H3~H5砂层组,地层厚度约250 m,以浅水三角洲环境下的砂泥岩沉积为主,有利储层为分流河道型储层。Y气田ODP实施后地质储量大幅减少,持续高效开发面临巨大挑战。

    目前,探井Y3井在H3、H5层均已钻遇有利相带储层,砂地比较高,而H4层尚未钻遇有利相带储层,已钻井整体为泥包砂埋藏相对较深,主要位于3000~3600 m的深度范围,储层物性为中低孔中低渗特征,该层单砂体厚度为2~20 m不等,均在调谐厚度(32 m)之内,受储层厚度影响H4层的潜力目标评价存在不确定性,因而高效开发对薄层AVO分析提出了更高的要求[22]

    图7为过Y气田花港组水道叠后地震剖面,其中Y3井H4层钻遇9.5 m薄层砂岩、H5层钻遇41 m厚层砂岩,其中H4砂岩厚度小于调谐厚度,H5砂岩厚度大于调谐厚度,但不论薄厚含气砂岩都表现为近似相同振幅能量地震响应特征。

    图  7  过Y气田花港组水道叠后地震剖面
    Figure  7.  Post-stack seismic section of the Huagang Formation over a Y gas field

    从叠前道集上看(图8),由于小于调谐厚度,Y3井H4层远道存在明显能量反转特征,A5 S井H4薄层含气砂岩也表现为类似的叠前AVO特征。而对于厚层含气砂岩,Y3井H5层AVO特征相对比较稳定。

    图  8  Y气田已钻井叠前道集AVO分析(位置见图7所示)
    Figure  8.  AVO analysis of pre-stack angle gather in a Y gas field (location shown in Figure 7)

    针对不稳定的叠前道集进行分频能量扫描分析(图9),发现叠前道集振幅能量在0~10 Hz、10~20 Hz、20~30 Hz、30~40 Hz、40~50 Hz、50~60 Hz多个频段都有有效信号反射,但在50~60 Hz以上AVO道集开始出现远角度能量反转现象,叠前道集多解性开始增强,为此选择AVO曲线稳定变化优势频带信息进行AVO曲线重构。

    图  9  叠前道集分频能量扫描分析(位置①)
    Figure  9.  Pre-stack angle gather crossover energy scan analysis (location shown in position ①)

    图10分别为全频带地震PG属性和优势频带信息重构PG属性剖面。相比全频带地震PG属性剖面,优势信息重构PG属性烃类检测异常整体明显改善,高低部位受地层调谐厚度影响较小,都呈明显的含气性异常,其中H4 c和H4 d水道内部叠置交切关系更为明显。结合H4层优势信息重构PG属性异常,在低部位河道主体部位部署了A11和A13两口开发评价井,在H4层分别钻遇18.5 m、20 m优质含气砂岩,孔隙度18%,基本探明了H4层有利相带储层含气区。通过该技术进步打开了Y气田花港组开发生产新局面,有效推动了东海Y气田的开发进程。

    图  10  常规地震PG属性和优势信息重构PG属性剖面
    Figure  10.  Conventional seismic PG attributes and advantageous information reconstructed PG attribute profiles

    本文主要从理论推导和正演模拟出发,以薄层地震反射特征为研究对象,系统开展了薄储层AVO变化规律及其分频AVO特征研究,主要获得以下几方面认识。

    (1)结合楔状模型正演,明确不同地震带宽下薄层叠前AVO响应的差异,有效信息频带宽度决定着砂岩储层的分辨能力及叠前道集AVO曲线的稳定程度;频带越宽,对于储层的分辨能力越强,叠前道集AVO曲线特征越稳定。

    (2)利用Wigner-Ville这种高精度时频分析方法,有利于精细分析不同频带下薄储层地震数据的叠前道集的AVO曲线特征,以选择优势频带进行叠加,保证薄储层AVO分析的稳定性。

    (3)实际地震资料表明,利用高精度时频分析方法对叠前地震数据进行分频扫描,并选择优势频带信息进行AVO曲线重构,能较好地削弱AVO分析不确定性,在东海Y气田获得较好的应用成效,更有利于对于中深层薄储层的AVO刻画分析。

    (4)在实际资料处理过程中,随机噪声、叠前地震道集近远道能量不一致以及频谱差异性等也会影响薄储层的AVO响应规律研究,因此需要在对叠前道集进行精细去噪和一致性处理后,再开展分频AVO分析。

  • 图  1   感兴趣区的绘制

    注:ROI1是门静脉主干,ROI2是肝实质,测量中使用工作站上的复制和粘贴功能。

    Figure  1.   Placement of the regions of interest (ROIs)

    图  2   主观评价得分图例

    注:所有图像中最低得分是2分,图像是同一患者同一层面,层厚是0.75 mm,窗宽和窗位均一致的条件下。

    Figure  2.   Subjective scoring legend

    图  3   VMI+ 40 ~ 80 keV的9组轴位图像的对比

    注:所有图像是同一患者同一层面,层厚是0.75 mm,窗宽和窗位均一致的条件下。

    Figure  3.   Comparison of images of VMI+ 40 ~ 80 keV groups

    图  4   VMI+40 keV组,70 keV组和对照组M-0.6图像的对比

    注:所有轴位图像是同一患者同一层面,层厚是0.75 mm,窗宽和窗位均一致的条件下;所有MIP图像是同一患者,层厚、窗宽和窗位均一 致的条件下。

    Figure  4.   Comparison of images of VMI+ 40 keV, 70 keV, and M-0.6 control groups.

    表  1   VMI+(40 ~ 80 keV)组与M-0.6组MPV图像的客观评价$(\bar x\pm s) $

    Table  1   Objective evaluation of MPV images of VMI+ (40–80 kev) group and M-0.6 group $(\bar x\pm s) $

    重建方法单能量组M-0.6统计检验
    40 keV45 keV50 keV55 keV60 keV65 keV70 keV75 keV80 keVt/FP
    MPVCT值/HU674±110b562±58b442±52b380±32b292±34b241±29b218±26191±21b163±22b217±2683.84<0.05a
    SNR11.6±1.3b11.2±0.9b10.8±0.8b10.1±0.9b9.8±1.0b9.5±0.9b8.5±0.88.0±0.7b7.8±0.7b8.5±0.922.71<0.05a
    CNR12.8±2.3b10.8±1.4b10.2±1.1b9.8±1.3b7.9±1.1b7.5±0.9b6.2±0.95.6±0.8b4.8±0.7b6.2±0.835.74<0.05a
    肝实质CT值/HU201±19b188±26b140±10b129±18b120±17b119±15b108±13101±12b96±10b106±1461.53<0.05a
    SD41.8±6.1b35.6±3.4b31.5±2.8b26.7±2.1b23.6±2.2b20.1±1.9b19.6±1.818.8±1.7b16.2±1.8b19.4±1.9b45.29<0.05a
    注:事后比较:a单能量各组比较,P<0.05;b与M-0.6组比较,P<0.05。
    下载: 导出CSV

    表  2   MPV及其侧枝血管的直径 $(\bar x\pm s) $

    Table  2   Diameters of MPV and its collateral vessels $(\bar x\pm s) $

    重建方法单能量组M-0.6统计检验
    40 keV45 keV50 keV55 keV60 keV65 keV70 keV75 keV80 keVFP
    MPV直径1.68±0.551.70±0.511.69±0.531.70±0.531.68±0.561.70±0.521.71±0.521.69±0.541.69±0.551.71±0.530.780.539
    IHLPV0.84±0.350.85±0.330.83±0.350.86±0.330.83±0.340.85±0.340.84±0.330.82±0.370.86±0.320.83±0.350.6110.551
    IHRPV0.97±0.381.00±0.350.98±0.360.97±0.370.99±0.360.98±0.350.99±0.350.97±0.390.99±0.380.98±0.360.2310.875
    SPV直径1.32±0.421.35±0.401.33±0.451.36±0.391.32±0.441.34±0.411.35±0.411.32±0.441.33±0.431.36±0.410.6170.604
    下载: 导出CSV

    表  3   门静脉高压伴发病变血管的检出例数(例)

    Table  3   Number of detected cases of associated abnormal blood vessels in portal hypertension (number of cases)

    重建方法40 keV45 keV50 keV55 keV60 keV65 keV70 keV75 keV80 keVM-0.6
    食管、胃底静脉曲张18181818181818181818
    脐周静脉曲张7777777777
    腹膜后交通开放6666666666
    下载: 导出CSV

    表  4   两名副主任医师的主观评价得分 (x±s )

    Table  4   Subjective evaluation scores by two deputy chief physicians (x±s )

    重建方法 单能量组 M-0.6 统计检验
    40 keV 45 keV 50 keV 55 keV 60 keV 65 keV 70 keV 75 keV 80 keV Z P
    主观评分1 4.53±0.50b 4.51±0.50b 4.52±0.45b 4.23±0.56b 4.10±0.61b 4.01±0.56b 3.67±0.62 3.28±0.64b 3.05±0.45b 3.66±0.65 40.03 <0.05a
    主观评分2 4.55±0.49b 4.48±0.51 4.50±0.46b 4.28±0.56b 4.12±0.62b 4.03±0.55b 3.70±0.64 3.30±0.67b 2.99±0.57b 3.57±0.68 39.41 <0.05a
    注:事后比较:a单能量各组比较,P<0.05;b与M-0.6组比较,P<0.05。
    下载: 导出CSV
  • [1] 刘腾, 卢俊宇, 何松. 肝硬化门静脉高压诊断方法的研究进展[J]. 重庆医学, 2023, 52(21): 3303−3307. DOI: 10.3969/j.issn.1671-8348.2023.21.019.

    LIU T, LU J Y, H S. Research advances in diagnostic methods of liver cirrhosis portal hypertension[J]. Chongqing Medical Journal, 2023, 52(21): 3303−3307. DOI: 10.3969/j.issn.1671-8348.2023.21.019. (in Chinese).

    [2]

    LAURA T, GUADALUPE G T. Portal Hypertension: Pathogenesis and Diagnosis[J]. Science Direct., 2019, 23(4): 573−587.

    [3] 高德军, 罗鹰, 陈首名, 等. 双源CT门静脉成像在肝硬化门静脉侧支循环的临床应用[J]. 中国CT和MRI杂志, 2019, 17(8): 100−102,106. DOI: 10.3969/j.issn.1672-5131.2019.08.030.

    GAO D J, Luo Y, CHEN S M, et a1. Clinical application of Dual source CT portographv in portal collateral circulation in hepatocirrhosis[J]. Chinese Journal of CT and MRI, 2019, 17(8): 100−102,106. DOI: 10.3969/j.issn.1672-5131.2019.08.030. (in Chinese).

    [4]

    RENGO M , BELLINI D , De CECCO C N , et al . The optimal contrast media policy in CT of the liver. Part I: Technical notes[J]. Acta Radiologica, 2011, 52(5): 467−472

    [5] 陈琳钰, 方姝, 陈勇, 等. 虚拟单能图像与常规CT图像的图像质量对比[J]. CT理论与应用研究, 2022, 31(2): 219−226.

    CHEN L Y, FANG S, CHEN Y, et al. Comparison of image quality between virtual monochromatic images and conventional CT images[J]. CT Theory and Applications, 2022, 31(2): 219−226. (in Chinese).

    [6] 胡斌, 陈静静, 徐文坚, 等. 双能CT改进的虚拟单能成像重建算法的门静脉图像质量评估[J]. 首都医科大学学报, 2019, 40(6): 943−947. DOI: 10.3969/j.issn.1006-7795.2019.06.023.

    HU B, CHEN J J, XU W J, et al. Assessent of portal vein based on a novel virtual monoenergetic algorithm in dual-energy CT[J]. Journal of Capital Medical University, 2019, 40(6): 943−947. DOI: 10.3969/j.issn.1006-7795.2019.06.023. (in Chinese).

    [7]

    ZENG Y W, GENG D Y, ZHANG J. Noise-optimized virtual monoenergetic imaging technology of the third-generation dual-source computed tomography and its clinical applications[J]. Quantitative Imaging in Medicine and Surgery, 2021, 11(11): 4627−4643. DOI: 10.21037/qims-20-1196.

    [8] 杨桐立, 王照谦, 杨志强, 等. 双能量CT Mono+技术临床应用的研究进展[J]. 临床放射学杂志, 2023, 42(2): 340−344.
    [9] 周彦宇, 林雁玲. 双源CT虚拟单能量成像技术不同keV值对门静脉成像质量提高的研究[J]. 现代医用影像学, 2022, 31(2): 199−203. DOI: 10.3969/j.issn.1006-7035.2022.2.xdyyyxx202202001.

    ZHOU Y Y, LIN Y L. Study on the improvement of portal vein imaging quality with different keV values of Dual-energy CT virtual single energy imaging technology[J]. Modern Medical Imageology, 2022, 31(2): 199−203. DOI: 10.3969/j.issn.1006-7035.2022.2.xdyyyxx202202001. (in Chinese).

    [10] 黄亚博, 杨莹, 王晋君. CT门静脉成像参数评估乙型肝炎肝硬化患者肝功能分级及EVB风险的价值[J]. 肝脏, 2024, 29(2): 162−165, 173.

    HUANG Y B, YANG Y, WANG J J. The value of portal vein imaging parameters of computer tomography in evaluating liver function grading and esophageal variceal bleeding risk in hepatitis B patients with cirrhosis[J]. Chinese Hepatology, 2024, 29(2): 162−165, 173. (in Chinese).

    [11] 莫彩建, 翁建辉, 陈金水, 等. MSCT 联合 MSCTP 技术在肝硬化门静脉高压症术前诊断中的应用价值[J]. 海南医学, 2022, 33(1): 86−89.

    MO C J, WONG J H, CHEN J S, et al. Application value of multi-slice spiral CT combined with multi-slice spiral CT portography in preoperative diagnosis of cirrhotic portal hypertension[J]. Hainan Medical Journal, 2022, 33(1): 86−89. (in Chinese).

    [12] 杨克柽, 陈伟. 肝脏门静脉期CT增强的影响因素与对策[J]. 国际医学放射学杂志, 2015, 38(6): 562−566,571. DOI: 10.3874/j.issn.1674-1897.2015.06.Z0612.

    YANG K G, CHEN W. Contrast enhancement of abdominal CT at portal venous phase: Considerations and solutions[J]. International Journal of Medical Radiology, 2015, 38(6): 562−566,571. DOI: 10.3874/j.issn.1674-1897.2015.06.Z0612. (in Chinese).

    [13]

    PARAKH A, MACRI F, SAHANI D. Dual-energy computed tomography: Dose reduction, series reduction, and contrast load reduction in dual-energy computed tomography[J]. Radiologic Clinics of North America, 2018, 56: 601−624. DOI: 10.1016/j.rcl.2018.03.002.

    [14] 王敏, 郭美玲, 李剑. CT低管电压对胸痛三联征排查 图像质量及辐射剂量的影响[J]. 中国医学影像学杂志, 2013, 21(9): 715−717,719. DOI: 10.3969/j.issn.1005-5185.2013.09.023.

    WANG M, GUO M L, LI J. Influence of CT low tube voltage on image quality and radiation dose in patients with chest pain[J]. Chinese Journal of Medical Imaging, 2013, 21(9): 715−717,719. DOI: 10.3969/j.issn.1005-5185.2013.09.023. (in Chinese).

    [15]

    D’HGELO T, CICEM G, MZIOTTI S, et a1. Dual energy computed tomography virual monoenergetic imaging: Technique clinical applications[J]. British Journal of Radiology, 2019, 92: 20180546.

    [16] 王珍, 张声旺, 顾潜彪, 等. CT影像组学评估门静脉高压及预测曲张静脉出血[J]. 中国医学影像技术, 2020, 36(12): 1896−1899.

    WANG Z, ZHANG S W, GU Q B, et al. CT radiomics for assessing portal hypertension and predicting variceal hemorrhage[J]. Chinese Journal of Medical Imaging Technology, 2020, 36(12): 1896−1899. (in Chinese).

    [17] 崔楚坤, 李家欢, 曾鹏程, 等. 多层螺旋CT血管造影在冠状动脉病变程度诊断中的临床价值[J]. 海南医学, 2019, 30(7): 875−877. DOI: 10.3969/j.issn.1003-6350.2019.07.018.

    CUI C K, LI J H, ZENG P C, et al. Clinical value of multi-slice spiral CT angiography in the diagnosis of coronary artery lesions[J]. Hainan Medical Journal, 2019, 30(7): 875−877. DOI: 10.3969/j.issn.1003-6350.2019.07.018. (in Chinese).

    [18] 王磊. 多层螺旋CT门静脉血管成像对食管胃底静脉曲张出血的预测价值[J]. 实用医技杂志, 2018, 25(8): 861−863. DOI: 10.19522/j.cnki.1671-5098.2018.08.026.
图(4)  /  表(4)
计量
  • 文章访问数:  53
  • HTML全文浏览量:  34
  • PDF下载量:  5
  • 被引次数: 0
出版历程
  • 收稿日期:  2024-06-04
  • 修回日期:  2024-09-17
  • 录用日期:  2024-09-24
  • 网络出版日期:  2024-11-04

目录

/

返回文章
返回
x 关闭 永久关闭