CN109146175B - 一种利用卡尔曼滤波的输电线路覆冰实时估计与预测方法 - Google Patents
一种利用卡尔曼滤波的输电线路覆冰实时估计与预测方法 Download PDFInfo
- Publication number
- CN109146175B CN109146175B CN201810955377.0A CN201810955377A CN109146175B CN 109146175 B CN109146175 B CN 109146175B CN 201810955377 A CN201810955377 A CN 201810955377A CN 109146175 B CN109146175 B CN 109146175B
- Authority
- CN
- China
- Prior art keywords
- moment
- value
- equivalent
- icing
- icing thickness
- 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
Links
- 238000001914 filtration Methods 0.000 title claims abstract description 28
- 230000005540 biological transmission Effects 0.000 title claims abstract description 26
- 238000000034 method Methods 0.000 title claims abstract description 26
- 239000011248 coating agent Substances 0.000 claims description 14
- 238000000576 coating method Methods 0.000 claims description 14
- 238000001556 precipitation Methods 0.000 claims description 11
- 238000011002 quantification Methods 0.000 claims description 7
- 238000005259 measurement Methods 0.000 claims description 4
- 238000013139 quantization Methods 0.000 claims description 4
- 238000005070 sampling Methods 0.000 claims description 3
- 238000012937 correction Methods 0.000 abstract description 2
- 238000011161 development Methods 0.000 description 2
- 230000005611 electricity Effects 0.000 description 2
- 238000012806 monitoring device Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000002159 abnormal effect Effects 0.000 description 1
- 230000004075 alteration Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000000611 regression analysis Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000010792 warming Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Business, Economics & Management (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- Human Resources & Organizations (AREA)
- Data Mining & Analysis (AREA)
- Marketing (AREA)
- Mathematical Analysis (AREA)
- Operations Research (AREA)
- Health & Medical Sciences (AREA)
- Pure & Applied Mathematics (AREA)
- General Business, Economics & Management (AREA)
- Mathematical Physics (AREA)
- Tourism & Hospitality (AREA)
- Mathematical Optimization (AREA)
- Computational Mathematics (AREA)
- Software Systems (AREA)
- Quality & Reliability (AREA)
- Development Economics (AREA)
- General Engineering & Computer Science (AREA)
- Game Theory and Decision Science (AREA)
- Bioinformatics & Computational Biology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Entrepreneurship & Innovation (AREA)
- Life Sciences & Earth Sciences (AREA)
- Evolutionary Biology (AREA)
- Databases & Information Systems (AREA)
- Algebra (AREA)
- Probability & Statistics with Applications (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Wind Motors (AREA)
Abstract
本发明涉及一种利用卡尔曼滤波的输电线路覆冰实时估计与预测方法,本发明利用K时刻的覆冰与微气象历史数据,建立基于等值覆冰厚度增量与微气象因子的多元非线性回归模型,结合预报的微气象数据,预测K+1时刻的等值覆冰厚度增量。基于卡尔曼滤波,结合K时刻的等值覆冰厚度最优估计值与K+1时刻的等值覆冰厚度增量预测值建立卡尔曼滤波预测方程,计算K+1时刻的等值覆冰厚度预测值,当获取K+1时刻的等值覆冰厚度观测值后建立卡尔曼滤波观测方程,进而获取K+1时刻的等值覆冰厚度最优估计值。本发明可消除覆冰观测序列中的噪声与粗差,通过覆冰模型预测值与覆冰观测值的相互校正,实现覆冰观测序列的实时估计与优化预测。
Description
技术领域
本发明属于架空输电线路在线监测技术领域,尤指一种利用卡尔曼滤波的输电线路覆冰实时估计与预测方法。
背景技术
国民经济增长对电力具有高度的依赖性,同时,经济增长也会促进用电量需求的增大。随着电网大规模的建设与发展,远距离、大容量的输电线路不可避免的要经过高寒、高海拔、高湿度、强降水地区或盆地、分水岭等特殊微地形造成的微气象地区,极易引发输电线路覆冰灾害,且全球气候变暖导致的极端异常气候频发,更加重了输电线路覆冰灾害,人们对输电线路覆冰问题的关注也越来越多。为了保证电力系统的安全运行,预防与减小覆冰灾害对输电线路造成的影响,研究输电线路覆冰实时估计与预测方法,并让管理人员实时了解线路覆冰现场状况与覆冰发展趋势显得尤其重要。
目前,国内外已开展了输电线路覆冰实时估计与预测方法的相关研究,由于覆冰预测模型易受覆冰观测值的影响,而覆冰观测值易受观测环境与观测噪声的干扰,使覆冰的实时估计与预测精度低,主要体现在以下几个方面:
(1)远距离、大容量的输电线路不可避免的要穿过山地、高原、盆地等特殊复杂地貌造成的微地形、微气象区域与高湿、高寒、高海拔、强降雨降雪地区,因此,覆冰在线监测装置的观测环境极其恶劣,进而影响其观测精度。
(2)覆冰实时估计值是线路管理人员实时了解覆冰现场状况的依据,而目前此方面的研究是直接将覆冰观测值作为覆冰实时估计值来评价当前时刻的覆冰状态,受观测噪声等因素的影响,覆冰观测值并不能准确的表达当前的覆冰状态。
(3)覆冰预测方法通常是建立覆冰预测模型,然后通过观测的覆冰与微气象历史数据的最小二乘回归拟合得到模型系数,进而计算覆冰预测值,而覆冰观测值不可避免的存在误差,必将影响覆冰预测精度。
发明内容
为了解决上述技术问题,本发明提供了一种利用卡尔曼滤波的输电线路覆冰实时估计与预测方法。
本发明所采用的技术方案是:
一种利用卡尔曼滤波的输电线路覆冰实时估计与预测方法,其特征在于,包括以下步骤:
步骤1:利用K时刻的覆冰与微气象历史数据,建立基于等值覆冰厚度增量与微气象因子的多元非线性回归模型,基于
ΔR=a2T2+b2Vn 2+c2P2+a1T+b1Vn+c1P+d1Hr+e0
其中,ΔR为等值覆冰厚度增量,T为温度,P为降水,Vn为有效风速,Hr为相对湿度;ai、bi、ci、di与e0均为模型系数,可通过覆冰与微气象历史数据的最小二乘回归求得;
步骤2:结合多元非线性回归模型与预报的微气象数据,预测K+1时刻的等值覆冰厚度增量,具体是将量化的预报微气象数据代入步骤1中的多元非线性回归模型,计算K+1时刻的等值覆冰厚度增量;
步骤3:基于卡尔曼滤波,计算K+1时刻的等值覆冰厚度预测值,并在获取K+1时刻的等值覆冰厚度观测值后,计算K+1时刻的等值覆冰厚度最优估计值。
在上述的一种利用卡尔曼滤波的输电线路覆冰实时估计与预测方法,步骤1中,将K时刻的覆冰与微气象历史数据进行量化,具体是:温度量化为单位为℃的温度数值;湿度量化为单位为%的相对湿度数值;降水量化为单位为mm/h的降水量数值;风速与风向量化为有效风速Vn=Vsinθ,单位为m/s,其中V为风速,θ为风向与导线方向之间的夹角;等值覆冰厚度量化为等值覆冰厚度增量ΔRK=RK-RK-1,单位为mm,其中RK为K时刻的等值覆冰厚度。
在上述的一种利用卡尔曼滤波的输电线路覆冰实时估计与预测方法,步骤2中所述预测K+1时刻的等值覆冰厚度增量,首先将K+1时刻预报的微气象数据进行量化,具体是:温度量化为单位为℃的温度数值;湿度量化为单位为%的相对湿度数值;降水量化为单位为mm/h的降水量数值;风速与风向量化为有效风速Vn=Vsinθ,单位为m/s,其中V为风速,θ为风向与导线方向之间的夹角。
在上述的一种利用卡尔曼滤波的输电线路覆冰实时估计与预测方法,步骤3中所述利用卡尔曼滤波进行实时估计与预测,其具体实现包括以下子步骤:
步骤3.1:结合K时刻的等值覆冰厚度最优估计值与K+1时刻的等值覆冰厚度增量预测值建立卡尔曼滤波预测方程(如下式所示),计算K+1时刻的等值覆冰厚度预测值:
其中,Δt=tk-tk-1为覆冰观测序列的采样间隔,Xk|k为K时刻的等值覆冰厚度最优估计值,Δxk|k为K+1时刻的等值覆冰厚度增量预测值,Wk为K时刻的动态噪声,Xk+1|k为K+1时刻的等值覆冰厚度预测值;
步骤3.2:当获取K+1时刻的等值覆冰厚度观测值后基于下式建立卡尔曼滤波观测方程,进而计算K+1时刻的等值覆冰厚度最优估计值:
其中,Vk+1为K+1时刻的观测噪声,Zk+1为K+1时刻的等值覆冰厚度观测值,Xk+1|k+1为K+1时刻的等值覆冰厚度最优估计值。
本发明可消除覆冰观测序列中的噪声与粗差,通过等值覆冰厚度的模型预测值与观测值相互校正,将覆冰观测序列进行优化,得到更加精确的等值覆冰厚度预测值,为输电线路覆冰实时估计与预测提供了新方法。
附图说明
图1是本发明实施例的流程图。
图2是本发明实施例的等值覆冰厚度增量预测流程图。
具体实施方式
为了便于本领域普通技术人员理解和实施本发明,下面结合附图及实施例对本发明作进一步的详细描述,应当理解,此处所描述的实施示例仅用于说明和解释本发明,并不用于限定本发明。
请见图1和图2,本发明提供的一种利用卡尔曼滤波的输电线路覆冰实时估计与预测方法,包括以下步骤:
步骤1:将K时刻的覆冰与微气象历史数据进行量化,利用量化后的覆冰与微气象历史数据,建立基于等值覆冰厚度增量与微气象因子的多元非线性回归模型,基于
ΔR=a2T2+b2Vn 2+c2P2+a1T+b1Vn+c1P+d1Hr+e0
其中,ΔR为等值覆冰厚度增量,T为温度,P为降水,Vn为有效风速,Hr为相对湿度;ai、bi、ci、di与e0均为模型系数,可通过覆冰与微气象历史数据的最小二乘回归求得;
步骤2:将K+1时刻预报的微气象数据进行量化,代入多元非线性回归模型,预测K+1时刻的等值覆冰厚度增量;
步骤3:根据卡尔曼滤波的状态预报方程,结合K时刻的等值覆冰厚度最优估计值与K+1时刻的等值覆冰厚度增量预测值,计算K+1时刻的等值覆冰厚度预测值,基于
其中,Δt=tk-tk-1为覆冰观测序列的采样间隔,Xk|k为K时刻的等值覆冰厚度最优估计值,Δxk|k为K+1时刻的等值覆冰厚度增量预测值,Wk为K时刻的动态噪声,Xk+1|k为K+1时刻的等值覆冰厚度预测值;
步骤4:在获取K+1时刻的等值覆冰厚度观测值后,结合K+1时刻的等值覆冰厚度预测值,建立卡尔曼滤波观测方程,进而计算K+1时刻的等值覆冰厚度最优估计值,基于
其中,Vk+1为K+1时刻的观测噪声,Zk+1为K+1时刻的等值覆冰厚度观测值,Xk+1|k+1为K+1时刻的等值覆冰厚度最优估计值。
本实施例的温度的量化,是将温度统一成单位为℃的温度数值。
本实施例的湿度的量化,是将湿度统一成单位为%的相对湿度数值。
本实施例的风速与风向的量化,是把风向的影响作用在风速上,将风速与风向转化为有效风速,单位依然为m/s。
本实施例的等值覆冰厚度的量化,是将相邻时刻的等值覆冰厚度进行差分,计算等值覆冰厚度增量,单位为mm。
本实施例的预报微气象数据的量化方法与微气象历史数据的量化方法相同。
本实施例的多元非线性回归模型系数,是基于量化后的覆冰与微气象历史数据,通过最小二乘回归分析计算得到的。
本实施例的卡尔曼滤波过程中,K时刻的等值覆冰厚度最优估计值与K+1时刻的等值覆冰厚度增量预测值组成了卡尔曼滤波的状态预报方程,K+1时刻的等值覆冰厚度观测值组成卡尔曼滤波观测方程,通过滤波优化处理,得到K+1时刻的等值覆冰厚度预测值与最优估计值。
本发明将卡尔曼滤波引入输电线路覆冰实时估计与预测过程,消除了覆冰观测序列中的噪声与粗差,通过等值覆冰厚度的模型预测值与观测值相互校正,将覆冰观测序列进行优化,得到更加精确的等值覆冰厚度预测值,为输电线路覆冰实时估计与预测提供了新方法。
应当理解的是,本说明书未详细阐述的部分均属于现有技术。
应当理解的是,上述针对较佳实施例的描述较为详细,并不能因此而认为是对本发明专利保护范围的限制,本领域的普通技术人员在本发明的启示下,在不脱离本发明权利要求所保护的范围情况下,还可以做出替换或变形,均落入本发明的保护范围之内,本发明的请求保护范围应以所附权利要求为准。
Claims (3)
1.一种利用卡尔曼滤波的输电线路覆冰实时估计与预测方法,其特征在于,包括以下步骤:
步骤1:利用K时刻的覆冰与微气象历史数据,建立基于等值覆冰厚度增量与微气象因子的多元非线性回归模型,基于
ΔR=a2T2+b2Vn 2+c2P2+a1T+b1Vn+c1P+d1Hr+e0
其中,ΔR为等值覆冰厚度增量,T为温度,P为降水,Vn为有效风速,Hr为相对湿度;ai、bi、ci、di与e0均为模型系数,可通过覆冰与微气象历史数据的最小二乘回归求得;
步骤2:结合多元非线性回归模型与预报的微气象数据,预测K+1时刻的等值覆冰厚度增量,具体是将量化的预报微气象数据代入步骤1中的多元非线性回归模型,计算K+1时刻的等值覆冰厚度增量;
步骤3:基于卡尔曼滤波,计算K+1时刻的等值覆冰厚度预测值,并在获取K+1时刻的等值覆冰厚度观测值后,计算K+1时刻的等值覆冰厚度最优估计值;
具体实现包括以下子步骤:
步骤3.1:结合K时刻的等值覆冰厚度最优估计值与K+1时刻的等值覆冰厚度增量预测值建立卡尔曼滤波预测方程(如下式所示),计算K+1时刻的等值覆冰厚度预测值:
其中,Δt=tk-tk-1为覆冰观测序列的采样间隔,Xk|k为K时刻的等值覆冰厚度最优估计值,Δxk|k为K+1时刻的等值覆冰厚度增量预测值,Wk为K时刻的动态噪声,Xk+1|k为K+1时刻的等值覆冰厚度预测值;
步骤3.2:当获取K+1时刻的等值覆冰厚度观测值后基于下式建立卡尔曼滤波观测方程,进而计算K+1时刻的等值覆冰厚度最优估计值:
其中,Vk+1为K+1时刻的观测噪声,Zk+1为K+1时刻的等值覆冰厚度观测值,Xk+1|k+1为K+1时刻的等值覆冰厚度最优估计值。
2.根据权利要求1所述的一种利用卡尔曼滤波的输电线路覆冰实时估计与预测方法,其特征在于,步骤1中,将K时刻的覆冰与微气象历史数据进行量化,具体是:温度量化为单位为℃的温度数值;湿度量化为单位为%的相对湿度数值;降水量化为单位为mm/h的降水量数值;风速与风向量化为有效风速Vn=Vsinθ,单位为m/s,其中V为风速,θ为风向与导线方向之间的夹角;等值覆冰厚度量化为等值覆冰厚度增量ΔRK=RK-RK-1,单位为mm,其中RK为K时刻的等值覆冰厚度。
3.根据权利要求1所述的一种利用卡尔曼滤波的输电线路覆冰实时估计与预测方法,其特征在于,步骤2中所述预测K+1时刻的等值覆冰厚度增量,首先将K+1时刻预报的微气象数据进行量化,具体是:温度量化为单位为℃的温度数值;湿度量化为单位为%的相对湿度数值;降水量化为单位为mm/h的降水量数值;风速与风向量化为有效风速Vn=Vsinθ,单位为m/s,其中V为风速,θ为风向与导线方向之间的夹角。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810955377.0A CN109146175B (zh) | 2018-08-21 | 2018-08-21 | 一种利用卡尔曼滤波的输电线路覆冰实时估计与预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810955377.0A CN109146175B (zh) | 2018-08-21 | 2018-08-21 | 一种利用卡尔曼滤波的输电线路覆冰实时估计与预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109146175A CN109146175A (zh) | 2019-01-04 |
CN109146175B true CN109146175B (zh) | 2022-04-01 |
Family
ID=64790793
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810955377.0A Active CN109146175B (zh) | 2018-08-21 | 2018-08-21 | 一种利用卡尔曼滤波的输电线路覆冰实时估计与预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109146175B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112949920B (zh) * | 2021-02-26 | 2023-04-07 | 中国电力工程顾问集团西南电力设计院有限公司 | 一种基于观冰代表站数据的区域覆冰预测及预警方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4812995A (en) * | 1987-05-19 | 1989-03-14 | Girgis Adly A | Adaptive Kalman Filtering in fault classification |
CN102297674A (zh) * | 2011-04-27 | 2011-12-28 | 中国电力工程顾问集团西南电力设计院 | 一种利用能见度模型预测输电线路冰厚的方法 |
CN102789447A (zh) * | 2012-07-09 | 2012-11-21 | 贵州电网公司输电运行检修分公司 | 基于灰多元线性回归的覆冰与气象关系的分析方法 |
CN105160599A (zh) * | 2015-09-17 | 2015-12-16 | 上海电力学院 | 基于滤波和支持向量机的电力线路覆冰评估方法 |
CN105654189A (zh) * | 2015-12-23 | 2016-06-08 | 西安工程大学 | 基于时间序列分析与卡尔曼滤波算法的覆冰短期预测方法 |
-
2018
- 2018-08-21 CN CN201810955377.0A patent/CN109146175B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4812995A (en) * | 1987-05-19 | 1989-03-14 | Girgis Adly A | Adaptive Kalman Filtering in fault classification |
CN102297674A (zh) * | 2011-04-27 | 2011-12-28 | 中国电力工程顾问集团西南电力设计院 | 一种利用能见度模型预测输电线路冰厚的方法 |
CN102789447A (zh) * | 2012-07-09 | 2012-11-21 | 贵州电网公司输电运行检修分公司 | 基于灰多元线性回归的覆冰与气象关系的分析方法 |
CN105160599A (zh) * | 2015-09-17 | 2015-12-16 | 上海电力学院 | 基于滤波和支持向量机的电力线路覆冰评估方法 |
CN105654189A (zh) * | 2015-12-23 | 2016-06-08 | 西安工程大学 | 基于时间序列分析与卡尔曼滤波算法的覆冰短期预测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109146175A (zh) | 2019-01-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107316095B (zh) | 一种耦合多源数据的区域气象干旱等级预测方法 | |
CN107423546B (zh) | 基于Copula函数的多变量水文不确定性处理方法 | |
CN110471024B (zh) | 一种基于量测数据分析的智能电表在线远程校验方法 | |
CN112213727B (zh) | 一种基于主被动微波联合探测的星载雷达的降水订正方法 | |
CN112070286B (zh) | 复杂地形流域的降水预报预警系统 | |
CN107169645B (zh) | 一种计及暴雨灾害影响的输电线路故障概率在线评估方法 | |
CN113159714B (zh) | 一种面向电网的气象数据订正方法 | |
CN103090831A (zh) | 覆冰区域输电线路覆冰厚度的判断方法 | |
CN116908941A (zh) | 一种基于多元线性回归模型的城市近地面臭氧浓度的预测方法 | |
CN111856621B (zh) | 一种基于模式与观测数据融合的集成演变svd转换方法 | |
Yu et al. | Load forecasting based on smart meter data and gradient boosting decision tree | |
CN109146175B (zh) | 一种利用卡尔曼滤波的输电线路覆冰实时估计与预测方法 | |
CN112686426A (zh) | 一种基于水电站流域关键点的来水量预警方法及系统 | |
CN114896872A (zh) | 一种高压输电线路覆冰状态综合评估方法 | |
CN115774834A (zh) | 一种基于多元线性回归模型的城市近地面臭氧浓度的预测方法 | |
Szolgayová et al. | Wavelet based deseasonalization for modelling and forecasting of daily discharge series considering long range dependence | |
CN108663727B (zh) | 利用蒸发率在世界海域范围内估算蒸发波导高度的方法 | |
CN110288117B (zh) | 一种电离层参数临界频率的区域重构方法 | |
CN108920655B (zh) | 一种道路天气信息系统时空覆盖范围量化方法及装置 | |
CN104751003A (zh) | 一种覆冰输电导线倾角预测计算方法 | |
CN104794304A (zh) | 工程海域海平面上升的快速推算方法 | |
CN114398423A (zh) | 一种基于多源数据的河涌水质时空预测方法及系统 | |
CN114564487A (zh) | 预报预测相结合的气象栅格数据更新方法 | |
CN112540165A (zh) | 水土流失预警系统及方法 | |
CN109740111A (zh) | 对地电场值预测算法 |
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 |