Combined Application of High-density Electrical Method and Transient Electromagnetic Method in Gobi Desert Area
-
摘要: 为准确掌握地下水资源的动态变化，有效支撑生态环境的可持续发展，在新疆哈密某矿集区利用高密度电法与瞬变电磁法联合进行物探勘查，对地下含水层与隔水层进行预判。经验证，高密度电法与瞬变电磁法的组合，有效指导水文钻探施工，为同类地质条件下物探手段的选取和应用提供参考。Abstract: This study probed the existence of groundwater in a mineral concentration area, using high-density electrical and transient electromagnetic methods. Our objective was to understand the dynamic changes in groundwater resources accurately and support the development of a sustainable ecological environment in the Gobi Desert area. The study demonstrates that the combination of the high-density electrical and transient electromagnetic methods guides the hydrological drilling construction in this area effectively. It also provides a reference for the selection and application of geophysical exploration methods under similar geological conditions for future studies.
-
伦琴于1895年发现X射线，而后英国工程师Hounsfield基于X射线发明出计算机断层成像（computed tomography，CT）机。自CT机问世几十年来，无论是在硬件性能（X光机、探测器）、扫描方式，还是重建算法上都获得了长足发展，在医疗、工业和安全检查等领域发挥着重要作用。更短的扫描和重建时间，以及更好的成像质量成为该领域研究人员追求的目标[1]。
双能CT是由Alvarez[2]在1976年首次提出，可以获取物体的等效原子序数和电子密度信息，从而更加准确地识别物质。多能谱CT是在双能CT的基础上，增加了参与扫描和重建的能区数量，引入了更多能谱信息，重建结果更加准确。
随着高计数率电子元器件的迅速发展，以光子计数探测器为核心的多能谱X射线成像技术成为射线检测领域的重要发展方向[3-4]。与传统能量积分探测器相比，光子计数探测器可以通过设置能量阈值来选择不同能区的信号，同时它将每个光子当作独立事件进行统计，一次扫描便能得到多段能谱，减小了辐射剂量[5]。但在实际的应用上也存在多能段光子数减少、信噪比下降，探测器性能不足造成的噪声等问题[6]，对影响光子计数X射线CT成像精度的因素及其规律进行基础性研究十分必要。
本文对系统能谱的获取、模体的扫描和重建过程进行仿真实验，以重建得到的等效原子序数和电子密度的标准差作为成像质量的评价标准，得出探测器能量分辨率、能谱阈值和能谱范围等因素对成像精度影响的有关结论，为后续相关技术的工程应用化提供参考。
1. 双能重建基本原理
本文主要研究探测器能量分辨率、能谱阈值和能谱范围等对重建图像精度的影响规律。以双能重建为例，根据光子计数探测原理[7]，进行双能成像的模拟。双能CT成像的重点研究内容之一是双能分解方法，常见的双能分解方法包含投影域前处理、图像域后处理和迭代方法[8]。
考虑到未来的工程化需求，本文选择的双能重建方法为基于基材料模型的双能投影预处理分解方法，基材料选取的是碳和铝。
材料模型公式：
$$ \mu \left(E\right)={{b}}_{1}{\mu }_{1}\left(E\right)+{{b}}_{2}{\mu }_{2}\left(E\right) \text{，} $$ (1) $ {\mu }_{1} $ (E)、$ {\mu }_{2} $ (E)分别为两种基材料的线性衰减系数；$ {b}_{1}$ 、$ {b}_{2}$ 对应两种基材料的分解系数，对于某一固定的物质，$ {b}_{1}$ 、$ {b}_{2}$ 是两个常数。双能分解即求解如下的双能投影积分方程组：
$$ \left\{\begin{aligned} {P}_{L}=\;& -{\ln}\bigg(\int \mathit{{S}_{L}} (E)\exp \Big({-{B}}_{1}{\mu }_{1}({E})-\\& {{B}}_{2}{\mu }_{2}({E})\Big){\rm d}{E}\bigg)+{\ln}\int \mathit{{S}_{L}}(E){\rm d}E\\ {P}_{H}=\;&-{\ln}\bigg(\int \mathit{{S}_{H}}(E) \exp \Big(-{{B}}_{1}{\mu }_{1}({E})-\\ & {{B}}_{2}{\mu }_{2}({E})\Big){\rm d}{E}\bigg)+{\ln}\int \mathit{{S}_{H}}(E){\rm d}E\end{aligned}\right. \text{，} $$ (2) 其中，SH、SL分别为高低能能谱，纵坐标记录的是光子数，PH、PL分别为高低能投影，
$ {{B}}_{1}=\displaystyle \int {{b}}_{1}\mathrm{d}{l}， {{B}}_{2}=\displaystyle \int {{b}}_{2}\mathrm{d}{l} $ 。求解
$ {B}_{1}$ 、$ {B}_{2}$ 的这个过程被称为投影分解过程。由于$ {B}_{1}$ 、$ {B}_{2}$ 为$ {b}_{1}$ 、$ {b}_{2}$ 的线积分投影值，求解出$ {B}_{1}$ 、$ {B}_{2}$ 后，根据CT重建的原理，利用滤波反投影重建算法[9]，便可计算出$ {b}_{1}$ 和$ {b}_{2}$ ，由此可以计算材质的等效原子序数和电子密度信息，以完成材料的探测识别。计算公式：$$ {Z}_{\rm{eff}}={\left(\frac{{b}_{1}{\rho }_{e1}{Z}_{1}^{n}+{b}_{2}{\rho }_{e2}{Z}_{2}^{n}}{{b}_{1}{\rho }_{e1}+{b}_{2}{\rho }_{e2}}\right)}^{\tfrac{1}{n}} \text{，} $$ (3) $$ {\rho }_{e}={{b}_{1}\rho }_{e1}+{{b}_{2}\rho }_{e2} ， $$ (4) 式（3）和式（4）中Z1、Z2分别为两种基材料的原子序数，ρe1、ρe2分别为两种基材料的电子密度，n一般取在3~4之间。
2. 光子计数CT成像系统能谱仿真
在Ubuntu系统下利用Geant 4软件，模拟高速电子撞击靶材射出X射线的过程。真空管内靶材为1 cm厚的钨，靶面倾角为25°，射线源射出能量为160 keV的电子，滤波片为2.5 mm厚的铝，射线穿过区域皆为真空，光子穿过虚拟探测器后被记录能量，得到射线源入射能谱。图1为能谱发射的几何示意，图2为得到的入射能谱。
通过对入射能谱与基于能量分辨率高斯函数的卷积来模拟探测器能量分辨率的影响，得到系统能谱[10]。后续实验所需的高低能能谱均基于系统能谱得到。探测器能量分辨率R通常用半高全宽（full width at half maximum，FWHM）W表示，定义：
$$ R=\frac{W}{E}\text{，} $$ (5) 其中，E是相对应的能量。W可以通过标准差σ与高斯分布的关系计算得到：
$$ W=2.355\;\sigma 。 $$ (6) 因此，标准差可以表示为：
$$ \sigma =\frac{R\times E}{2.355} 。 $$ (7) 根据探测器指定能量下的标准能量分辨率R0和能量E0，可由式（8）获得其他能量下的能量分辨率。由式（7）获得每个能量点下的标准差σ，并使用该标准差对原始入射能谱进行高斯模糊[11]。
$$ R={R}_{0}\sqrt{\frac{{E}_{0}}{E}} 。 $$ (8) 3. 投影数据加噪仿真
由于光子计数的特性遵循泊松分布[12-13]，所以本文为投影数据添加泊松噪声，用于模拟真实环境。泊松分布的概率函数见式（9），该函数表示在给定时间间隔或空间区域内发生k次事件的概率，其中，
$ \lambda $ 是泊松分布的均值和方差。$$ P\Big(X=k\Big)=\frac{{\lambda }^{k}{\exp}{(-\lambda )}}{k!} 。 $$ (9) 具体加噪过程。
（1）光子穿过模体后的衰减规律：
$$ {n}=\int S\left(E\right){\exp}{\left(\int -l\mu \left(E\right)\mathrm{d}l\right)}\mathrm{d}E \text{，} $$ (10) 其中，n为入射后探测器接收到的光子数，S(E)为系统能谱，exp为自然常数，l表示射线穿过模体的厚度，μ(E)为水模的线性衰减系数。
将各个单能（E1、E2、E3、···）下的光子数ni=
$ S\left({E}_{\mathrm{i}}\right){\exp}{\left(\int -l\mu \left({E}_{\mathrm{i}}\right)\mathrm{d}l\right)} $ 分别代入式（9）中的$ \lambda $ ，将对应泊松分布下生成的随机数作为该能量下加噪后的光子数ni_noise。（2）将能谱各单能段加噪后的光子数ni_noise累加，并与入射前探测器接收到的光子数n0按下式进行变换，便可得到加噪后的投影Pnoise：
$$ {P}_{\mathrm{n}\mathrm{o}\mathrm{i}\mathrm{s}\mathrm{e}}=\mathrm{l}\mathrm{n}\left(\frac{{n}_{0}}{{n}_{\mathrm{n}\mathrm{o}\mathrm{i}\mathrm{s}\mathrm{e}}}\right) \text{，} $$ (11) 其中，入射前探测器接收到的光子数
$ {n}_{0}= \displaystyle \int S\left(E\right)\mathrm{d}E $ 。4. 实验仿真和结果分析
4.1 实验模体和成像几何
设置的仿真模体为一圆柱状物体，材料为钙（原子序数20，密度1.55 g/cm3），半径为3 cm，模体截面如图3所示。
进行扫描仿真实验的探测系统几何结构如图4所示。射线源为扇形束，扇束张角为22.6°；探测器线性等距排列，共256个探测单元，总长512 mm；射线源、模体、探测器中心位于同一条直线；射线源到模体中心和探测器中心的距离分别为640 mm和
1280 mm；每隔1° 扫描一次模体，每次扫描发射的光子数为105。4.2 能量分辨率的影响
市场上性能相对较好的光子计数探测器材料为碲锌镉（CZT），其在60~100 keV下的能量分辨率大约是10%。因此参考碲锌镉材料的分辨率作为最优值，从10% @ 100 keV到40% @ 100 keV依次递减能量分辨率，间隔10%取4组能量分辨率进行实验，并配以理想能量分辨率（0%）作为对比，高低能能量阈值选为80 keV，各实验组参数见表1。各实验组高斯模糊化之后的高低能能谱见图5。
表 1 能量分辨率实验对照表Table 1. Experimental comparison table of energy resolution编号 能量分辨率 低能区光子数 高能区光子数 1 0% @ 100 keV 82978 17286 2 10% @ 100 keV 82964 17271 3 20% @ 100 keV 82964 17218 4 30% @ 100 keV 82963 17132 5 40% @ 100 keV 82932 17005 穿过模体的光子由所设能量阈值被划入不同的能量通道，从而得到高低能投影数据。为投影数据添加泊松噪声后，进行双能重建[14-15]。
上述实验组的成像结果见表2。不同能量分辨率下的重建结果标准差曲线如图6所示。
表 2 能量分辨率实验成像结果表Table 2. Imaging results of energy resolution experiment编
号能量
分辨率等效原子
序数等效原子
序数标准差电子密度 电子密
度标准差1 0% @
100 keV19.538647 0.268186 1.518764 0.039321 2 10% @
100 keV19.549584 0.291602 1.517544 0.042370 3 20% @
100 keV19.544331 0.295689 1.518037 0.042150 4 30% @
100 keV19.544904 0.313434 1.518503 0.044175 5 40% @
100 keV19.570808 0.350212 1.515818 0.048496 ![]()
图 6 不同能量分辨率下的成像结果对比Figure 6. Comparison of imaging results at different energy resolutions上述结果表明，随着探测器能量分辨率的提高，成像结果受噪声影响程度下降，等效原子序数相较电子密度而言，受噪声影响程度更大。不难分析，光子计数探测器能量分辨率越低，表示其对目标能区光子的统计区域越大，这会使得能谱重叠范围增加，不利于能量分解。但研究发现，能量分辨率在30%甚至40%的情况下，密度和原子序数的标准差变化并不十分明显。
图7是不同探测器能量分辨率的重建结果对比图，电子密度图的灰度显示窗口为[1.44 1.60]，等效原子序数图的灰度显示窗口为[18.3 20.7]。图8为同一排像素灰度曲线比较。
![]()
图 7 不同能量分辨率的重建结果对比Figure 7. Comparison of reconstruction results at different energy resolutions4.3 能量阈值及能量范围的影响
本部分考察的是高低能能谱重合程度对于成像的影响，在探测器能量分辨率相同的条件下，通过选取不同能量阈值，改变能量范围，可以得到重合度不一的高低能能谱，重合度的定义为：
$$ {\text{重合度}} = \frac{\text{能谱重叠部分面积}}{\text{全能谱范围面积}}。 $$ (12) 探测器能量分辨率选用30% @ 100 keV，具体对照实验见表3。各实验组能谱对比见图9。
表 3 能量范围实验对照表Table 3. Experimental comparison table of energy range编号 低能范围EL 高能范围EH 重合度Overlap_ratio 低能段光子数 高能段光子数 1 0~70 keV 110~160 keV 0.0025 74978 4419 2 0~80 keV 100~160 keV 0.0134 82323 7438 3 0~80 keV 80~160 keV 0.0648 82963 17132 4 0~100 keV 60~160 keV 0.1526 62561 37533 实验结果见表4，图10为重建结果的均值变化。图11为重建切片图，电子密度图的灰度显示窗口为[1.41 1.61]，等效原子序数图的灰度显示窗口为[18.6 20.5]。图12为同一排像素的灰度比较。
表 4 能量范围实验结果Table 4. Results of energy range experiment编号 重合度Overlap_ratio 等效原子序数 等效原子序数标准差 电子密度 电子密度标准差 1 0.0025 19.550129 0.316059 1.517544 0.049246 2 0.0134 19.551358 0.298083 1.517176 0.043675 3 0.0648 19.544904 0.313434 1.518503 0.044175 4 0.1526 19.554059 0.409670 1.517345 0.055390 ![]()
图 10 不同能谱重合度下的成像结果对比Figure 10. Comparison of imaging results at different energy overlap-ratios![]()
图 11 不同能量范围的重建结果对比Figure 11. Comparison of reconstruction results within different energy ranges从结果可以看出，设置能谱阈值改变能谱范围，使得高低能能谱重合度发生改变，会对成像质量造成影响。盲目减小高低能能谱重合度，并不一定使重建结果更好，因为当光子数减少时，噪声会增大。如何平衡能谱重合度和噪声是下一步的研究工作。
5. 小结
光子计数探测器作为多能谱成像的关键器件，成为研究能谱成像的一个重要内容。本文通过实验，模拟光子计数探测器在不同能量分辨率、不同能谱阈值和能量范围下的成像实验。通过比较重建的等效原子序数图和电子密度图的标准差，评价成像结果的好坏，得到上述因素对光子计数X射线CT成像精度的影响规律。
（1）针对安检CT中常用的160 kV射线源能量，光子计数探测器能量分辨率越高，成像效果越好，但是能量分辨率对于材料的分辨差异并不显著。
（2）在探测器能量分辨率确定的情况下，恰当选取能谱阈值对于成像很有帮助。通过设置合适的阈值，使高低能能谱重合度减小，往往可以得到更好的成像结果。但是随着能谱变窄、重合度下降，噪声也会加大，需要平衡这两种因素，否则更低的能谱重合度也有可能带来更差的成像结果。
本文的研究成果将有助于能谱CT系统的总体设计，如从双能成像材料分辨的角度，并不一定过高追求能量分辨率，这对于能量成像CT系统的硬件选型具有积极的指导意义。
-
[1] 王瑞丰, 温来福, 程久龙, 等. 高密度电法与瞬变电磁法联合勘查河北承德地区基岩裂隙水[J]. 地球科学与环境学报, 2020,42(6): 784−790. DOI: 10.19814/j.jese.2020.06013. WANG R F, WEN L F, CHENG J L, et al. Joint detection of bedrock fissure water using high-density electrical method and transient electromagnetic method in Chengde area of Hebei, China[J]. Journal of Earth Science and Environment, 2020, 42(6): 784−790. DOI: 10.19814/j.jese.2020.06013. (in Chinese).
[2] 张保祥, 刘春华. 瞬变电磁法在地下水勘查中的应用综述[J]. 地球物理学进展, 2004,19(3): 537−542. doi: <10.3969/j.issn.1004-2903.2004.03.008 ZHANG B X, LIU C H. Summarization on applications to groundwater exploration by using transient electromagnetic methods[J]. Progress in Geophysics, 2004, 19(3): 537−542. (in Chinese). doi: <10.3969/j.issn.1004-2903.2004.03.008
[3] 李国占, 孙银行. 地下水地球物理勘查技术模式[J]. 物探与化探, 2010,34(2): 202−204. LI G Z, SUN Y H. A tentative discussion on a technological mode for groundwater geophysical exploration[J]. Geophysical and Geochemical Exploration, 2010, 34(2): 202−204. (in Chinese).
[4] 王星明, 郭栋, 李嘉. 水资源勘查中综合电法勘探方法技术与应用[J]. 物探与化探, 2011,35(1): 65−69. WANG M X, GUO D, LI J. The integrated electrical prospecting technology and its application to water resource exploration[J]. Geophysical and Geochemical Exploration, 2011, 35(1): 65−69. (in Chinese).
[5] 宋希利, 宋鹏, 田明阳, 等. 物探方法在侵入岩地区抗旱找水定井中的应用[J]. 地球物理学进展, 2012,27(3): 1280−1286. DOI: 10.6038/j.issn.1004-2903.2012.03.057. SONG X L, SONG P, TIAN M Y, et al. Geophysical prospecting method in intrusive rocks area fight a drought to find water wells set[J]. Progress in Geophysics, 2012, 27(3): 1280−1286. DOI: 10.6038/j.issn.1004-2903.2012.03.057. (in Chinese).
[6] 苏永军, 马震, 孟利山, 等. 高密度电阻率法和激发极化法在抗旱找水定井位中的应用[J]. 现代地质, 2015,29(2): 265−271. doi: 10.3969/j.issn.1000-8527.2015.02.007 SU Y J, MA Z, MENG L S, et al. Application of high-density resistivity method and induced polarization method to determine a good well location in groundwater prospecting[J]. Geoscience, 2015, 29(2): 265−271. (in Chinese). doi: 10.3969/j.issn.1000-8527.2015.02.007
[7] 孙中任, 杨殿臣, 赵雪娟. 综合物探方法寻找深部地下水[J]. 物探与化探, 2017,41(1): 52−57. DOI: 10.11720/wtyht.2017.1.08. SUN Z R, YANG D C, ZHAO X J. The application of integrated geophysical methods to the prospecting for deep geothermal resource[J]. Geophysical and Geochemical Exploration, 2017, 41(1): 52−57. DOI: 10.11720/wtyht.2017.1.08. (in Chinese).
[8] 康方平, 蒋建良, 彭杰, 等. 综合物探方法在湖南某贫水板岩地区找水的应用研究[J]. 工程地球物理学报, 2020,17(2): 258−264. DOI: 10.3969/j.issn.1672-7940.2020.02.018. KANG F P, JIANG J L, PENG J, et al. Application of integrated geophysical method to water search in a poor slate region of Hunan province[J]. Chinese Journal of Engineering Geophysics, 2020, 17(2): 258−264. DOI: 10.3969/j.issn.1672-7940.2020.02.018. (in Chinese).
[9] 潘剑伟, 张成丽, 鲁恺, 等. SNMR联合SP用于滑动带水体赋存状态的探测研究[J]. CT理论与应用研究, 2021,30(1): 23−34. DOI: 10.15953/j.1004-4140.2021.30.01.03. PAN J W, ZHANG C L, LU K, et al. The combined use of surface nuclear magnetic resonance and self-potential in the research of groundwater state in the landslide[J]. CT Theory and Applications, 2021, 30(1): 23−34. DOI: 10.15953/j.1004-4140.2021.30.01.03. (in Chinese).
[10] 黄国民, 李世平, 陶毅, 等. 广西碎屑岩地区电法找水实例[J]. 物探与化探, 2019,43(1): 77−83. DOI: 10.11720/wtyht.2019.2134. HUANG G M, LI S P, TAO Y, et al. A case study of water prospecting by electrical method in clastic rock area of Guangxi[J]. Geophysical and Geochemical Exploration, 2019, 43(1): 77−83. DOI: 10.11720/wtyht.2019.2134. (in Chinese).
[11] 马吉静. 高密度电阻率法的异常识别和推断-以溶洞探测和寻找地下水为例[J]. 地球物理学进展, 2019,34(4): 1489−1498. DOI: 10.6038/pg2019CC0054. MA J J. Anomaly identification and inference of high density resistivity method: Take Karst cave exploration and groundwater exploration as an example[J]. Progress in Geophysics, 2019, 34(4): 1489−1498. DOI: 10.6038/pg2019CC0054. (in Chinese).
[12] 底青云, 倪大来, 王若. 高密度电阻率成像[J]. 地球物理学进展, 2003,18(2): 323−326. doi: 10.3969/j.issn.1004-2903.2003.02.023 DI Q Y, NI D L, WANG R. High-density resistivity image[J]. Progress in Geophysics, 2003, 18(2): 323−326. (in Chinese). doi: 10.3969/j.issn.1004-2903.2003.02.023
[13] 董浩斌, 王传雷. 高密度电法的发展与应用[J]. 地学前缘, 2003,10(1): 171−176. doi: 10.3321/j.issn:1005-2321.2003.01.020 DONG H B, WANG C L. Development and application of 2D resistivity imaging surveys[J]. Earth Science Frontiers, 2003, 10(1): 171−176. (in Chinese). doi: 10.3321/j.issn:1005-2321.2003.01.020
[14] 严加永, 孟贵祥, 吕庆田, 等. 高密度电法的进展与展望[J]. 物探与化探, 2012,36(4): 576−584. YAN J Y, MENG G X, LV Q T, et al. The progress and prospect of the electrical resistivity imaging survey[J]. Geophysical and Geochemical Exploration, 2012, 36(4): 576−584. (in Chinese).
[15] 陈松, 余绍文, 刘怀庆, 等. 高密度电法在水文地质调查中的应用研究−以江平圩幅为例[J]. 地球物理学进展, 2017,32(2): 0849−0855. DOI: 10.6038/pg20170254. CHEN S, YU S W, LIU H Q, et al. Application and research high density electrical method in hydrogeological prospceting: A case study on Jiangping town map[J]. Progress in Geophysics, 2017, 32(2): 0849−0855. DOI: 10.6038/pg20170254. (in Chinese).
[16] 刘艳秋, 徐洪苗, 胡俊杰. 综合物探方法在水库堤坝隐患探测中的应用[J]. 工程地球物理学报, 2019,16(4): 546−551. DOI: 10.3969/j.issn.1672-7940.2019.04.018. LIU Y Q, XU H M, HU J J. Application of comprehensive geophysical exploration technique to detecting hidden defects of reservoir dams[J]. Chinese Journal of Engineering Geophysics, 2019, 16(4): 546−551. DOI: 10.3969/j.issn.1672-7940.2019.04.018. (in Chinese).
[17] 薛国强, 李貅, 底青云. 瞬变电磁法理论与应用研究进展[J]. 地球物理学进展, 2007,22(4): 1195−1200. doi: 10.3969/j.issn.1004-2903.2007.04.026 XUE G Q, LI X, DI Q Y. The progress of TEM in theory and application[J]. Progress in Geophysics, 2007, 22(4): 1195−1200. (in Chinese). doi: 10.3969/j.issn.1004-2903.2007.04.026
[18] 严良俊, 徐世浙, 胡文宝, 等. 中心回线瞬变电磁测深法快速电阻率成像方法及应用[J]. 煤田地质与勘探, 2002,30(6): 58−60. doi: 10.3969/j.issn.1001-1986.2002.06.020 YAN L J, XU S Z, HU W B, et al. A rapid resistivity imaging method for central loop transient electromagnetic sounding and its application[J]. Coal Geology & Exploration, 2002, 30(6): 58−60. (in Chinese). doi: 10.3969/j.issn.1001-1986.2002.06.020
[19] 杨文钦, 胡东祥. 断层富水性的综合探测技术与应用[J]. 煤田地质与勘探, 2002,30(2): 51−53. doi: 10.3969/j.issn.1001-1986.2002.02.019 YANG W Q, HU D X. The technology and it’s application survey for water-bearing condition of fault[J]. Coal Geology & Exploration, 2002, 30(2): 51−53. (in Chinese). doi: 10.3969/j.issn.1001-1986.2002.02.019
[20] 王庆, 石磊, 兰云飞, 等. 微小盆地倾斜煤层采空区瞬变电磁法探测[J]. CT理论与应用研究, 2020,29(5): 576−583. DOI: 10.15953/j.1004-4140.2020.29.05.08. WANG Q, SHI L, LAN Y F, et al. Detection of mine tunnels in inclined coal seam of small coal bearing basin using TEM[J]. CT Theory and Applications, 2020, 29(5): 576−583. DOI: 10.15953/j.1004-4140.2020.29.05.08. (in Chinese).
[21] XUE G Q, ZHANG L B, ZHOU N N, et al. Development measurements of TEM sounding in China[J]. Geological Journal, 2019: 1−8. DOI: 10.1002/gj.3544.
[22] 彭赟, 李燕, 杨淮. 瞬变电磁法在积水采空区探测中的应用−以贵州某煤矿为例[J]. 工程地球物理学报, 2019,16(6): 849−855. doi: 10.3969/j.issn.1672-7940.2019.06.010 PENG Y, LI Y, YANG H. Detection of water-accumulating goaf based on transient electromagnetic method: A case study of a coal mine in guizhou[J]. Chinese Journal of Engineering Geophysics, 2019, 16(6): 849−855. (in Chinese). doi: 10.3969/j.issn.1672-7940.2019.06.010
[23] 日丹诺夫 M S. 地球物理反演理论与应用[M]. 底青云, 薛国强, 李貅, 等, 译. 北京: 科学出版社, 2018. -
期刊类型引用(2)
1. 李梦雨，段诗苗，张雷，周咏春. 提高肺SBRT精度的多窗位动态组合诊疗手段. 中国CT和MRI杂志. 2024(04): 63-65 . 
百度学术
2. 李玉辉，刘龙进，徐乐意，刘远高. 胸部CT不同图像算法对人工智能辅助诊断软件肺结节检出效果的影响研究. 中国医学装备. 2023(12): 10-14 . 百度学术
其他类型引用(3)
下载:
计量
- 文章访问数: 425
- HTML全文浏览量: 139
- PDF下载量: 61
- 被引次数: 5
