1995年日本神户地震建筑物倒塌与S波和盆地边缘次生面波干涉的研究

用地震波在不均匀介质传播中的波动方程的数值模拟方法, 分析了两维剖面SH以及P-SV地震波传播过程中波的叠加干涉与1995年日本神户地震中建筑物倒塌率分布之间的关系. 模拟中使用了错格实数Fourier变换的快速、高精度的微分算子. 利用地震波在不均匀地下介质传播的瞬间分布图分析了地震波在空间叠加干涉的过程. 结果表明, 模拟所得加速度波形的峰值很可能是由基岩表面传入沉积层的次生面波和沉积层中体波在地面表层中的干涉所致. 模拟所得地面加速度振幅的分布与建筑物倒塌率的分布极为一致. 远离地震断层的建筑物峰值倒塌率也和加速度的峰值吻合. 干涉过程的分析表示了建筑物倒塌率的分布与地质构造和地震波速度结构密切相关. 地震体波与盆地边缘次生面波的干涉很可能是形成倒塌极大值的原因. 由此可见用数值模拟方法研究地震波动传播过程以及其干涉过程对于地震灾害研究是十分必要的.     

张北-尚义地震同震形变场雷达差分干涉测量

用合成孔径雷达差分干涉测量技术获取1998年1月10日张北-尚义地震(Ms = 6.2级)的同震形变图, 利用地震前后欧洲遥感卫星1和2号合成孔径雷达(ERS-1/2 SAR)的三景数据做重复轨道差分处理得到震前震后的干涉纹图. 由视向形变图可知, 此次地震造成隆起形变, 形变中心位于114°20′E, 40°57′N, 最大视向位移量达25 cm, 形变集中在300 km2范围内, 根据地震同震形变场的空间分布特征分析了震源机理和孕震构造.     

地震地面运动数值模拟中的土层增幅和波的多重干涉效应

阐述了在研究复杂地震灾害问题时不均匀介质中土层增幅率和地震波多重干涉效应对于用数值模拟方法研究地震地面运动的影响. 在数值模拟中, 首先分析了区域地质构造, 使用了由地表微动矩阵观测结果而得来的高精度地震波速度结构. 又运用地表浅层40 m深的钻孔岩(土)芯结构的地质和地震波速度资料, 进一步分析了浅地表土层对于地震波地面运动增幅率分布的影响. 揭示了重灾区地表土层对地面运动增幅率的非线形关系特征. 在地质结构特征分析基础上, 用数值模拟方法计算了日本神户地震灾区地面运动SH波分量加速度波形. 波动方程数值模拟解法中使用了错格实数傅里叶微分算子的拟Fourier谱微分法. 地面运动数值模拟结果表明地表土层的增幅率和地下速度结构特征对地面加速度波形及其最大振幅分布有很大影响. 由数值模拟地面运动结果与该区的1995年日本神户地震建筑物倒塌率关系的分析可知, 模拟所得地面加速度波形最大振幅与建筑物倒塌率的分布极为一致. 远离地震断层的建筑物的两重峰值倒塌率也和加速度的峰值吻合. 利用模拟的瞬间地震波不均匀多层介质中空间传播快拍图像分析了两重地面运动峰值与地质构造的关系. 结果表明两重峰值的地面运动可能是由基岩表面传入沉积层表层的次生面波和沉积层中体波干涉形成了第1次地面运动峰值; 表层下层通道波次生的面波与体波干涉形成了第2次波动干涉, 导致了第2个地面运动峰值. 这表明了波的多重干涉与地下结构有关. 由此可知详细分析土层增幅率和与地质结构相关的波的多重干涉现象对于用地面运动数值模拟方法研究地震灾害及建筑物倒塌率分布特征是十分有益的. 地质构造和结构的调查研究, 地面运动数值模拟分析对于地震波传播和抗震研究都是十分重要的.     

2022年青海门源Mw 6.6地震的发震断层及孕震构造模式

