CN114089443B - 一种基于tec积分量及季节变化系数的uhf频段电离层闪烁事件预报方法 - Google Patents

一种基于tec积分量及季节变化系数的uhf频段电离层闪烁事件预报方法 Download PDF

Info

Publication number
CN114089443B
CN114089443B CN202111206590.XA CN202111206590A CN114089443B CN 114089443 B CN114089443 B CN 114089443B CN 202111206590 A CN202111206590 A CN 202111206590A CN 114089443 B CN114089443 B CN 114089443B
Authority
CN
China
Prior art keywords
itec
tec
frequency band
ionosphere
data
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202111206590.XA
Other languages
English (en)
Other versions
CN114089443A (zh
Inventor
张红波
王飞飞
盛冬生
刘玉梅
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Institute of Radio Wave Propagation CETC 22 Research Institute
Original Assignee
China Institute of Radio Wave Propagation CETC 22 Research Institute
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by China Institute of Radio Wave Propagation CETC 22 Research Institute filed Critical China Institute of Radio Wave Propagation CETC 22 Research Institute
Priority to CN202111206590.XA priority Critical patent/CN114089443B/zh
Publication of CN114089443A publication Critical patent/CN114089443A/zh
Application granted granted Critical
Publication of CN114089443B publication Critical patent/CN114089443B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01W—METEOROLOGY
    • G01W1/00—Meteorology
    • G01W1/10—Devices for predicting weather conditions
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/03—Cooperating elements; Interaction or communication between different cooperating elements or between cooperating elements and receivers
    • G01S19/07—Cooperating elements; Interaction or communication between different cooperating elements or between cooperating elements and receivers providing data for correcting measured positioning data, e.g. DGPS [differential GPS] or ionosphere corrections
    • G01S19/072—Ionosphere corrections
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Environmental & Geological Engineering (AREA)
  • Environmental Sciences (AREA)
  • Ecology (AREA)
  • Biodiversity & Conservation Biology (AREA)
  • Atmospheric Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics And Detection Of Objects (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明公开了一种基于TEC积分量及季节变化系数的UHF频段电离层闪烁事件预报方法，包括如下步骤：步骤1，计算海口地区十年间每天的季节变化系数以及提取十年间每天的夜间UHF频段电离层闪烁发生情况；步骤2，通过决策树算法形成一张电离层闪烁事件发生与否的1、0统计预报表；步骤3，应用时，在本地19:30时刻计算当天的季节变化系数和积分量；步骤4，每年年底增加当年数据后，重复步骤2的工作以更新统计预报表，用于未来一整年的海口地区UHF频段电离层闪烁预报工作。本发明所公开预报方法，能够为海口等低纬地区提供预报时效为0—6小时的夜间UHF频段电离层闪烁事件预报信息，有效解决电离层闪烁事件预报问题。

Description

一种基于TEC积分量及季节变化系数的UHF频段电离层闪烁事件预报方法
技术领域
本发明属于卫星通信和雷达领域，特别涉及该领域中的一种基于TEC积分量及季节变化系数的UHF频段电离层闪烁事件预报方法，服务于UHF或P频段天基SAR、远程目标监视雷达系统、卫星通信的低纬地区电离层闪烁影响预报。
背景技术
在电离层正常的电子密度背景上漂浮着密度不等、尺度不等的电离体，称为电离层不均匀体。当卫星信号在电离层不均匀体中传播时，会引起传播路径和传播时间的改变，使得信号的振幅、相位以及到达角发生快速起伏变化，导致信号衰落、频移等，被称为电离层闪烁。电离层闪烁效应会直接影响卫星通信、雷达、测控等电子信息系统的误码率、通信质量、定位精度、作用距离等性能指标，因此，通过监测、建模、告警与预报预警等空间天气保障工作，可减缓电离层闪烁对UHF、P频段天基SAR、远程目标监视雷达系统、卫星通信等电子信息系统的不利影响，具有重要应用价值。
在现有技术中，尚未发现预报时效为0—6小时且采用本专利所用技术路线的与UHF频段电离层闪烁事件预报相关的学术论文或者专利。类似技术中，美国军方研发了一个UHF频段卫星通信电离层闪烁异常短期预报模型，该模型通过全球中低纬地区部署的UHF频段电离层闪烁观测站的实时数据驱动，利用观测到的电离层闪烁发生区域，综合电离层闪烁漂移运动规律，可预报观测区域以东地区在未来3小时内能否出现电离层闪烁事件情况，该模式预报时效为0—3小时。
发明内容
本发明所要解决的技术问题就是提供一种UHF频段电离层闪烁事件预报方法，基于TEC积分量及季节变化系数与海口等低纬地区UHF频段电离层闪烁发生之间的相关性，解决了预报时效为0—6小时的海口等低纬地区夜间UHF频段电离层闪烁事件预报问题。
本发明采用如下技术方案：
一种基于TEC积分量及季节变化系数的UHF频段电离层闪烁事件预报方法，其改进之处在于，包括如下步骤：
步骤1，基于中国低纬地区电离层总电子含量TEC图历史数据计算十年间每天的东经120°南北驼峰附近日间TEC积分量ITECy，d；计算海口地区十年间每天的季节变化系数Cy，d以及提取十年间每天的夜间UHF频段电离层闪烁发生情况Sy，d；
步骤2，构建以Cy，d为横轴，ITECy，d为纵轴的二维坐标空间，将Cy，d，ITECy，d历史数据归到分辨率为0.05*2.5的相应网格空间，通过决策树算法形成一张电离层闪烁事件发生与否的1、0统计预报表TSC(i,j)，i＝0,1,......,20，j＝0,1,......,96；
步骤3，应用时，在本地19:30时刻计算当天的Cy，d和ITECy，d，查找该Cy，d，ITECy，d数据在统计预报表TSC(i,j)中对应的单元格数值，如果是1，则预报海口地区夜间会发生UHF频段电离层闪烁事件；
步骤4，每年年底增加当年全部Cy，d，ITECy，d，Sy，d数据后，重复步骤2的工作以更新统计预报表TSC(i,j)，用于未来一整年的海口地区UHF频段电离层闪烁预报工作。
进一步的，所述步骤1具体为：
A，基于中国低纬地区TEC图历史数据计算十年间每天的东经120°南北驼峰附近日间TEC积分量ITECy，d，单位为无量纲，具体计算方法如下：
假定当前年份为y0，从国际全球导航卫星系统GNSS服务网站下载前十年间[y0-1,y0-2,…,y0-10]时间分辨率为15分钟的全球TEC图历史数据，每天有96张全球TEC图，每张图的经度范围为东经-180°至180°，经度分辨率为5°；纬度范围为北纬87.5°至-87.5°，纬度分辨率为2.5°；时间范围为世界时UT0点至UT24点；
首先，从每一天的UT0:00点至UT11:15点，本地时为8:00点至19:15点之间46张全球TEC图中，提取东经120°线上北纬40°至-30°之间的网格点上的全部数据vTEClati，t，lati＝1，2，......，29，t＝1，......，46，纬度下标lati所对应的纬度为42.5-lati*2.5，单位为度；时间下标t所对应的世界时为t*0.25-0.25，单位为小时；
其次，搜索出十年间全部vTEClati，t数据中的最大值vTECmax，将全部vTEClati，t进行归一化处理，即：
normTEClati，t＝vTEClati，t/vTECmax
再次，针对每一天的normTEClati，t数据进行以下处理：搜索地磁场线0度以北的网格点normTEClati，t，lati＝1，2，......，13，t＝1，......，46，数据中最大值，及其所在的纬度下标latnmax；搜索地磁场线0度以南的网格点normTEClati，t，lati＝14，15，......，29，t＝1，......，46，数据中最大值，及其所在的纬度下标latsmax；
最后，计算每一天的南北驼峰附近日间TEC积分量ITECy，d，单位为无量纲，具体计算公式如下：
Figure BDA0003307161740000031
B，计算海口地区十年间每天的季节变化系数Cy，d，单位为无量纲，具体计算公式如下：
Figure BDA0003307161740000032
式中φangle为海口地区地磁子午面与本地电离层F2层日落线之间夹角，单位为度；Ws取值为13.0；
C，提取十年间每天的夜间UHF频段电离层闪烁发生情况Sy，d，具体计算方法为：
首先，计算海口地区本地日落时间tutsunset，tutsunset为海口地区本地电离层F2层日落的世界UT时间，即本地电离层F2层处太阳天顶角为107°时的时间，单位为小时，tutsunset默认值为11.50；
其次，基于海口观测站UHF频段电离层闪烁历史观测数据进行分析，如果在tutsunset至tutsunset+6.0时间段内发生了UHF频段电离层闪烁事件，则Sy，d＝1，否则Sy，d＝0；
D，最后将ITECy，d，Cy，d和Sy，d数据按天组织成一维矢量ITECy，d，Cy，d，Sy，d。
进一步的，所述步骤2具体为：
A，构建以Cy，d为横轴，ITECy，d为纵轴的二维坐标空间D(x,y)，Cy，d取值范围从0至1，ITECy，d取值范围从0至240，将该坐标空间以横轴分辨率0.05，纵轴分辨率2.5进行网格化处理D(xi,yj)，则横轴网格点坐标xi＝i*0.05，i＝0,1,......,20，纵轴网格点坐标yj＝j*2.5，j＝0,1,......,96；同时构建一张对应的二维表TSC(i,j)，i＝0,1,......,20，j＝0,1,......,96，初始值全部为0；
B，通过决策树算法形成电离层闪烁事件发生与否的1、0统计预报表TSC(i,j)：
首先，针对统计预报表TSC(i,j)中的每一个网格i,j，从十年间Cy，d，ITECy，d，Sy，d数据中筛选出落到网格空间D(xi,yj)的全部数据，筛选规则如下：
i＝＝floor(Cy,d/0.05)and j＝＝floor(ITECy,d/2.5)
式中floor(x)函数为取小于x的最大整数；
其次，采用决策树二元分类算法进行分析，如果网格空间D(xi,yj内Sy，d为1的Cy，d，ITECy，d数据百分比超过0.32，则TSC(i,j)赋值为1，其余情况则赋值为0，最终形成一张电离层闪烁事件发生与否的1、0统计预报表TSC(i,j)。
进一步的，所述步骤3具体为：
A，应用时，在UT11:30，本地19:30时刻计算当天的Cy，d和ITECy，d，计算其在统计预报表TSC(i,j)的单元格坐标i,j，计算公式如下：
i＝floor(Cy,d/0.05)
j＝floor(ITECy,d/2.5)
如果统计预报表单元格i,j中的数值是1，则预报海口地区夜间会发生UHF频段电离层闪烁事件。
进一步的，所述步骤4具体为：
在当前积累的ITECy，d，Cy，d和Sy，d历史数据基础上，年底增加y0年的全部Cy，d，ITECy，d，Sy，d数据后，重复步骤2的工作以更新统计预报表TSC(i,j)，用于未来一整年的海口地区UHF频段电离层闪烁预报工作。
本发明的有益效果是：
本发明所公开基于TEC积分量及季节变化系数的UHF频段电离层闪烁事件预报方法，能够为海口等低纬地区提供预报时效为0—6小时的夜间UHF频段电离层闪烁事件预报信息，有效解决电离层闪烁事件预报问题。基于海口地区2011年至2020年十年数据对本发明方法的测试验证结果为：UHF频段电离层闪烁事件预报准确度为90％，UHF频段电离层闪烁事件预报精度为65％，UHF频段电离层闪烁事件召回率为80％。因此，本发明方法是解决电离层闪烁事件短时预报问题的一种有效技术途径，有助于提升UHF或P频段天基SAR、远程目标监视雷达系统、卫星通信等系统应对电离层闪烁事件不利影响的能力，具有非常好的应用价值。
此外，本申请发明人张红波等已经申请的预报时效为0—30分钟的一种未来半小时内单站UHF频段电离层闪烁事件发生时长的预测方法(申请号：202010746590.8)，解决的是在已知正在发生的电离层闪烁事件后预测其后续发生时长的临近预报问题，而本申请解决的则是在电离层闪烁事件发生前预报夜间电离层闪烁是否能发生的短时预报问题。
附图说明
图1是本发明方法的流程示意图；
图2是海口地区2011年至2020年十年间的全部(Cy，d，ITECy，d)数据；
图3是将统计预报表TSC(i,j)中单元格为1所对应的(Cy，d，ITECy，d)数据进行绘制的结果。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图和实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。
实施例1，本实施例公开了一种基于TEC积分量及季节变化系数的UHF频段电离层闪烁事件预报方法，如图1所示，包括如下步骤：
步骤1，基于中国低纬地区电离层总电子含量TEC图历史数据计算十年间每天的东经120°南北驼峰附近日间TEC积分量ITECy，d；计算海口地区十年间每天的季节变化系数Cy，d以及提取十年间每天的夜间UHF频段电离层闪烁发生情况Sy，d，具体为：
A，基于中国低纬地区TEC图历史数据计算十年间每天的东经120°南北驼峰附近日间TEC积分量ITECy，d，单位为无量纲，具体计算方法如下：
假定当前年份为y0，从国际全球导航卫星系统GNSS服务网站下载前十年间[y0-1,y0-2,…,y0-10]时间分辨率为15分钟的全球TEC图历史数据，每天有96张全球TEC图，每张图的经度范围为东经-180°至180°，经度分辨率为5°；纬度范围为北纬87.5°至-87.5°，纬度分辨率为2.5°；时间范围为世界时UT0点至UT24点；
首先，从每一天的UT0:00点至UT11:15点(本地时为8:00点至19:15点)之间46张全球TEC图中，提取东经120°线上北纬40°至-30°之间的网格点上的全部数据vTEClati，t，lati＝1，2，......，29，t＝1，......，46，纬度下标lati所对应的纬度为42.5-lati*2.5，单位为度；时间下标t所对应的世界时为t*0.25-0.25，单位为小时；
其次，搜索出十年间全部vTEClati，t数据中的最大值vTECmax，将全部vTEClati，t进行归一化处理，即：
normTEClati，t＝vTEClati，t/vTECmax
再次，针对每一天的normTEClati，t数据进行以下处理：搜索地磁场线0度以北的网格点normTEClati，t(lati＝1，2，......，13，t＝1，......，46)数据中最大值，及其所在的纬度下标latnmax；搜索地磁场线0度以南的网格点normTEClati，t(lati＝14，15，......，29，t＝1，......，46)数据中最大值，及其所在的纬度下标latsmax；
最后，计算每一天的南北驼峰附近日间TEC积分量ITECy，d，单位为无量纲，具体计算公式如下：
Figure BDA0003307161740000061
B，计算海口地区十年间每天的季节变化系数Cy，d，单位为无量纲，具体计算公式如下：
Figure BDA0003307161740000062
式中φangle为海口地区地磁子午面与本地电离层F2层(350千米高度)日落线之间夹角，单位为度；Ws取值为13.0；
C，提取十年间每天的夜间UHF频段电离层闪烁发生情况Sy，d，具体计算方法为：
首先，计算海口地区本地日落时间tutsunset，tutsunset为海口地区本地电离层F2层(350千米高度)日落的世界UT时间，即本地电离层F2层处太阳天顶角为107°时的时间，单位为小时，tutsunset默认值为11.50；
其次，基于海口观测站UHF频段电离层闪烁历史观测数据进行分析，如果在tutsunset至tutsunset+6.0时间段内发生了UHF频段电离层闪烁事件，则Sy，d＝1，否则Sy，d＝0；
D，最后将ITECy，d，Cy，d和Sy，d数据按天组织成一维矢量(ITECy，d，Cy，d，Sy，d)。
步骤2，构建以Cy，d为横轴，ITECy，d为纵轴的二维坐标空间，将(Cy，d，ITECy，d)历史数据归到分辨率为0.05*2.5的相应网格空间，通过决策树算法形成一张电离层闪烁事件发生与否的1/0统计预报表TSC(i,j)，i＝0,1,......,20，j＝0,1,......,96，具体为：
A，构建以Cy，d为横轴，ITECy，d为纵轴的二维坐标空间D(x,y)，Cy，d取值范围从0至1，ITECy，d取值范围从0至240，将该坐标空间以横轴分辨率0.05，纵轴分辨率2.5进行网格化处理D(xi,yj)，则横轴网格点坐标xi＝i*0.05，i＝0,1,......,20，纵轴网格点坐标yj＝j*2.5，j＝0,1,......,96；同时构建一张对应的二维表TSC(i,j)，i＝0,1,......,20，j＝0,1,......,96，初始值全部为0；
B，通过决策树算法形成电离层闪烁事件发生与否的1/0统计预报表TSC(i,j)：
首先，针对统计预报表TSC(i,j)中的每一个网格(i,j)，从十年间(Cy，d，ITECy，d，Sy，d)数据中筛选出落到网格空间D(xi,yj)的全部数据，筛选规则如下：
i＝＝floor(Cy,d/0.05)andj＝＝floor(ITECy,d/2.5)
式中floor(x)函数为取小于x的最大整数操作；
其次，采用决策树二元分类算法进行分析，如果网格空间D(xi,yj内Sy，d为1的(Cy，d，ITECy，d)数据百分比超过0.32，则TSC(i,j)赋值为1，其余情况则赋值为0，最终形成一张电离层闪烁事件发生与否的1/0统计预报表TSC(i,j)。
步骤3，应用时，在本地19:30时刻计算当天的Cy，d和ITECy，d，查找该(Cy，d，ITECy，d)数据在统计预报表TSC(i,j)中对应的单元格数值，如果是1，则预报海口地区夜间会发生UHF频段电离层闪烁事件，具体为：
A，应用时，在UT11:30(本地19:30)时刻计算当天的Cy，d和ITECy，d，计算其在统计预报表TSC(i,j)的单元格坐标(i,j)，计算公式如下：
i＝floor(Cy,d/0.05)
j＝floor(ITECy,d/2.5)
如果统计预报表单元格(i,j)中的数值是1，则预报海口地区夜间会发生UHF频段电离层闪烁事件。
步骤4，每年年底增加当年全部(Cy，d，ITECy，d，Sy，d)数据后，重复步骤2的工作以更新统计预报表TSC(i,j)，其可应用于未来一整年的海口地区UHF频段电离层闪烁预报工作，具体为：
在当前积累的ITECy，d，Cy，d和Sy，d历史数据基础上，年底增加y0年的全部(Cy，d，ITECy，d，Sy，d)数据后，重复步骤2的工作以更新统计预报表TSC(i,j)，其可应用于未来一整年的海口地区UHF频段电离层闪烁预报工作。
图2为海口地区2011年至2020年十年间的全部(Cy，d，ITECy，d)数据，圆圈为夜间有UHF频段电离层闪烁事件的(Cy，d，ITECy，d)数据，x为夜间无UHF频段电离层闪烁事件的(Cy，d，ITECy，d)数据。
图3为将统计预报表TSC(i,j)中单元格为1所对应的(Cy，d，ITECy，d)数据进行绘制的结果。圆圈为夜间有UHF频段电离层闪烁事件的(Cy，d，ITECy，d)数据，占十年间全部夜间有UHF频段电离层闪烁事件的(Cy，d，ITECy，d)数据的80％，即召回率指标。
综上所述，本发明方法可有效解决预报时效为0—6小时的海口地区夜间UHF频段电离层闪烁事件预报问题。基于海口地区2011年至2020年的十年数据，对本发明方法的测试验证结果为UHF频段电离层闪烁事件预报准确度为90％，UHF频段电离层闪烁事件预报精度为65％，UHF频段电离层闪烁事件召回率为80％。因此，本发明方法是解决UHF/P频段天基SAR/远程目标监视雷达系统/卫星通信等系统电离层闪烁事件短时预报问题的一种有效技术途径，具有非常好的应用价值。

Claims (5)

1.一种基于TEC积分量及季节变化系数的UHF频段电离层闪烁事件预报方法，其特征在于，包括如下步骤：
步骤1，基于中国低纬地区电离层总电子含量TEC图历史数据计算十年间每天的东经120°南北驼峰附近日间TEC积分量ITECy，d；计算海口地区十年间每天的季节变化系数Cy，d以及提取十年间每天的夜间UHF频段电离层闪烁发生情况Sy，d；
步骤2，构建以Cy，d为横轴，ITECy，d为纵轴的二维坐标空间，将Cy，d，ITECy，d历史数据归到分辨率为0.05*2.5的相应网格空间，通过决策树算法形成一张电离层闪烁事件发生与否的1、0统计预报表TSC(i，j)，i＝0，1，......，20，j＝0，1，......，96；
步骤3，应用时，在本地19:30时刻计算当天的Cy，d和ITECy，d，查找该Cy，d，ITECy，d数据在统计预报表TSC(i，j)中对应的单元格数值，如果是1，则预报海口地区夜间会发生UHF频段电离层闪烁事件；
步骤4，每年年底增加当年全部Cy，d，ITECy，d，Sy，d数据后，重复步骤2的工作以更新统计预报表TSC(i，j)，用于未来一整年的海口地区UHF频段电离层闪烁预报工作。
2.根据权利要求1所述基于TEC积分量及季节变化系数的UHF频段电离层闪烁事件预报方法，其特征在于，所述步骤1具体为：
A，基于中国低纬地区TEC图历史数据计算十年间每天的东经120°南北驼峰附近日间TEC积分量ITECy，d，单位为无量纲，具体计算方法如下：
假定当前年份为y0，从国际全球导航卫星系统GNSS服务网站下载前十年间[y0-1，y0-2，......，y0-10]时间分辨率为15分钟的全球TEC图历史数据，每天有96张全球TEC图，每张图的经度范围为东经-180°至180°，经度分辨率为5°；纬度范围为北纬87.5°至-87.5°，纬度分辨率为2.5°；时间范围为世界时UT0点至UT24点；
首先，从每一天的UT0:00点至UT11:15点，本地时为8:00点至19:15点之间46张全球TEC图中，提取东经120°线上北纬40°至-30°之间的网格点上的全部数据VTEClati，t，lati＝1，2，......，29，t＝1，......，46，纬度下标lati所对应的纬度为42.5-lati*2.5，单位为度；时间下标t所对应的世界时为t*0.25-0.25，单位为小时；
其次，搜索出十年间全部vTEClati，t数据中的最大值vTECmax，将全部vTEClati，t进行归一化处理，即：
normTEClati，t＝vTEClati，t/vTECmax
再次，针对每一天的normTEClati，t数据进行以下处理：搜索地磁场线0度以北的网格点normTEClati，t，lati＝1，2，......，13，t＝1，......，46，数据中最大值，及其所在的纬度下标latnmax；搜索地磁场线0度以南的网格点normTEClati，t，lati＝14，15，......，29，t＝1，......，46，数据中最大值，及其所在的纬度下标latsmax；
最后，计算每一天的南北驼峰附近日间TEC积分量ITECy，d，单位为无量纲，具体计算公式如下：
Figure FDA0003307161730000021
B，计算海口地区十年间每天的季节变化系数Cy，d，单位为无量纲，具体计算公式如下：
Figure FDA0003307161730000022
式中φangle为海口地区地磁子午面与本地电离层F2层日落线之间夹角，单位为度；Ws取值为13.0；
C，提取十年间每天的夜间UHF频段电离层闪烁发生情况Sy，d，具体计算方法为：
首先，计算海口地区本地日落时间tutsunset，tutsunset为海口地区本地电离层F2层日落的世界UT时间，即本地电离层F2层处太阳天顶角为107°时的时间，单位为小时，tutsunset默认值为11.50；
其次，基于海口观测站UHF频段电离层闪烁历史观测数据进行分析，如果在tutsunset至tutsunset+6.0时间段内发生了UHF频段电离层闪烁事件，则Sy，d＝1，否则Sy，d＝0；
D，最后将ITECy，d，Cy，d和Sy，d数据按天组织成一维矢量ITECy，d，Cy，d，Sy，d。
3.根据权利要求1所述基于TEC积分量及季节变化系数的UHF频段电离层闪烁事件预报方法，其特征在于，所述步骤2具体为：
A，构建以Cy，d为横轴，ITECy，d为纵轴的二维坐标空间D(x，y)，Cy，d取值范围从0至1，ITECy，d取值范围从0至240，将该坐标空间以横轴分辨率0.05，纵轴分辨率2.5进行网格化处理D(xi，yj)，则横轴网格点坐标xi＝i*0.05，i＝0，1，......，20，纵轴网格点坐标yj＝j*2.5，j＝0，1，.....，96；同时构建一张对应的二维表TSC(i，j)，i＝0，1，......，20，j＝0，1，......，96，初始值全部为0；
B，通过决策树算法形成电离层闪烁事件发生与否的1、0统计预报表TSC(i，j)：
首先，针对统计预报表TSC(i，j)中的每一个网格i，j，从十年间Cy，d，ITECy，d，Sy，d数据中筛选出落到网格空间D(xi，yj)的全部数据，筛选规则如下：
i＝＝floor(Cy，d/0.05)and j＝floor(ITECy，d/2.5)
式中floor(x)函数为取小于x的最大整数；
其次，采用决策树二元分类算法进行分析，如果网格空间D(xi，yj内Sy，d为1的Cy，d，ITECy，d数据百分比超过0.32，则TSC(i，j)赋值为1，其余情况则赋值为0，最终形成一张电离层闪烁事件发生与否的1、0统计预报表TSC(i，j)。
4.根据权利要求1所述基于TEC积分量及季节变化系数的UHF频段电离层闪烁事件预报方法，其特征在于，所述步骤3具体为：
A，应用时，在UT11:30，本地19:30时刻计算当天的Cy，d和ITECy，d，计算其在统计预报表TSC(i，j)的单元格坐标i，j，计算公式如下：
i＝floor(Cy，d/0.05)
j＝floor(ITECy，d/2.5)
如果统计预报表单元格i，j中的数值是1，则预报海口地区夜间会发生UHF频段电离层闪烁事件。
5.根据权利要求2所述基于TEC积分量及季节变化系数的UHF频段电离层闪烁事件预报方法，其特征在于，所述步骤4具体为：
在当前积累的ITECy，d，Cy，d和Sy，d历史数据基础上，年底增加y0年的全部Cy，d，ITECy，d，Sy，d数据后，重复步骤2的工作以更新统计预报表TSC(i，j)，用于未来一整年的海口地区UHF频段电离层闪烁预报工作。
CN202111206590.XA 2021-10-17 2021-10-17 一种基于tec积分量及季节变化系数的uhf频段电离层闪烁事件预报方法 Active CN114089443B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111206590.XA CN114089443B (zh) 2021-10-17 2021-10-17 一种基于tec积分量及季节变化系数的uhf频段电离层闪烁事件预报方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111206590.XA CN114089443B (zh) 2021-10-17 2021-10-17 一种基于tec积分量及季节变化系数的uhf频段电离层闪烁事件预报方法

Publications (2)

Publication Number Publication Date
CN114089443A CN114089443A (zh) 2022-02-25
CN114089443B true CN114089443B (zh) 2023-04-25

Family

ID=80297105

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111206590.XA Active CN114089443B (zh) 2021-10-17 2021-10-17 一种基于tec积分量及季节变化系数的uhf频段电离层闪烁事件预报方法

Country Status (1)

Country Link
CN (1) CN114089443B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114897010A (zh) * 2022-04-29 2022-08-12 东南大学 基于重采样改进的梯度提升树的强电离层闪烁检测方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2015017824A2 (en) * 2013-08-02 2015-02-05 Atmospheric & Space Technology Research Associates Llc Real-time autonomous weather and space weather monitoring
CN105116469A (zh) * 2015-08-05 2015-12-02 中国电子科技集团公司第二十二研究所 一种中国低纬地区电离层闪烁发生概率预报方法
WO2016034252A1 (en) * 2014-09-05 2016-03-10 Fugro Intersite B.V. Ionospheric scintillation prediction
CN111830596A (zh) * 2020-07-30 2020-10-27 中国电波传播研究所（中国电子科技集团公司第二十二研究所） 一种未来半小时内单站uhf频段电离层闪烁事件发生时长的预测方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2015017824A2 (en) * 2013-08-02 2015-02-05 Atmospheric & Space Technology Research Associates Llc Real-time autonomous weather and space weather monitoring
WO2016034252A1 (en) * 2014-09-05 2016-03-10 Fugro Intersite B.V. Ionospheric scintillation prediction
CN105116469A (zh) * 2015-08-05 2015-12-02 中国电子科技集团公司第二十二研究所 一种中国低纬地区电离层闪烁发生概率预报方法
CN111830596A (zh) * 2020-07-30 2020-10-27 中国电波传播研究所（中国电子科技集团公司第二十二研究所） 一种未来半小时内单站uhf频段电离层闪烁事件发生时长的预测方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
张红波 ; 王飞飞 ; 徐良 ; 盛冬生 ; 刘玉梅 ; .一种考虑春秋分不对称性的电离层闪烁发生概率季节变化模型.电波科学学报.2019,(第02期),全文. *

Also Published As

Publication number Publication date
CN114089443A (zh) 2022-02-25

Similar Documents

Publication Publication Date Title
Alvarez et al. Comparison of different wind products and buoy wind data with seasonality and interannual climate variability in the southern Bay of Biscay (2000–2009)
CN103698750A (zh) 海洋二号卫星散射计海面风场反演方法和装置
CN105116469B (zh) 一种中国低纬地区电离层闪烁发生概率预报方法
CN112649899B (zh) 一种全球电离层数据同化和预报方法
CN111401602A (zh) 基于神经网络的卫星以及地面降水测量值同化方法
CN112946784B (zh) 一种基于深度学习的SuperDARN雷达对流图短期预报方法
CN104392113B (zh) 一种近海海面冷空气大风风速的估算方法
CN114089443B (zh) 一种基于tec积分量及季节变化系数的uhf频段电离层闪烁事件预报方法
Desai et al. Impacts of intense geomagnetic storms on NavIC/IRNSS system
CN111932017A (zh) 一种适合于单站电离层tec的短期预报方法
CN111830596B (zh) 一种未来半小时内单站uhf频段电离层闪烁事件发生时长的预测方法
Wu et al. Sea surface wind speed retrieval based on empirical orthogonal function analysis using 2019–2020 CYGNSS data
Wang et al. A rapid identification and warning method for severe weather via Doppler radar based on an improved TITAN algorithm
CN116029162B (zh) 利用星载gnss-r数据的洪涝灾害淹没范围监测方法和系统
Fang et al. Multi-scheme corrected dynamic–analogue prediction of summer precipitation in northeastern China based on BCC_CSM
Guha et al. Characteristics of severe thunderstorms studied with the aid of VLF atmospherics over North–East India
CN114355411B (zh) 一种基于北斗或gps载噪比观测值的洪水探测方法
Priyadarshi et al. Machine Learning-based ionospheric modelling performance during high ionospheric activity
CN114384556B (zh) 一种区域高分辨率电离层tec地图重构方法
Curtis et al. An analysis of the lightning jump algorithm using geostationary lightning mapper flashes
CN110727719B (zh) 一种基于动力松弛逼近的闪电定位资料同化方法
CN111008361A (zh) 一种电离层参数重构方法
Wang et al. Coastal Sea Surface Temperature Inversion from Microwave Radiometer using Radial Basis Function Neural Network
CN111323807B (zh) 一种夏季区域Es变化趋势的预测方法
Frishfelds et al. Aggregating sea surface hydrodynamic forecasts from multi-models for European seas

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant