Experiment and Forward Modeling Analysis of Microgravity Detection of Urban Underground Space
-
摘要: 随着国家城市快速发展，城市地下空间探测与开发利用需求加大，传统物探方法在人文活动区域受干扰因素较多，无法获取真实准确的探测数据。微重力方法相对而言受干扰因素较小，对于城市建筑物和人类活动遗迹干扰可以通过模型正演校正的方法消除，从而获取高精度重力数据，进而通过有效的反演方法，可以获得城市地下空间的隧道、采空区、空洞、塌陷区、管廊等空间位置信息。本文用地面移动式高精度重力测量仪器，进行城市探测中受影响的因素进行试验性探测与分析，并结合正演模型校正研究，微重力方法在城市地下空间探测上有不错的效果。Abstract: With the rapid development of national cities, the demand for urban underground space exploration development and utilization has increased. Due to the interference factors in human activities areas, traditional geophysical methods can't obtain true and accurate detection data. The microgravity method is relatively less affected by interference factors. The interference from urban buildings and human activities can be eliminated by the method of model forward correction, so as to obtain high-precision gravity collection data, and then the spatial location information of tunnel, goaf, cavity, collapse area and pipe gallery in urban underground space can be obtained through effective inversion method. In this paper, through experimental detection and analysis of the theoretically affected factors in the urban detection carried out by the ground mobile high-precision gravity measuring instrument, combined with the forward model correction research, the microgravity method shows good effect in the detection of urban underground space.
-
Keywords:
- forward model. /
- urban space /
- microgravity /
- gravity detection
-
急性肺栓塞（pulmonary embolism，PE）具有较高的发病率及病死率，是常见的三大致死性心血管疾病之一。CT肺动脉成像（computed tomography pulmonary arteriography，CTPA）是确诊PE影像学检查的首选方法[1]，具有很高的敏感性和特异性。64排以上CT扫描速度快，CTPA的数据采集时间都可以控制在3 s以内，准确把握肺动脉增强峰值时间是检查成功的关键。CTPA扫描延迟时间的选择目前常用的方法有小剂量团注测试技术和对比剂团注跟踪技术，小剂量团注测试技术操作步骤多，受操作者个人经验影响较大，对比剂团注跟踪技术受限于机器的过渡延迟时间（即ROI到达阈值切换到正式扫描之间的时间），易导致肺静脉过度强化影响肺动脉的观察。
本研究在两种技术的基础上，通过改良对比剂注射方案，即两段式对比剂同速注射的方法，探讨一种获得优质CTPA成像效果的扫描技术。
1. 材料与方法
1.1 资料与分组
本研究通过首都医科大学附属北京同仁医院伦理委员会的批准（TREC2022-KY037）。收集2022年7月至12月在首都医科大学附属北京同仁医院因怀疑肺栓塞行肺动脉CT增强检查的30例患者作为试验组，其中男17例、女13例，年龄31～87岁，平均（64±16）岁，体重53～84 kg，平均（66.0±8.9）kg，身体质量指数（body mass index，BMI）为（24.03±2.61）kg/m2，增强扫描时采用两段式对比剂同速注射结合团注跟踪技术的方法。
收集2021年1月至2021年12月间30例患者为对照组，男16例、女14例，年龄34～85岁，平均（61±14）岁，体重50～85 kg，平均（65.0±10.5）kg，BMI为（23.51±2.56）kg/m2，采用小剂量团注测试技术的方法。
排除标准：肝肾功能不全及严重心功能不全；既往碘对比剂过敏者，妊娠、甲状腺亢进、癫痫及严重哮喘患者；图像中存在金属、呼吸运动等伪影明显影响诊断；体重大于85 kg的患者未纳入本研究。
1.2 检查方法及扫描参数
使用Philips IQon双层探测器光谱CT机和拜耳双筒高压注射器（StellantD-CE超级版），对比剂为碘普罗胺注射液（拜耳医药保健有限公司），碘浓度370 mg/mL。患者仰卧于CT扫描床上，双手上举，放于头部上方，根据患者情况选择左则或右侧肘静脉注射对比剂。
试验组患者的监测兴趣区（region of interest，ROI）置于肺动脉主干，设定阈值为100 HU。对比剂和生理盐水注射顺序：①经 A注射筒注射对比剂10 mL；②经 B筒注射生理盐水30 mL；③经 A筒注射对比剂20 mL；④经 B筒注射生理盐水30 mL。上述注射步骤均按设定顺序由高压注射器自动无间隔完成，对比剂和生理盐水注射流率均为5 mL/s。在开始注射第1段对比剂时同时跟踪肺动脉主干CT值，达到设定阈值后延迟10 s开始扫描（图1）。
![]()
图 1 两段式对比剂注射方法，注射流率均为5 mL/s，两段对比剂开始注射时间间隔8 s，第1段对比剂达到设定阈值后10 s相当于第2段对比剂达到阈值后2 sFigure 1. In the two-stage injection method, the injection flow rate was 5 mL/s, and the interval between the start of the two-stage contrast injection was 8 s; 10 s after the first stage of contrast reaches the set threshold is equivalent to 2 s after the second stage of contrast reaches the threshold对照组患者使用小剂量团注测试技术。监测ROI同样置于肺动脉主干，注射流率5 mL/s，先经A筒注射对比剂10 mL，再经B筒注射生理盐水30 mL，注射开始同时不间断同层连续扫描，每1.5 s扫描1次，共扫描15次，获得肺动脉主干ROI的时间-密度曲线，测量增强峰值时间，本研究中测得的肺动脉主干对比剂达峰时间为7～13 s，平均（10.8±0.8）s。以此达峰时间+1 s作为延迟时间，然后经A筒注射20 mL对比剂，经B筒注射30 mL生理盐水，以此延迟时间进行扫描。
CTPA扫描参数：使用螺旋扫描方式，管电压120 kV，旋转时间0.27 s，使用Philips IQon CT自带管电流自动调制技术（DoseRight），DoseRight Index选择22；参考有效管电流量129 mAs/层，准直宽度为64×0.625 mm，螺距1.015，使用迭代重建（iDose），迭代等级选择3；滤过核选择Standard（B），重建图像层厚1 mm，层间距1 mm，视野300～350 mm，依患者胸部大小调整，矩阵512×512，窗宽窗位使用600/100。扫描范围从肺尖至膈顶，所有患者均在扫描前嘱平静呼吸下呼气后屏气扫描。
1.3 图像后处理方法
将原始薄层图像传至Philips星云工作站，进行多平面重组（multi planar reformation，MPR），获得横断面、冠状面和矢状面图像，同时进行最大密度投影（maximum intensity projection，MIP）。使用容积再现（volume rendering，VR），去骨，结合MIP同时显示肺动脉各段分支。
1.4 图像质量评价
1.4.1 客观评价
两组影像均测量以下血管的CT值：肺动脉主干、右肺动脉、左肺动脉，肺动脉周围分支（右肺上叶动脉、右肺中叶动脉、右肺下叶动脉、左肺上叶动脉、左肺下叶动脉）、肺静脉（右上肺静脉、右下肺静脉、左上肺静脉、左下肺静脉）、注射对比剂侧的锁骨下静脉及升主动脉。如果存在肺栓塞，仅测量未受影响的血管。ROI大小为0.05～0.6 cm2，记录CT值，并计算同侧肺动脉与肺静脉的CT值差，公式为肺动脉CT值-（同侧上肺静脉CT值+下肺静脉CT值）/2。
1.4.2 主观评价
由2名具有10年以上诊断经验的副主任医师以盲法分别对试验组和对照组图像质量进行评价。
根据既往临床诊断经验以4分制标准评价肺动脉图像质量：4分为肺动脉层面CT值≥250 HU，肺动脉及分支血管均显示良好，但是肺动脉与肺静脉CT差值>150 HU，肺静脉略强化，不影响肺动脉分支的鉴别；3分为肺动脉层面CT值≥250 HU，肺动脉及分支血管均显示良好，但是肺动脉与肺静脉CT差值≤150 HU，肺静脉明显强化，影响肺动脉分支的鉴别；2分为150 HU<肺动脉层面CT值<250 HU，肺动脉二、三级血管显影不佳；1分为肺动脉层面 CT值≤150 HU，说明肺动脉强化效果弱，无法明确诊断分支血管。
采用4分制评分法评价上腔静脉硬化伪影：4分，无硬化伪影；3分，硬化伪影可忽略；2分，硬化伪影明显，但肺动脉可显示，尚能诊断；1分，硬化伪影严重影响肺动脉的显示和诊断。
1.5 统计学处理
采用SPSS 23.0软件进行统计分析，采用Shapiro-Wilk检验对数据行正态性分析，符合正态分布的计量资料用
$(\bar{x}\pm s)$ 表示，不符合正态分布用${M}({Q_1}, {Q_3})$ 表示。试验组与对照组患者性别比较采用卡方检验；年龄、BMI、血管CT值的比较采用独立样本t检验；肺动脉图像质量评分、上腔静脉硬化伪影评分的比较采用非参数Mann-Whitney U检验。两名医师间图像质量主观评价的一致性采用Kappa检验，Kappa值<0.40提示两者一致性差，0.41～0.60为一致性一般，0.61～0.80为一致性好，>0.80为一致性非常好。P<0.05为差异有统计学意义。
2. 结果
2.1 一般资料分析结果
试验组和对照组患者性别、年龄、BMI差异均无统计学意义。试验组和对照组通过CTPA检查诊断为肺栓塞的分别为3例和2例。
试验组30例患者均获得良好的肺动脉强化效果，3例患者出现肺静脉强化，干扰肺动脉的观察，而对照组有8例患者肺静脉强化明显。试验组肺静脉强化原因是跟踪ROI时因为呼吸运动ROI有的期相不在跟踪层面，导致ROI达到阈值时间延迟。对照组肺静脉强化原因主要是测试达峰时间误差或启动扫描延迟所致（图2）。
![]()
图 2 对照组患者使用小剂量团注测试技术，试验组患者使用两段式对比剂注射结合团注跟踪技术（a）和（b）为对照组组患者，（a）横轴位厚层MIP图像，（b）冠状位厚层MIP图像，因扫描启动过晚导致肺静脉强化，肺动脉与肺静脉CT值差值小于150 HU，图像质量评分为3分；（c）和（d）为试验组患者，（c）横轴位厚层MIP图像，（d）冠状位厚层MIP图像，肺静脉略显影，肺动脉处于强化峰值，与肺静脉有良好的对比，CT值差值大于150 HU，图像质量评分为4分。Figure 2. The control groups, using the test bolus technique, and the experimental groups, using a two-stage injection of contrast agent combined with bolus tracking technique2.2 客观评价
试验组和对照组所测量血管CT值见表1，其中左肺动脉、右肺上叶动脉、右肺中叶动脉、右肺下叶动脉、左肺上叶动脉、升主动脉CT值差异有统计学意义，试验组CT值大于对照组；肺动脉主干、右肺动脉、左肺下叶动脉、右上肺静脉、右下肺静脉、左上肺静脉、左下肺静脉、锁骨下静脉、右侧动静脉CT值差值、左侧动静脉CT值差值差异无统计学意义。
表 1 试验组和对照组测量血管的平均值、最小值和最大值的比较Table 1. Comparison of mean, minimum, and maximum values of measured vessels between the experimental and control groups测量血管 CT值/HU 统计检验 试验组$\bar{x}\pm s$（min～max） 对照组$\bar{x}\pm s$（min～max） t P 肺动脉主干 340±50（264～431） 312±65（183～483） 1.86 0.067 右肺动脉 339±51（251～435） 313±58（180～463） 1.86 0.068 左肺动脉 345±54（252～431） 314±57（188～466） 2.21 0.031 右肺上叶动脉 367±62（281～498） 329±61（185～497） 2.40 0.020 右肺中叶动脉 359±61（269～468） 322±58（202～472） 2.39 0.020 右肺下叶动脉 363±58（268～511） 329±61（206～469） 2.19 0.033 左肺上叶动脉 362±54（258～441） 330±59（197～467） 2.18 0.033 左肺下叶动脉 362±57（259～457） 332±63（195～474） 1.89 0.064 右上肺静脉 137±39（85～225） 128±57（40～263） 0.74 0.462 右下肺静脉 132±30（86～192） 122±53（45～225） 0.86 0.392 左上肺静脉 137±37（89～226） 134±56（42～265） 0.20 0.841 左下肺静脉 135±35（91～228） 124±52（40～244） 0.90 0.375 锁骨下静脉 935±250（469～1728） 840±303（317～1645） 1.32 0.193 升主动脉 120±21（84～165） 71±19（44～120） 9.50 0.000 右侧动静脉CT值差值 204±57（86～289） 188±89（-17～374） 0.87 0.389 左侧动静脉CT值差值 210±61（70～301） 184±86（6～376） 1.31 0.195 注：$\bar{x}\pm s$为平均值±标准差，min为最小值，max为最大值。 2.3 主观评价
两名医师对试验组和对照组肺动脉影像质量、上腔静脉硬化伪影评分结果见表2。两名医师间评分结果一致性好（Kappa值均>0.61，表2）。试验组和对照组肺动脉影像质量评分差异无统计学意义；试验组和对照组上腔静脉硬化伪影评分差异无统计学意义。
表 2 两名医师对试验组和对照组肺动脉影像质量、上腔静脉伪影评分分布和一致性分析结果Table 2. Distribution and consistency of pulmonary artery image quality and superior vena sclerotic artifact scores analyzed by two physicians in the experimental and control groups组别 医师1/例 医师2/例 评分
$({M}({Q}_{1}, t{Q}_{3} ))$一致性分析 4分 3分 2分 1分 4分 3分 2分 1分 Kappa值 P 肺动脉图像质量评分 试验组 22 8 0 0 23 7 0 0 4.0（3.5,4.0） 0.73 ＜0.001 对照组 15 13 2 0 19 11 0 0 3.75（3.0,4.0） 0.62 ＜0.001 上腔静脉硬化伪影评分 试验组 4 22 4 0 2 24 4 0 3.0（3.0,3.0） 0.66 ＜0.001 对照组 2 24 4 0 2 19 9 0 3.0（2.5,3.0） 0.63 ＜0.001 3. 讨论
优质的CTPA图像是肺动脉及分支充分强化，肺静脉略强化，并与肺动脉有良好的对比，以区分肺动脉和肺静脉，同时上腔静脉没有对比剂硬化伪影或硬化伪影极少，不影响肺动脉的观察[2]，这决定了CTPA要在肺动脉达到峰值时，上腔静脉内对比剂峰值已过，同时要在肺静脉充分强化前完成扫描，由于肺动脉循环很快，肺动脉到肺静脉的循环时间仅2～3 s，可利用的增强时间窗非常窄，这就要求对比剂的注射时间短，同时扫描速度要快，而64排CT基本可以将肺动脉扫描时间控制在3 s，因此对比剂用量不需要太多，根据既往研究[2]，按照注射流率5 mL/s，20 mL对比剂完全可以满足肺动脉增强要求，但是较短的对比剂注射时间对肺动脉增强峰值的把握提出了较高的要求，准确确定CTPA扫描延迟时间异常关键。
延迟时间的选择目前常用的方法有小剂量团注测试技术和对比剂团注跟踪技术，小剂量团注测试技术能够有效地反映肺循环信息，有利于把握准确的扫描时机[2-4]，但是操作步骤多，且受操作环境及人为因素影响较大，对比剂启动注射与CT启动扫描的一致性很难做到绝对一致，由此产生的误差影响了预实验测试肺动脉增强峰值时间的准确性，并且有的机器并不具备测试达峰时间的程序。对比剂团注跟踪技术作为一种简单有效地方法，普遍应用于各种CT血管成像[5-6]，相较于小剂量团注测试技术，对比剂团注跟踪技术更具个性化[7]，受患者个体循环差异影响小。
过渡延迟时间是对比剂团注跟踪技术中一个重要的参数，指的是ROI到达阈值切换到正式扫描之间的时间，这段时间内设备需要完成扫描前准备工作（扫描床移动到扫描起始位置）以及向患者发出检查所需要的呼吸指令，此时间也与设备本身有关，一般在4～7 s。多数研究将ROI置于肺动脉主干[8-9]，由于此时间大于肺循环时间，这就导致肺静脉强化，影响肺动脉的观察，同时需要更多的对比剂注射时间以及对比剂用量。
针对此问题，有研究将上腔静脉及头臂干静脉设为ROI[10-11]，虽然肺静脉不成像但是受心功能影响较大，容易造成检查失败。也有学者使用3期团注法对比剂注射方案[12]或多期双流混合注射方案[13]，以肺动脉主干为ROI，第Ⅰ期注射混合液的对比剂浓度低，起到达阈值触发扫描的作用，同时注射混合液所用时间抵消触发后启动扫描的时间，第Ⅱ期对比剂到达时开始扫描，第Ⅰ期混合液进入肺静脉强化程度低，对肺动脉的观察不会造成影响，提高了检查的成功率。此方案虽然第一期混合液与第Ⅱ期对比剂注射流率一致，但是对比剂浓度不一致，两期对比剂的相对注射流率不一致，其达峰时间不一样，混合液达到设定阈值的时间过晚仍会导致检查失败发生。
本研究借鉴CT血管成像中小剂量团注测试技术的一些研究[2-4,7,14-15]，这些研究测试达峰时间的小剂量对比剂与进行正式扫描的对比剂都使用了相同的浓度和相同的注射流率。基于此，本研究跟踪肺动脉主干的CT值，将测试达峰时间的小剂量对比剂和正式扫描的对比剂以相同注射流率（5 mL/s）分两段注射，中间注射相同流率生理盐水以抵消过渡延迟时间，注射方法是先注射2 s的对比剂，然后注射6 s生理盐水，再注射4 s对比剂，然后注射6 s生理盐水，过渡延迟时间设为10 s。由于两段对比剂的开始注射时间间隔是8 s，且注射流率相同，第1段对比剂达到设定阈值后8 s，第2段对比剂也将达到设定阈值，在此基础上再延迟2 s扫描，相当于第2段对比剂达到设定阈值后2 s开始扫描，此时第2段对比剂刚好完成一次肺循环，不会对肺静脉产生明显强化，而第1段对比剂达到设定阈值后10 s，大部分已经循环出肺静脉。此方法间接实现了将过渡延迟时间设定为2 s，解决了多数设备过渡延迟时间过长的问题，节省了对比剂用量，提高了检查的成功率。
本研究中，试验组部分肺动脉血管CT值高于对照组，其他肺动脉血管CT值及锁骨下静脉CT值、肺静脉CT值、动静脉CT值差值均与对照组无差异，且两组图像在肺动脉图像质量和上腔静脉硬化伪影评分方面均无差异。本研究中，试验组升主动脉的CT值高于对照组，这与第1段对比剂流入胸主动脉导致胸主动脉部分强化有关，本研究中升主动脉CT值均＜165 HU，不影响肺动脉的观察。同时本研究在扫描前嘱患者平静呼吸呼气末屏气，以减少因为过度吸气导致下腔静脉血大量回流入右心房，稀释右心腔对比剂产生肺动脉强化不佳的情况[16-18]。
综上，本研究在进行肺动脉CTA检查时使用两段式对比剂同速注射结合团注跟踪技术，可获得稳定的图像质量，其图像质量不低于传统的小剂量团注测试技术，且团注跟踪技术操作简单易行，受操作者个人影响较小。同时本研究通过在两段对比剂中间注射生理盐水抵消过渡延迟时间的方法，将过渡延迟时间设定为10 s，适用于大多数机器，值得在临床上推广。
本研究的不足：①第一段对比剂在 10 s后大部分流入胸主动脉，导致胸主动脉强化，更有部分心功能不佳患者，第一段对比剂在10 s仍残留于肺静脉，影响肺动脉观察，这种情况要增加两段对比剂之间的生理盐水注射时间，同时增加过渡延迟时间，使第1段对比剂有充足的时间循环流出，而过渡延迟时间只需设定为两段对比剂开始注射间隔时间+2 s即可；②本研究第 2段对比剂用量为 20 mL，注射时间为4 s，要求设备扫描快，适用于将扫描时间控制在2～3 s的设备，如果设备扫描速度过慢，要适当增加第2段对比剂的用量；③本研究未将超重人群纳入，对于超重人群，20 mL对比剂是否适用需要进一步研究。
-
![]()
图 1 人文活动对重力观测影响结果曲线图
Figure 1. Curve of the influence of human activities on the results of gravity observation
![]()
图 2 多层建筑对重力观测影响结果曲线图
Figure 2. Curve of influence of multi-storey buildings on gravity observation
![]()
图 3 高层建筑对重力观测影响结果曲线图
Figure 3. Curve of influence of high-rise buildings on gravity observation
![]()
图 4 交通道路对重力观测影响结果曲线图
Figure 4. Curve of influence of traffic roads on gravity observation
![]()
图 10 未经模型数据校正的布格重力异常图
Figure 10. Bouguer gravity anomaly map uncorrected by model data
![]()
图 11 剩余布格重力异常与空洞位置图
Figure 11. Remaining Bouguer gravity anomaly and location map of voids
表 1 人文活动对重力观测影响试验重力值
Table 1 Gravity value in the experiment of the influence of human activities on the gravity observation
距离/m 重力值/（10-5m/s2） LCR-D-88# CG-5-40769# 差值 10 -0.234 -0.233 0.001 20 -0.200 -0.189 0.011 30 -0.183 -0.180 0.003 40 -0.173 -0.175 0.002 50 -0.157 -0.152 0.005 
下载: 导出CSV
表 2 多层建筑对重力观测影响试验重力值
Table 2 Gravity value of multi-storey building's influence test on gravity observation
距离/m 重力值/（10-5m/s2） LCR-D-88# CG-5-40769# 差值 10 -0.034 -0.023 0.011 20 -0.004 -0.006 -0.002 30 0.032 0.035 0.003 40 0.033 0.022 -0.011 50 0.051 0.054 0.003
下载: 导出CSV
表 3 高层建筑对重力观测影响试验重力值
Table 3 Gravity value of high-rise building's influence test on gravity observation
距离/m 重力值/（10-5m/s2） LCR-D-88# CG-5-40769# 差值 10 -0.043 -0.051 -0.008 20 -0.003 -0.023 -0.020 30 0.029 0.013 -0.016 40 0.033 0.028 -0.005 50 0.054 0.052 -0.002
下载: 导出CSV
表 4 交通道路对重力观测影响试验观测数据
Table 4 Observational data of the impact of traffic roads on gravity observation
距离/m 重力值/（10-5m/s2） LCR-D-88# 1 2 3 4 5 6 平均值 30 -0.229 -0.228 -0.228 -0.228 -0.229 -0.230 -0.229 50 -0.155 -0.155 -0.156 -0.156 -0.158 -0.161 -0.157 70 -0.123 -0.124 -0.123 -0.124 -0.124 -0.125 -0.124 90 -0.027 -0.028 -0.028 -0.029 -0.030 -0.031 -0.029 110 0.024 0.023 0.022 0.022 0.023 0.023 0.023 距离/m CG-5-40769# 1 2 3 4 5 6 平均值 30 -0.220 -0.220 -0.217 -0.218 -0.224 -0.223 -0.220 50 -0.138 -0.132 -0.128 -0.131 -0.130 -0.128 -0.131 70 -0.084 -0.083 -0.074 -0.085 -0.078 -0.082 -0.081 90 -0.035 -0.035 -0.039 -0.039 -0.044 -0.045 -0.040 110 0.038 0.035 0.038 0.036 0.038 0.039 0.037
下载: 导出CSV
表 5 交通道路对重力观测影响试验重力值
Table 5 Gravity value in the test of the impact of traffic road on gravity observation
距离/m 重力值/（10-5m/s2） LCR-D-88#均值 CG-5-40769#均值 两台仪器均值 两台仪器均值互差 30 -0.229 -0.220 -0.225 0.008 50 -0.157 -0.131 -0.144 0.026 70 -0.124 -0.081 -0.102 0.043 90 -0.029 -0.040 -0.034 0.011 110 0.023 0.037 0.030 0.014
下载: 导出CSV
-
[1] 路利春, 周明霞, 李小龙, 等. CG-5重力仪外业工作中常见问题及解决方法[J]. 地质装备, 2018,19(4): 28−31. doi: 10.3969/j.issn.1009-282X.2018.04.012 LU L C, ZHOU M X, LI X L, et al. Common problems and solutions in field work of CG-5 gravimeter[J]. Geological Equipment, 2018, 19(4): 28−31. (in Chinese). doi: 10.3969/j.issn.1009-282X.2018.04.012
[2] 焦新华, 吴燕冈. 重力与磁法勘探[M]. 北京: 地质出版社, 2009. [3] 孟玲顺, 杜晓娟. 勘探重力学与地磁学[M]. 北京: 地质出版社, 2008. [4] 罗孝宽, 郭绍雍. 应用地球物理教程—重力磁法[M]. 北京: 地质出版社, 1991. [5] 王谦身, 安玉林, 张赤军, 等. 重力学[M]. 北京: 地震出版社, 2003. [6] 陈善. 重力勘探[M]. 北京: 地质出版社, 1986. [7] 方俊. 重力测量与地球形状学[M]. 北京: 科学出版社, 1975. [8] 何绍基. 重力测量学[M]. 北京: 测绘出版社, 1957. [9] 许厚泽, 王谦身, 陈益惠. 中国重力测量与研究的进展[J]. 地球物理学报, 1994,(S1): 339−352. XU H Z, WANG Q S, CHEN Y H. Progress of gravity survey and research in China[J]. Chinese Journal of Geophysics, 1994, (S1): 339−352. (in Chinese).
[10] 王艺霖, 刘宽厚, 刘晓兰. 大型地面重力仪设备的现状比较及优化配置[J]. 地质装备, 2017,18(2): 22−24. doi: 10.3969/j.issn.1009-282X.2017.02.004 WANG Y L, LIU K H, LIU X L. Status comparison and optimal configuration of large ground gravimeter equipment[J]. Geological Equipment, 2017, 18(2): 22−24. (in Chinese). doi: 10.3969/j.issn.1009-282X.2017.02.004
[11] 陈明. 微重力测量观测精度提高方法探讨[J]. 华南地震, 2020,40(3): 76−83. DOI: 10.13512/j.hndz.2020.03.011. CHEN M. Discussion on improving the observation accuracy of microgravity measurement[J]. South China Earthquake, 2020, 40(3): 76−83. DOI: 10.13512/j.hndz.2020.03.011. (in Chinese).
[12] 李玉君, 任芳祥, 杨立强, 等. 稠油注蒸汽开采蒸汽腔扩展形态4D微重力测量技术[J]. 石油勘探与开发, 2013,40(3): 381−384. doi: 10.11698/PED.2013.03.19 LI Y J, REN F X, YANG L Q, et al. 4D microgravity measurement technology of steam cavity expansion in steam injection recovery of heavy oil[J]. Petroleum Exploration and Development, 2013, 40(3): 381−384. (in Chinese). doi: 10.11698/PED.2013.03.19
[13] 卢鹏羽. 二维微重力测量与目标发现率研究[J]. 吉林大学学报(地球科学版), 2010,40(S1): 1−5. LU P Y. Two-dimensional microgravity measurement and target detection rate research[J]. Journal of Jilin University (Earth Science Edition), 2010, 40(S1): 1−5. (in Chinese).
[14] 贾民育. 微重力测量技术的应用[J]. 地震研究, 2000,23(4): 452−456. doi: 10.3969/j.issn.1000-0666.2000.04.014 JIA M Y. Application of microgravity measurement technology[J]. Seismological Research, 2000, 23(4): 452−456. (in Chinese). doi: 10.3969/j.issn.1000-0666.2000.04.014
[15] 岳建华, 刘树才, 于景村. 煤矿井下微重力测量方法与应用[J]. 中国煤田地质, 1995,(4): 93−96. YUE J H, LIU S C, YU J C. Measurement method and application of microgravity in coal mine[J]. China Coalfield Geology, 1995, (4): 93−96. (in Chinese).
[16] 王谦身, 周文虎, 武传真, 等. 微重力方法在考古工程中的应用──明茂陵地下陵殿探查[J]. 地球物理学进展, 1995,(2): 85−94. WANG Q S, ZHOU W H, WU C Z, et al. Application of microgravity method in archaeological engineering: Exploration of Ming Mausoleum[J]. Progress in Geophysics, 1995, (2): 85−94. (in Chinese).
[17] 张赤军. 微重力测量的应用及其改善[J]. 地球物理学进展, 1988,(4): 1−6. ZHANG C J. Application and improvement of microgravity measurement[J]. Advances in Geophysics, 1988, (4): 1−6. (in Chinese).
[18] 刘平利, 乔天荣, 张鸿祥. 老旧防空洞探测方法与分析[J]. 测绘与空间地理信息, 2021,44(4): 189−191. doi: 10.3969/j.issn.1672-5867.2021.04.051 LIU P L, QIAO T R, ZHANG H X. Detection methods and analysis of old air-raid shelters[J]. Surveying and Mapping and Spatial Geographic Information, 2021, 44(4): 189−191. (in Chinese). doi: 10.3969/j.issn.1672-5867.2021.04.051
[19] 卢进延. 微重力测量及浅层地震在岩溶勘查中的应用[J]. 广东化工, 2020,47(20): 124−126. doi: 10.3969/j.issn.1007-1865.2020.20.056 LU J Y. Application of microgravity measurement and shallow earthquake in Karst exploration[J]. Guangdong Chemical Engineering, 2020, 47(20): 124−126. (in Chinese). doi: 10.3969/j.issn.1007-1865.2020.20.056
[20] 胡强, 伍吉仓, 郑二龙, 等. 利用微重力测量探测城市地下孔洞[J]. 工程勘察, 2015,43(11): 74−78. HU Q, WU J C, ZHENG E L, et al. Detecting urban underground holes by microgravity measurement[J]. Engineering Investigation, 2015, 43(11): 74−78. (in Chinese).
[21] 曹金国, 王来鹏, 翟广卿, 等. CG-5重力仪及应用[M]. 北京: 解放军出版社, 2007. [22] 吴天彪. 国外新型重磁仪器述评[J]. 地质装备, 2002,(3): 3−7. doi: 10.3969/j.issn.1009-282X.2002.03.001 WU T B. Review of foreign new gravity and magnetic instruments[J]. Geological Equipment, 2002, (3): 3−7. (in Chinese). doi: 10.3969/j.issn.1009-282X.2002.03.001
[23] 耿启立. 重力仪器国外代表产品及国内研发最新进展[J]. 地质装备, 2016,17(1): 27−30. doi: 10.3969/j.issn.1009-282X.2016.01.006 GENG Q L. Foreign representative products of gravity instruments and latest research and developmentprogress in China[J]. Geological Equipment, 2016, 17(1): 27−30. (in Chinese). doi: 10.3969/j.issn.1009-282X.2016.01.006
[24] 王延涛, 潘瑞林. 微重力法在采空区勘查中的应用[J]. 物探与化探, 2012,36(S1): 61−64. WANG Y T, PAN R L. Application of microgravity method in goaf exploration[J]. Geophysical and Geochemical Exploration, 2012, 36(S1): 61−64. (in Chinese).
[25] 陈贻祥. 地面微重力方法在地质灾害调查中的应用效果[J]. 中国岩溶, 1995,(2): 176−185. CHEN Y X. Application effect of ground microgravity method in geological disaster investigation[J]. Karst in China, 1995, (2): 176−185. (in Chinese).
[26] 陈贻祥. 地面微重力测量在工程地质勘察中的应用效果[J]. 水文地质工程地质, 1995,(4): 43−46. CHEN Y X. Application effect of ground microgravity survey in engineering geological survey[J]. Hydrogeology and Engineering Geology, 1995, (4): 43−46. (in Chinese).
[27] 高好林. 微重力测量在地裂缝探测中的应用[J]. 物探与化探, 2005,(5): 414−417. GAO H L. Application of microgravity measurement in ground fissure detection[J]. Geophysical and Geochemical Exploration, 2005, (5): 414−417. (in Chinese).
[28] 路利春, 胡登攀, 张冲, 等. 概率成像技术在重力位场中的应用研究[J]. 矿产与地质, 2018,32(5): 888−894. doi: 10.3969/j.issn.1001-5663.2018.05.015 LU L C, HU D P, ZHANG C, et al. Application of probabilistic imaging technology in gravity potential field[J]. Minerals and Geology, 2018, 32(5): 888−894. (in Chinese). doi: 10.3969/j.issn.1001-5663.2018.05.015
[29] 路利春, 文博, 程斌, 等. 成像技术在矿区重力勘探中的应用研究[J]. 地质与资源, 2018,27(1): 99−102,88. doi: 10.3969/j.issn.1671-1947.2018.01.014 LU L C, WEN B, CHENG B, et al. Study on the application of imaging technology in gravity exploration in mining areas[J]. Geology and Resources, 2018, 27(1): 99−102,88. (in Chinese). doi: 10.3969/j.issn.1671-1947.2018.01.014
[30] 路利春, 赵炳坤, 周明霞, 等. 高精度重磁测量在渭河盆地氦气调查中的应用研究[J]. 陕西地质, 2017,35(2): 52−59. doi: 10.3969/j.issn.1001-6996.2017.02.009 LU L C, ZHAO B K, ZHOU M X, et al. Application of high-precision gravity and magnetic survey in helium survey in Weihe Basin[J]. Shaanxi Geology, 2017, 35(2): 52−59. (in Chinese). doi: 10.3969/j.issn.1001-6996.2017.02.009
[31] 王园, 陈丽森, 戴山岭, 等. 重力异常模型的建立和精度分析[J]. 大地测量与地球动力学, 2017,37(2): 160−162. WANG Y, CHEN L S, DAI S L, et al. Establishment of gravity anomaly model and accuracy analysis[J]. Geodesy and Geodynamics, 2017, 37(2): 160−162. (in Chinese).
[32] 冯兰天. 台阶(断层)模型重力异常正演计算的一种新算法[J]. 世界地质, 2014,33(2): 477−483. doi: 10.3969/j.issn.1004-5589.2014.02.026 FENG L T. A new algorithm for forward calculation of gravity anomaly in step (fault) model[J]. World Geology, 2014, 33(2): 477−483. (in Chinese). doi: 10.3969/j.issn.1004-5589.2014.02.026
[33] 王芃, 张忠杰, 张晰, 等. 基于Matlab平台实现二维复杂地质模型多边形网格建模及重力异常正演计算[J]. 桂林理工大学学报, 2014,34(2): 254−259. doi: 10.3969/j.issn.1674-9057.2014.02.006 WANG P, ZHANG Z J, ZHANG X, et al. Realization of polygon mesh modeling and gravity anomaly forward calculation of 2D complex geological model based on Matlab platform[J]. Journal of Guilin University of Technology, 2014, 34(2): 254−259. (in Chinese). doi: 10.3969/j.issn.1674-9057.2014.02.006
[34] GIANI L, MIRANDA T, PIATTELLA O F. Cosmology and Newtonian limit in a model of gravity with nonlocally interacting metrics[J]. Physics of the Dark Universe, 2019, 26: 100357. doi: 10.1016/j.dark.2019.100357
[35] GHOMSI F, SEVERIN N, MANDAL A, et al. Cameroon's crustal configuration from global gravity and topographic models and seismic data[J]. Journal of African Earth Sciences, 2020, 161(Jan.): 103657.1−103657.13.
[36] 张瑞芳, 贾全山, 路利春. 高精度质子磁力仪性能校验与干扰因素分析[J]. 地质装备, 2020,21(5): 23−28. doi: 10.3969/j.issn.1009-282X.2020.05.005 ZHANG R F, JIA Q S, LU L C. Performance check of high-precision proton precession magnetometer and analysis of interference factors[J]. Geological Equipment, 2020, 21(5): 23−28. (in Chinese). doi: 10.3969/j.issn.1009-282X.2020.05.005
[37] 王立发, 雷晓东, 何祎. 平原区深层隐伏岩溶的重力正演模拟[J]. 城市地质, 2018,13(4): 19−23. doi: 10.3969/j.issn.1007-1903.2018.04.004 WANG L F, LEI X D, HE Y. Gravity forward modeling of deep hidden Karst in plain area[J]. Urban Geology, 2018, 13(4): 19−23. (in Chinese). doi: 10.3969/j.issn.1007-1903.2018.04.004
[38] 贾敦新. 基于重力模型的重庆及周边城市群关系研究[J]. 地理空间信息, 2019,17(2): 48−50,10. doi: 10.3969/j.issn.1672-4623.2019.02.013 JIA D X. Study on the relationship between Chongqing and its surrounding urban agglomerations based on gravity model[J]. Geospatial Information, 2019, 17(2): 48−50,10. (in Chinese). doi: 10.3969/j.issn.1672-4623.2019.02.013
[39] 张洪波, 赵珞成, 邓洪涛. 相对重力测量虚拟仿真系统的建立与模型研究[J]. 地理空间信息, 2018,16(10): 51−53,66. doi: 10.3969/j.issn.1672-4623.2018.10.015 ZHANG H B, ZHAO L C, DENG H T. Establishment of virtual simulation system for relative gravity measurement and model research[J]. Geospatial Information, 2018, 16(10): 51−53,66. (in Chinese). doi: 10.3969/j.issn.1672-4623.2018.10.015
[40] 刘璇, 刘艳芳, 张梦珂, 等. 基于改进重力模型的镇域网络化规划研究[J]. 地理空间信息, 2017,15(3): 8−11,133. LIU X, LIU Y F, ZHANG M K, et al. Research on town network planning based on improved gravity model[J]. Geospatial Information, 2017, 15(3): 8−11,133. (in Chinese).
-
期刊类型引用(3)
1. 杨文杰，梁烨，杨柯. 双流法肺动脉技术在肺动脉扫描中的应用评价. 影像诊断与介入放射学. 2024(04): 278-282 . 
百度学术
2. 胡嘉诚，刘云福，王新艳，王倩，马梓轩，张永县，李伟，刘荣，牛延涛. 对比剂分次团注联合能量成像获得CTU和CTA联合成像的应用研究. CT理论与应用研究. 2024(06): 692-700 . 本站查看
3. 肖爱兰. 基于肺动脉CTA对肺栓塞检出流程优化的探讨. 现代医用影像学. 2024(09): 1644-1646+1656 . 百度学术
其他类型引用(0)
计量
- 文章访问数: 651
- HTML全文浏览量: 245
- PDF下载量: 86
- 被引次数: 3
