ISSN 1004-4140
CN 11-3017/P

磁法和电法联合探测在浅覆盖区寻找铜镍多金属矿的应用——以新疆东天山镜儿泉地区为例

贾世豪, 李勇, 龚胜平, 赵克强, 郭红宇

贾世豪, 李勇, 龚胜平, 等. 磁法和电法联合探测在浅覆盖区寻找铜镍多金属矿的应用——以新疆东天山镜儿泉地区为例[J]. CT理论与应用研究(中英文), 2025, 34(3): 392-400. DOI: 10.15953/j.ctta.2024.038.
引用本文: 贾世豪, 李勇, 龚胜平, 等. 磁法和电法联合探测在浅覆盖区寻找铜镍多金属矿的应用——以新疆东天山镜儿泉地区为例[J]. CT理论与应用研究(中英文), 2025, 34(3): 392-400. DOI: 10.15953/j.ctta.2024.038.
JIA S H, LI Y, GONG S P, et al. Application of Comprehensive Electromagnetic Methods in the Search for Copper-Nickel Polymetallic Ores in Shallow-Covered Areas: A Case Study of The Jing'erquan Area in the East Tianshan Mountains of Xinjiang[J]. CT Theory and Applications, 2025, 34(3): 392-400. DOI: 10.15953/j.ctta.2024.038. (in Chinese).
Citation: JIA S H, LI Y, GONG S P, et al. Application of Comprehensive Electromagnetic Methods in the Search for Copper-Nickel Polymetallic Ores in Shallow-Covered Areas: A Case Study of The Jing'erquan Area in the East Tianshan Mountains of Xinjiang[J]. CT Theory and Applications, 2025, 34(3): 392-400. DOI: 10.15953/j.ctta.2024.038. (in Chinese).

磁法和电法联合探测在浅覆盖区寻找铜镍多金属矿的应用——以新疆东天山镜儿泉地区为例

基金项目: 

国家重点研发计划(深部矿产资源勘查评价技术联合研究(2018YFE0208300));中国地质调查局(新疆东天山等重点成矿区带浅覆盖区铜镍多金属矿调查(DD20230309))。

详细信息
    作者简介:

    贾世豪,男,地球探测与信息技术专业硕士研究生,主要从事勘探地球物理电磁法应用研究,E-mail:jsh180121@126.com

    通讯作者:

    李勇✉,男,正高级工程师、硕士生导师,主要从事地球物理地面电磁系统设计、电磁法数值模拟与反演成像、方法技术应用研究,E-mail:cgslyong@mail.cgs.gov.cn

  • 中图分类号: P 631;P 618.2