2022年1月8日青海门源盆地北缘发生Mw 6.6地震,震源机制反演表明此次地震属于左旋走滑事件.震后10 d内,近600个余震被检测到,最大余震为M 5.1级.此次地震发生在祁连-海原左旋走滑断裂系统的冷龙岭段,该断裂段全长127 km,由古地震研究确定的特征地震大小在Mw 7.3～7.5.为了更为全面理解此次地震的震源机制以及当地孕震模式,我们分析了地震波形,获取了主震和17个Ms≥3.0余震的震源机制与矩心深度.利用升、降轨道SAR数据获取的像元偏移数据和同震干涉相位(interferometric synthetic aperture radar, InSAR)确定了两条地表破裂带的位置,并利用InSAR数据反演了主震的滑动模型.研究发现,此次地震破裂带对应于冷龙岭断裂西段和托莱山断裂的阶区,发震断层存在3个形变中心,最大滑动量约为4 m,出现在冷龙岭断裂上,形变中心深度为4 km.滑动模型显示释放了累计能量～1.58×10¹⁹Nm,约合矩震级Mw 6.68,与本文利用地震学方法得到的Mw 6.58接近.结合区域活动构造特征、1986和2016年两次门源地...     

关于《张北震区附近形变的成因》的意见

《张北震区附近形变的成因》一文（简称“成因”）针对2000年12月《科学通报》王超等人《张北-尚义地震同震形变场雷达差分干涉测量》（简称“形变”）和2001年11月《科学通报》张红等人《基于SAR差分干涉测量的张北-尚义地震震形参数反演》（简称“反演”）两篇报道，对1998年1月10日张北地震的雷达差分干涉测量（INSAR）形变场的成因提出了不同的看法，认为基于SAR差分干涉测量得到的形变不是震源断层剪切错动在地面的表现，提出了形变区下方的热物质上涌并导致火山活化的假设。     

2022年泸定MS6.8地震同震形变特征及周边强震危险性

2022年四川泸定6.8级地震发生于鲜水河断裂磨西段,为揭示此次地震的同震形变特征,收集了震中200 km范围内的连续全球定位系统(global positioning system, GPS)测站和震中50 km范围内的强震动数据,进行了高精度处理以提取测站的同震水平位移.此外,还收集了覆盖震中区域的Sentinel和ALOS-2降轨数据,通过干涉差分技术处理并获得了卫星视线向同震形变场.结果显示:水平同震形变呈四象限分布,表明泸定地震的震源机制为左旋走滑;距震中16 km的强震动台站记录到了明显的位移波形,其峰值位移达14 cm、永久变形约为12 cm;合成孔径雷达干涉测量(interferometric synthetic aperture radar, InSAR)显示断层近场30 km×30 km区域具有明显的同震变形,断层西侧卫星视线向最大位移达15 cm,挖角乡附近的局部形变超过15 cm;野外调查推测在“二台子-爱国村”间的断层段发育地表破裂,但从InSAR形变场上难以判定,同震是否破裂到地表仍需详细的野外考察予以确定.反演了川滇块体北东边界主干断裂在泸定地震前的闭锁分...     

一种提高积雪深度估算精度的D-InSAR方法

积雪深度(雪深)是反映积雪时空变化规律的重要参数,是全球及区域气候变化与水文循环等研究不可或缺的观测变量.差分合成孔径雷达干涉(differential interferometric synthetic aperture radar, D-In SAR)技术利用降雪前后微波穿透积雪层形成的差分干涉相位与雪深建立几何函数关系,被广泛地应用于区域小尺度雪深估算研究.然而,其估算精度受干涉像对相干性、局地地形和积雪介电常数等因素的影响.本研究基于高分辨率Sentinel-1 SAR数据,通过引入高相干系数区域、Sentinel-2光学影像、Google Earth影像和地表覆盖类型选取地面控制点优化雪深差分干涉相位解缠精度,利用卫星局地入射角和实测积雪密度降低斜距-相位关系模型经验误差,估算得到青藏高原东北部八宝河流域2021年消融期雪深时空分布,同时根据卫星同步积雪实测资料详细探讨了雪深估算误差来源.122个地面雪深实测数据(气象站点+野外测量)验证结果表明,优化后的D-In SAR差分干涉处理能提高雪深估算精度, RMSE为3.9 cm, MAPE为20.03%, R²     

