申请/专利权人：中国空间技术研究院

申请日：2020-04-30

公开（公告）日：2024-04-12

公开（公告）号：CN111797491B

主分类号：G06F30/20

分类号：G06F30/20;G06Q10/04

优先权：

专利状态码：有效-授权

法律状态：2024.04.12#授权;2020.11.06#实质审查的生效;2020.10.20#公开

摘要：本发明公开了一种华北平原地壳垂向位移季节性和时空变化分析方法和系统，该方法包括：获取地壳形变序列、GRACE卫星同期RL06球谐系数产品和降雨数据产品；解算得到地壳负荷形变；分解得到两个模态分量；确定两个模态分量各自的季节项分量和趋势项；根据两个模态分量各自的季节项分量和趋势项进行重构，确定季节项分和趋势项的重构结果；根据重构结果，确定华北平原地区地壳形变的时空变化关系，以实现对当地水资源的有效管理。本发明分析了华北平原地壳垂向位移与降雨量的时空变化关系，对当地水资源的有效管理以及居民生活具有重要意义。

主权项：1.一种华北平原地壳垂向位移季节性和时空变化分析方法，其特征在于，包括：获取地壳形变序列yGPS、GRACE卫星同期RL06球谐系数产品和降雨数据产品；根据负荷弹性形变理论、滑动窗拟合多项式法、高斯滤波及GAC，对GRACE卫星同期RL06球谐系数产品进行改正，得到地壳负荷形变采用变分模态分解VMD方法分解地壳形变序列yGPS和地壳负荷形变得到yGPS和各自的模态分量IMFGPS和IMFGRACE；根据能量谱值和L-S谱分析法，确定IMFGPS的季节项分量SEAGPS和趋势项TRGPS、以及IMFGRACE的季节项分量SEAGRACE和趋势项TRGRACE；分别对IMFGPS的季节项分量SEAGPS和趋势项TRGPS、以及IMFGRACE的季节项分量SEAGRACE和趋势项TRGRACE进行重构，确定季节项分和趋势项的重构结果；根据重构结果，利用加权均方根法进行分析，确定SEAGPS与SEAGRACE的季节项关系，并结合TRGPS、TRGRACE和降雨数据产品确定华北平原地区地壳形变的时空变化关系；获取地壳形变序列yGPS，包括：获取中国环境监测网络CMONOC提供的CORS站坐标时间序列，将CORS站坐标时间序列作为GPS卫星观测数据Ut；剔除GPS卫星观测数据Ut中大于三倍标准差的异常值，得到地壳形变序列yGPS；根据负荷弹性形变理论、滑动窗拟合多项式法、高斯滤波及GAC，对GRACE卫星同期RL06球谐系数产品进行改正，得到地壳负荷形变包括：采用滑动窗拟合多项式法、高斯滤波及GAC，对GRACE卫星同期RL06球谐系数产品进行改正，得到地球重力场球谐系数C的变化量ΔClm和地球重力场球谐系数S的变化量ΔSlm；根据负荷弹性形变理论，确定负荷弹性形变公式： 通过负荷弹性形变公式，解算得到地壳负荷形变其中，θ和分别表示余纬度和经度，A表示地球半径，hl表示l阶引力潮作用下的径向勒夫参数，kl表示l阶引力潮作用下的水平勒夫参数；Wl表示高斯平滑核函数，表示完全规格化缔合勒让德函数，m表示固定球谐系数的阶次，l表示勒夫参数的阶次； 其中，r表示高斯滤波平滑半径，a表示地球的平均半径，γ表示球面上任意两点与之间角度；采用滑动窗拟合多项式法、高斯滤波及GAC，对GRACE卫星同期RL06球谐系数产品进行改正，得到地球重力场球谐系数C的变化量ΔClm和地球重力场球谐系数S的变化量ΔSlm，包括：确定GRACE卫星同期RL06球谐系数产品中的C20项，将GRACE卫星同期RL06球谐系数产品中的C20项采用CSR机构提供的激光数据文件进行替换，并读取得到地球重力场球谐系数C的原始变化量ΔC和地球重力场球谐系数S的原始变化量ΔS；采用固定球谐系数的阶次m、大小为w的滑动窗口，对各阶球谐系数Cl-2i,m，…，Cl-2,m，Cl,m，Cl+2,m，…，Cl+2i,m；Sl-2i,m，…，Sl-2,m，Sl,m，Sl+2,m，…，Sl+2i,m进行二阶多项式拟合，得到对应阶次的第一拟合值C和第二拟合值S，将第一拟合值C和第二拟合值S作为改正值；其中，下标p，m表示阶数，p＝l-2i、…l-2、l、l+2、…、l+2i；获取GAC球谐系数C的变化量ΔCGAC和球谐系数S的变化量ΔSGAC，并依次通过如下公式解算得到地球重力场球谐系数C的变化量ΔClm和地球重力场球谐系数S的变化量ΔSlm：ΔC0＝ΔC-C，ΔS0＝ΔS-SΔC′＝ΔC0+ΔCGAC，ΔS′＝ΔS0+ΔSGACΔClm＝ΔC′+ΔC1，ΔSlm＝ΔS′+ΔS1其中，ΔC0表示去相关后的地球重力场球谐系数C的变化量，ΔS0表示去相关后的地球重力场球谐系数S的变化量，ΔC1表示地球重力场球谐系数C的一阶改正，ΔS1表示地球重力场球谐系数S的一阶改正；采用变分模态分解VMD方法分解地壳形变序列yGPS和地壳负荷形变得到yGPS和各自的模态分量IMFGPS和IMFGRACE，包括：构造变分方程： 或其中，yGPS和均可作为待分解的序列，yk表示模态函数，wk表示IMF1-～IMFk两侧对称频率，k表示设定的模态分量个数，δt表示均段脉冲函数，表示泛函数对时间t的一阶偏导，j为虚数单位；确定拉格朗日增广表达式： 其中，α表示惩罚因子，λt表示拉格朗日乘子；利用交替方向乘子法进行求解，得到变分解 其中，和表示yt、ykt和λt的傅里叶变换结果；将滤波后的信号进行傅里叶逆变换处理，当时，得到yGPS时域下的模态分量IMFGPS，当时，得到时域下的模态分量IMFGRACE；其中，IMFGPS＝IMF1+IMF2+…+IMFn或IMFGRACE＝IMF1+IMF2+…+IMFn；根据能量谱值和L-S谱分析法，确定IMFGPS的季节项分量SEAGPS和趋势项TRGPS、以及IMFGRACE的季节项分量SEAGRACE和趋势项TRGRACE，包括：求解各IMF分量IMF1、IMF2、…、IMFn的能量谱指数值，筛选得到第一个能量谱指数的极小值点，将所述第一个能量谱指数的极小值点之前的IMF分量确定为残差序列；根据确定的残差序列，通过L-S谱分析法，判断IMF分量中季节项部分并对季节项部分进行重构，将IMF分量中除季节项部分之外的其余部作为即趋势项部分；其中，当IMFGPS＝IMF1+IMF2+…+IMFn时，得到IMFGPS的季节项分量SEAGPS和趋势项TRGPS；当IMFGRACE＝IMF1+IMF2+…+IMFn时，得到IMFGRACE的季节项分量SEAGRACE和趋势项TRGRACE。

全文数据：

权利要求：

百度查询： 中国空间技术研究院 华北平原地壳垂向位移季节性和时空变化分析方法和系统

免责声明
1、本报告根据公开、合法渠道获得相关数据和信息，力求客观、公正，但并不保证数据的最终完整性和准确性。
2、报告中的分析和结论仅反映本公司于发布本报告当日的职业理解，仅供参考使用，不能作为本公司承担任何法律责任的依据或者凭证。

阅读全文 双屏查看 官方信息 专利公告 收藏专利 下载PDF 下载WORD