Application of Comprehensive Electromagnetic Methods in the Search for Copper-Nickel Polymetallic Ores in Shallow-Covered Areas: A Case Study of The Jing'erquan Area in the East Tianshan Mountains of Xinjiang

  • 摘要:

    新疆东天山镜儿泉地区属于浅覆盖区,传统地表地质勘查手段难以准确判断矿体的存在。为克服这一难题,研究采纳磁法、大功率时间域激电法(TDIP)和可控源电磁测深法(CSEM)等先进技术手段进行地下勘查。勘探工作首先通过磁法勘查确定异常区域,进而利用TDIP法评估矿化范围和构造分布,最终通过CSEM法的电阻率断面反演技术,揭示地下构造特征和矿体的空间分布情况。基于综合地球物理资料处理解释,成功圈定2处磁异常区和2处激电异常区。这些成果不仅为镜儿泉地区后续的钻探工作提供科学依据,而且为其他具有类似地质特征地区的铜镍多金属矿床勘查提供宝贵的方法和经验参考。

    Abstract:

    The Jing'erquan area in the East Tianshan Mountains of Xinjiang is characterized by shallow coverage, posing challenges to surface geological exploration for ore bodies. This study utilized a combination of geophysical methods to overcome these challenges. Magnetic anomaly mapping was used initially to delineate areas of interest, followed by the application of the high-power time-domain induced polarization (TDIP) method to estimate mineralization extents and structural distributions within the mining area. Subsequently, the controllable source electromagnetic sounding method was employed to generate resistivity section maps, revealing subsurface structural characteristics and ore body occurrences. Comprehensive data processing and interpretation identified two magnetic anomalies and four IP anomalies. The results of these applications provide foundational support for future drilling projects within the mining area and offer valuable insights for workers in similar regions seeking copper–nickel polymetallic deposits.

  • 间质性肺异常(interstitial lung anomaly,ILA)是在轻度症状或轻微肺功能损害患者中偶然发现的,提示早期间质性肺疾病(interstitial lung disease,ILD)[1]。Fleischner协会提出,ILA是在无疑似ILD患者中发现与ILD相符的特异性CT表现,并建议对其进行评估[2]。研究表明,ILA很常见[3],且与死亡率增加、肺癌风险和癌症治疗相关并发症显著相关[4-5]。Cedars-Sinai等[6]通过对类风湿关节炎(rheumatoid arthritis,RA)患者进行纵向随访,发现基线时存在ILA的患者1年后出现影像学进展。因此,有必要对ILA进行早期诊断和适当管理。

    ILD是结缔组织病(connective tissue disease,CTD)最常见的肺部并发症[7],通常与不良结局和早期死亡率有关[8]。结缔组织疾病相关间质性肺疾病(connective tissue disease-associated interstitial lung disease,CTD-ILD)常表现为咳嗽等非特异性症状[9],通常在病程晚期被诊断[10],且与其他疾病重叠,导致ILD诊断常被延迟和低估[11],使患者预后恶化[12]。而抗纤维化治疗已被证明可降低ILD患者的死亡率和急性加重的风险[13],延缓肺纤维化进展[14]。因此,保持对ILD的高度怀疑,并尽早识别和治疗有助于改善患者的临床结局[15-16]

    目前,高分辨率计算机断层扫描(high-resolution CT,HRCT)是诊断ILD的金标准[17],也是筛查ILA的重要工具[2]。然而,肉眼识别肺区5%异常具有挑战性[18],且早期ILA影像表现与其他肺部病变(如重力依赖型肺不张、局灶胸椎旁纤维化等[19])相似。而定量CT可克服这些局限性,在ILA评估中发挥着重要作用[20]。该技术不仅能量化不同影像表现(如磨玻璃影、纤维化等)的占比,还能评估肺部整体损伤(如平均肺密度等)情况,并能早期发现病变,为疾病评估和早期诊断提供重要依据。

    然而,目前针对CTD患者的ILA的研究仍较少。本研究纳入CTD-ILA及CTD-ILD患者,通过定量CT评估患者肺损伤情况,探讨定量CT指标对于CTD-ILA、CTD-ILD的预测作用及诊断效能,建立基于定量CT的CTD患者筛查方法,为早期诊断与个性化管理提供参考依据。

    选取延安大学附属医院2018年1月至2024年11月CTD-ILA患者33例、CTD-ILD患者140例,对照组109例。

    CTD-ILA/ILD纳入标准:类风湿关节炎(RA)、系统性硬化症(systemic sclerosis,SSc)、皮肌炎(dermatomyositis,DM)、系统性红斑狼疮(systemic lupus erythematosus,SLE)、抗中性粒细胞胞浆抗体(anti-neutrophil cytoplasmic antibody,ANCA)相关血管炎的诊断分别符合2010年[21]、2013年[22]、2017年[23]、2019年[24]、2022年[25]美国风湿病学会指南,干燥综合征(sjögren’s syndrome,SS)、弥漫性结缔组织病(mixed connective tissue disease,MCTD)诊断符合2016年[26]、2019年[27]美欧共识小组指南,ILA的定义符合2020年Fleischner协会指南[2],ILD的定义符合2013年英国胸科协会指南[28]。ILA和ILD示例图如图1所示。

    图  1  CTD-ILA(a)和CTD-ILD(b)HRCT示例图
    注:(a)一名59岁系统性红斑狼疮女性患者,按照指南诊断ILA,患者尚无呼吸系统症状。(b)一名66岁类风湿关节炎男性患者,按照指南诊断ILD,患者存在呼吸系统症状。
    Figure  1.  Representative HRCT images of CTD-ILA (a) and CTD-ILD (b)

    CTD-ILA/ILD排除标准:无薄层高分辨率CT;胸部CT存在伪影;3D slicer无法分割;存在胸部肿瘤、感染[29-30]、肺气肿、胸腔积液;胸部手术史。对照组的纳入标准:符合CTD诊断但未患ILD的人群,其排除标准与CTD-ILD/ILD患者排除标准一致。

    使用128层螺旋CT扫描仪(上海联影UCT-760)进行检查,受试者取仰卧位,于深吸气末扫描,扫描范围为肺尖至肺底膈面。扫描参数:管电压120 kV、管电流自动调节,扫描层厚5 mm,使用1 mm标准算法重建。

    由两名经验丰富的放射科主任医生判断患者胸部CT质量并对其进行影像分型。定量CT指标使用3D-Slicer(5.6.2版本,http://www.slicer.org)获得。使用密度直方图和密度阈值法获得定量CT指标。具体步骤如图2所示。

    图  2  定量CT分析流程图
    Figure  2.  Flow diagram of quantitative CT analysis protocol

    密度阈值法的阈值设置:全肺像素阈值设定为 − 1024 HU至 − 200 HU,进一步细分为正常肺区(− 950 HU至 − 700 HU)、磨玻璃密度区(− 700 HU至 − 500 HU)和纤维化区(− 500 HU至 − 200 HU)。计算这些区域像素与全肺像素比例,得到正常肺区域(normal lung,NL)、磨玻璃密度区域(ground-glass opacity,GGO)和纤维化区域(fibrosis,F)百分比,以及异常病变区域(abnormal area,AA)的百分比(AA%为GGO%与F%之和),定量分割结果如图3所示。

    图  3  密度阈值分割法示例图
    注:(a)为CT横断位图,(b)为定量CT分割结果图。蓝色代表正常肺组织的百分比(NL%),橙色表示磨玻璃样混浊区域的百分比(GGO%),粉红色表示纤维化区域的百分比(F%),红色表示血管。磨玻璃样混浊区和纤维化区的百分比之和被设置为异常病变区的百分比(AA%)。
    Figure  3.  Schematic illustration of density threshold segmentation methodology

    密度直方图法通过统计不同密度值出现的频率生成灰度直方图,得到高衰减区(high attenuation area,HAA)、平均肺密度(mean lung density,MLD)、标准差(standard deviation,SD)、峰度(Kurtosis)和偏度(Skewness)参数。

    采用Shapiro-Wilk W进行正态检验。正态资料以(均数±标准差)表示,非正态资料以中位数(四分位数间距)表示,分类变量以频率(百分比(%))表示。组间差异用单因素方差分析或Kruskal-Wallis H检验分析,采用Bonferroni法进行事后比较。

    应用受试者操作特征曲线(receiver operator characteristic curve,ROC)评估肺定量CT指标在区分对照组与CTD-ILA、对照组与CTD-ILD、CTD-ILA与CTD-ILD的诊断效能,使用多元逐步Logistic回归筛选预测CTD-ILA及CTD-ILD的定量CT指标。

    使用SPSS 25.0进行统计分析,以P<0.05为差异具有统计学意义。

    本研究共纳入282名患者,包括对照组109例,CTD-ILD组140例,CTD-ILA组33例。3组患者在性别、年龄和CTD类型方面存在差异,而病程、BMI以及CT分型方面无显著差异(表1)。

    表  1  患者临床资料表
    Table  1.  Demographic and clinical characteristics of study participants
    项目 组别 统计检验
    对照组(n=109) CTD-ILA(n=33) CTD-ILD(n=140) F/H P
    性别(%) 14.52 0.001
      男 19(17.3) 14(42.4) 53(37.9)
      女 90(81.8) 19(57.6) 87(62.1)
    年龄/岁 40(19) 64(11) 62.5(14) 104.52 <0.001
    BMI/(kg/m2 22.54(2.81) 21.89(3.08) 22.34(4.39) 0.29 0.864
    病程 3(7) 5(11.25) 4(10.75) 3.25 0.187
    CTD类型(%) 9.66 0.008
      类风湿关节炎 63(57.3) 22(66.7) 67(47.9)
      系统性红斑狼疮 24(21.8) 5(15.2) 10(7.1)
      系统性硬化症 2(1.8) 2(6.1) 15(10.7)
      干燥综合症 9(8.2) 2(6.1) 16(11.4)
      皮肌炎 2(6.1) 8(5.7)
      弥漫性结缔组织病 7(6.4) −  15(10.7)
      ANCA相关血管炎 4(3.6) 2(6.1) 9(6.4)
    CT分型(%) 1.12 0.291
      普通型间质性肺炎 7(21.2) 41(29.3)
      非特异性间质性肺炎 18(54.5) 71(50.7)
      淋巴细胞性间质性肺炎 4(12.1) 22(15.7)
      机化性肺炎 4(12.1) 6(4.3)
    下载: 导出CSV 
    | 显示表格

    定量CT指标分析显示(表2),组间定量CT指标均存在统计学差异(P<0.05)。在密度阈值法中,从无ILD、ILA到ILD,NL%逐渐下降,而GGO%、F%、AA% 则逐渐上升。

    表  2  各组间定量CT指标差异
    Table  2.  Intergroup differences in quantitative CT metrics among controls, CTD-ILA, and CTD-ILD cohorts
    项目 组别 统计检验
    对照组(n=109) CTD-ILA(n=33) CTD-ILD(n=140) F/H P
    NL% 74(7) 68.0(7.5)* 65.5(9)* 89.82 <0.001
    GGO% 5.7(3.7) 11.3(8.75)* 12.65(9.1)* 76.00 <0.001
    F% 2.9(1.25) 5.2(3.4)* 6.6(4.98)* 113.72 <0.001
    AA% 8.8(4.9) 16.40(11.80)* 19.35(12.67)* 89.51 <0.001
    HAA 4.04(2.18) 9.25±5.26* 9.44(6.92)* 114.05 <0.001
    MLD −830.16(49.19) −777.08(59.02)* −768.19(93.69)* 70.18 <0.001
    SD 179.94(16.38) 205.52(25.83)* 220.19(41.43)* 139.60 <0.001
    Kurtosis 13.62±4.20 6.77(4.91)* 5.28(5.64)* 122.59 <0.001
    Skewness 3.27(0.71) 2.34(0.74)* 2.19±0.60* 115.41 <0.001
    注:*表示与对照组相比P<0.05。NL%为正常肺区域百分比;GGO%为磨玻璃密度区域百分比;F%为纤维化区域百分比;AA%为异常病变区域百分比;HAA为高衰减区;MLD为平均肺密度;SD为标准差;Kurtosis为峰值;Skewness为偏度。
    下载: 导出CSV 
    | 显示表格

    密度直方图法中,MLD、HAA、SD呈上升趋势,而Kurtosis和Skewness则逐渐减小。

    ROC曲线分析结果显示,SD、HAA、Kurtosis和F%在区分对照组与CTD-ILA(表3图4)、对照组与CTD-ILD(表4图5)方面表现优异。其中,SD的曲线下面积(AUC)最大,在区分对照组与CTD-ILA(AUC=0.862)和对照组与CTD-ILD(AUC=0.923)时均表现出较高的诊断效能。在密度阈值法中,F%区分两者AUC值最大。

    表  3  定量CT指标区分对照组与CTD-ILA的ROC曲线分析结果
    Table  3.  ROC curve analysis of quantitative CT metrics for discriminating control groups from CTD-ILA patients
    项目 AUC P 最佳截断值 灵敏度 特异度 约登指数
    NL% 0.770 <0.001 < 71.50 0.758 0.743 0.501
    GGO% 0.774 <0.001 > 8.000 0.788 0.734 0.522
    F% 0.814 <0.001 > 4.050 0.758 0.807 0.565
    AA% 0.785 <0.001 > 12.20 0.788 0.762 0.549
    HAA 0.830 <0.001 > 6.121 0.788 0.835 0.623
    MLD 0.763 <0.001 > −800.9 0.697 0.798 0.495
    SD 0.862 <0.001 > 191.2 0.849 0.817 0.665
    Kurtosis 0.821 <0.001 < 9.615 0.727 0.844 0.571
    Skewness 0.816 <0.001 < 2.795 0.788 0.817 0.604
    注:NL% 为正常肺区域百分比;GGO%为磨玻璃密度区域百分比;F%为纤维化区域百分比;AA%为异常病变区域百分比;HAA为高衰减区;MLD为平均肺密度;SD为标准差; Kurtosis为峰值; Skewness为偏度。
    下载: 导出CSV 
    | 显示表格
    图  4  定量CT指标区分对照组与CTD-ILA ROC曲线
    Figure  4.  Receiver operating characteristic (ROC) curves of quantitative CT metrics for discriminating between control group and CTD-ILA patients
    图  5  定量CT指标区分对照组与CTD-ILD ROC曲线
    Figure  5.  Receiver operating characteristic (ROC) curves of quantitative CT metrics for differentiating between control group from CTD-ILD patients

    在区分CTD-ILA与CTD-ILD的ROC曲线分析(表5图6),SD的诊断价值最佳(AUC=0.649),其次为Kurtosis(AUC=0.638),密度阈值法中F%的AUC最大为0.617。

    表  4  定量CT指标区分对照组与CTD-ILD的ROC曲线分析结果
    Table  4.  ROC curve analysis of quantitative CT metrics for differentiating control groups from CTD-ILD patients
    项目 AUC P 最佳截断值 灵敏度 特异度 约登指数
    NL% 0.844 <0.001 < 71.50 0.829 0.743 0.572
    GGO% 0.815 < 0.001 > 8.350 0.771 0.752 0.524
    F% 0.886 < 0.001 > 4.050 0.727 0.807 0.535
    AA% 0.843 < 0.001 > 11.30 0.857 0.716 0.573
    HAA 0.885 < 0.001 > 6.205 0.727 0.844 0.571
    MLD 0.801 < 0.001 > −798.2 0.636 0.826 0.462
    SD 0.923 < 0.001 > 191.2 0.864 0.817 0.680
    Kurtosis 0.896 < 0.001 < 11.74 0.818 0.706 0.525
    Skewness 0.889 < 0.001 < 2.795 0.727 0.817 0.544
    注:NL%为正常肺区域百分比;GGO%为磨玻璃密度区域百分比;F%为纤维化区域百分比;AA%为异常病变区域百分比;HAA为高衰减区;MLD为平均肺密度;SD为标准差; Kurtosis为峰值; Skewness为偏度。
    下载: 导出CSV 
    | 显示表格
    表  5  定量CT指标区分CTD-ILA与CTD-ILD的ROC曲线分析结果
    Table  5.  ROC curve analysis of quantitative CT metrics for differentiating CTD-ILA from CTD-ILD
    项目 AUC P 最佳截断值 灵敏度 特异度 约登指数
    NL% 0.600 0.074 < 63.50 0.379 0.788 0.167
    GGO% 0.544 0.435 > 12.10 0.536 0.576 0.112
    F% 0.617 0.037 > 5.450 0.650 0.576 0.226
    AA% 0.574 0.187 > 16.85 0.600 0.576 0.176
    HAA 0.585 0.128 > 10.06 0.486 0.727 0.213
    MLD 0.568 0.224 > −758.5 0.429 0.758 0.186
    SD 0.649 0.008 > 219.1 0.529 0.758 0.286
    Kurtosis 0.638 0.014 < 5.095 0.479 0.849 0.327
    Skewness 0.614 0.042 < 1.932 0.364 0.879 0.243
    注:NL%为正常肺区域百分比;GGO%为磨玻璃密度区域百分比;F%为纤维化区域百分比;AA%为异常病变区域百分比;HAA为高衰减区;MLD为平均肺密度;SD为标准差; Kurtosis为峰值; Skewness为偏度。
    下载: 导出CSV 
    | 显示表格
    图  6  定量CT指标区分CTD-ILA与CTD-ILD ROC曲线
    Figure  6.  Receiver operating characteristic (ROC) curves of quantitative CT metrics for discriminating between CTD-ILA and CTD-ILD

    将NL%、GGO%、F%、AA%、HAA、MLD、SD、Kurtosis、Skewness作为自变量,分组结果(ILA/ILD)为因变量,进行多元逐步logistic回归。

    结果显示(表6),在区分对照组与ILD及ILA时,F%(β=0.620,P=0.006,OR=1.895;β=0.854,P=0.000,OR=2.349)、GGO%(β=−0.302,P=0.002,OR=0.739;β=−0.454,P=0.000,OR=0.635)、SD(β=0.026,P=0.009,OR=1.016;β=0.016,P=0.000,OR=1.026)、Kurtosis(β=−0.370,P=0.000,OR=0.691;β=−0.439,P=0.000,OR=0.645)具有统计学差异。

    表  6  定量CT指标预测ILA与ILD多因素多元logistic回归结果
    Table  6.  Multivariate logistic regression analysis of quantitative CT metrics in predicting ILA vs. ILD
    项目 β BE wald P OR 95% CI
    下限 上限
    ILA F% 0.620 0.225 7.629 0.006 1.859 1.197 2.887
    GGO% −0.302 0.098 9.546 0.002 0.739 0.610 0.895
    SD 0.015 0.006 6.784 0.009 1.016 1.004 1.027
    Kurtosis −0.370 0.069 28.727 0.000 0.691 0.603 0.791
    ILD F% 0.854 0.209 16.666 0.000 2.349 1.559 3.540
    GGO% −0.454 0.092 24.385 0.000 0.635 0.530 0.760
    SD 0.026 0.005 23.855 0.000 1.026 1.016 1.037
    Kurtosis −0.439 0.062 50.962 0.000 0.645 0.571 0.727
    注:GGO%为磨玻璃密度区域百分比;F%为纤维化区域百分比;SD为标准差; Kurtosis为峰值。
    下载: 导出CSV 
    | 显示表格

    ILA及ILD均与CTD患者不良预后和早期死亡率密切相关,早期识别和管理对改善患者结局至关重要。但目前临床筛查仍依赖视觉评估,易导致早期病变漏诊。与既往研究相比,本研究系统量化了正常肺组织、到ILA、再到ILD的连续演变规律,为疾病发展提供客观证据;其次,构建基于定量CT的CTD患者筛查方法,为建立风险分层提供依据;最后,使用开源分析软件,免费、省时,有利于技术推广与普及。

    从无ILD的CTD对照组、CTD-ILA到CTD-ILD,NL%、峰度(Kurtosis)和偏度(Skewness)逐渐下降,而GGO%、F%、AA%、MLD、HAA以及SD则逐渐上升,与Ahn等[31]及杨凯等[32]的研究一致,表明ILA是ILD的早期阶段。

    定量CT指标的ROC曲线和Logistic回归分析表明,F%、GGO%、SD和Kurtosis是区分对照组与CTD-ILA/ILD的敏感指标。在密度阈值法,F%在区分对照组与CTD-ILA(AUC=0.814)和CTD-ILD(AUC=0.886)时表现最佳,提示其在早期诊断ILA和ILD方面具有重要意义。与Zhang等[33]研究结果一致,F%同样被证实为识别健康对照与特发性肺间质纤维化(idiopathic pulmonary fibrosis,IPF)的最佳指标(AUC=0.962)。但本研究F%的AUC值低于其在IPF中的表现,这提示炎症主导的CTD-ILD与纤维化驱动的IPF存在差异[34],表明定量指标筛查应结合患者病因。

    此外,磨玻璃密度区域百分比(GGO%)在诊断ILA及ILD方面同样具有重要价值(β1=−0.302,P=0.002,OR=0.739;β2=−0.454,P=0.000,OR=0.635)。这与马震忠等[35]发现皮肌炎/多发性肌炎相关ILD(DM/PM-ILD)患者HRCT表现以磨玻璃影为主(87.2%,123/141)一致,徐光兴等[36]证实GGO%在鉴别DM/PM-ILD与健康对照组时具有较高的诊断效能(AUC=0.91),但DM/PM-ILD患者的影像表现以非特异性间质性肺炎(NSIP)为主[9],而本研究还纳入了约30%以普通型间质性肺炎(usual interstitial pneumonia,UIP)表现为主的患者,这表明GGO%在鉴别不同亚型ILD时具有广泛适用性。

    密度直方图法中,标准差(SD)和峰度(Kurtosis)是区分正常组和CTD-ILA/ILD有意义的指标。SD的曲线下面积最大(AUC=0.862,AUC=0.923),Kurtosis区分ILD(AUC=0.896)仅次于SD(AUC=0.923),与Guisado-Vasco等[37]研究一致。研究表明SD还是区分局限性和弥漫性SSc-ILD和SS-ILD的最佳鉴别参数[38-39],进一步证明SD对于判断疾病严重程度也具有重要价值。

    而区分ILA与ILD的ROC曲线显示SD的AUC最大(AUC=0.649),SD>219.1为区分ILA和ILD的最佳截断值。在密度阈值法,F%的AUC最大(AUC=0.617),F%>5.450为区分ILA和ILD的最佳截断值。然而,基于定量CT区分ILD和ILA的结果仅供参考,实际应用还需结合患者临床症状、肺功能检查及随访结果进行综合评估。

    基于当前证据,我们提出CTD肺部评估的大致路径:初筛阶段采用SD、F%、GGO%、和Kurtosis排除高危人群;对于异常患者采用SD和F%判断患者是否患ILA或ILD。针对ILA患者,定期随访监测病情变化;而ILD患者则评估病情程度,尽早治疗改善预后。

    尽管高衰减区(HAA)未被纳入回归模型,但其在区分对照组与ILA患者时表现较好(AUC=0.830)。既往研究也表明,HAA在无症状类风湿关节炎(RA)识别ILD风险的能力最强[40],不同肺区域HAA与肺重塑的生物标志物、ILA的风险和全因死亡率相关[41],呼气HAA被认为是所有直方图分析指标中ILD严重程度的最佳预测因子[42],均提示HAA对于早期诊断的价值。此外,Shiraishi等[43]纵向评估COPD患者ILA,发现ILDvol%(磨玻璃影、网状影和蜂窝影的体积之和)与ILA 的发生相关,是识别和监测COPD患者ILA的可重复方法。这与我们的AA%指标类似,但AA%区分对照组与CTD-ILA/ILD AUC值较小,可能是3D slicer无法区分蜂窝影与肺气肿的局限性,未来需通过深度学习算法优化分割精度。

    本研究的局限性。①样本量较小,尤其是CTD-ILA组。由于前期研究重点是CTD-ILD,对ILA人群关注较少,但本研究仍为该领域提供了有价值的依据,未来将增加样本量,优化研究结果。②基线资料年龄和性别存在差异,但将其纳入回归分析后模型参数并未显著优化,因此未将其纳入最终分析。③3D slicer目前无法区分肺气肿与蜂窝影,因此纳入人群时已排除合并肺气肿的患者,可能存在一定的选择偏移。④研究主要聚焦于定量CT诊断价值,忽略了血清学指标、肺功能等的诊断意义,后续研究将纳入临床指标进行分析。

    综上所述,本研究不仅验证了定量CT在CTD-ILD/ILA早期诊断中的价值,还提出基于定量CT的CTD患者筛查方式,为患者精准管理提供有力支持。未来需在大样本中验证其适用性,并探索其与临床指标的协同作用。

  • 图  1   镜儿泉某测区地质图与工作布置

    注:1-第四纪冲积物;2-镜儿泉岩组;3-晚石炭世黑云母闪长岩;4-晚泥盆世二长花岗岩;5-晚泥盆糜棱岩化黑云母斜长花岗岩;6-辉长岩;7-花岗岩;8-地质界线;9-断裂;10-韧性剪切带断裂;11-片理及产状;12-磁法测点;13-激电测点;14-CSEM测点。

    Figure  1.   Geological map and work arrangement of the survey area in Jing'erquan

    图  2   磁异常化极与延拓等值线平面图

    Figure  2.   Contour maps of magnetic anomaly reduction to the pole and magnetic anomaly continuation

    图  3   激电中梯视电阻率与视充电率等值线平面图

    Figure  3.   Contour maps of the apparent resistivity and apparent chargeability of the TDIP survey

    图  4   90线综合电磁法勘查剖面及推断解释图

    Figure  4.   Integrated electrical survey profile and inferred interpretation map of line 90

    图  5   91线综合电磁法勘查剖面及推断解释图

    Figure  5.   Integrated electrical survey profile and inferred interpretation map of line 91

    表  1   岩石物性统计表

    Table  1   Petrophysical properties of the rock

    岩石名称 磁化率$ {K}/({10}^{-6}\mathrm{S}\mathrm{I}) $ 电阻率$ \rho / (\mathrm{\Omega }\cdot {\mathrm{m}}) $ 极化率$ \eta / \% $
    标本数 变化范围 常见值 标本数 变化范围 常见值 变化范围 常见值
    砂岩 50 295~3773 972 0.13~2.60 1.03
    凝灰质砂岩 186 0~3900 550
    灰岩 108
    凝灰岩 923 0~21800 950 35 552~4728 2210 0.51~2.14 1.30
    安山岩 3 459~859 630 928 352~1741 928 0.84~2.23 1.37
    火山角砾岩 20 0~2600 950 71 131~221 191
    二长花岗岩 39 13~522 186 25 175~1146 655 0.47~2.04 1.27
    花岗闪长岩 104 0~1191 50
    闪长岩 40 22~1504 207 46 269~4559 1560 0.49~2.72 1.54
    辉长石 127 20~2929 542 22 605~1279 941 0.82~2.10 1.35
    铜镍矿化及矿体 20 560~14200 10500 9 3~7250 420 2.3~35 6.00
    下载: 导出CSV
  • [1]

    QIN K Z. Overview of major Au, Cu, Ni and Fe deposits and metallogenic evolution of the eastern Tianshan Mountains, Northwestern China[J]. Tectonic Evolution and Metallogeny of the Chinese Altay and Tianshan, 2003, 227-248.

    [2]

    SONG X Y, CHEN L M, DENG Y F, et al. Syncollisional tholeiitic magmatism induced by asthenosphere upwelling owing to slab detachment at the southern margin of the Central Asian Orogenic Belt[J]. Journal of the Geological Society, 2013, 170(6): 941-950. DOI: 10.1144/jgs2012-130.

    [3] 崔敏利, 张宝林, 梁光河, 等. 黄土覆盖区钼矿综合地球物理找矿技术组合: 以沙坡岭钼矿为例[J]. 地球物理学进展, 2010, 25(2): 602-611. DOI: 10.3969/j.issn.1004-2903.2010.02.033.

    CUI M L, ZHANG B L, LIANG G H, et al. The technical combination of comprehensive geophysical prospecting in the Molybdenum mines with loess-covered: A case study at the Shapoling molybdenum mine[J]. Progress in Geophysics, 2010, 25(2): 602-611. DOI: 10.3969/j.issn.1004-2903.2010.02.033. (in Chinese).

    [4] 王钦军, 魏永明, 陈玉, 等. 低植被覆盖区斑岩铜矿遥感找矿模型及其应用——以环巴尔喀什–西准噶尔成矿带为例[J]. 地质学报, 2017, 91(2): 400-410.

    WANG Q J, WEI Y M, CHEN Y, et al. Remote sensing porphyry copper deposit exploration model and its application in low vegetated areas: An example from the taken Balkhash-western Junggar metallogenic zones[J]. Acta Geologica Sinica, 2017, 91(2): 400-410. (in Chinese).

    [5] 吴新刚, 刘福胜. 综合电磁法在内蒙某铁锌矿床勘查中的应用[J]. 中国煤炭地质, 2015, 27(12): 71-75.

    WU X G, LIU F S. Application of integrated electromagnetic method in an Iron-Zinc deposit prospecting in inner Mongolia[J]. Coal Geology of China, 2015, 27(12): 71-75. (in Chinese).

    [6] 吴新刚, 陆桂福, 杨亚斌. 铁锌矿床物探异常特征及勘查实例[J]. 物探化探计算技术, 2016, 38(3): 347-352.

    WU X G, LU G F, YANG Y B. The anomaly characteristics of the geophysical prospecting and exploration examples for iron and zinc ore bed[J]. Journal Computing Techniques For Geophysical and Geochemical, 2016, 38(3): 347-352. (in Chinese).

    [7] 王振亮, 邓友茂, 孟银生, 等. 综合物探方法在维拉斯托铜多金属矿床北侧寻找隐伏矿体的应用[J]. 物探与化探, 2019, 43(5): 958-965.

    WANG Z L, DENG Y M, MENG Y S, et al. The application of integrated geophysical prospecting method to the prospecting for concealed orebodies in the northern area of the Weilasituo copper polymetallic deposit[J]. Geophysical and Geochemical Exploration, 2019, 43(5): 958-965. (in Chinese).

    [8] 王振亮, 邓友茂, 林天亮, 等. 综合物探方法在内蒙古巴彦乌拉铜多金属矿中的应用[J]. 矿产勘查, 2018, 9(11): 2151-2158.

    WANG Z L, DENG Y M, LIN T L, et al. Application of comprehensive geophysical method for exploration on Bayanwulacopper polymetallic depositin Inner Mongolia[J]. Mineral Exploration, 2018, 9(11): 2151-2158. (in Chinese).

    [9] 严加永, 孟贵祥, 吕庆田, 等. 综合地球物理在荒漠覆盖区隐伏矿床预测与定位中的应用: 以新疆拉伊克勒克铜多金属矿床为例[J]. 地球物理学报, 2021, 64(11): 4117-4133. DOI: 10.6038/cjg2021P0054.

    YAN J Y, MENG G X, LV Q T, et al. Prediction and location of concealed deposits in desertgobi coverage areas using integrated geophysics: An example of the Layikeleke copper polymetallic deposit in Xinjiang, Northwest China[J]. Chinese Journal of Geophysics, 2021, 64(11): 4117-4133. DOI: 10.6038/cjg2021P0054. (in Chinese).

    [10] 龚胜平, 陆桂福, 席明杰, 等. 干旱荒漠区综合物化探方法寻找铜多金属矿[J]. 物探与化探, 2021, 45(1): 1-10.

    GONG S P, LU G F, XI M J, et al. The application of integrated geophysical and geochemical methods to the prospecting of copper polymetallic deposits in the arid desert area[J]. Geophysical and Geochemical Exploration, 2021, 45(1): 1-10. (in Chinese).

    [11] 林方丽, 王光杰, 杨晓勇. 综合电磁法在矿区深部成矿机制中的应用研究——以皖南乌溪多金属矿区为例[J]. 地球物理学报, 2016, 59(11): 4323-4337. DOI: 10.6038/cjg20161132.

    LIN F L, WANG G J, YANG X Y. Application of comprehensive electromagnetic study in deep mineralization mechanism: A case study of the Wuxi polymetallic ore deposit, south Anhui[J]. Chinese Journal of Geophysics, 2016, 59(11): 4323-4337. DOI: 10.6038/cjg20161132. (in Chinese).

    [12] 李建华, 林品荣, 杜炳锐, 等. 综合电法在隐伏矿产勘查中的应用研究[J]. 物探化探计算技术, 2022, 44(4): 468-476. DOI: 10.3969/j.issn.1001-1749.2022.04.09.

    LI J H, LING P R, DU B R, et al. The application study of integrated electrical method to the prospecting for concealed ore deposits[J]. Computing Techniques for Geophysical and Geochemical Exploration, 2022, 44(4): 468-476. DOI: 10.3969/j.issn.1001-1749.2022.04.09. (in Chinese).

    [13] 赵克强, 王振亮, 杨剑洲, 等. 黄山−镜儿泉成矿带东段铜镍矿成矿模式及浅覆盖区找矿预测初探[J]. 新疆地质, 2023, 41(S1): 53.
    [14] 邓宇峰, 宋谢炎, 颉炜, 等. 黄山–镜儿泉铜镍成矿带地层时代的厘定及其地质意义探讨[J]. 地质学报, 2021, 95(2): 362-376.

    DENG Y F, SONG X Y, JIE W, et al. Determination of sedimentary ages of strata in the Huangshan-Jingerquan mineralization belt and its geological significance[J]. Acta Geologica Sinica, 2021, 95(2): 362-376. (in Chinese).

    [15] 三金柱. 黄山−镜儿泉铜镍矿带区域成矿规律探讨——以图拉尔根铜镍矿为例[J]. 西北地质, 2012, 45(4): 176-184.

    SAN J Z. An approach to regional metallogenic regularities of Cu-Ni ore belt in Huangshan-Jinger’quan in Xinjiang, China: A case for Tulaergen Cu-Ni deposits[J]. Northwestern Geology, 2012, 45(4): 176-184. (in Chinese).

    [16] 郝泽江, 张强, 龚胜平, 等. 东天山镜儿泉地区地壳电性结构特征及其地质含义[J]. 地质学报, 2025, 99(2): 417-427. DOI: 10.19762/j.cnki.dizhixuebao.2023186.

    HAO Z J, ZHANG Q, GONG S P, et al. Electrical structure of the crust in the Jingerquan region, eastern Tianshan orogenic belt and its geological implications[J]. Acta Geologica Sinica, 2025, 99(2): 417-427. DOI:10.19762/j.cnki.dizhixuebao.2023186. (in Chinese).

    [17] 宋谢炎, 邓宇峰, 颉炜, 等. 新疆黄山–镜儿泉铜镍硫化物成矿带岩浆通道成矿特征及其找矿意义[J]. 矿床地质, 2022, 41(6): 1108-1123. DOI: 10.16111/j.0258-7106.2022.06.002.

    SONG X Y, DENG Y F, XIE W, et al. Ore-forming processes in magma plumbing systems and significances for prospecting of Huangshan-Jinger’quan Ni-Cu sulfide metallogenetic belt, Xinjiang, NW China[J]. Mineral Deposits, 2022, 41(6): 1108-1123. DOI: 10.16111/j.0258-7106.2022.06.002. (in Chinese).

    [18] 娄德波, 刘欢, 张长青, 等. 新疆黄山–镜儿泉地区与铜镍硫化物矿床有关的元素组合特征及找矿方向[J]. 大地构造与成矿学, 2017, 41(1): 133-144. DOI: 10.16539/j.ddgzyckx.2017.01.011.

    LOU D B, LIU H, ZHANG C Q, et al. Element assemblage characteristics of magmatic Cu-Ni sulfide deposits in Huangshan-Jing’erquan Area, Xinjiang and its prospecting significance[J]. Geotectonica et Metallogenia, 2017, 41(1): 133-144. DOI: 10.16539/j.ddgzyckx.2017.01.011. (in Chinese).

    [19] 赵云, 杨永强, 柯君君. 含铜镍岩浆起源及硫饱和机制: 以新疆黄山南岩浆铜镍硫化物矿床Sr-Nd-Pb-S同位素和元素地球化学研究为例[J]. 岩石学报, 2016, 32(7): 2086-2098.

    ZHAO Y, YANG Y Q, KE J J. Origin of Cu- and Ni-bearing magma and sulfide saturation mechanism: A case study of Sr-Nd-Pb-S isotopic composition and element geochemistry on the Huangshannan magmatic Ni-Cu sulfide deposit, Xinjiang[J]. Acta Petrologica Sinica, 2016, 32(7): 2086-2098. (in Chinese).

    [20] 娄德波, 肖克炎, 丁建华. 黄山−镜儿泉镍铜成矿带地球物理、地球化学特征及找矿意义[J]. 地质学刊, 2013, 37(3): 372-377.

    LOU D B, XIAO K Y, DING J J. On geophysical and geochemical characteristics of Huangshan-Jing'erquan Ni-Cu belt and its ore prospecting significance[J]. Journal of Geology, 2013, 37(3): 372-377. (in Chinese).

图(5)  /  表(1)
计量
  • 文章访问数:  129
  • HTML全文浏览量:  27
  • PDF下载量:  22
  • 被引次数: 0
出版历程
  • 收稿日期:  2024-03-07
  • 修回日期:  2024-05-21
  • 录用日期:  2024-07-01
  • 网络出版日期:  2024-08-14
  • 刊出日期:  2025-05-04

目录

/

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