基于SAR差分干涉测量的张北-尚义地震震源参数反演

基于欧洲遥感卫星（ERS）合成孔径雷达差分干涉测量技术（D-InSAR）获取1998年1月10日张北-尚义地震形变图,通过弹性半空间位移模型模拟形变场,由最佳参数解得到张北-尚义地震的震源参数.模拟结果表明,张北-尚义地震的发震构造为走向95°,倾向西南,倾角30°的具右族性质逆断层,断层面长12 km,宽14 km,震源（断层面中心）深度7.5km,震中位置为北纬40°58＇,东经114°21＇.断层滑动角105.95°,滑动方向为北西I3.26°,位移量为0.728m,地震矩M0为2.69×1018N·m.根据震源参数和区域构造背景分析,断裂活动可能与张家口-蓬莱断裂北西端的活动有关.     

结合时间相关特征改进中国区域二氧化碳卫星观测数据的地统计分析

通过卫星从空间观测全球大气二氧化碳(CO₂)浓度为全球碳循环研究提供了新的数据源. 虽然卫星可以进行高密度点观测, 但是由于云和卫星观测模式等影响, 卫星观测点数据在空间上呈不规则分布且存在大量无观测值的空白区域. 精确填补这些空白区域将有利于研究全球和区域大气CO₂浓度特征. 地统计学分析方法利用卫星观测的大气CO₂数据固有的自相关结构进行插值, 可以解决上述问题. 该分析方法包括两个方面, 通过变异函数探讨数据的相关结构和对空白区域进行最优线性估计, 即克里格插值. 过去的研究通常采用仅考虑了空间相关特征的空间地统计学方法. 然而, 这种空间地统计学方法对存在于CO₂数据中的时间相关特征并没有充分考虑.     

多基线距DInSAR技术反演城市地表缓慢形变

利用多基线距DInSAR技术反演城市地表缓慢形变已经成为当前DInSAR研究的热点. 该技术的基本流程包括干涉图像对组合, 差分干涉图生成, 高相干点提取, Delaunay三角网建立, 网络增量计算和集成, 残余相位解缠和定标, 大气相位和非线性形变相位估计等. 其中网络增量计算是线性形变反演的关键, 残余相位解缠和定标是非线性形变反演的关键. 针对前人在这两方面研究的不足, 提出基于改进的模型相关系数方程的增量计算方法, 以及基于三角形“环”处理方法进行离散点相位解缠和定标的算法. 在上述算法的基础上, 利用24景ERS数据, 反演了苏州市在1993-02~2000-12之间的地面沉降. 结果表明, 在此期间内苏州市大部分城区平均沉降速率在20~30 mm/a之间, 而且沉降速率呈逐年递减趋势, 并在2000年之后趋于稳定. 水准数据证明了结果的准确性和算法的可行性.     

利用脉冲星观测检测和研究板块运动和地震

1 板块运动与地震甚长基线干涉测量(VLBI)、卫星激光测距(SLR)和全球定位系统(GPS)等空间新技术的发展开辟了监测和研究地球运动的新途径。利用这些空间技术高精度地测量地壳运动(板块运动和区域性地壳形变),进一步研究这些运动的力学机制以及与其他地球物理现象和地球内部结构的关系,是当今地球科学研究的前沿。     

GPS测定的尼泊尔M_w7.9和M_w7.3级地震同震形变场

2015年尼泊尔M w7.9和M w7.3级地震致灾范围包括尼泊尔、印度北部、巴基斯坦、孟加拉和中国藏南地区,地震应变应力调整将对震区和邻区的地震活动产生不同程度的影响.尼泊尔境内的GPS连续观测数据、"中国大陆构造环境监测网络"GPS基准站和地震应急流动站观测数据计算结果揭示了尼泊尔地震的静态和动态同震形变场.对M w7.9级地震,尼泊尔境内近场静态同震GPS水平位移最大为1.89 m,距震中100~400 km的中国藏南地区观测到几毫米到几十厘米的静态水平位移,毗邻尼泊尔的聂拉木县最大形变54.0 cm.M w7.3级地震静态同震形变最大为2 cm,局限在震中附近200 km范围内.距震中约2000 km内的GPS基准站均记录到M w7.9大震明显的动态形变信号,高频GPS动态形变幅度与地震破裂方向有关,位于破裂方向上的测站动态形变幅度明显大于其他方向的测站.用弹性半空间位错模型正演模拟了震区和青藏高原南部格网点上的同震形变,并分析了地震应变影响,认为尼泊尔地震对中国藏南地区产生一定程度的拉张型应力变化,需要持续关注.     

顶得住地震的新建筑

大多数摩天楼在地震时会倒塌。因此一些大学正就关于如何减缓地震扭转效应的新建议进行研究。这些建议包括:用缆索缚住建筑物的一侧,从而使之能在压力中相对拉紧;或者使空气迸裂,就象喷气发动机那样,以沿着建筑物的侧壁施加一个力同地震破坏力相抵销。     

张北震区附近形变的成因--与王超、张红等人商榷

<科学通报>2000年第23期王超等人[1]的<张北-尚义地震同震形变场雷达差分干涉测量>一文(后面简称"形变"),介绍了1998年1月10日张北地震附近的同震隆起形变,形变中心位于114°20'E,40°57'N,最大视向位移量达25 cm,形变集中在300km2范围内.     

嫦娥一号干涉成像光谱仪数据TiO2 反演初步结果

反演月表TiO₂含量的分布对于月球科学研究和未来月球资源开采与利用具有重要意义. 利用干涉成像光谱仪数据对月表TiO₂ 反演进行了初步尝试, 获得了基于月球着陆点数据和月面真值的TiO₂ 反演公式. 将该公式应用于嫦娥一号干涉成像光谱仪数据, 获得了TiO₂ 在Apollo 16 登陆点附近高地、月球澄海MS2 地区的低钛玄武岩和高钛玄武岩等地区的局域分布特征. 经与Clementine UVVIS 数据反演结果进行比较, 嫦娥一号干涉成像光谱仪数据的TiO₂ 含量初步模型对于所研究的高地地区预测偏高0.7 wt.%左右, 对于月海低钛地区的预测偏高1.5 wt.%左右, 对于月海高钛地区的含量预测偏低0.8 wt.%左右.     

青藏高原东缘地震各向异性、应力及汶川地震影响

使用国家地震台网固定台站和区域流动地震台站资料,通过对最长达18年记录开展的剪切波分裂分析,得到青藏高原东缘构造域(包括龙门山断裂带)的上地壳各向异性空间分布和时间变化特征.通过快剪切波(快波)偏振和慢剪切波(慢波)时间延迟参数,获得了快波偏振的分区分布特征和汶川地震前后不同分区统计特征的时间变化.剪切波分裂参数受到应力场和断裂构造的影响,块体边界或断裂附近的快波偏振特征比块体内部更为复杂,导致一些分区显示两个优势方向.在汶川MS8.0地震前后,可以观察到一些台站或区域的剪切波分裂参数的变化,既有快波偏振方向的改变,也有慢波时间延迟的变化.地震前后有明显变化的台站,基本符合两种情况:一是震中距较近;二是处于大的断裂或构造边界附近.研究发现,汶川地震前后,龙门山断裂带区域的慢波时间延迟的降幅显著大于周边区域;龙门山断裂带北段的慢波时间延迟降幅,大于龙门山断裂带中段,大于龙门山断裂带南端与鲜水河断裂、安宁河断裂的交汇区,这反映了地壳应力及介质物性状态的变化.通过分析剪切波分裂参数,推断出龙门山断裂带域及周边区域的主压应力分布.本文同时还分析了芦山地震前后剪切波分裂特征的变化.研究认为,利用地震剪切波特性监测应力变化,进而可应用于地震应力预测研究.需要注意是,不同的构造部位对应力变化有不同的响应.     

