Characteristics of Geological Anomalies Shown on Seismic Time Sections
-
摘要:
三维地震勘探相关技术随着近年来的快速发展,取得了比较丰富的成果。利用三维地震勘探技术在解释地下煤层赋存形态及其构造方面具有较为明显的优势,尤其是在解释地质异常区以及圈定其范围方面取得较好的效果。本文通过地震时间剖面上不同的特征表现,对采空区、陷落柱、煤层露头、岩浆岩侵入、火烧区、煤层冲刷带以及煤层分叉等较为常见的地质异常体进行解释、圈定,得到比较理想的地质成果。
Abstract:With the rapid development of three-dimensional (3D) seismic exploration technology in recent years, its application has demonstrated significant advantages in interpreting the form and structure of underground coal seams, particularly in identifying geological anomalies and delineating their extent. Common geological anomalies such as goafs, collapse columns, coal seam outcrops, magmatic rock intrusions, fire zones, coal seam scouring zones, and coal seam bifurcations were identified through the distinct characteristics observed in seismic time sections, yielding highly accurate geological interpretations.
-
前列腺癌(prostate cancer,PCa)是男性最常见的恶性肿瘤之一,已成为了全球范围内影响男性健康问题的疾病。近年来,我国前列腺癌的发病率有逐年上升的趋势,并且已证实超过了欧美等国家[1]。血清前列腺特异性抗原(prostate specific antigen,PSA)是PCa常用的筛查工具。由于PSA为4~10 ng/mL时,PCa与前列腺增生(benign prostate hyperplasia,BPH)存在很大程度上的重叠,有文献报道[2]此区间 PCa的穿刺检出率仅为20%~30%,在鉴别诊断上存在诸多挑战。因此,如何无创性的提高PSA灰区(PSA值处于4~10 ng/mL时)前列腺癌的检出率,避免不必要的穿刺活检,已成为临床关注热点。
2019年欧洲泌尿生殖协会联合美国放射协会在对第二版前列腺影像报告与数据系统(prostate imaging-reporting and data system version 2.0,PI-RADS v2.0)进一步细化修订的基础上,推出了PI-RADS v2.1,提高了多参数磁共振(multiparametric magnetic resonance imaging,mpMRI)对PCa的诊断准确性[3]。但是对于PI-RADS 3分的病灶是否进行穿刺活检仍然没有明确的指导意见,临床决策陷入两难境地,有研究报道[4],在前列腺PI-RADS 3分病灶穿刺活检的阴性率达到了67%。前列腺特异性抗原密度(prostate specific antigen density,PSAD)被认为可以较好的为评估PI-RADS 3分患者的患癌风险提供参考依据[5]。
表观扩散系数(apparent diffusion coefficient,ADC)通过对水分子扩散的扩散程度进行定量分析,在前列腺癌的诊断中得到了较好的应用。由于ADC值会因扫描条件的不同而改变,近年来,有研究者提出了标准化ADC值(normalized ADC,ADCn)的概念[6],认为可以避免患者本身及扫描技术条件(如b值)对ADC绝对值产生的影响。鉴于PSA灰区时PI-RADS v2.1评分3分病灶诊断的不确定性,本研究旨在探讨ADCn、PSAD对良、恶性病变鉴别诊断价值。
1. 资料与方法
1.1 一般资料
回顾性收集2018年1月至2022年6月经连云港市第二人民院病理证实的58例前列腺疾病患者临床和MRI影像资料。
纳入标准:①患者术前均行完善的 MRI检查,检查序列包括磁共振常规序列T1WI、T2WI、动态对比增强磁共振成像(dynamic contrast-enhanced magnetic resonance imaging,DCE-MRI)、DWI,PI-RADS v2.1评分为3分;②病理结果在 MRI检查后2周内获得,并且具备明确的Gleason评分分级;③患者在 MRI检查前均未接受任何相关治疗;④PSA值为4~10 ng/mL。排除标准:①患者MRI检查不完善或者图像质量不佳无法进行PI-RADS评分及测量ADC值;②前列腺增生明显导致外周带显示不清;③患者检查前接受相关治疗;④临床、病理资料不完善。最终入组患者58例,年龄48~86岁,平均(67.8±8.6)岁。
患者及家属均知情同意,且签署知情同意书;本研究经连云港市第二人民医院医学伦理委员会批准,伦理审查号LW2023001。
1.2 检查设备和方法
患者术前所有MRI检查在3.0 T MRI机器上进行(SIEMENS syngo MRI E11)。检查前嘱患者尽量清洁肠道内容物,以避免肠腔内容物影响图像质量,适量饮水充盈膀胱。取仰卧位,选择腹部线圈。扫描的上下范围包括前列腺、精囊腺和盆腔结构。均行轴位、矢状位、冠状位高分辨T1WI、T2WI、高b值 DWI(b值取 0,1500 s/mm2)。
表1为扫描序列及参数。DCE动态扫描时肘静脉注射钆喷酸葡胺(Gd-DTPA)20 mL,流率3.0 mL/s,采用3D-VIBE序列。
表 1 前列腺MRI扫描具体参数Table 1. Specific parameters of prostate MRI scans扫描序列 TR/ms TE/ms Fov/mm 层厚/mm 矩阵 层间距/mm 横断位T1WI 850 20 250×220 4 224×380 1 横断位T2WI 4000 110 250×220 4 224×380 1 矢状位T2WI 3000 110 240×220 4 240×160 1 横断位DWI 6000 76 290×240 4 128×128 1 DCE-MRI 5.28 1.35 220×220 4 240×330 1 1.3 图像处理和数据分析
由两名10年以上工作经验的腹部MRI诊断副主任医师在不知道患者临床信息和病理结果的情况下分别对病灶进行独立诊断和定位,给出PI-RADS评分,若结果不一致时由两位医师分析讨论达成共同意见后再进行评分。参照高分辨T1WI、T2WI、DWI图像,在工作站自动生成的ADC图像上手动勾画感兴趣区(region of interest,ROI),测量ADC值。具体方法如下:
①在病变最可疑区域的横轴位最大层面选择感兴趣区(小于病变的范围,且避开出血、钙化、囊变区),并测量其ADC值,定义为ADC病灶。②选取同侧或对侧正常外周带组织设置感兴趣区,需避开病灶侵犯及中央带与外周带结合区域,测量ADC值,定义为ADC外周带。根据所得结果计算出标准化ADC值,即ADCn=ADC病灶/ADC外周带;按照上述方法分别测量3次,取其平均值。
以上所选ROI均分别取病灶或前列腺外周带3个不同部位,范围大致保持一致。在高分辨率T2WI图像上测量前列腺体积,范围为前列腺底部至尖部,轴位图像上测量最大左右径LR和最大前后径AP,矢状位图像上测量前列腺最长上下径CC,根据长椭圆形公式[7]得出相应的前列腺体积V(cm3),即V=CC×AP×LR×0.52,然后根据患者的血清tPSA值,计算PSAD=tPSA/V。
1.4 病理学资料分析
所有患者病理标本经前列腺癌根治术或穿刺获得,由经验丰富的病理医师进行常规病理检查和免疫组织化学检查,进行病理诊断并记录病理结果,最终由病理科主任医师对结果进行审核。
1.5 统计学分析
应用SPSS 25.0统计分析软件。采用K-S检验各观察指标是否符合正态分布。采用独立样本t检验与秩和检验分别对符合正态分布和非正态分布的计量资料进行分析,比较各观察指标在PCa和非PCa组间差异是否具有统计学意义(P<0.05为差异有统计学意义)。采用MedCalc软件绘制ADCn、PSAD及两者联合应用鉴别前列腺癌的受试者工作特征曲线(ROC),获取曲线下面积(AUC),在最佳诊断阈值下,计算相应的诊断灵敏度、特异度;采用Z检验比较各AUC之间差异。
2. 结果
2.1 入组观察指标及单因素统计分析结果
入组患者经病理结果确定为PCa 20例,非PCa 38例。年龄、PSA符合正态分布,ADCn、PSAD不符合正态分布。独立样本t检验和秩和检验结果显示ADCn、PSAD在PCa组和非PCa组间差异具有统计学意义,ADCn值:PCa组<非PCa组;PSAD值:PCa组>非PCa组;两组间年龄和PSA差异无统计学意义(表2)。
表 2 各变量在PCa组和非PCa组间差异比较Table 2. Comparison of the differences in variables between the PCa and non-PCA groups项目 组别 统计检验 PCa组 非PCa组 t/Z P 例数 20 38 − − 年龄/岁 67.60±10.78 67.84±7.37 -0.101 0.021 PSA/(ng/mL) 7.15±2.17 6.39±1.75 1.45 0.052 PSAD/(ng/mL·cm3) 0.28(0.23,0.31) 0.18(0.14,0.24) -3.89 <0.001 ADCn 0.52(0.47,0.57) 0.69(0.64,0.78) -4.36 <0.001 注:表中括号内数字为第 25 百分位数和第 75 百分位数。 2.2 ADCn、PSAD对PI-RADS 3分病灶PCa的诊断效能结果
分别绘制ADCn、PSAD单独及联合使用诊断PCa的ROC曲线(图1),并进行分析。当ADCn取最佳诊断阈值0.588时,AUC、敏感度、特异度为0.849(95%CI 0.730~0.929)、85.2%、81.6%;PSAD取最佳诊断阈值0.211 ng/mL·cm3时,AUC、敏感度、特异度为0.813(95%CI 0.688~0.903)、85.0%、78.4%,两者间AUC差异无统计学意义。ADCn+PSAD联合使用,最佳诊断阈值为0.373,诊断PCa的敏感度、特异度为90.0%和89.5%,AUC达到了0.962(95%CI 0.876~0.995),高于两者单独使用,且差异具有统计学意义(表3和表4)。ADCn、PSAD对两例PI-RADS v2.1评分3分病灶的诊断与病理结果对照(图2和图3)。
表 3 ADCn、PSAD及两者联合的灵敏度、特异度、最佳阈值及AUC的比较Table 3. Comparison of sensitivity, specificity, optimal threshold, and area under the AUC of ADCn, PSAD, and their combination观察变量 AUC 95%CI 灵敏度/% 特异度/% 约登指数 阈值 标准误差 P ADCn 0.849 0.730,0.929 85.2 81.6 0.666 0.588 0.051 <0.001 PSAD 0.813 0.688,0.903 85.0 78.4 0.534 0.211 0.055 <0.001 ADCn+PSAD 0.962 0.876,0.995 90.0 89.5 0.795 0.373 0.021 <0.001 表 4 ADCn、PSAD及两者联合的受试者工作特征曲线(ROC)比较Table 4. Comparison of ROCs of ADCn, PSAD, and their combination观察变量 ADCn PSAD ADCn+PSAD ADCn − 0.664/0.434 0.012/2.524 PSAD 0.664/0.434 − 0.004/2.863 ADCn+PSAD 0.012/2.524 0.004/2.863 − 注:表中数字为不同变量之间ROC比较P/Z的结果,前者为P值,后者为Z值。如0.012/2.524即为ADCn与ADCn+PSAD之间的AUC比较,P值为0.012,Z值为2.524。 3. 讨论
2019年PI-RADS委员会发布了多参数MRI的前列腺影像报告和数据系统版本2.1(prostate imaging reporting and data system,version 2.1,PI-RADS V2.1),对PI-RADS v2.0版的评分细则进行了修改,提高了敏感性和特异性,被诸多研究者用于前列腺癌的检测和危险分层[8-9]。
PI-RADS v2.1认为1~2分病灶患癌的可能性极低,建议随访观察;4~5分病灶则建议穿刺活检,以达到尽早的明确诊断。然而,对于3分病灶的良恶性判定仍然存在不确定性,并没有给出明确的指导意见,临床决策进退两难。PI-RADS v2.1评分为3分的病灶影像学表现多不典型,往往是在排除了其他评分后归为此类,国内外学者对于此类病灶中csPCa的检出率的报道也不尽相同。
Osses等[10]研究认为,PI-RADS 3分病灶的穿刺阳性率为32%,其中23% 为临床有意义PCa(clinically significant PCa,csPCa)。Schlenker等[11]研究发现,PI-RADS 3分病灶PCa及csPCa穿刺阳性率分别为26.8% 和14.6%。这也就提示我们,若对PI-RADS 3分病灶均采取穿刺活检,虽可避免一部分漏诊,但带来的是阳性率较低,以及穿刺后出血、感染等并发症的风险。因此,目前对于PI-RADS 3分尤其是处于PSA灰区的患者,如何发掘更多有应用潜力的鉴别诊断指标,减少过渡诊断和治疗,将会是临床持续关注的热点。
MRI扩散加权成像(diffusion weighted imaging,DWI)是目前唯一能够无创性反映活体组织水分子扩散情况的功能性成像技术手段,可以通过测量ADC值对肿瘤组织微观结构变化进行量化评估,消除了T2透射效应影响,从而反映组织或病变性质。
前列腺癌时,由于组织细胞异质性造成细胞排列紧密,细胞外间隙扭曲、变形,间质成分减少,水分子扩散受限,在DWI图像表现为高信号,其ADC值减低。而良性的前列腺增生结节多由丰富的肥大腺体细胞和间质组成,并且细胞分化良好,水分子自由扩散不受限,DWI图像多呈等或低信号,与正常前列腺信号差异并不明显。
根据v2.1的标准,PI-RADS评分3分的病灶定义为ADC图表现为显著低信号或高b值的 DWI上表现为显著高信号,二者不能兼有;T2WI上表现为边界不清、不均匀低信号区亦或边界不清、信号不均的结节。我们分析入组的PI-RADS 3分诊断为PCa患者的MRI参数特征,认为对于DWI表现出来的高信号需要结合ADC值综合考虑分析,对病灶给与定性诊断。已有研究证实,PCa患者在MRI图像中具有更高的DWI信号和更低的ADC值[12]。
前列腺的ADC值容易受MRI设备、场强、b值选择以及年龄和种族的影响[13],前列腺炎、前列腺增生等疾病状态的生理变化也使其产生偏差,因此需要一个相对标准化的指标来平衡绝对ADC值的差异性。Hermie等[14]研究认为,前列腺癌与对侧正常组织ADC比值具备预测PI-RADS 3分csPCa的潜力,敏感性为59%,特异性为88%。我们认为其选择对侧正常的前列腺组织作为参考是有困难的,因为前列腺增生结节在ADC图像上亦可表现为不均匀信号或者稍低信号,会给测量结果带来不确定性,又由于前列腺增生常发生于中央腺体区,因此我们决定选择外周带正常前列腺组织作为参考较为妥当(即ADCn=ADC病灶/ADC外周带)。
本研究结果显示,PI-RADS 3分病灶中前列腺癌组的ADCn小于非前列腺癌组,差异具有统计学意义,当取最佳诊断阈值利用0.588时,绘制的ROC曲线下面积AUC达到0.849,敏感度85.2% 和特异度为81.6%,显示出了较高的预测价值。这也说明ADCn在一定程度上消除了个体差异影响,能够更加真实的量化PI-RADS 3分病灶组织扩散程度,纠正不同疾病组织ADC值的重叠影响,增加不同个体间的可比性。
PSA作为PCa重要的筛查标志物,当大于10 ng/mL时常提示患癌风险,由于其具有器官特异性而非肿瘤特异性的特点,表现出来的敏感性和特异性并不满意。临床中一部分患者PSA处于“灰区”,即4~10 ng/mL时,很难通过检验指标鉴别PCa和BPH,从而避免不必要的穿刺活检[15]。PSAD作为校正后的检测参数,体现了PSA与前列腺体积之间的相关性,纠正了病灶体积对PSA的影响,对PCa具有较高的诊断价值,也被认为是PSA灰区较为理想的参考指标。Song等[16]研究认为,PSAD对预测PSA灰区PCa表现出了较高的应用价值,同时对于PI-RADS 3分患者的临床决策也具有良好的参考意义[17]。
本文将PSAD纳入研究中,结果显示PCa组的PSAD值明显高于非PCa组(0.28 ng/mL·cm3 vs. 0.18 ng/mL·cm3),这也提示在PSA灰区PI-RADS 3分病灶中PSAD值越高,患癌的风险越大,其诊断PCa的AUC为0.813,敏感度、特异度分别为85%和78.4%,其中特异度略低于国内学者吴静等[18]的研究结果(81.6%),分析原因可能为本研究纳入的患者PSA水平均位于4~10 ng/mL区间,PCa患病率相对较低,且样本量小容易放大个体化差异,一定程度上降低PSAD的诊断特异性,后续研究将会加大样本量进一步验证。有学者研究认为,PSAD能够提高PI-RADS 3分病灶PCa的检出率,建议最佳诊断阈值为0.15 ng/mL·cm3,低于本研究结果0.211 ng/mL·cm3。我们认为原因可能在于该研究组前列腺体积较大,以及就医不及时、PCa分期晚等造成PSA值偏高等有关[19]。
本研究结果提示,ADCn、PSAD对PI-RADS 3分病灶中PCa的检出具有较好的作用。目前,临床上为减少模棱两可的诊断结果,建议不应单一的采取某特定的标准,推荐应用多指标联合进行全面系统评估。本研究联合应用ADCn与PSAD,诊断的敏感度和特异度为90.0%和89.5%,均优于单一指标,其AUC更是达到了0.962,具有明显的优越性。我们认为,对于PSA灰区PI-RADS 3分病灶,联合应用ADCn和PSAD可以有效提高PCa的预测能力,可大幅提高诊断准确性,避免不必要的穿刺活检和过度治疗,具有极高的参考和应用价值。
本研究存在的不足:①本研究为回顾性分析设计,存在一定的选择偏倚风险;②本研究的入组样本量较少,可能会降低预测价值,后期需要增加样本量进一步论证;③本研究部分病理结果来自穿刺活检,同样存在漏诊风险。
综上所述,在对PSA灰区PI-RADS 3分病灶的临床决策仍然存在困难,单纯依靠常规影像学表现来诊断PCa是有局限性的,ADCn、PSAD可以有效提高对于此类患者的患癌风险评估,两者联合鉴别诊断效果达到了很好的水平,为术前无创性准确检出PCa提供了新的方法,值得临床推广应用。
-
-
[1] 陆基孟. 地震勘探原理(上册)[M]. 东营: 中国石油大学出版社, 2006. LU J M. Principles of seismic exploration (vol. 1)[M] . Dongying: China University of Petroleum Press, 2006. (in Chinese).
[2] 张向鹏, 聂荣花. 地震勘探技术在解释地质异常体的应用分析[J]. 工程地球物理学报, 2013, 10(4): 465-471. DOI: 10.3969/j.issn.1672-7940.2013.04.007. ZHANG X P, NIE R H. Application analysis of seismic exploration technology in interpreting geological anomaly bodies[J]. Engineering Geophysics, 2013, 10(4): 465-471. DOI: 10.3969/j.issn.1672-7940.2013.04.007. (in Chinese).
[3] 张向鹏. 特殊地质异常体在地震时间剖面的特征[J]. 科技视界, 2014, (2): 113-115. ZHANG X P. Characteristics of special geological anomalies in seismic time section[J]. Technology Horizons, 2014, (2): 113-115. (in Chinese).
[4] 衡亮. 地震属性技术的研究和应用[J]. 化工管理, 2017, (6): 131. DOI: 10.3969/j.issn.1008-4800.2017.06.118. HENG L. Research and application of seismic attribute technology[J]. Chemical Engineering Management, 2017, (6): 131. DOI: 10.3969/j.issn.1008-4800.2017.06.118. (in Chinese).
[5] 李振春, 张军华. 地震数据处理方法[M]. 东营: 中国石油大学出版社, 2004. LI Z C, ZHANG J H. Seismic data processing method[M]. Dongying: China University of Petroleum Press, 2004. (in Chinese).
[6] 耿晓兵. 三维地震勘探技术在解释采空区的应用分析[J]. 中国科技信息, 2022, (19): 98-100. DOI: 10.3969/j.issn.1001-8972.2022.19.zgkjxx202219032. GENG X B. Application analysis of 3D seismic exploration technology for interpreting goaf areas[J]. China Science and Technology Information, 2022, (19): 98-100. DOI: 10.3969/j.issn.1001-8972.2022.19.zgkjxx202219032. (in Chinese).
[7] 王宪峰. 煤田陷落柱的特征浅析[J]. 山西煤炭, 2011, 31(12): 69-70. DOI: 10.3969/j.issn.1672-5050.2011.12.026. WANG X F. Analysis of the characteristics of coal field collapse column[J]. Shanxi Coal, 2011, 31(12): 69-70. DOI: 10.3969/j.issn.1672-5050.2011.12.026. (in Chinese).
[8] 李劲松, 左杰海, 封云杰, 等. 三维地震勘探在赵庄煤矿陷落柱探查中的应用[J]. 煤矿开采, 2015, 20(5): 15-18. LI J S, ZUO J H, FENG Y J. Application of 3-d seismic exploration in detecting collapse column in Zhaozhuang Colliery[J]. Coal Mining Technology, 2015, 20(5): 15-18. (in Chinese).
[9] 张宏, 谢文伟. 地震勘探中煤层隐伏露头的解释精度[J]. 中国煤炭地质, 2016, 28(3): 65-69. DOI: 10.3969/j.issn.1674-1803.2016.03.13. ZHANG H, XIE W W. Interpretation accuracy of coal seam hidden outcrops in seismic exploration[J]. China Coal Geology, 2016, 28(3): 65-69. DOI: 10.3969/j.issn.1674-1803.2016.03.13. (in Chinese).
[10] 王树威. 三维地震勘探技术在新疆某矿区煤层露头预测中的应用[J]. 现代矿业, 2019, 35(8): 24-26. DOI: 10.3969/j.issn.1674-6082.2019.08.008. WANG S W. Application of 3D seismic exploration technology in coal seam outcrop prediction in a mining area in Xinjiang[J]. Modern Mining, 2019, 35(8): 24-26. DOI: 10.3969/j.issn.1674-6082.2019.08.008. (in Chinese).
[11] 陈鸣祺, 宋春华. 岩浆岩侵入对地震资料解释的影响分析[J]. 内蒙古煤炭经济, 2018, (6): 159-160. DOI: 10.3969/j.issn.1008-0155.2018.06.084. CHEN M Q, SONG C H. Analysis of the influence of magmatic intrusion on seismic data interpretation[J]. Inner Mongolia Coal Economy, 2018, (6): 159-160. DOI: 10.3969/j.issn.1008-0155.2018.06.084. (in Chinese).
[12] 单蕊, 张广忠, 王千遥, 等. 高密度三维地震数据驱动的煤层岩浆岩侵入区综合解释方法与应用[J]. 煤田地质与勘探, 2020, 48(6): 72-79. DOI: 10.3969/j.issn.1001-1986.2020.06.010. SHAN R, ZHANG G Z, WANG Q Y, et al. High-density 3D seismic data-driven comprehensive interpretation of magmatic intrusion zones in coal seams and their applications[J]. Coal Geology & Exploration, 2020, 48(6): 72-79. DOI: 10.3969/j.issn.1001-1986.2020.06.010. (in Chinese).
[13] 苏彦丁. 煤层火烧区探测中物探技术的应用与探索[J]. 勘察科学技术, 2016, (4): 62-64. DOI: 10.3969/j.issn.1001-3946.2016.04.017. SU Y D. Application and exploration of geophysical exploration technology in coal fire area detection[J]. Site Investigation Science and Technology, 2016, (4): 62-64. DOI: 10.3969/j.issn.1001-3946.2016.04.017. (in Chinese).
[14] 韩少明. 地震勘探资料解释煤层冲刷带的应用效果[J]. 华北国土资源, 2017, (3): 39-41. DOI: 10.3969/j.issn.1672-7487.2017.03.014. HAN S M. Application of seismic exploration data to interpretation of coalbed erosion zone[J]. Huabei Natural Resources, 2017, (3): 39-41. DOI: 10.3969/j.issn.1672-7487.2017.03.014. (in Chinese).
[15] 赵国辉. 张双楼矿区煤层分叉合并规律分析[J]. 化工管理, 2014, (18): 113. DOI: 10.3969/j.issn.1008-4800.2014.18.101. ZHAO G H. Analysis of coal seam bifurcation and merger law in the Zhangshuanglou mining area[J]. Chemical Engineering Management, 2014, (18): 113. DOI: 10.3969/j.issn.1008-4800.2014.18.101. (in Chinese).
[16] 刘磊, 李伟, 杜玉山, 等. 基于Stacking集成学习的分频地震属性融合储层预测方法[J]. 石油地球物理勘探, 2025, 59(1): 12-22. DOI: 10.13810/j.cnki.issn.1000-7210.2024.01.002. LIU L, LI W, DU Y S, et al. Frequency Division Seismic Attribute Fusion Reservoir Prediction Method Based on Stacking Ensemble Learning[J]. Oil Geophysical Prospecting, 2024, 59(1): 12-22. DOI: 10.13810/j.cnki.issn.1000-7210.2024.01.002. (in Chinese).