1.本发明涉及预测目的的数据处理
技术领域：
：，特别是涉及一种地震后物理参数变化的预测方法及系统。
背景技术：
：：2.本部分的陈述仅仅是提到了与本发明相关的
背景技术：
：，并不必然构成现有技术。3.地震作为普遍的自然灾害之一，具有突发性强、破坏性大、连锁强的特点，给人们生产生活以及社会经济发展都带来了不容小觑的影响。而震后极容易引发泥石流、滑坡、暴雨等与降水有关的次生灾害。而在这些与降水有关的极端天气和灾害中，大气可降水量与降水密切相关。4.同时，地震过程中会引发各种物理和化学变化，尤其为能量释放造成的震中地区地表温度的显著变化。大气可降水量主要与温度有关，随温度的变化而变化，温度升高，空气容纳水分的能力以克劳修斯-克拉佩伦速率增加。5.目前，针对地震过程中大气可降水量的研究甚少，大多数是利用地表温度探测地震前的热异常来研究地震预测，亦或者是仅对地震过程中水汽和地表温度两者之一的研究，并没有将两者互相结合来探测震后气象参量异常。6.然而，缺少结合地表温度热异常对震后水汽的研究，忽视了震后是在震中位置极大可能发生次生灾害的时间点，这非常不利于灾后救援与重建工作。技术实现要素：7.为了解决现有技术的不足，本发明提供了一种地震后物理参数变化的预测方法及系统；采用北斗卫星导航系统反演地震后震中及附近地区大气可降水量，并进一步分析大气可降水量与地表温度热异常关系。8.第一方面，本发明提供了一种地震后物理参数变化的预测方法；一种地震后物理参数变化的预测方法，包括：获取历史时间段内，目标监测区域地震发生后的大气可降水量；对历史时间段内采集的大气可降水量，进行精度验证；获取历史时间段内，目标监测区域地震发生后的地表温度热异常数据；基于历史地表温度热异常数据，和通过精度验证的历史大气可降水量，确定大气可降水量峰值与地表温度热异常峰值之间的时间差；当目标监测区域再次发生地震时，检测目标监测区域的地表温度热异常峰值对应的时间点；将时间点加上时间差得到大气可降水量峰值的预测时间点，将预测时间点作为次生灾害发生的时间点输出。9.第二方面，本发明提供了一种地震后物理参数变化的预测系统；一种地震后物理参数变化的预测系统，包括：第一获取模块，其被配置为：获取历史时间段内，目标监测区域地震发生后的大气可降水量；精度验证模块，其被配置为：对历史时间段内采集的大气可降水量，进行精度验证；第二获取模块，其被配置为：获取历史时间段内，目标监测区域地震发生后的地表温度热异常数据；时间差确定模块，其被配置为：基于历史地表温度热异常数据，和通过精度验证的历史大气可降水量，确定大气可降水量峰值与地表温度热异常峰值之间的时间差；预测输出模块，其被配置为：当目标监测区域再次发生地震时，检测目标监测区域的地表温度热异常峰值对应的时间点；将时间点加上时间差得到大气可降水量峰值的预测时间点，将预测时间点作为次生灾害发生的时间点输出。10.与现有技术相比，本发明的有益效果是：本发明提出了一种利用北斗卫星导航系统(bds)技术反演地震后震中及附近地区大气可降水量(pwv)并进一步分析其与地表温度热异常关系的方法。该方法可以真正获得bds反演的地震后震中地区大气可降水量，并与无线电探空仪(rs)和再分析资料(era5)数据对比验证精度，证明北斗卫星导航系统反演水汽技术的精度，从而为采用bds反演的大气可降水量(bds-pwv)来体现震后水汽变化具有信服性和参考价值。其次，结合地表温度产品，从两个时间序列尺度提取地震引发的热异常，一方面为在地震过程中地表温度对比，另一方面为基于10年背景场数据的多年同期同地区地表温度对比。最后，结合采用bds反演的大气可降水量bds-pwv变化与地震后震中位置地表温度热异常变化，获取大气可降水量与热异常的延迟间隔，从而为将来该区域地震发生后的次生灾害预测和灾后重建工作提供参考价值。附图说明11.构成本发明的一部分的说明书附图用来提供对本发明的进一步理解，本发明的示意性实施例及其说明用于解释本发明，并不构成对本发明的不当限定。12.图1为实施例一的方法流程图；图2(a)-图2(b)为实施例一的rtklib获取的ztd值与igs发布的ztd产品拟合图；图3(a)-图3(c)为实施例一的三次地震多源数据的pwv时间序列图；图4(a)-图4(f)为实施例一的反演获取的pwv与rs-pwv及era5-pwv精度验证。具体实施方式13.应该指出，以下详细说明都是示例性的，旨在对本发明提供进一步的说明。除非另有指明，本文使用的所有技术和科学术语具有与本发明所属
技术领域：
：的普通技术人员通常理解的相同含义。14.需要注意的是，这里所使用的术语仅是为了描述具体实施方式，而非意图限制根据本发明的示例性实施方式。如在这里所使用的，除非上下文另外明确指出，否则单数形式也意图包括复数形式，此外，还应当理解的是，术语“包括”和“具有”以及他们的任何变形，意图在于覆盖不排他的包含，例如，包含了一系列步骤或单元的过程、方法、系统、产品或设备不必限于清楚地列出的那些步骤或单元，而是可包括没有清楚地列出的或对于这些过程、方法、产品或设备固有的其它步骤或单元。15.在不冲突的情况下，本发明中的实施例及实施例中的特征可以相互组合。16.本实施例所有数据的获取都在符合法律法规和用户同意的基础上，对数据的合法应用。17.地震是由于地壳在地球构造力的作用下发生突然破裂和错位而引起的振动，导致能量迅速释放，引起许多物理和化学变化。震后又极容易发生泥石流、滑坡、暴雨等与降水有关的次生灾害，对灾后救援和重建工作产生很大影响。而在这些与降水有关的极端天气和灾害中，水汽发挥着必不可少的作用。水汽在生物圈中通过降水、地表蒸发、地表径流、水汽输送等参与地表系统与大气之间的水循环，在物质交换和能量传递过程中起着重要作用。作为一种中性气体，水汽虽然在大气中含量较少，但变化迅速，极为活跃，存在固、液、气三种相态，与众多天气变化密切相关。同时，在震中区域地表温度的变化是经常探测到的，水蒸气容易受到温度影响。因此探究震后大气可降水量变化对次生灾害应对和震后重建有非常重要的意义。大气可降水量是指测站上空整个空气柱中的水分全部凝结成雨、雪降落后形成的降水量。北斗卫星导航系统(beidounavigationsatellitesystem)简称bds。18.常用的大气可降水量观测技术有无线电探空技术、星载辐射计微波辐射计、无线电掩星技术，雷达和地基gnss(globalnavigationsatellitesystem)观测(包括bds)。gps气象学的提出为气象研究提供了一种新的高精度技术。bds于2020年6月完成全球组网，成功标志着中国卫星导航系统走向独立自主。高性能、高精度、高可靠的bds建设，既是对北斗系统的完善，更是为世界卫星导航发展贡献了中国智慧。北斗气象学是最早应用于bds的重要领域之一，现已成为北斗全面应用的新局面。其反演大气可降水量具有高精度、全天候、低成本等优点。卫星信号在穿过大气层时会发生信号折射，其中就包含对流层延迟，即卫星信号在穿过高度为50km以下的未被电离的中性大气层时所发生的延迟，使传播路径发生弯曲、传播速度发生变化。天顶对流层延迟(totalzenithdelay,ztd)由天顶干延迟(zenithhydrostaticdelay,zhd)和湿延迟(zenithwetdelay,zwd)组成，干延迟虽在对流层总延迟中占比大但相对稳定，湿延迟虽占比小，但由于水分含量多且非常活跃，使得其具有变化迅速的特点，可从中提取出大气可降水量。天顶对流层延迟ztd通过bds系统处理获得，zhd利用对流层经典模型获得，两者相减可获得湿延迟zwd，结合水汽转换系数可求得大气可降水量，结果为bds大气可降水量，记为bds-pwv。19.结合当前的研究，大部分都是将重点放在地震前的地震预测与监测，少有研究对地震发生后的物理参量和现象进行研究。据统计，地震引发的山体崩塌、滑坡和泥石流是目前造成人员伤亡最多、最严重的地震后次生灾害，5.0级左右的地震就可以诱发滑坡和泥石流，震级越大，诱发面积也越大。同时，震后极容易发生强降水，地壳积聚的能量不断释放，使得地表温度升高，促进水分蒸发进入到空气中，从而在地面横截面形成了很强的上升气流，空气柱的水汽越来越多，凝结后形成降水。这些与降水有关的次生灾害极大地影响了灾后救援和重建工作。目前，大多数研究是对地震前水汽的研究，亦或是对震前热异常的研究，少有对地震后水汽和地表温度的探讨，也缺乏将大气可降水量和基于不同时间序列的地表温度热异常结合起来进一步分析前后关系。虽然gnss技术在地震相关水汽中有较多的研究工作，但利用bds研究地震后的水汽变化特性还鲜有研究。20.术语解释：pwv，大气可降水量，precipitablewatervapor；bds,北斗卫星导航系统，beidounavigationsatellitesystem；lst，地表温度，landsurfacetemperature；ztd，天顶对流层延迟，zenithtroposphericdelay；rs，无线电探空仪，radiosonde；era5，第五代欧洲中期天气预报中心再分析资料，thefifthgenerationecmwfatmosphericreanalysisoftheglobalclimate；本实施例提供了一种地震后物理参数变化的预测方法；如图1所示，一种地震后物理参数变化的预测方法，包括：s101：获取历史时间段内，目标监测区域地震发生后的大气可降水量；s102：对历史时间段内采集的大气可降水量，进行精度验证；s103：获取历史时间段内，目标监测区域地震发生后的地表温度热异常数据；s104：基于历史地表温度热异常数据，和通过精度验证的历史大气可降水量，确定大气可降水量峰值与地表温度热异常峰值之间的时间差；s105：当目标监测区域再次发生地震时，检测目标监测区域的地表温度热异常峰值对应的时间点；将时间点加上时间差得到大气可降水量峰值的预测时间点，将预测时间点作为次生灾害发生的时间点输出。21.进一步地，s101：获取历史时间段内，目标监测区域地震发生后的大气可降水量；具体包括：s1011：采用北斗卫星导航系统，获取地震发生后目标监测区域的大气可降水量bds-pwv；s1012：基于无线电探空方式，获取地震发生后目标监测区域的大气可降水量rs-pwv；s1013：基于era5再分析资料，获取地震发生后目标监测区域的大气可降水量era5-pwv。22.进一步地，所述s1011：采用北斗卫星导航系统，获取地震发生后目标监测区域的大气可降水量bds-pwv；具体包括：s1011-1：根据北斗卫星导航系统观测数据，获取地震发生后目标监测区域的天顶对流层延迟ztd；s1011-2：利用saastamoinen对流层经典改正模型和姚宜斌大气加权平均温度模型，最终获得大气可降水量bds-pwv。23.进一步地，所述利用saastamoinen对流层经典改正模型和姚宜斌大气加权平均温度模型，最终获得大气可降水量bds-pwv；具体包括：ztd由zhd和zwd组成，zhd采用saastamoinen经典对流层模型，输入纬度、气压和高度计算得出：ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ(1)(2)其中，是测站气压；是和的函数，是测站弧度制纬度，是大地高。24.zwd含有大量水汽而变化迅速，用间接方法获取，如下式：ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ(3)pwv定义为测站上方空气柱中水汽全部凝结成雨、雪降落后形成的降水量，表示为zwd的函数：ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ(4)其中，为水汽转换系数，由公式(5)计算得出：ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ(5)其中，；是水汽比常数；和是大气折射率常数，分别为和；是大气加权平均温度。25.采用姚宜斌提出的全球纬度相关的大气加权平均温度模型，选择研究区所在纬度范围内的模型公式：ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ(6)其中，是地表温度。26.最终获取的pwv记为bds-pwv，其中bds表示北斗卫星导航系统。27.图2(a)-图2(b)为实施例一的rtklib获取的ztd值与igs发布的ztd产品拟合图。28.进一步地，所述s1012：基于无线电探空方式，获取地震发生后目标监测区域的大气可降水量rs-pwv；具体包括：无线电探空资料是通过在每日的0点和12点释放探空气球获取各种气象数据值。29.虽然精度高，但该方法结果的时间分辨率低。探空资料获取的pwv记为rs-pwv，rs表示无线电探空。30.进一步地，所述s1013：基于era5再分析资料，获取地震发生后目标监测区域的大气可降水量era5-pwv；具体包括：采用era5再分析资料的空气柱水汽(totalcolumnwatervapor,tcwv)，提取震中目标区域的水汽含量。31.获取的pwv记为era5-pwv，era5表示第五代欧洲中期天气预报中心再分析资料。32.进一步地，s102：对历史时间段内采集的大气可降水量，进行精度验证；具体包括：s1021：对采用北斗卫星导航系统获取的大气可降水量bds-pwv、基于无线电探空方式获取的大气可降水量rs-pwv和基于era5再分析资料获取的大气可降水量era5-pwv进行高程统一；s1022：计算采用北斗卫星导航系统获取的大气可降水量bds-pwv与基于无线电探空方式获取的大气可降水量rs-pwv的第一差值；s1023：计算采用北斗卫星导航系统获取的大气可降水量bds-pwv与基于era5再分析资料获取的大气可降水量era5-pwv的第二差值；s1024：如果第一差值和第二差值均在设定的范围内，则表示采用北斗卫星导航系统获取的大气可降水量反演的大气可降水量bds-pwv精度验证通过；否则表示精度验证不通过。33.进一步地，所述s1021：对采用北斗卫星导航系统获取的大气可降水量bds-pwv、基于无线电探空方式获取的大气可降水量rs-pwv和基于era5再分析资料获取的大气可降水量era5-pwv进行高程统一；具体包括：s1021-1：将era5数据的坐标系统由位势统一到bds坐标系统大地高下；具体包含：s1021-11：将位势转换为位势高，公式如下：其中，gph为位势高；gh为位势；为重力加速度；s1021-12：将位势高转换为正高，公式如下：12：将位势高转换为正高，公式如下：12：将位势高转换为正高，公式如下：其中，是正高；是纬度为时地球的有效半径；=9.806665；为旋转椭球表面的正常重力值；s1021-13：采用egm2008模型计算大地水准面差距n。在matlab中使用geoidegm2008grid.mat功能，输入纬度、经度，选择egm2008即可获得大地水准面差距n。34.s1021-14：将正高系统转换到大地高系统下，公式如下：其中，为大地高。35.s1021-2：将rs数据的坐标系统由位势高统一到bds坐标系统大地高下；具体包含：s1021-21：将位势高转换为正高，公式如下：21：将位势高转换为正高，公式如下：其中，是正高；是纬度为时地球的有效半径；=9.806665；为旋转椭球表面的正常重力值；s1021-22：采用egm2008模型计算大地水准面差距n。在matlab中使用geoidegm2008grid.mat功能，输入纬度、经度，选择egm2008即可获得大地水准面差距n。36.s1021-23：将正高系统转换到大地高系统下，公式如下：其中，为大地高。37.s1021-3：多源数据(rs和era5)的pwv高程统一，获取到皆高程统一到大地高系统下的pwv，公式如下：其中，为高程统一后pwv；为原始pwv；为站点gnss的大地高；为rs或era5对应的大地高；最终获取高程统一的bds-pwv、rs-pwv和era5-pwv。38.进一步地，所述s103：获取历史时间段内，目标监测区域地震发生后的地表温度热异常数据；具体是采用modis(moderateresolutionimagingspectroradiometer)传感器提取地震发生后目标监测区域的地表温度热异常数据。39.进一步地，所述s103：获取历史时间段内，目标监测区域地震发生后的地表温度热异常数据；具体包括：s1031：对地震过程中的地表温度值进行前后时间序列对比，得到第一对比结果；s1032：将地震发生后地震当天的地表温度与过去十年同一时间点同一地理位置的地表温度值进行对比，得到第二对比结果；s1033：将第一对比结果和第二对比结果取并集，得出地表温度热异常结果。40.应理解地，基于modis11a1(modis/terralst/edailyl3global1kmsingrid)单日地表温度l3(三级产品)全球网格数据，利用地表温度在目标区域从两个时间序列来探测震后热异常，一方面为在地震过程中的地表温度对比，另一方面为基于十年背景场的同期同地的地表温度对比。41.进一步地，所述s1031：对地震过程中的地表温度值进行前后时间序列对比，得到第一对比结果；具体包括：针对地震发生过程中前后n天内的地表温度前后对比，从而获取震后地表温度较震前和发震时的差异；得出设定时间范围的地表温度热异常结果，n为正整数。42.进一步地，所述s1032：将地震发生后地震当天的地表温度与过去十年同一时间点同一地理位置的地表温度值进行对比，得到第二对比结果；具体包括：将地震发生后地震当天的地表温度，与对过去十年的地表温度值建立背景场进行多年同期对比，将alice指标作为热异常指标，得到地表温度热异常结果。43.例如第i年5月6日地震，那么我们以十年每年5月6日的地表温度值建立背景场，再以此来对比第i年5月6日这一天的异常；具体包括：将地震发生后的地表温度和基于10年背景场与地震同期同地的地表温度进行对比，突出时空域上的异常，提取空间热异常时间点，将alice(absolutelyllocalindexofchangeoftheenvironment)指数作为热异常指标：ꢀꢀꢀꢀꢀꢀꢀꢀꢀ(7)其中，是指地震时期的lst值；为多年平均lst值；是多年地表温度的标准差。44.进一步地，所述s104：基于历史地表温度热异常数据，和通过精度验证的历史大气可降水量，确定大气可降水量峰值与地表温度热异常峰值之间的时间差；具体包括：s1041：获取每一年大气可降水量峰值与地表温度热异常峰值之间的时间差；s1042：对所有年份的时间差求平均值，得到平均时间差，将平均时间差作为大气可降水量峰值与地表温度热异常峰值之间的时间延迟。45.bds-pwv具有升高变化时，有极大可能发生降水相关的次生灾害。而在pwv变化之前往往伴随有lst热异常出现，因此可以通过先捕获lst热异常现象，再根据pwv峰值与热异常出现两者之间的时间延迟判定pwv突增至峰值的时间点，从而采取相关措施来预防次生灾害的发生，从而达到预测的目的。bds-pwv峰值变化和热异常出现的前后时间延迟。热异常出现分别为：1.lst在地震发生后的增加现象；2.以alice指数作为指标的高值。46.结合北斗卫星导航系统监测的大气可降水量的变化与地震后震中位置地表温度热异常的变化，获取大气可降水量与热异常的延迟信息，从而为将来该区域地震发生后的次生灾害预测和灾后重建工作提供参考价值。47.本发明基于北斗卫星导航系统观测数据和地表温度产品，提出了一种利用北斗卫星导航系统技术反演地震后震中及附近地区大气可降水量并进一步分析其与地表温度热异常关系的方法。48.本发明基于陆态网北斗卫星导航系统观测数据和mod11a1单日地表温度格网产品，针对在某地区发生的三次地震提出了一种利用北斗卫星导航系统技术反演震后震中及附近地区大气可降水量并进一步分析其与地表温度热异常关系的方法。其中，将北斗卫星导航系统技术反演的大气可降水量与rs及era5数据进行比对，验证了bds-pwv精度良好，可用来描述地震后水汽变化。本发明属于北斗卫星导航系统与遥感交叉的大气水汽研究范畴。49.随机选择三个发生过地震的区域，分别利用距离三个地震震中位置最近的北斗卫星导航系统的观测数据，按照反演流程获取bds-pwv。50.根据图3(a)-图3(c)，第一区域、第二区域和第三区域地震后的所有大气可降水量数据都出现了异常增加的情况。在第一区域地震爆发时，bds-pwv、rs-pwv和era5-pwv均处于最小值，在此之后便出现了pwv上升的变化，呈上升趋势，并在大地震发生后3天由于5月7日的余震从而达到峰值（bds-pwv约14mm，rs和era5约为11mm）。余震进一步促进能量释放，水汽不断碰撞和聚集，形成可降水。bds-pwv在第二区域地震发生三天后出现了急剧的上升，最高可达6mm，是整个地震过程中的最高值，同时，相比于rs-pwv和era5-pwv，bds-pwv的变化明显要明显很多。对于第三区域地震，三个数据源的pwv都在12月18日达到最大值（地震发生在12月13日23时，持续到第二天凌晨2时），其中bds和era5的pwv约为8.2mm，rs探测结果基本达到10mm，然后又急剧下降。51.上述分析表明，地震后存在大气可降水量异常，主要表现为先增加随后减少的趋势。地震导致岩石在构造应力作用下孔隙裂缝开闭，释放出大量地热能，如热气和热水。热效应的产生，加速震中上空水分的蒸发和水汽的运动，这也是地震过程中大气可降水量异常变化的原因。52.bds-pwv的坐标系统为大地高，探空资料的坐标系统为位势高，而era5再分析资料则为位势。高程统一到大地高后，进行精度验证。53.相关系数和误差如图4(a)-图4(f)所示。第一区域地震的bds-pwv与rs-pwv和era5-pwv非常接近，rmse分别为1.28mm和1.78mm。在第二区域地震中，bds-pwv与rs-pwv和era5-pwv的rmse分别为1.06mm和0.7mm。第三区域地震的rmse则为1.34mm和0.97mm。从以上分析可以看出，北斗卫星导航系统反演大气可降水量与其他两个数据源结果的rmse均小于2mm，其反演的准确性和可靠性非常高，能够为地震过程的气象研究提供有效大气可降水量参考。54.通过分析三个地震所对应的地表温度在地震过程中的变化，可以看出，第一区域地震爆发当天（st呈现大面积最大值为-6~-2，比地震前三天高约12~16。对于第二区域地震，震中位置lst有向西递增的趋势，稳定在-10~6，但在震后5天之内又开始逐渐下降，随后直到2月20日，lst恢复到震前水平。第三区域地震中，lst在震中及周边地区上升并持续了3天。从整个过程来看，在地震发生的当天或次日，lst都比前几天高。三次地震后的lst都比震前存在一个明显上升，并在高值处维持了三天左右，然后又下降到正常水平，其原因可能为地震过程中释放的能量在震后几天逐渐消散。55.其次，本发明又采用rst(robustsatellitetechnique，卫星鲁棒分析技术)算法，基于十年背景场数据来突出多年同地震时期、同地震位置的空间热异常。将alice指数作为热异常指标。56.第一区域地震过程的alice指数显示，某月6日和7日震中及周边地区出现严重的热异常，最大值为1.5，基于十年背景场出现了地表温度热异常；虽然8日至10日有云层覆盖缺失数据，但从11日开始，异常强度逐渐减弱。57.第二区域地震过程中的alice指数在某月13日至15日均为1以上，随后恢复正常（0±0.5），说明某年发生该地震时的lst比过去十年同范围区域的lst都高，存在热异常现象。58.第三区域地震后，震中及附近地区的alice指数上升至1，震后lst高于多年lst，基于10年数据的热异常相当明显。59.总之，地震后alice指数有所上升，地表温度热异常高于往年同期。此后几乎没有任何变化。60.以上分别从地震发生过程中lst前后对比和基于十年背景数据的多年同期同地lst对比两个方面提取了地震热异常，发现地震能够引起地表温度热异常。61.第一区域震中区最大lst在某月6日，同时pwv开始增加，次日达到第一个峰值，达到最大pwv峰值是在某月9日。62.对于第二区域地震，最高lst是在发震当天测得，bds-pwv在某月16日达到最大。63.第三区域地震于某月13日23:00发生，余震一直持续到次日凌晨2时，lst在震中及周边地区在15日最高，而在某月18日bds-pwv达到峰值，比lst最大值晚3天。64.综上所述，实验所选取的3次地震发生后lst和bds-pwv变化之间存在延迟，即大气可降水量变化紧随地表温度变化之后。bds-pwv变化与地表温度热异常存在一定的关系，大气可降水量峰值出现晚于地表温度升高3天。65.参考背景数据发现，震中地表温度热异常越严重，bds-pwv越高；否则，地表温度异常较弱或没有异常，则在震中没有bds-pwv的突出异常变化。66.地震诱发地表和地下热异常，引发更多的热能释放反应，促进水汽活跃，特别是水汽向上运动过程中的蒸发、输送和积聚，大量pwv沉积在大气中，超过一定程度可以产生降水。通过北斗卫星导航系统技术反演发现，震后pwv会增加，与地表温度热异常存在延迟关系。基于这一发现，本发明认为这对震后与降水相关的次生灾害的预测预防和震后重建具有重要意义。67.对待研究区域，就按照上述流程结合当前待研究区域的历史地震数据，得出bds-pwv峰值和热异常出现的结果，最终获得延迟时间；基于延迟时间，为以后当前待研究区域再次发生地震时通过提前检测到热异常，在加以历史延迟时间计算出pwv峰值时间点，预测次生灾害。68.本发明利用北斗卫星导航系统反演震后震中位置的pwv，并结合地表温度热异常进一步探求其与pwv关系。69.本发明北斗卫星导航系统反演获得震后震中及附近地区pwv并进一步分析其与地表温度热异常关系，为震后与降水相关的次生灾害的预测监测和震后重建工作提供参考。70.实施例二本实施例提供了一种地震后物理参数变化的预测系统；一种地震后物理参数变化的预测系统，包括：第一获取模块，其被配置为：获取历史时间段内，目标监测区域地震发生后的大气可降水量；精度验证模块，其被配置为：对历史时间段内采集的大气可降水量，进行精度验证；第二获取模块，其被配置为：获取历史时间段内，目标监测区域地震发生后的地表温度热异常数据；时间差确定模块，其被配置为：基于历史地表温度热异常数据，和通过精度验证的历史大气可降水量，确定大气可降水量峰值与地表温度热异常峰值之间的时间差；预测输出模块，其被配置为：当目标监测区域再次发生地震时，检测目标监测区域的地表温度热异常峰值对应的时间点；将时间点加上时间差得到大气可降水量峰值的预测时间点，将预测时间点作为次生灾害发生的时间点输出。71.此处需要说明的是，上述第一获取模块、精度验证模块、第二获取模块、时间差确定模块和预测输出模块对应于实施例一中的步骤s101至s105，上述模块与对应的步骤所实现的示例和应用场景相同，但不限于上述实施例一所公开的内容。需要说明的是，上述模块作为系统的一部分可以在诸如一组计算机可执行指令的计算机系统中执行。72.以上所述仅为本发明的优选实施例而已，并不用于限制本发明，对于本领域的技术人员来说，本发明可以有各种更改和变化。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。当前第1页12当前第1页12
技术领域：
：，特别是涉及一种地震后物理参数变化的预测方法及系统。
背景技术：
：：2.本部分的陈述仅仅是提到了与本发明相关的
背景技术：
：，并不必然构成现有技术。3.地震作为普遍的自然灾害之一，具有突发性强、破坏性大、连锁强的特点，给人们生产生活以及社会经济发展都带来了不容小觑的影响。而震后极容易引发泥石流、滑坡、暴雨等与降水有关的次生灾害。而在这些与降水有关的极端天气和灾害中，大气可降水量与降水密切相关。4.同时，地震过程中会引发各种物理和化学变化，尤其为能量释放造成的震中地区地表温度的显著变化。大气可降水量主要与温度有关，随温度的变化而变化，温度升高，空气容纳水分的能力以克劳修斯-克拉佩伦速率增加。5.目前，针对地震过程中大气可降水量的研究甚少，大多数是利用地表温度探测地震前的热异常来研究地震预测，亦或者是仅对地震过程中水汽和地表温度两者之一的研究，并没有将两者互相结合来探测震后气象参量异常。6.然而，缺少结合地表温度热异常对震后水汽的研究，忽视了震后是在震中位置极大可能发生次生灾害的时间点，这非常不利于灾后救援与重建工作。技术实现要素：7.为了解决现有技术的不足，本发明提供了一种地震后物理参数变化的预测方法及系统；采用北斗卫星导航系统反演地震后震中及附近地区大气可降水量，并进一步分析大气可降水量与地表温度热异常关系。8.第一方面，本发明提供了一种地震后物理参数变化的预测方法；一种地震后物理参数变化的预测方法，包括：获取历史时间段内，目标监测区域地震发生后的大气可降水量；对历史时间段内采集的大气可降水量，进行精度验证；获取历史时间段内，目标监测区域地震发生后的地表温度热异常数据；基于历史地表温度热异常数据，和通过精度验证的历史大气可降水量，确定大气可降水量峰值与地表温度热异常峰值之间的时间差；当目标监测区域再次发生地震时，检测目标监测区域的地表温度热异常峰值对应的时间点；将时间点加上时间差得到大气可降水量峰值的预测时间点，将预测时间点作为次生灾害发生的时间点输出。9.第二方面，本发明提供了一种地震后物理参数变化的预测系统；一种地震后物理参数变化的预测系统，包括：第一获取模块，其被配置为：获取历史时间段内，目标监测区域地震发生后的大气可降水量；精度验证模块，其被配置为：对历史时间段内采集的大气可降水量，进行精度验证；第二获取模块，其被配置为：获取历史时间段内，目标监测区域地震发生后的地表温度热异常数据；时间差确定模块，其被配置为：基于历史地表温度热异常数据，和通过精度验证的历史大气可降水量，确定大气可降水量峰值与地表温度热异常峰值之间的时间差；预测输出模块，其被配置为：当目标监测区域再次发生地震时，检测目标监测区域的地表温度热异常峰值对应的时间点；将时间点加上时间差得到大气可降水量峰值的预测时间点，将预测时间点作为次生灾害发生的时间点输出。10.与现有技术相比，本发明的有益效果是：本发明提出了一种利用北斗卫星导航系统(bds)技术反演地震后震中及附近地区大气可降水量(pwv)并进一步分析其与地表温度热异常关系的方法。该方法可以真正获得bds反演的地震后震中地区大气可降水量，并与无线电探空仪(rs)和再分析资料(era5)数据对比验证精度，证明北斗卫星导航系统反演水汽技术的精度，从而为采用bds反演的大气可降水量(bds-pwv)来体现震后水汽变化具有信服性和参考价值。其次，结合地表温度产品，从两个时间序列尺度提取地震引发的热异常，一方面为在地震过程中地表温度对比，另一方面为基于10年背景场数据的多年同期同地区地表温度对比。最后，结合采用bds反演的大气可降水量bds-pwv变化与地震后震中位置地表温度热异常变化，获取大气可降水量与热异常的延迟间隔，从而为将来该区域地震发生后的次生灾害预测和灾后重建工作提供参考价值。附图说明11.构成本发明的一部分的说明书附图用来提供对本发明的进一步理解，本发明的示意性实施例及其说明用于解释本发明，并不构成对本发明的不当限定。12.图1为实施例一的方法流程图；图2(a)-图2(b)为实施例一的rtklib获取的ztd值与igs发布的ztd产品拟合图；图3(a)-图3(c)为实施例一的三次地震多源数据的pwv时间序列图；图4(a)-图4(f)为实施例一的反演获取的pwv与rs-pwv及era5-pwv精度验证。具体实施方式13.应该指出，以下详细说明都是示例性的，旨在对本发明提供进一步的说明。除非另有指明，本文使用的所有技术和科学术语具有与本发明所属
技术领域：
：的普通技术人员通常理解的相同含义。14.需要注意的是，这里所使用的术语仅是为了描述具体实施方式，而非意图限制根据本发明的示例性实施方式。如在这里所使用的，除非上下文另外明确指出，否则单数形式也意图包括复数形式，此外，还应当理解的是，术语“包括”和“具有”以及他们的任何变形，意图在于覆盖不排他的包含，例如，包含了一系列步骤或单元的过程、方法、系统、产品或设备不必限于清楚地列出的那些步骤或单元，而是可包括没有清楚地列出的或对于这些过程、方法、产品或设备固有的其它步骤或单元。15.在不冲突的情况下，本发明中的实施例及实施例中的特征可以相互组合。16.本实施例所有数据的获取都在符合法律法规和用户同意的基础上，对数据的合法应用。17.地震是由于地壳在地球构造力的作用下发生突然破裂和错位而引起的振动，导致能量迅速释放，引起许多物理和化学变化。震后又极容易发生泥石流、滑坡、暴雨等与降水有关的次生灾害，对灾后救援和重建工作产生很大影响。而在这些与降水有关的极端天气和灾害中，水汽发挥着必不可少的作用。水汽在生物圈中通过降水、地表蒸发、地表径流、水汽输送等参与地表系统与大气之间的水循环，在物质交换和能量传递过程中起着重要作用。作为一种中性气体，水汽虽然在大气中含量较少，但变化迅速，极为活跃，存在固、液、气三种相态，与众多天气变化密切相关。同时，在震中区域地表温度的变化是经常探测到的，水蒸气容易受到温度影响。因此探究震后大气可降水量变化对次生灾害应对和震后重建有非常重要的意义。大气可降水量是指测站上空整个空气柱中的水分全部凝结成雨、雪降落后形成的降水量。北斗卫星导航系统(beidounavigationsatellitesystem)简称bds。18.常用的大气可降水量观测技术有无线电探空技术、星载辐射计微波辐射计、无线电掩星技术，雷达和地基gnss(globalnavigationsatellitesystem)观测(包括bds)。gps气象学的提出为气象研究提供了一种新的高精度技术。bds于2020年6月完成全球组网，成功标志着中国卫星导航系统走向独立自主。高性能、高精度、高可靠的bds建设，既是对北斗系统的完善，更是为世界卫星导航发展贡献了中国智慧。北斗气象学是最早应用于bds的重要领域之一，现已成为北斗全面应用的新局面。其反演大气可降水量具有高精度、全天候、低成本等优点。卫星信号在穿过大气层时会发生信号折射，其中就包含对流层延迟，即卫星信号在穿过高度为50km以下的未被电离的中性大气层时所发生的延迟，使传播路径发生弯曲、传播速度发生变化。天顶对流层延迟(totalzenithdelay,ztd)由天顶干延迟(zenithhydrostaticdelay,zhd)和湿延迟(zenithwetdelay,zwd)组成，干延迟虽在对流层总延迟中占比大但相对稳定，湿延迟虽占比小，但由于水分含量多且非常活跃，使得其具有变化迅速的特点，可从中提取出大气可降水量。天顶对流层延迟ztd通过bds系统处理获得，zhd利用对流层经典模型获得，两者相减可获得湿延迟zwd，结合水汽转换系数可求得大气可降水量，结果为bds大气可降水量，记为bds-pwv。19.结合当前的研究，大部分都是将重点放在地震前的地震预测与监测，少有研究对地震发生后的物理参量和现象进行研究。据统计，地震引发的山体崩塌、滑坡和泥石流是目前造成人员伤亡最多、最严重的地震后次生灾害，5.0级左右的地震就可以诱发滑坡和泥石流，震级越大，诱发面积也越大。同时，震后极容易发生强降水，地壳积聚的能量不断释放，使得地表温度升高，促进水分蒸发进入到空气中，从而在地面横截面形成了很强的上升气流，空气柱的水汽越来越多，凝结后形成降水。这些与降水有关的次生灾害极大地影响了灾后救援和重建工作。目前，大多数研究是对地震前水汽的研究，亦或是对震前热异常的研究，少有对地震后水汽和地表温度的探讨，也缺乏将大气可降水量和基于不同时间序列的地表温度热异常结合起来进一步分析前后关系。虽然gnss技术在地震相关水汽中有较多的研究工作，但利用bds研究地震后的水汽变化特性还鲜有研究。20.术语解释：pwv，大气可降水量，precipitablewatervapor；bds,北斗卫星导航系统，beidounavigationsatellitesystem；lst，地表温度，landsurfacetemperature；ztd，天顶对流层延迟，zenithtroposphericdelay；rs，无线电探空仪，radiosonde；era5，第五代欧洲中期天气预报中心再分析资料，thefifthgenerationecmwfatmosphericreanalysisoftheglobalclimate；本实施例提供了一种地震后物理参数变化的预测方法；如图1所示，一种地震后物理参数变化的预测方法，包括：s101：获取历史时间段内，目标监测区域地震发生后的大气可降水量；s102：对历史时间段内采集的大气可降水量，进行精度验证；s103：获取历史时间段内，目标监测区域地震发生后的地表温度热异常数据；s104：基于历史地表温度热异常数据，和通过精度验证的历史大气可降水量，确定大气可降水量峰值与地表温度热异常峰值之间的时间差；s105：当目标监测区域再次发生地震时，检测目标监测区域的地表温度热异常峰值对应的时间点；将时间点加上时间差得到大气可降水量峰值的预测时间点，将预测时间点作为次生灾害发生的时间点输出。21.进一步地，s101：获取历史时间段内，目标监测区域地震发生后的大气可降水量；具体包括：s1011：采用北斗卫星导航系统，获取地震发生后目标监测区域的大气可降水量bds-pwv；s1012：基于无线电探空方式，获取地震发生后目标监测区域的大气可降水量rs-pwv；s1013：基于era5再分析资料，获取地震发生后目标监测区域的大气可降水量era5-pwv。22.进一步地，所述s1011：采用北斗卫星导航系统，获取地震发生后目标监测区域的大气可降水量bds-pwv；具体包括：s1011-1：根据北斗卫星导航系统观测数据，获取地震发生后目标监测区域的天顶对流层延迟ztd；s1011-2：利用saastamoinen对流层经典改正模型和姚宜斌大气加权平均温度模型，最终获得大气可降水量bds-pwv。23.进一步地，所述利用saastamoinen对流层经典改正模型和姚宜斌大气加权平均温度模型，最终获得大气可降水量bds-pwv；具体包括：ztd由zhd和zwd组成，zhd采用saastamoinen经典对流层模型，输入纬度、气压和高度计算得出：ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ(1)(2)其中，是测站气压；是和的函数，是测站弧度制纬度，是大地高。24.zwd含有大量水汽而变化迅速，用间接方法获取，如下式：ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ(3)pwv定义为测站上方空气柱中水汽全部凝结成雨、雪降落后形成的降水量，表示为zwd的函数：ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ(4)其中，为水汽转换系数，由公式(5)计算得出：ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ(5)其中，；是水汽比常数；和是大气折射率常数，分别为和；是大气加权平均温度。25.采用姚宜斌提出的全球纬度相关的大气加权平均温度模型，选择研究区所在纬度范围内的模型公式：ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ(6)其中，是地表温度。26.最终获取的pwv记为bds-pwv，其中bds表示北斗卫星导航系统。27.图2(a)-图2(b)为实施例一的rtklib获取的ztd值与igs发布的ztd产品拟合图。28.进一步地，所述s1012：基于无线电探空方式，获取地震发生后目标监测区域的大气可降水量rs-pwv；具体包括：无线电探空资料是通过在每日的0点和12点释放探空气球获取各种气象数据值。29.虽然精度高，但该方法结果的时间分辨率低。探空资料获取的pwv记为rs-pwv，rs表示无线电探空。30.进一步地，所述s1013：基于era5再分析资料，获取地震发生后目标监测区域的大气可降水量era5-pwv；具体包括：采用era5再分析资料的空气柱水汽(totalcolumnwatervapor,tcwv)，提取震中目标区域的水汽含量。31.获取的pwv记为era5-pwv，era5表示第五代欧洲中期天气预报中心再分析资料。32.进一步地，s102：对历史时间段内采集的大气可降水量，进行精度验证；具体包括：s1021：对采用北斗卫星导航系统获取的大气可降水量bds-pwv、基于无线电探空方式获取的大气可降水量rs-pwv和基于era5再分析资料获取的大气可降水量era5-pwv进行高程统一；s1022：计算采用北斗卫星导航系统获取的大气可降水量bds-pwv与基于无线电探空方式获取的大气可降水量rs-pwv的第一差值；s1023：计算采用北斗卫星导航系统获取的大气可降水量bds-pwv与基于era5再分析资料获取的大气可降水量era5-pwv的第二差值；s1024：如果第一差值和第二差值均在设定的范围内，则表示采用北斗卫星导航系统获取的大气可降水量反演的大气可降水量bds-pwv精度验证通过；否则表示精度验证不通过。33.进一步地，所述s1021：对采用北斗卫星导航系统获取的大气可降水量bds-pwv、基于无线电探空方式获取的大气可降水量rs-pwv和基于era5再分析资料获取的大气可降水量era5-pwv进行高程统一；具体包括：s1021-1：将era5数据的坐标系统由位势统一到bds坐标系统大地高下；具体包含：s1021-11：将位势转换为位势高，公式如下：其中，gph为位势高；gh为位势；为重力加速度；s1021-12：将位势高转换为正高，公式如下：12：将位势高转换为正高，公式如下：12：将位势高转换为正高，公式如下：其中，是正高；是纬度为时地球的有效半径；=9.806665；为旋转椭球表面的正常重力值；s1021-13：采用egm2008模型计算大地水准面差距n。在matlab中使用geoidegm2008grid.mat功能，输入纬度、经度，选择egm2008即可获得大地水准面差距n。34.s1021-14：将正高系统转换到大地高系统下，公式如下：其中，为大地高。35.s1021-2：将rs数据的坐标系统由位势高统一到bds坐标系统大地高下；具体包含：s1021-21：将位势高转换为正高，公式如下：21：将位势高转换为正高，公式如下：其中，是正高；是纬度为时地球的有效半径；=9.806665；为旋转椭球表面的正常重力值；s1021-22：采用egm2008模型计算大地水准面差距n。在matlab中使用geoidegm2008grid.mat功能，输入纬度、经度，选择egm2008即可获得大地水准面差距n。36.s1021-23：将正高系统转换到大地高系统下，公式如下：其中，为大地高。37.s1021-3：多源数据(rs和era5)的pwv高程统一，获取到皆高程统一到大地高系统下的pwv，公式如下：其中，为高程统一后pwv；为原始pwv；为站点gnss的大地高；为rs或era5对应的大地高；最终获取高程统一的bds-pwv、rs-pwv和era5-pwv。38.进一步地，所述s103：获取历史时间段内，目标监测区域地震发生后的地表温度热异常数据；具体是采用modis(moderateresolutionimagingspectroradiometer)传感器提取地震发生后目标监测区域的地表温度热异常数据。39.进一步地，所述s103：获取历史时间段内，目标监测区域地震发生后的地表温度热异常数据；具体包括：s1031：对地震过程中的地表温度值进行前后时间序列对比，得到第一对比结果；s1032：将地震发生后地震当天的地表温度与过去十年同一时间点同一地理位置的地表温度值进行对比，得到第二对比结果；s1033：将第一对比结果和第二对比结果取并集，得出地表温度热异常结果。40.应理解地，基于modis11a1(modis/terralst/edailyl3global1kmsingrid)单日地表温度l3(三级产品)全球网格数据，利用地表温度在目标区域从两个时间序列来探测震后热异常，一方面为在地震过程中的地表温度对比，另一方面为基于十年背景场的同期同地的地表温度对比。41.进一步地，所述s1031：对地震过程中的地表温度值进行前后时间序列对比，得到第一对比结果；具体包括：针对地震发生过程中前后n天内的地表温度前后对比，从而获取震后地表温度较震前和发震时的差异；得出设定时间范围的地表温度热异常结果，n为正整数。42.进一步地，所述s1032：将地震发生后地震当天的地表温度与过去十年同一时间点同一地理位置的地表温度值进行对比，得到第二对比结果；具体包括：将地震发生后地震当天的地表温度，与对过去十年的地表温度值建立背景场进行多年同期对比，将alice指标作为热异常指标，得到地表温度热异常结果。43.例如第i年5月6日地震，那么我们以十年每年5月6日的地表温度值建立背景场，再以此来对比第i年5月6日这一天的异常；具体包括：将地震发生后的地表温度和基于10年背景场与地震同期同地的地表温度进行对比，突出时空域上的异常，提取空间热异常时间点，将alice(absolutelyllocalindexofchangeoftheenvironment)指数作为热异常指标：ꢀꢀꢀꢀꢀꢀꢀꢀꢀ(7)其中，是指地震时期的lst值；为多年平均lst值；是多年地表温度的标准差。44.进一步地，所述s104：基于历史地表温度热异常数据，和通过精度验证的历史大气可降水量，确定大气可降水量峰值与地表温度热异常峰值之间的时间差；具体包括：s1041：获取每一年大气可降水量峰值与地表温度热异常峰值之间的时间差；s1042：对所有年份的时间差求平均值，得到平均时间差，将平均时间差作为大气可降水量峰值与地表温度热异常峰值之间的时间延迟。45.bds-pwv具有升高变化时，有极大可能发生降水相关的次生灾害。而在pwv变化之前往往伴随有lst热异常出现，因此可以通过先捕获lst热异常现象，再根据pwv峰值与热异常出现两者之间的时间延迟判定pwv突增至峰值的时间点，从而采取相关措施来预防次生灾害的发生，从而达到预测的目的。bds-pwv峰值变化和热异常出现的前后时间延迟。热异常出现分别为：1.lst在地震发生后的增加现象；2.以alice指数作为指标的高值。46.结合北斗卫星导航系统监测的大气可降水量的变化与地震后震中位置地表温度热异常的变化，获取大气可降水量与热异常的延迟信息，从而为将来该区域地震发生后的次生灾害预测和灾后重建工作提供参考价值。47.本发明基于北斗卫星导航系统观测数据和地表温度产品，提出了一种利用北斗卫星导航系统技术反演地震后震中及附近地区大气可降水量并进一步分析其与地表温度热异常关系的方法。48.本发明基于陆态网北斗卫星导航系统观测数据和mod11a1单日地表温度格网产品，针对在某地区发生的三次地震提出了一种利用北斗卫星导航系统技术反演震后震中及附近地区大气可降水量并进一步分析其与地表温度热异常关系的方法。其中，将北斗卫星导航系统技术反演的大气可降水量与rs及era5数据进行比对，验证了bds-pwv精度良好，可用来描述地震后水汽变化。本发明属于北斗卫星导航系统与遥感交叉的大气水汽研究范畴。49.随机选择三个发生过地震的区域，分别利用距离三个地震震中位置最近的北斗卫星导航系统的观测数据，按照反演流程获取bds-pwv。50.根据图3(a)-图3(c)，第一区域、第二区域和第三区域地震后的所有大气可降水量数据都出现了异常增加的情况。在第一区域地震爆发时，bds-pwv、rs-pwv和era5-pwv均处于最小值，在此之后便出现了pwv上升的变化，呈上升趋势，并在大地震发生后3天由于5月7日的余震从而达到峰值（bds-pwv约14mm，rs和era5约为11mm）。余震进一步促进能量释放，水汽不断碰撞和聚集，形成可降水。bds-pwv在第二区域地震发生三天后出现了急剧的上升，最高可达6mm，是整个地震过程中的最高值，同时，相比于rs-pwv和era5-pwv，bds-pwv的变化明显要明显很多。对于第三区域地震，三个数据源的pwv都在12月18日达到最大值（地震发生在12月13日23时，持续到第二天凌晨2时），其中bds和era5的pwv约为8.2mm，rs探测结果基本达到10mm，然后又急剧下降。51.上述分析表明，地震后存在大气可降水量异常，主要表现为先增加随后减少的趋势。地震导致岩石在构造应力作用下孔隙裂缝开闭，释放出大量地热能，如热气和热水。热效应的产生，加速震中上空水分的蒸发和水汽的运动，这也是地震过程中大气可降水量异常变化的原因。52.bds-pwv的坐标系统为大地高，探空资料的坐标系统为位势高，而era5再分析资料则为位势。高程统一到大地高后，进行精度验证。53.相关系数和误差如图4(a)-图4(f)所示。第一区域地震的bds-pwv与rs-pwv和era5-pwv非常接近，rmse分别为1.28mm和1.78mm。在第二区域地震中，bds-pwv与rs-pwv和era5-pwv的rmse分别为1.06mm和0.7mm。第三区域地震的rmse则为1.34mm和0.97mm。从以上分析可以看出，北斗卫星导航系统反演大气可降水量与其他两个数据源结果的rmse均小于2mm，其反演的准确性和可靠性非常高，能够为地震过程的气象研究提供有效大气可降水量参考。54.通过分析三个地震所对应的地表温度在地震过程中的变化，可以看出，第一区域地震爆发当天（st呈现大面积最大值为-6~-2，比地震前三天高约12~16。对于第二区域地震，震中位置lst有向西递增的趋势，稳定在-10~6，但在震后5天之内又开始逐渐下降，随后直到2月20日，lst恢复到震前水平。第三区域地震中，lst在震中及周边地区上升并持续了3天。从整个过程来看，在地震发生的当天或次日，lst都比前几天高。三次地震后的lst都比震前存在一个明显上升，并在高值处维持了三天左右，然后又下降到正常水平，其原因可能为地震过程中释放的能量在震后几天逐渐消散。55.其次，本发明又采用rst(robustsatellitetechnique，卫星鲁棒分析技术)算法，基于十年背景场数据来突出多年同地震时期、同地震位置的空间热异常。将alice指数作为热异常指标。56.第一区域地震过程的alice指数显示，某月6日和7日震中及周边地区出现严重的热异常，最大值为1.5，基于十年背景场出现了地表温度热异常；虽然8日至10日有云层覆盖缺失数据，但从11日开始，异常强度逐渐减弱。57.第二区域地震过程中的alice指数在某月13日至15日均为1以上，随后恢复正常（0±0.5），说明某年发生该地震时的lst比过去十年同范围区域的lst都高，存在热异常现象。58.第三区域地震后，震中及附近地区的alice指数上升至1，震后lst高于多年lst，基于10年数据的热异常相当明显。59.总之，地震后alice指数有所上升，地表温度热异常高于往年同期。此后几乎没有任何变化。60.以上分别从地震发生过程中lst前后对比和基于十年背景数据的多年同期同地lst对比两个方面提取了地震热异常，发现地震能够引起地表温度热异常。61.第一区域震中区最大lst在某月6日，同时pwv开始增加，次日达到第一个峰值，达到最大pwv峰值是在某月9日。62.对于第二区域地震，最高lst是在发震当天测得，bds-pwv在某月16日达到最大。63.第三区域地震于某月13日23:00发生，余震一直持续到次日凌晨2时，lst在震中及周边地区在15日最高，而在某月18日bds-pwv达到峰值，比lst最大值晚3天。64.综上所述，实验所选取的3次地震发生后lst和bds-pwv变化之间存在延迟，即大气可降水量变化紧随地表温度变化之后。bds-pwv变化与地表温度热异常存在一定的关系，大气可降水量峰值出现晚于地表温度升高3天。65.参考背景数据发现，震中地表温度热异常越严重，bds-pwv越高；否则，地表温度异常较弱或没有异常，则在震中没有bds-pwv的突出异常变化。66.地震诱发地表和地下热异常，引发更多的热能释放反应，促进水汽活跃，特别是水汽向上运动过程中的蒸发、输送和积聚，大量pwv沉积在大气中，超过一定程度可以产生降水。通过北斗卫星导航系统技术反演发现，震后pwv会增加，与地表温度热异常存在延迟关系。基于这一发现，本发明认为这对震后与降水相关的次生灾害的预测预防和震后重建具有重要意义。67.对待研究区域，就按照上述流程结合当前待研究区域的历史地震数据，得出bds-pwv峰值和热异常出现的结果，最终获得延迟时间；基于延迟时间，为以后当前待研究区域再次发生地震时通过提前检测到热异常，在加以历史延迟时间计算出pwv峰值时间点，预测次生灾害。68.本发明利用北斗卫星导航系统反演震后震中位置的pwv，并结合地表温度热异常进一步探求其与pwv关系。69.本发明北斗卫星导航系统反演获得震后震中及附近地区pwv并进一步分析其与地表温度热异常关系，为震后与降水相关的次生灾害的预测监测和震后重建工作提供参考。70.实施例二本实施例提供了一种地震后物理参数变化的预测系统；一种地震后物理参数变化的预测系统，包括：第一获取模块，其被配置为：获取历史时间段内，目标监测区域地震发生后的大气可降水量；精度验证模块，其被配置为：对历史时间段内采集的大气可降水量，进行精度验证；第二获取模块，其被配置为：获取历史时间段内，目标监测区域地震发生后的地表温度热异常数据；时间差确定模块，其被配置为：基于历史地表温度热异常数据，和通过精度验证的历史大气可降水量，确定大气可降水量峰值与地表温度热异常峰值之间的时间差；预测输出模块，其被配置为：当目标监测区域再次发生地震时，检测目标监测区域的地表温度热异常峰值对应的时间点；将时间点加上时间差得到大气可降水量峰值的预测时间点，将预测时间点作为次生灾害发生的时间点输出。71.此处需要说明的是，上述第一获取模块、精度验证模块、第二获取模块、时间差确定模块和预测输出模块对应于实施例一中的步骤s101至s105，上述模块与对应的步骤所实现的示例和应用场景相同，但不限于上述实施例一所公开的内容。需要说明的是，上述模块作为系统的一部分可以在诸如一组计算机可执行指令的计算机系统中执行。72.以上所述仅为本发明的优选实施例而已，并不用于限制本发明，对于本领域的技术人员来说，本发明可以有各种更改和变化。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。当前第1页12当前第1页12
再多了解一些
本文用于企业家、创业者技术爱好者查询，结果仅供参考。