“嫦娥一号”干涉成像光谱仪数据FeO反演初步结果

月表FeO含量的分布对于未来月球资源利用和科学研究具有重要意义.在Apollo和Luna取样点的干涉成像光谱数据与月表FeO含量真值的基础上,得到了利用干涉成像光谱仪数据进行月表FeO填图的初步数学模型.将该公式应用于"嫦娥一号"干涉成像光谱仪数据,成功得到了月表FeO的局域地区分布.经过与Clementine UVVIS数据进行比较,嫦娥干涉成像光谱仪数据的FeO含量初步模型对于月海地区的预测偏低(2%),对于高地地区的预测偏高(3%).     

2011年3月11日日本地震震前电离层异常变化分析

以2011年3月11日日本本州东海岸附近海域发生的Mw9.0级地震为背景,利用震中附近GNSS观测站数据计算得到的VTEC和电离层测高仪站的foF2数据,分析了地震前后的电离层异常.发现3月8日存在明显的正异常,且明显的异常持续了6h.进一步分析电离层异常的空间分布后,认为3月8日中午至日落前后发生的电离层正异常很可能是由于孕育地震而引起的.通过给出的3月8日VTEC二维异常分布图发现异常的区域在经度方向跨度有50°,在纬度方向接近20°,长度大致为3:1,电离层VTEC异常峰值出现在15:00~19:00LT之间,峰值点并不位于震中的正上方,而是位于震中偏向赤道方向,在赤道共轭区也存在异常,而在全球其他区域没有显著的异常.为深入分析该天电离层异常的变化情况,采用电离层层析技术重构了震区附近电离层电子密度时空分布,确认3月8日的电离层异常可看作是日本311地震的前兆.     

星载合成孔径雷达高精度快速定位方法

星载合成孔径雷达(SAR)定位是其影像应用的基础, 研究高效能的定位方法对促进SAR 影像应用具有重要意义. 从SAR 成像机理出发, 分析揭示了SAR 成像几何的两个基本特性, 并对两个特性的合理性进行了论证. 基于此, 提出了一种高精度快速定位解算方法. 针对目前国际上最为典型的四种SAR 卫星, 德国TerraSAR-X、意大利COSMOSkyMed、日本ALOS-PalSAR 和加拿大Radarsat-2 进行了精度分析, 结果表明本文方法的精度完全能够满足这四种卫星的实际需求. 最后, 利用TerraSAR-X 聚束模式SAR 实测数据进行快速定位处理, 性能评估结果表明在极大降低运算量的情况下保持了极高的定位精度, 验证了方法的高效性和正确性.     

一种改进的地震灾害倒塌房屋遥感信息自动识别方法

地震灾害倒塌房屋的识别是地震灾害损失评价的重要环节, 也是判断地震强度的指标之一. 5.12汶川大地震发生后, 中欧地震灾害评估项目组使用自动识别技术处理0.5 m分辨率的ADS40航空遥感数据, 需要克服航空遥感数据中未倒塌混凝土桥梁、滑坡体滚石与倒塌房屋反射光谱的相似性. 经过多次实验改进后的技术, 一次性识别的倒塌房屋面积覆盖率优于95%. 介绍了这次实验过程和结果, 实验表明, 该方法也可以用于其他灾害造成的倒塌房屋的自动识别.