Urban Buildings MIMO Downward-Looking Array SAR 3D Simulation Under Non-ideal Trajectory
-
摘要: 在缺少真实下视阵列合成孔径雷达（synthetic aperture radar，SAR）系统而带来数据获取和成像算法研究等困难的当下，下视阵列SAR三维仿真的研究具有重要意义。为了还原载机平台真实的飞行航迹，验证阵列SAR技术在高层建筑物密集城区的地形测绘能力，摒弃传统仿真研究基于匀速直线运动的假设，以高度骤变的城市建筑物为对象，分析构建了多输入多输出（multiple input multiple output，MIMO）下视阵列SAR非理想航迹运动误差模型，提出了非理想航迹城市建筑群MIMO下视阵列SAR三维距离多谱勒成像算法。利用航空平台的航迹与姿态建模仿真技术，以及快速高效的回波仿真技术对成像算法进行了仿真实验，验证了方法的正确性和有效性。
-
关键词:
- 多输入多输出(MIMO) /
- 下视阵列SAR /
- 城市建筑群 /
- 非理想航迹
Abstract: In the era of absence of realistic downward-looking array synthetic aperture radar (SAR) system, the difficulties in data acquisition and further imaging algorithms research have been brought about. So the research of downward-looking array SAR three-dimensional (3D) simulation has great significance. In order to restore the aircraft's real flight path and to verify the topographical mapping ability of array SAR technology in high-rise buildings, the traditional simulation research based on the assumption of uniform linear motion is abandoned. And the non-ideal trajectory motion error model of MIMO(multiple input multiple output) downward-looking array SAR is analyzed and constructed with high-abrupt urban buildings. The 3D range-Doppler (RD) imaging algorithm of urban buildings MIMO downward-looking array SAR under non-ideal trajectory is proposed. Then the flight path and attitude modeling simulation technique of the aeronautical platform, and a fast and efficient echo simulation technique are used to carry out simulation experiments. The correctness and effectiveness of the imaging algorithm are verified by the imaging results.-
Keywords:
- MIMO /
- downward-looking array SAR /
- urban buildings /
- non-ideal trajectory
-
地球重力场是反映物质迁移的基本物理场，直接反映地球内部构造运动、地表质量迁移的本质和过程。地震、火山活动等地球动力学过程都会引起一定程度的地球重力场变化[1-7]。重力场对质量运动或变化信号敏感，适合定量研究各种地球系统过程的时空特性。高精度的地表重力场观测由于距离地壳内部场源近、观测位置固定重复、观测仪器精度高等特点，有利于发现与地壳内部场源直接相关的物理信号。
地震研究不仅依赖相关基础理论的发展，更依赖地震前后时空完整的高精度观测资料。流动重力观测获得的重力时变信号能较好地反映深部物质运移与地壳密度变化等构造活动信息，强震前区域重力场可能会观测到显著的重力异常变化特征，通过对重力时变信号的深入分析，能捕获到与孕震区物性变化有关的物理信息，利用高精度重力测量数据分析震前重力异常变化的时空分布特征，有助于判定未来大震的高风险区域，是利用重力手段开展地震预测的基本出发点。
1 流动重力监测网布局
1.1 重力网布局思路
流动重力监测网是进行地震重力测量的基础，是直接为地震预报服务的。地震监测预报业务能力的提升离不开高精度、高时空分辨率的重力场变化监测资料，而建设布局科学合理的台网是获取此类监测资料的前提。根据历史地震资料，地震学家在中国大陆划出了几十条地震带，涉及地域极广。在这些地震带上都布设密集的地震重力监测网是不切实际的，不仅自然条件和监测能力不允许，实际上也是不必要的。因此，按照地震工作的方针，只在有一定经济建设规模、人口比较稠密、社会影响显著的主要地震带布设相对密集的重力监测网，如南北地震带、天山地震带、华北地震带等。重力监测网的布设始终坚持服务于中强以上地震预测，尤其是强震/大震预测这一主题。
流动重力测量是对已经建立的重力点按观测程序进行联测，定期进行重复观测，形成重力场变化观测网络，并辅以绝对重力观测加以控制，获得高精度重力场的空间和时间变化。
1.2 重力网布局概况
地震重力监测网是全国地震前兆监测系统的一个重要组成部分和主要监测手段，主要用于精确监视区域重力场的动态演化，并用于地震孕育、发展、发生过程的追踪，为地震预报和相关的地球动力学研究服务。
中国大陆的流动重力测量工作始于1966年邢台地震后，1981年后得到长足发展。1981年前主要使用石英弹簧重力仪，测量精度在30~40 µGal，之后引进LaCoste-Romberg高精度重力仪，精度在10 µGal内。1990年前后，中国大陆的主要活动构造带或地震危险地区陆续建立了30多个区域重力监测网，形成了由近4 000个点组成的定期复测重力监测网，区域网各自独立，呈网状或条状，平均范围小于300 km×300 km[8-9]。1998年，中国地震局与原国家测绘局、原总参谋部测绘局和中国科学院等四部委共建了中国地壳运动观测网络工程（下简称网络工程），该工程建成了中国大陆统一的全国地震重力基本网[10-11]。
2008年汶川8.0级地震后，中国地震局高度重视流动重力监测与预报工作，认为各区域地震重力监测网覆盖区域较小，捕捉地震前兆的能力有限，大部分监测网虽然分布在中国大陆的一些主要活动构造带上，但强度不够、监测能力不足，如川滇地区的许多地震都发生在监测网外[8-9]。另外，各区域重力监测网彼此独立、自成体系，缺乏有效的联测和绝对重力基准控制，降低了观测资料的内在精度，不能得到区域重力场变化的绝对量，严重影响了观测资料的有效应用，使其不能在地震预报中发挥更大作用。因此，有必要将区域重力监测网连接成整体，建成统一重力观测基准，使地震重力监测网建设和观测技术得到快速发展[12-13]。在华北强震震情监视与跟踪、中国大陆地球物理场综合观测、中国大陆构造环境监测网络等重大项目资助下，对零散分布的地震重力监测网进行了多次整体结构优化改造，在重力网布局思路上不断创新，为中国大陆成场分布的地表重力场监测预报工作奠定了重要基础。2010年，由相对重力联测网和绝对重力控制网组成的中国大陆整体重力测量网[13]建成，每年定期开展观测。
21世纪以来，中国大陆发生了8次7级以上大震，包括2001年昆仑山口西8.1级、2008年汶川8.0级和于田7.3级、2010年玉树7.1级、2013年芦山7.0级、2014年于田7.3级、2017年九寨沟7.0级和2021年玛多7.4级地震，除了发生在青藏高原腹地的2001年昆仑山口西8.1级震中附近没有重力测点外，其他7次大震的震中附近均有一定的重力测点，尤其是2008年汶川8.0级、2013年芦山7.0级和2017年九寨沟7.0级地震，震中四周均有重力监测点，这为可靠地提取强震前的重力变化前兆信息打下了良好基础[9-11]。
2 重力数据处理
利用重力时变资料进行地震分析预测研究需要进行重力测量数据处理。只有通过必要的数据处理，从观测资料中提取可靠的重力变化信息，才能实现观测资料的价值。数据处理的目的是获得真实的重力测量值，描述重力场演变过程，从观测资料中获取反映地震孕育发生过程中的重力变化信息。重力观测资料的质量既依赖监测成果质量，也依赖处理资料方法，不同的资料处理方法可能会得到差异较大的结果，这在资料分析中尤为重要。
流动重力测量一般成网状布设，它包括基准点上的绝对重力测量以及联测点间的相对重力联测。相对重力联测具有机动灵活、省时省力、经济成本低等优点，能快速获取地震重点监视区的区域重力场变化。但相对重力联测获得的观测数据是相邻测点间的重力段差，需进行数据处理，又称平差。依据最小二乘原理进行平差，可以合理消除网中各种不符值，削弱偶然误差，消除系统误差，评定测量精度，同时获得流动重力测网中各点的重力值等。在平差之前需要进行预处理，包括对观测值进行固体潮、气压、一次项、仪器高改正等。
进行重力网平差，首先需要选择适当的平差基准，重力网平差基准的选择将直接影响重力场变化结果分析的可靠性[14]。重力网平差基准的选择可以有不同方式，包括经典平差、自由网平差和拟稳平差[15-16]。（1）当网内有2个以上绝对重力点时，一般采用经典平差，选取绝对重力点值作为平差的起算基准；（2）当网中无绝对重力控制、且多期测网测点或网型存在一定变化、但测网内存在部分稳定的重力测点时，通常采用拟稳平差；（3）当网中无绝对重力点、前后两期测网的测点和网型无大变化时，通常采用自由网平差。
1997年以前，中国地震局各施测单位独立的重力测网基本没有绝对重力点，地震重力的数据处理与分析工作由各单位对其测区进行，重力网数据处理无法满足经典平差方法的需求，因此，区域重力网的数据处理通常采用自由网平差与拟稳平差。1998—2009年，网络工程的重力网采用绝对重力控制的经典平差，其他区域网由于缺少绝对重力观测点仍采用自由网平差与拟稳平差。2010年以来，中国大陆流动重力监测网是由相对重力联测网和绝对重力控制网共同组成的，将中国大陆构造环境监测网络重力网与各施测单位的区域重力测网连接成整体，形成中国大陆流动重力统一观测网。因此，重力数据的处理从各省区地震重力网的独立计算，发展到相邻省区联合处理的整体平差计算，既扩大了研究区范围，又能有效提取各省区交界地带的重力场变化信息[17]。
综合利用地面绝对重力、相对重力资料的整体处理分析方法，获得测区内以绝对重力为基准的区域重力场变化。数据处理的重点是将绝对重力观测资料与同期的相对重力观测资料结合，其中绝对重力点构成高精度控制网，相对重力观测视为与该网的定期联测，形成具有绝对基准的区域动态监测网。具体数据处理步骤为：（1）使用中国地震局LGADJ软件对重力观测数据进行统一平差计算[18]；（2）采用稳健估计法对少数存在误差较大的观测段差实行降权处理和粗差剔除，利用计算程序自动优化、合理确定各台仪器的先验方差后进行整体平差计算[19-20]；（3）重力网平差采用弱基准，若基准点有误差、甚至异常，易从平差结果中发现，从而降低误差影响，且可通过平差使基准点重力值得到改善[19-20]；（4）实测标定仪器的一次项[21]，在野外作业中使用基线场的标定值，但随着测区的不同和观测时间的变化，重力仪的格值系数会发生变化，尤其是CG-5/6型重力仪格值系数会产生较大变化，应利用测网中的绝对重力点控制解算出区域适定解，精确标定仪器的格值系数，获取可靠重力变化[22]；（5）采用贝叶斯方法提供的非线性零漂计算模型对仪器非线性漂移的漂移率变化进行估算，消除仪器零漂的影响[23]；（6）用最小二乘配置等方法对重力观测数据进行拟合推估和滤波[9，11]，以便突出显示构造因素的重力效应。
3 重力变化与构造运动及地震活动
根据流动重力监测情况和重力数据处理结果，可以获得不同时空尺度的区域重力场，作为分析研究重力场动态变化图像的依据，利用最小二乘配置、小波分析、双三次样条分析等方法，提取可靠重力场变化信息，系统分析研究多时空尺度的区域重力场动态变化图像演化特征与规律。
3.1 重力场时空演化特征
由于地壳介质的非均匀性和构造运动过程的复杂性，区域重力场的时空动态变化具有明显的不均匀性[11-13]，主要表现为：
1）区域重力监测网范围不同，在不同时空尺度上的重力场动态变化图像存在明显差异。中国大陆整体重力网的观测资料能反映强震活动的主体区，高密度的地震重点监视区重力网的观测资料能反映出发震构造活动引起的局部重力异常。
2）不同时段的重力场在空间分布上表现出明显的差异性。在无显著变化时段，重力场空间分布比较零散，重力变化没有明显的规律性；在显著变化时段，区域重力场空间分布则表现出大范围的有序性变化和局部异常的相对集中性。
3）重力场变化与所处的地质构造环境等因素有关，沿块体构造活动带呈分区变化现象，一般表现为山区重力负值变化、盆地重力正值变化的继承性构造运动。
3.2 重力变化与构造活动的关系
地表重力场动态变化主要由地表观测点的位置变化、变形运动以及地球内部因构造块体变形运动的密度变化效应叠加引起，包含了丰富的构造运动信息[24-28]。
1）重力场变化能够较好地反映测区内主要构造断裂的构造运动过程。由于构造形变沿活动断裂带分布，当块体构造单元运动速率发生改变时，易在其边缘地区出现与地质构造线一致的重力变化梯度带，反映了块体构造单元的相对运动。
2）重力场变化与地壳深部活动断裂构造密切相关，区域重力场变化较显著的梯度带走向与深部构造活跃断裂带走向基本一致，构造单元的边缘往往容易出现重力等值线形态的转折和密集，形成重力变化高梯度带。
3）重力变化能较好地反映与地震孕育发生有关的断层构造活动，在震中区附近出现不同形态、不同规模的局部重力异常区及其伴生的重力变化梯度带。
总之，区域重力场动态变化与构造活动存在密切联系，重力资料能较好地反映在区域构造应力场和震源应力场作用下断层活动引起的重力效应。
3.3 地震孕育、发生过程的重力场变化特征
1966年邢台地震后，中国开始重力场时变的观测研究。20世纪70年代观测到1975年海城7.3级、1976年唐山7.8级等几次大震前重力变化，证实了震前重力异常的存在[1]；20世纪90年代观测到1995年中缅交界7.3级地震前出现约110 µGal重力异常[28]，1996年丽江7.0级地震前震中附近出现约120 µGal重力异常[8]；近年来观测到2008年汶川8.0级和于田7.3级、2010年玉树7.1级、2013年芦山7.0级、2014年于田7.3级、2017年九寨沟7.0级和2021年玛多7.4级等大震前后重力变化[25-32]。
根据中国大陆不同时空尺度的区域重力场演化特征与规律[28-37]，有学者提出以时空区域重力场作为分析重力场动态变化图像的依据，利用不同时空尺度的重力场变化，如相邻两期重力变化、年尺度重力变化、累积重力变化等，选取重力场变化趋势、重力异常形态、重力异常区的范围、持续时间、异常量级等作为判定指标，综合进行地震孕育发生的预测，具体见表 1。
表 1 重力异常与地震孕育发生的关系Table 1. Relation Between Gravity Anomaly and Earthquake Inoculation震级 潜在发震地点的重力场特征 趋势特征 异常范围/km 持续时间 异常量级/µGal 5 重力变化高梯度带 出现局部异常区 ≥100 0.5~1 a ≥50 6 与构造活动有关的重力变化高梯度带转弯部位，四象限分布特征中心附近 出现局部异常区 ≥200 1~3 a或更长 ≥80 7 与块体边界有关的重力变化高梯度带转弯部位，四象限分布特征中心附近 出现正、负局部异常区 ≥400 3~5 a或更长 ≥100 8 与一级块体边界有关的重力变化高梯度带转弯部位，四象限分布特征中心附近 出现正、负局部异常区 ≥600 5 a或更长 ≥120 由表 1可知，要了解地震孕育、发生过程的重力场变化特征，应从重力异常变化特征、空间范围、异常持续时间和异常量级4个因素来考虑。
1）强震发生之前，区域重力场出现大范围的有序性变化，震源区附近产生与地震孕育发生有关的局部重力异常区，并沿区域主要发震构造断裂带出现显著重力变化梯度带。由于孕震构造环境的差异，不同地震的重力变化响应具有明显差异。
2）强震易发生在与构造活动有关联的重力变化正、负异常区过渡的高梯度带上，重力变化等值线的拐弯部位或四象限中心附近。强震前震中区及其附近观测到明显的区域性重力异常及重力变化高梯度带，可能是地震孕育过程中观测到的重力前兆信息。
3）地震震级与重力异常变化的范围、持续时间和变化量级密切相关[17]。一般来讲，观测资料积累的时间越长越有利于判断强震发震震级，重力变化的空间分布特征与发震地点有关；震前重力异常变化覆盖范围越大，则震级越大；震前重力异常变化量级越大，则震级越大；震前2~3 a的观测对发震时间预测非常关键，6级以上地震相邻两期重力变化显著，2 a以上的累计重力变化更为显著，5级地震不存在这种现象。
大量研究成果表明[28-37]，区域构造活动引起的重力时变异常是与地震孕育发生相关的前兆异常之一，特别是对于6级以上地震的孕震过程。
4 重力时变与强震预测
获取可靠的地震前兆信息是实现地震预测的前提。实践证明，通过对重力时变信号的深入分析，能捕获与孕震区物性变化有关的物理信息，是利用重力手段开展地震预测的基本出发点。
在地震预测实践方面，中国地震局于20世纪90年代就对滇西和河西地区利用流动重力观测资料进行了地震预测研究探索，观测到滇西地区1992年永胜5.1级和5.4级、1993年大姚5.3级地震前的重力场变化，1996年丽江7.0级地震前的重力异常变化，并对这几次地震曾进行了一定程度的预测[5，8]。1995年甘肃永登5.8级和2000年甘肃景泰5.9级地震前，河西地区也观测到明显的重力异常变化，并对这两次地震均进行了较准确的中期预测[38-39]。21世纪以来，李辉等[11]利用网络工程重力观测资料获得的1998年以来中国大陆2~3 a尺度重力场动态变化图像，较好地反映了中国大陆地壳构造运动和主要强震活动的基本轮廓。祝意青等[24]利用1998-2008年期间网络工程的重力观测数据，获得中国大陆地区以绝对重力为统一起算基准的重力场动态变化图像，认为中国大陆重力场变化与构造活动断裂和与强震孕育发展有着密切联系，强震易发生在沿构造活动断裂出现的重力变化正、负异常区过渡的高梯度带上以及重力变化等值线的拐弯部位。此外，祝意青等[25]还利用网络工程1998-2005年期间的重力资料对2008年汶川8.0级和于田7.3级等大震做出了准确的中期预测。
预测的主要依据：1998-2005年期间的累积重力变化表明，川滇菱形块体内部出现大范围的重力正值变化异常区，川滇菱形块体北部的四川地区出现较大范围的重力负值变化异常区，两异常区的升降差异变化超过100 µGal，并在马尔康-汶川-成都一带形成重力变化高梯度带（图 1）。
![]()
图 1 汶川Ms 8.0地震前（1998-2005年）重力场变化Figure 1. Gravity Changes Before the Wenchuan Ms 8.0 Earthquake (1998-2005)预测情况：发震时间：2007-2008年；发震地点：（103.7°E，31.6°N）为中心，半径200 km；震级：6~7级。
结果：2008年5月12日，汶川8.0级地震（103.4°E，31.0°N）发生在预测区内，预测震中位置距离汶川8.0级震中72 km，在北川与映秀两个极震区之间（图 2），时间、地点正确，震级提到7级。
![]()
预测依据：1998-2005年较长时期的网络重力变化表明，新疆与西藏交界地区出现了显著重力异常变化，最大变化差异量超过90 µGal，并在新疆于田-和田一带沿构造块体边界带形成显著的重力变化梯度带。
预测情况：发震时间：2007-2008年；发震地点：（80.0°E，36.0°N）为中心，半径200 km；震级：6~7级。
结果：2008年3月21日，于田7.3级地震（81.6°E，35.6°N）发生在预测区内，预测震中位置距离于田7.3级震中149 km，这在地震监测能力比较弱的地区是很准确的（新藏交界地区测点间距100~200 km，测点较稀）。
2008年汶川地震后，中国地震局重力学科管理部加强了重力监测预报研究工作，积极探索利用重力观测资料进行地震年度危险性预测，充分利用日益丰富的重力观测资料，在精细处理区域重力网观测资料的基础上，深入细致地研究区域重力场的动态演化特征与规律，从复杂多变的重力异常现象中寻找构造活动和地震前兆信息，开展震例回溯研究，跟踪分析区域性重力场的异常变化及其可能的中短期异常信息，提出了多种在强震前有预报意义的不同时空尺度的重力场变化特征模式[2，13，35-36]。在强震预测方法和判别指标的基础上，利用重力场动态变化资料对重点监视区的地震危险性进行年度预测和震后趋势判定，在全国地震重点危险区强震地点判定中发挥了重要作用，取得了一定程度的减灾实效[17，40-41]，对近年来发生的2013年芦山7.0级和岷县6.6级、2014年鲁甸6.5级和于田7.3级、2016年门源6.4级、2017年九寨沟7.0级、2020年伽师6.4级、2021年漾濞6.4级和2022年门源6.9级地震等均进行了年度预测[29-35]，尤其是强震发震地点的判定。表 2列出了近年来，利用重力观测资料对中国2013年以来10多次强震预测的情况。
表 2 2013年以来6级以上强震年度预测情况Table 2. Midterm Prediction Information of More than 6 Strong Earthquakes Since 2013预测震中 预测震级 预测时窗 发震地点 实际震级/级 实际震中 发震时间 102.2°E, 30.2°N 6级左右 2013年 四川芦山 7.0 103.0°E, 30.3°N 2013-04-20 105.2°E, 35.5°N 6~7级 2011—2013年 甘肃岷县 6.6 104.2°E, 34.5°N 2013-07-22 81.6°E, 35.6°N 7级左右 2014年 新疆于田 7.3 82.5°E, 36.1°N 2014-02-12 102.1°E, 28.8°N 7级左右 2014年 云南鲁甸
四川康定6.5
6.3103.3°E, 27.1°N
101.7°E, 30.3°N2014-08-03
2014-11-22102.2°E, 37.5°N 6级左右 2016年 青海门源 6.4 101.62°E, 37.68°N 2016-01-21 75.7°E, 39.3°N 6~7级 2016年 新疆阿克陶 6.7 74.04°E, 39.27°N 2016-11-25 85.5°E, 43.5°N 6级左右 2016年 新疆呼图壁 6.2 86.35°E, 43.83°N 2016-12-08 101.9°E, 34.7°N 6~7级 2016—2017年 四川九寨沟 7.0 103.82°E, 33.20°N 2017-08-08 75.8°E, 39.6°N 6~7级 2020年 新疆伽师 6.4 77.21°E, 39.83°N 2020-01-19 100.5°E, 25.8°N 6级左右 2020年 云南漾濞 6.4 99.87°E, 25.67°N 2021-05-21 101.2°E, 37.8°N 6级左右 2021—2022年 青海门源 6.9 101.26°E, 37.77°N 2022-01-08 由表 2可以看出，近年来利用流动重力资料对西部有监测能力地区发生的6级以上强震/大震均进行了较准确的中期预测，但也存在强震漏报的情况，如2017-11-18西藏米林6.9级、2020-07-23西藏尼玛6.6级等地震均发生在青藏高原重力监测能力较弱的地区，不能进行有效监控。
总的来说，流动重力观测资料能较好地反映出地震孕育发生过程中的重力异常变化，尤其是6级以上的强震。通过对流动重力观测资料及时进行缜密分析与研究，有可能对未来强震/大震做出较准确的中期预测，尤其是强震地点的判定。如2022年青海门源6.9级地震地点的预测（见表 2），地震发生在预测的祁连—海原断裂带冷龙岭段，预测的震中位置距离门源6.9级震中不到10 km，基本上一致。
5 问题与展望
获取可靠的地震前兆信息是实现地震预测的前提。重力变化测量是当前大地测量中最精密和困难的工作，地震监测获取的是微伽级变化。观测技术很重要，数据处理也十分关键，由于目前主要以绝对重力控制下的相对重力测量为主，因此，数据处理显得格外重要，应加强数据处理与分析，确保第一手资料的可靠性。此外，不同环境及物理因素导致的重力变化也要考虑，如环境变化、地下水、气压、冰后回弹等。
5.1 问题
在各种自然灾害中，地震是对人类生活影响最严重的自然灾害之一。通过地震监测预报能预测大震何时何地发生，对于造福人类来讲是非常重要的。地震监测的最终目的是为地震预报服务，虽然重力时变监测目前已在中国大陆的地震监测与研究中取得了较为成功的经验，但距离地震预测预防等减灾需求尚有不足，主要存在以下问题：
1）重力变化监测网布局分布不均匀。重力监测点主要分布于南北地震带与大华北地区，而在青藏高原及其周缘地区分布较少，时空分辨率不够，不足以实现对块体边界带的有效监测，所得到的重力变化监测信息是残缺不全的，对6级强震仍不完全具备监测能力，不能有效地捕捉到强震孕育发生过程中出现的完整前兆信息。
2）绝对重力控制程度较低。绝对重力与相对重力没有进行准同步测量，重力监测网缺少足够的绝对重力控制以及对相对重力测量仪器的格值系数精确标定。
3）地面重力与卫星重力缺少有效融合。卫星重力测量主要依托GRACE/GRACE-FO卫星观测资料，虽填补了青藏高原腹地无人区等观测空白，但未能与陆地重力测量形成有效融合，且囿于权限，只能获取国外机构提供的预处理结果，不能获取第一手原始观测资料。
4）新技术应用不足。重力观测仪器设备的核心传感器大都是20世纪研发的弹簧传感器，核心技术已逐渐落后，不能满足获取微小重力变化的需求。近年来，超导重力仪观测技术已发展成熟，但在现有地震重力监测台网中的装备很少，地震监测应用不足；量子重力仪发展快速，相关产品已具备较高的可用性，但目前尚未在地震重力监测中得到应用。
5）时变重力数据分析和研究中，构造和水的影响是造成地面重力变化最大的因素，一些地下水、地表变形、局部近地表水储量变化等干扰因素难以被精确剔除，如何扣除掉地下水等对重力变化的影响和贡献，有效提取反映构造和动力学变化的信号是亟待加强的工作。
5.2 展望
地震预报一直是世界公认的科学难题，地震机理极其复杂，但并不是完全不可知的。地震的发生本是地下构造运动由慢变快的变形过程，问题在于是否能观测、记录到，即使观测到了，又能否识别[42]。从科学技术上讲，实现有效的地震预测预报，首先需要对地震机理、地震孕育发生的物理过程有全面完整的认识；然后需要能够有效地监视地震孕育、发生过程，以便根据孕震过程的发震阶段进行地震预测。这两点一方面是科学理论，一方面是技术能力。近年来，中国地震局注重地震重力监测预报队伍建设与人才培养，一方面加强重力变化机理研究，另一方面对中国大陆地震重力观测网进行整体优化完善，加强地震监测能力建设，流动重力在地震预测预报中发挥着越来越重要的作用，虽取得了多次中期预测成功案例，但目前仍存在地震孕育发生过程中的重力变化机理认识不清、重力观测技术需进一步提高等问题。从长远发展来看，应加强以下工作：
1）发展高精度绝对重力观测技术代替流动相对重力观测技术；发展低漂移（超导）或无漂移式定点连续观测技术代替弹簧式定点连续观测技术；发展低空卫星重力观测技术弥补地面重力观测技术的不足[35]。
2）加强协作，发挥多部门优势，联合大学、科学院所开展中国大陆重力场时变及强震孕育发生机理研究，深化强震发生机理的认识。
3）在中国地震科学实验场的高风险区，建立高密度的绝对重力观测网络或具有良好绝对重力控制的相对重力观测网络，开展密集的强化监测，这有可能揭示高风险区震中附近的地下结构分布特征，提取强震/大震孕育过程中震源变化伴随的绝对重力变化信号，获取高精度绝对重力动态变化图像，为地震孕育环境和孕育过程监测提供前所未有的强大支持。
未来，随着观测技术的进步、重力场时空监测分辨力的增强，强震孕育发生过程中的重力变化机理研究的深入，时变重力观测将进一步提高其在地震分析预报中的能力，服务于中国防震减灾事业的发展需要。
-
![]()
图 1 MIMO下视阵列SAR运动误差模型
Figure 1. Motion Compensation Model of MIMO Downward-Looking Array SAR
![]()
图 2 非理想航迹MIMO下视阵列SAR成像方法
Figure 2. Imaging Algorithm of MIMO Downward-Looking Array SAR Under Non-ideal Trajectory
![]()
图 15 飞行航迹(跨航向加入正弦曲线)
Figure 15. Flight Path (Join Sine Curve in the Cross-Track Direction)
![]()
图 22 飞行航迹(跨航向、高程向加入正弦曲线)
Figure 22. Flight Path (Join Sine Curve in Both Cross-Track Direction and Elevation Direction)
表 1 仿真参数
Table 1 mulation Parameters
发射阵
元数目接收阵
元数目脉冲重复
频率/Hz载频
/GHz采样频率
/MHz阵元间距
/m平台速度
/(m·s-1)脉冲宽度
/μs平台高度
/m带宽
/MHz30 40 200 37.5 180 0.004 100 0.1 600 150 
下载: 导出CSV
表 2 三维成像结果定量指标统计
Table 2 Quantitative Statistics Indicators of 3D Imaging Results
数据 消除阴影区域影响
后高程误差标准差/m消除阴影区域影响
后高程误差均值/m整个场景误差在半个
分辨率之内的概率/%单独建筑物区域误差在
半个分辨率之内的概率/%第1组数据 3.216 1 0.358 0 97.41 92.53 第2组数据 3.154 3 0.122 5 97.50 92.78 第3组数据 3.612 0 0.290 9 96.70 90.51
下载: 导出CSV
-
[1] 靳国旺, 徐青, 张红敏.合成孔径雷达干涉测量[M].北京:国防工业出版社, 2014 Jin Guowang, Xu Qing, Zhang Hongmin. Synthetic Aperture Radar Interferometry[M]. Beijing:National Defense Industry Press, 2014
[2] 靳国旺. InSAR获取高精度DEM关键处理技术研究[D].郑州: 信息工程大学, 2007 Jin Guowang. Research on Key Processing Technique for Deriving Accurate DEM from InSAR[D]. Zhengzhou: Information Engineering University, 2007
[3] 刘辉, 李葛爽.无人机MIMO下视阵列SAR天线构型设计[J].测绘通报, 2016(12):39-43 Liu Hui, Li Geshuang. Antenna Configuration Design of MIMO Downward-Looking Array SAR of Unmanned Aerial Vehicle[J]. Bulletin of Surve-ying and Mapping, 2016(12):39-43
[4] Liu Hui, Xu Qing, Jin Guowang, et al. MIMO Antenna Polynomial Weighted Average Design Method of Downward-Looking Array SAR[J]. Internatio-nal Journal of Antennas and Propagation, 2017(3):1-18
[5] 杜磊.阵列天线下视合成孔径雷达三维成像模型、方法与实验研究[D].北京: 中国科学院电子学研究所, 2010 Du Lei. Study on Model, Algorithm and Experiment for Downward-Looking Synthetic Aperture Radar Three-Dimensional Imaging Based on Linear Array Antennas[D]. Beijing: Institute of Electro-nics, Chinese Academy of Sciences, 2010
[6] 王银波.新型阵列三维SAR关键技术研究[D].成都: 电子科技大学, 2009 Wang Yinbo. Research on Key Technologies of New Three-Dimensional Array SAR[D]. Chengdu: University of Electronic Science and Technology of China, 2009
[7] Gierull C H. On a Concept for an Airborne Downward-Looking Imaging Radar[J]. International Journal of Electronics and Communications, 1999, 53(6):295-304
[8] Giret R, Jeuland H, Enert P. A Study of a 3D-SAR Concept for a Millimeter-Wave Imaging Radar Onboard an UAV[C]. European Radar Conference 2004, Amsterdam, Netherland, 2004
[9] Nouvel J, Jeuland H, Bonin G, et al. A Ka-band Imaging Radar: DRIVE On board ONERA Motorglider[C]. IEEE International Geoscience and Remote Sensing Symposium, Denver, USA, 2006
[10] Klare J, Brenner A R, Ender J H G. A New Airborne Radar for 3D Imaging: Image Formation Using the ARTINO Principle[C]. The 6th European Conference on Synthetic Aperture Radar, Dresden, Germany, 2006
[11] Klare J. A New Airborne Radar for 3D Imaging-Simulation Study of ARTINO[C]. The 6th Euro-pean Conference on Synthetic Aperture Radar, Dresden, Germany, 2006
[12] 李伟华.基于PCA原理的线阵三维SAR成像原理与算法研究[D].成都: 电子科技大学, 2009 Li Weihua. Linear 3D SAR Imaging Principle and Algorithm Research Based on PCA Principle[D]. Chengdu: University of Electronic Science and Technology of China, 2009
[13] 杜磊, 王彦平, 洪文, 等.基于俯仰角压缩的阵列天线合成孔径雷达三维成像算法研究[J].中国科学院研究生院学报, 2010, 27(6):800-808 Du Lei, Wang Yanping, Hong Wen, et al. Three-Dimensional Imaging Algorithm for Synthetic Aperture Radar with Linear Array Antennas Based on Elevation Angle Compression Principle[J]. Journal of the Graduate School of the Chinese Academy of Sciences, 2010, 27(6):800-808
[14] 彭学明, 王彦平, 谭维贤, 等.基于跨航向稀疏阵列的机载下视MIMO 3D-SAR三维成像算法[J].电子与信息学报, 2012, 34(4):943-949 Peng Xueming, Wang Yanping, Tan Weixian, et al. Airborne Downward-Looking MIMO 3D-SAR Imaging Algorithm Based on Cross-Track Thinned Array[J]. Journal of Eletronics and Information Technology, 2012, 34(4):943-949
[15] 滕秀敏, 李道京.机载交规稀疏阵列天线雷达的下视三维成像处理[J].电子与信息学报, 2012, 34(6):1311-1317 Teng Xiumin, Li Daojing. Downward-Looking 3D Imaging Processing for Airborne Cross-Track Sparse Array Radar[J]. Journal of Electronics and Information Technology, 2012, 34(6):1311-1317
[16] 彭学明, 王彦平, 谭维贤, 等.基于感兴趣区域搜寻的机载下视阵列3D SAR波数域快速成像方法[J].电子与信息学报, 2013, 35(7):1526-1531 Peng Xueming, Wang Yanping, Tan Weixian, et al. Fast Wavenumber Domain Imaging Algorithm for Airborne Downward-Looking Array 3D-SAR Based on Region of Interest Pick[J]. Journal of Eletronics and Information Technology, 2013, 35(7):1526-1531
[17] 李学仕, 孙光才, 徐刚, 等.基于压缩感知的下视三维成像新方法[J].电子与信息学报, 2012, 34(5):1017-1023 Li Xueshi, Sun Guangcai, Xu Gang, et al. A New Method of Downward-Looking 3D-SAR Based on Compressed Sensing[J]. Journal of Electronics and Information Technology, 2012, 34(5):1017-1023
[18] 罗煜川.基于压缩感知的阵列SAR三维成像方法研究[D].成都: 电子科技大学, 2015 Luo Yuchuan. Array SAR 3D Imaging Method Based on Compressed Sensing[D]. Chengdu: University of Electronic Science and Technology of China, 2015
[19] 师君, 张晓玲, 韦顺军, 等.基于变分模型的阵列三维SAR最优DEM重建方法[J].雷达学报, 2015, 4(1):20-28 Shi Jun, Zhang Xiaoling, Wei Shunjun, et al. An Optimal DEM Reconstruction Method for Linear Array Synthetic Aperture Radar Based on Variational Model[J]. Journal of Radars, 2015, 4(1):20-28
[20] 赵逸超, 朱宇涛, 杨猛, 等.基于压缩冗余采样的线阵三维SAR成像方法[J].现代电子技术, 2015, 38(16):76-80 Zhao Yichao, Zhu Yutao, Yang Meng, et al. Linear Array SAR Imaging Method Based on Redundant Compressive Sampling[J]. Modern Electronics Technique, 2015, 38(16):76-80
[21] 刘永坦.雷达成像技术[M].哈尔滨:哈尔滨工业大学出版社, 1999 Liu Yongtan. Radar Imaging Technology[M]. Harbin:Harbin Institute of Technology Press, 1999
[22] 保铮, 邢孟道, 王彤.雷达成像技术[M].北京:电子工业出版社, 2005 Bao Zheng, Xing Mengdao, Wang Tong. Radar Ima-ging Technology[M]. Beijing:Publishing House of Electronics Industry, 2005
-
期刊类型引用(16)
1. 祝意青，赵云峰，隗寿春，张国庆，刘芳. 2022年青海门源M_S6.9和四川泸定M_S6.8地震震中的精准判定及思考. 地震地质. 2025(01): 235-245 . 
百度学术
2. 王伟力，郑兵，郝洪涛，温军军，刘宏伟，马伶俐，廖明辉. 四川地区流动重力观测及其在地震预测中的应用. 大地测量与地球动力学. 2024(01): 89-94 . 百度学术
3. 杨雄，祝意青，赵云峰，隗寿春，操聪，吕宏斌，曾亚光. 2020年伽师M_S6.4地震前后重力变化. 大地测量与地球动力学. 2024(03): 277-282 . 百度学术
4. 艾力夏提·玉山，刘代芹，阿卜杜塔伊尔·亚森，陈丽，李杰，赵磊，李秉烨. 北天山流动重力测网场源分辨能力. 大地测量与地球动力学. 2024(03): 272-276 . 百度学术
5. 李树鹏，胡敏章，祝意青，朱成林，崔华伟，陆汉鹏，张刚，王锋吉. 不同类型地震震级与重力变化异常的统计特征及其机理分析. 地球物理学报. 2024(08): 3032-3059 . 百度学术
6. 李树鹏，胡敏章，祝意青，郝洪涛，殷海涛，贾媛，崔华伟，陆汉鹏，张刚，王锋吉，刘洪良. 2023年山东平原M_S5.5地震前重力变化. 地震地质. 2024(05): 1172-1191 . 百度学术
7. 梁霄，储飞，徐如刚，孙鸿博，肖伟鹏，王俊. 2014年霍山M_S4.3地震前后重磁场变化特征及机理分析. 地震地质. 2024(05): 1151-1171 . 百度学术
8. 孔祥魁，刘代芹，艾力夏提·玉山，李杰，陈丽，李瑞，陈荣柳. 乌鲁木齐及周边地区重力场及视密度变化特征. 地震地质. 2024(05): 1123-1150 . 百度学术
9. 李秉烨，艾力夏提·玉山，刘代芹，陈丽，阿卜杜塔伊尔·亚森，阿尔达克·卡德尔别克，赵磊. 乌什M_S7.1地震重力场变化特征分析研究. 内陆地震. 2024(04): 366-375 . 百度学术
10. 赵云峰，祝意青，隗寿春，郑兵，刘芳，杨雄，毛经伦，孙和平. 基于重力数据的2022年青海门源M_S 6.9及四川泸定M_S 6.8地震预测. 科学通报. 2023(16): 2116-2123 . 百度学术
11. 付广裕. 中国地震重力研究现状及其面临的挑战. 武汉大学学报(信息科学版). 2023(06): 858-869 . 百度学术
12. 郑家军，李铂，朱成林，黄金水，杨立涛，董敏，熊玮，韩博，刘海林，池国民，吴婧. 2022年青州M_L4.1地震前后的重力场变化分析. 地震. 2023(03): 102-119 . 百度学术
13. 杨雄，祝意青，赵云峰，隗寿春. 2021年漾濞M_S6.4地震前后重力场动态特征分析. 地震学报. 2023(05): 863-874 . 百度学术
14. Jin Wei，Huining Huang，Ying Jiang，Ziwei Liu. Abnormal characteristics and effectiveness evaluation of the micro-seismic signal before the Debao M_S4.8 earthquake. Geodesy and Geodynamics. 2023(06): 605-613 . 必应学术
15. 祝杰，李婉秋，陶志刚，刘洋洋，马未宇. 青海门源M_S6.9地震震前OLR异常与InSAR同震形变关联性分析. 中国地震. 2023(04): 857-867 . 百度学术
16. 许才军，王晓航，温扬茂，李伟. 地震大地测量确定凹凸体研究进展与展望. 武汉大学学报(信息科学版). 2022(10): 1701-1712 . 百度学术
其他类型引用(4)
计量
- 文章访问数: 1625
- HTML全文浏览量: 94
- PDF下载量: 206
- 被引次数: 20
