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
-
原发性眼部淋巴瘤是指病理证实为眼部淋巴瘤,同时没有伴随眼部以外其他部位的淋巴瘤,或既往没有淋巴瘤诊断与治疗史。原发性眼部淋巴瘤较少见,约占所有结外淋巴瘤的5%~10%,占所有眼部肿瘤的10%。病理上多为非何杰金式淋巴瘤,且以B细胞淋巴瘤为主[1-4]。由于较低的发病率,关于原发眼部淋巴瘤影像学特征,文献报道相对较少。
Tunlayadechanont等[5]通过11例眼部淋巴瘤的能谱CT研究,认为相较于眼部淋巴组织增生性疾病,淋巴瘤具有更高的碘密度值;Yuan等[6]通过对36例眼部淋巴瘤的MRI分析,发现MRI联合影像组学对眼部淋巴瘤的诊断具有一定的诊断价值。但影像组学的鲁棒性与实用性一直受到临床与影像医生的争议,目前尚未在临床开展应用。传统的形态学影像特征仍是临床诊疗的重要依据,早期诊断、规范治疗,是影响原发眼部淋巴瘤预后的重要因素之一。但因临床与影像医师对该病的认识尚不足,早期容易被误诊,不同医生、不同个体间的临床治疗存在一定差异性。
基于此,本研究尝试通过回顾性分析,总结该病的发病特点、影像与病理特征,以期为临床早期诊断与个体化治疗提供参考依据。
1. 资料与方法
1.1 研究对象
通过医院信息管理系统(Hospital Information System,HIS),回顾性分析2015年1月至2022年5月间在国内两家三甲医院诊断与治疗的原发性眼部淋巴瘤。入组标准:①经手术或活检病理证实为眼部淋巴瘤,并有明确的淋巴瘤病理亚分期;②无其他部位同时发生的淋巴瘤,且既往无淋巴瘤诊断与治疗史;③有完整的诊疗资料;④至少有 1种治疗前增强CT或MRI检查。最终共入组82例,其中男性49例,女性33例,年龄18~83岁,中位年龄53岁。
1.2 影像检查
MRI检查采用GE Signa Excite 1.5 T 或 GE Signa Excite HD 3.0 T扫描仪。相控阵头线圈,扫描序列包括平扫T1加权像(T1 weighted imaging,T1WI)及T2加权像(T2 weighted imaging,T2WI)。扫描矩阵384×256,层厚3~5 mm,间隔0.3~0.5 mm,视野18 cm×18 cm,T1WI TR/TE=340~400 ms/8.4~11.4 ms,T2WI TR/TE=2590~3460 ms/108~118 ms,激励次数1~2。
增强扫描包括横断面、冠状面、矢状面,TR/TE=440~700 ms/11.1~11.3 ms,矩阵512×512,层厚3 mm,激励次数2,团注钆对比剂总量15 mL,流率2~2.5 mL/s,动态扫描共6个时相,总扫描时间1 min 25 s。
CT增强检查采用GE Lightspeed 16层或64层螺旋扫描仪,高压注射器经肘前静脉团注非离子型对比剂(300 mgI/mL),总量100 mL,流率2.5 mL/s。管电压120 kV,管电流100 mAs。注射对比剂后30 s开始扫描,层厚3~5 mm。
1.3 病理检查
将手术切除标本或活检标本经10% 中性福尔马林浸泡24~48 h。然后多点取材、固定、脱水、透明与包埋,经HE染色制成常规病理石蜡切片。免疫组化采用免疫荧光染色法。
1.4 临床与影像特征
1.4.1 临床与病理特征
分析并记录本组原发性眼部淋巴瘤的临床特征,具体包括:眼球突出、眼睑肿胀、疼痛、结膜炎、眼球运动障碍、自觉触及肿物、视力或视野改变、眼睑下垂及眼肌麻痹等。当一位患者同时有多项症状时,需分别记录。
原发性眼部淋巴瘤的病理取材及病理诊断包括:将标本置于10% 中性福尔马林浸泡24~48 h。然后多点取材、固定、脱水、包埋,经HE染色制成常规病理石蜡切片。免疫组化采用免疫荧光染色法。病理诊断包括是否淋巴瘤及具体亚型、免疫组化等信息。
本组原发性眼部淋巴瘤治疗方法包括:手术、放疗或化疗。如为两种或两种以上治疗方法的联合,需分别记录:包括手术+化疗;手术+放疗;放疗+化疗;手术+放化疗。
1.4.2 影像特征
原发性眼部淋巴瘤的影像特征具体包括:单眼或双眼受累;单发或多发病灶(其中多发病灶包括单眼多发或双眼同时发生的病灶);肿瘤发生部位(将眼部划分为眼睑、结膜、泪腺、眼球、肌锥、肌锥内、肌锥外7个区,记录肿瘤主体所在部位,多发病灶时记录最大肿瘤的主体所在部位);生长方式(局限性或浸润性生长);MRI信号特征(以同层肌肉信号强度作为参照,T1WI及T2WI信号特征分为高、等、低信号3类);MRI强化程度(以同层肌肉作为参照,分为明显、中等或轻度强化);MRI强化均匀性(均匀或不均匀);CT强化程度(以同层肌肉作为参照,分为明显、中等或轻度强化);CT强化均匀性(均匀或不均匀);邻近骨质改变。
1.5 统计方法
运用SPSS 22.0软件进行数据分析,计量资料先进行正态分布性检验,对符合正态分布的数据,采用均值±标准差表示;不符合正态分布的数据,采用中位数(四分位数)描述;计数资料采用频数与百分比进行描述。
2. 结果
2.1 临床与病理特征
2.1.1 临床症状
本组原发性眼部淋巴瘤患者临床症状缺乏典型特异性,56例患者表现为多种症状并存。具体包括:眼睑肿胀52例;眼球突出43例;眼睛疼痛34例;结膜炎31例;眼球运动障碍29例,自觉触及肿物25例;视力或视野改变22例;眼睑下垂16例;眼肌麻痹13例(表1)。
表 1 原发性眼部淋巴瘤的临床与影像征象分析Table 1. Analysis of the clinical characteristics and imaging features of primary orbital lymphoma观察指标 例数 观察指标 例数 年龄/岁,中位数(四分位数) 53(32,73) 生长方式(n=82),例数(%) 性别(n=82),例数(%) 局限性 26(31.7) 男 49(59.8) 浸润性 56(68.3) 女 33(40.2) MRI-T1 WI信号(n=66),例数(%) 临床症状(n=82),例数(%) 等 60(90.9) 眼睑肿胀 52(63.4) 低 6(9.1) 眼球突出 43(52.4) MRI-T2 WI信号(n=66),例数(%) 眼睛疼痛 34(41.5) 高 51(77.3) 结膜炎 31(37.8) 等 15(22.7) 眼球运动障碍 29(35.4) 低 0(0) 自觉触及肿物 25(30.5) MRI强化程度(n=66),例数(%) 视力或视野改变 22(26.8) 明显 49(74.2) 眼睑下垂 16(19.5) 中等 13(19.7) 眼肌麻痹 13(15.9) 轻度 4(6.1) 病理亚型(n=82),例数(%) MRI强化方式(n=66),例数(%) 粘膜相关淋巴组织淋巴瘤 60(73.2) 均匀 59(89.4) 弥漫大B细胞淋巴瘤 13(15.9) 不均匀 7(10.6) 套细胞淋巴瘤 3(3.7) CT强化特征(n=61),例数(%) 滤泡性淋巴瘤 2(2.4) 明显强化 43(70.5) 淋巴浆细胞淋巴瘤 1(1.2) 中等强化 13(21.3) T细胞淋巴瘤 3(3.7) 轻度强化 5(8.2) 单眼或双眼受累(n=82),例数(%) CT强化方式(n=61),例数(%) 单眼 59(72.0) 均匀 51(83.6) 双眼 23(28.0) 不均匀 10(16.4) 单灶或多灶(n=82),例数(%) CT骨质破坏(n=61),例数(%) 单灶 47(57.3) 有 2(3.3) 多灶 35(42.7) 无 59(96.7) 肿瘤主体位置(n=82),例数(%) 治疗方法(n=82),例数(%) 结膜区 29(35.4) 手术 9(11.0) 眼睑区 16(19.5) 放疗 14(17.1) 肌锥区 11(13.4) 化疗 7(8.5) 肌锥内区 10(12.2) 手术+放疗 24(29.3) 泪腺区 9(11.0) 手术+化疗 17(20.7) 肌锥外区 5(6.1) 手术+放化疗 8(9.8) 眼球 2(2.4) 放疗+化疗 3(3.7) 2.1.2 病理类型
本组病例中,病理以B细胞淋巴瘤为主,共79例,包括黏膜相关淋巴组织淋巴瘤60例(73.2%),弥漫大B细胞淋巴瘤13例(15.9%),套细胞淋巴瘤3例(3.7%),滤泡性淋巴瘤2例(2.4%),淋巴浆细胞淋巴瘤1例(1.2%)。另有T细胞淋巴瘤3例,包括NK-T细胞淋巴瘤2例(2.4%),T-淋巴母细胞瘤1例(1.2%),具体病例分布见表1。
2.1.3 治疗方法
本组病例中,手术包括全切或局部切除;放疗采用总剂量30~46 Gy,单次剂量1.8~2 Gy,共15~23次;化疗采用CHOP或R-CHOP方案,根据个体具体情况进行4~8周期。其中单纯手术9例(11.0%),单纯放疗14例(17.1%),单纯化疗7例(8.5%),手术+放疗 24例(29.3%),手术+化疗17例(20.7%),手术+放化疗 8例(9.8%),放疗+化疗 3例(3.7%)。
2.2 影像特征
本组病例共有66例(66/82,80.5%)行增强MRI检查,61例(61/82,74.4%)行增强CT检查。其中45例(45/82,54.9%)同时行CT及MRI检查。具体影像特征结果分析见表1。
2.2.1 眼部受累
本组病例中,单眼受累59例(72.0%),双眼受累23例(28.0%);单发病灶47例(57.3%),多发病灶35例(42.7%)(表1)。
2.2.2 肿瘤主体部位
肿瘤主体位于结膜区29例(35.4%),眼睑区16例(19.5%),肌锥区11例(13.4%),肌锥内区10例(12.2%),泪腺区9例(11.0%),肌锥外区5例(6.1%),眼球2例(2.4%)(表1)。
2.2.3 生长方式
本组病例中,肿物呈局限性生长26例(31.7%),浸润性生长56例(68.3%)。
2.2.4 MRI信号特征
本组共有66例患者行MRI检查,其中TWI表现为等信号60例(90.9%),低信号6例(9.1%);T2 WI表现为高信号51例(77.3%),等信号15例(22.7%);明显强化49例(74.2%),中等强化13例(19.7%),轻度强化4例(6.1%)(图1);肿物呈均匀强化59例(89.4%),不均匀强化7例(10.6%)(表1)。
2.2.5 CT特征
本组共有61例行增强CT检查。肿物呈明显强化43例(70.5%),中等强化13例(21.3%),轻度强化5例(8.2%)(图2);均匀强化51例(83.6%),不均匀强化10例(16.4%);2例浸润性生长的病例可见邻近骨质破坏(表1)。
3. 讨论
原发性眼部淋巴瘤以B细胞淋巴瘤为主,文献报道占90% 以上[1-3]。本组病例中,B细胞淋巴瘤共79例,占96.3%,其中主要为粘膜相关淋巴组织淋巴瘤,其次为弥漫大B细胞淋巴瘤。眼部淋巴瘤男性略高于女性,文献报道为1.25∶1[4],本组为1.45∶1。病变多发生于中老年人,本组中位年龄为53岁,与文献报道基本一致[1-3]。
本组病例中,患者临床症状表现多样,主要症状包括眼睑肿胀、眼球突出、疼痛、结膜炎等,同一患者可表现为多种症状。有文献认为单眼突出可以作为淋巴瘤的特异性临床征象,需积极进行MRI或CT等影像检查除外眼部淋巴瘤[4]。淋巴瘤可累及眼眶任何部位,本组最多见于结膜,占35.4%,与文献报道相仿[1,7],其次为眼睑、肌锥周围等部位。
原发性眼部淋巴瘤以单眼受累多见,尤其是粘膜相关淋巴组织淋巴瘤及弥漫大B细胞淋巴瘤,而双眼受累多见于套细胞淋巴瘤[1]。本组研究中,有23例双眼受累,占28.0%,高于文献报道。分析其原因,可能与就诊该两家医院的病例多为较复杂病例有关。3例套细胞淋巴瘤中,有2例为双眼受累,与文献报道一致。Olsen等[1]认为双眼较单眼受累提示预后不良,这可能与肿瘤负荷重、较高的临床分期等因素有关。双眼受累以结膜部位最多见,本组共有14例双侧结膜受累,占双眼受累病例的60.9%。
本组原发性眼部淋巴瘤MRI TWI多表现为等信号,占90.9%,T2 WI多为高信号,占77.3%。增强扫描超过2/3病例呈明显强化,少数呈中等或轻度强化,且强化多较均匀,这与Lecler等[8]报道相仿。不均匀强化的肿瘤多因为生长速度快,血供不足导致肿瘤内部出现坏死,提示该肿瘤具有较高的恶性度[9-10]。本研究中,有7例病灶在增强MRI显示不均匀强化,分别为弥漫大B细胞淋巴瘤与套细胞淋巴瘤,与该两种肿瘤高恶性度相一致。
与增强MRI表现相仿,原发性眼部淋巴瘤在增强CT上仍以均匀明显强化为主,其中均匀强化占82.3%,明显强化占70.5%。与MRI相比,CT可以更直观、清楚观察肿物邻近骨质的改变。本组有2例病例由于肿物呈浸润性生长,邻近骨质受肿物压迫、侵袭而出现吸收破坏,因此提醒临床医生在手术时或肿瘤放疗时要注意手术或放疗野范围及放疗剂量的调整。本研究未进一步分析不同亚型间淋巴瘤的影像表现,Juniat等[11]认为通过影像学难以进一步区分淋巴瘤的病理亚型。
原发性眼部淋巴瘤的发生在近几年有增高的趋势,其发病原因尚不明确。有作者认为肿瘤发生可能与机体长期存在的慢性感染有关,尤其是粘膜相关淋巴组织淋巴瘤。这些慢性感染包括衣原体感染、幽门螺杆菌感染、慢性肝炎等[12-14]。机体免疫状态也被认为与眼部淋巴瘤的发生具有一定的相关性。Nutting等[15]报道桥本氏甲状腺炎引起的甲状腺眼病可能与粘膜相关淋巴组织淋巴瘤的发生有关。李开明等[16]也曾报道IgG4相关性眼病与眼附属器黏膜相关淋巴组织淋巴瘤的发生相关。
发生于眼眶的炎性假瘤由于具有与原发性眼部淋巴瘤的相似的症状,在临床工作中容易混淆,应注意二者鉴别,眼眶炎性假瘤是一种与免疫相关的疾病,以青少年与中年相对多见,伴有眼眶的炎性改变是其与原发性眼部淋巴瘤的重要鉴别点之一[17]。根据受累部位可分为眶前炎症、弥漫性眼眶炎症、眼眶肌炎或泪腺炎。二者有时临床鉴别有一定困难,必要时需组织活检,以免贻误诊断。
对于原发性眼部淋巴瘤,手术多在术前诊断不明确、肿瘤表浅或肿瘤较大需缓解压迫症状时进行,单一的手术治疗不足以改善患者的预后[1,7,13]。对于低度恶性的淋巴瘤如粘膜相关淋巴组织淋巴瘤,中等度的放射治疗剂量可以取得较好的治疗效果[7,16,18]。而对于恶性度较高的肿瘤,如弥漫大B细胞淋巴瘤,多需采用化疗或化疗联合放疗[1,10]。本组最多见的治疗方法为手术+放疗或化疗,占总例数的1/2。关于不同治疗方案或不同亚型淋巴瘤对患者预后的影响,我们将在未来的工作中进一步研究。
总之,原发性眼部淋巴瘤以B细胞淋巴瘤多见,尤其是粘膜相关淋巴组织淋巴瘤。单眼较双眼受累常见,可局限性或浸润性生长,病灶在MRI或CT上以均匀、明显强化为主,是其特征性表现之一。影像特征结合临床病史有助于原发性眼部淋巴瘤的早期诊断,为临床精准、个体化治疗提供重要的参考依据。
-
表 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 表 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 表 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 表 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 表 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 -
[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).