CN113328457A - 一种提高并网友好性的风电机组虚拟控制方法 - Google Patents

一种提高并网友好性的风电机组虚拟控制方法 Download PDF

Info

Publication number
CN113328457A
CN113328457A CN202011431695.0A CN202011431695A CN113328457A CN 113328457 A CN113328457 A CN 113328457A CN 202011431695 A CN202011431695 A CN 202011431695A CN 113328457 A CN113328457 A CN 113328457A
Authority
CN
China
Prior art keywords
virtual
converter
grid
virtual synchronous
power
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.)
Pending
Application number
CN202011431695.0A
Other languages
English (en)
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.)
Xiajin Power Supply Co Of State Grid Shandong Electric Power Co
Dezhou Power Supply Co of State Grid Shandong Electric Power Co Ltd
Original Assignee
Xiajin Power Supply Co Of State Grid Shandong Electric Power Co
Dezhou Power Supply Co of State Grid Shandong Electric Power Co Ltd
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 Xiajin Power Supply Co Of State Grid Shandong Electric Power Co, Dezhou Power Supply Co of State Grid Shandong Electric Power Co Ltd filed Critical Xiajin Power Supply Co Of State Grid Shandong Electric Power Co
Priority to CN202011431695.0A priority Critical patent/CN113328457A/zh
Publication of CN113328457A publication Critical patent/CN113328457A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/38Arrangements for parallely feeding a single network by two or more generators, converters or transformers
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/11Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2300/00Systems for supplying or distributing electric power characterised by decentralized, dispersed, or local generation
    • H02J2300/20The dispersed energy generation being of renewable origin
    • H02J2300/28The renewable source being wind energy
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E10/00Energy generation through renewable energy sources
    • Y02E10/70Wind energy
    • Y02E10/76Power conversion electric or electronic aspects

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Data Mining & Analysis (AREA)
  • Theoretical Computer Science (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Operations Research (AREA)
  • Algebra (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Control Of Eletrric Generators (AREA)

Abstract

本发明提供一种提高并网友好性的风电机组虚拟控制方法,具体为一种基于惯量阻尼系数自适应调整的风电机组变流器并网控制方法。本发明通过建立励磁绕组电磁动态方程和同步发电机的转子运动方程,得到适用于风电机组接入电网的变流器的虚拟同步发电机控制数学模型。本发明分析转动惯量和阻尼系数对系统稳定性的影响,建立虚拟同步闭环传递函数,优化风电变流器的励磁控制参数,求解得到系统平衡点,提高风电机组并网自适应虚拟同步控制性能,实现风电机组并网最大功率点跟踪及并网适应性,降低了风电机组功率波动性对电网不利影响,提高了高比例风电并网友好性。

Description

一种提高并网友好性的风电机组虚拟控制方法
技术领域
本发明涉及高比例新能源并网控制,更具体地,涉及一种提高并网友好性的风电机组虚拟控制方法。
背景技术
高比例风电通过电力电子器件比如整流器、变流器将电能并网或者输送到负荷,相比于运行成熟的传统发电设备(如同步发电机),电力电子变换器具有响应迅速,控制灵活等优势,但是也有缺乏惯性和阻尼等缺点,从而为电力系统的稳定性带来了一系列的挑战。虚拟同步机技术是一种使得电力电子器件模拟出类似同步电机所具有的旋转惯量和阻尼特性的新型控制技术,为解决电网电力电子化程度不断加深所面临的一系列挑战提供了新的思路。
随着风电新能源占比不断提高,将致使系统惯性能力减弱,给系统的安全稳定运行带来严重威胁。因此,为保证电力系统正常运行,风力发电系统应当具有与传统发电机组类似的惯性和一次调频的电网辅助功能。
发明内容
为应对上述传统方法的不足,本发明一种提高并网友好性的风电机组虚拟控制方法,以通过风电机组并网变流器改善电力系统暂态稳定性和动态品质。该控制器构造的李雅普诺夫函数采用了递归的设计方法,并充分利用了虚拟同步机本身的暂态能量函数。
如图3所示,本发明的具体技术解决方案如下:
S10测量变流器输入输出端口的电压电流频率特征量,计算得到虚拟同步发电机输入及输出功率;
S20结合最大功率点跟踪建立背靠背风电机组并网系统的变流器的虚拟同步发电机的转子运动方程,得到变流器PWM调制波相角;
所述虚拟转子角速度ω及调制波信号相角δ计算公式如下:
Figure BDA0002826743650000011
Figure BDA0002826743650000012
S30建立虚拟同步机的虚拟励磁绕组电磁动态方程;
S40建立具有励磁控制的虚拟同步发电机单机无穷大系统数学模型并定义其平衡点:
所述具有励磁控制的虚拟同步发电机单机无穷大系统数学模型为:
Figure BDA0002826743650000021
所述平衡点代数方程为
Figure BDA0002826743650000022
S50结合步骤四所得模型和其平衡点,递归设计直接补偿阻尼系数的虚拟励磁控制器;
所述设计虚拟励磁控制律为:
Figure BDA0002826743650000023
式中kDωr为直接补偿的虚拟同步机阻尼系数。
李雅普诺夫函数在该控制律条件下为正定,即说明所述系统在平衡点的一定邻域内渐进稳定,且该闭环控制系统在平衡点处是稳定的。并且
Figure BDA0002826743650000024
可推出ωr=0,e=0。又由e的定义可知ωr=0,e=0可得到ΔV′q=0,根据平衡点的方程得到Pset=a1V′qssinδ-a2sin2δ,进而得到δ最终收敛于δs,其收敛速度取决于kD、ke的大小。根据拉萨尔不变集原理可知,系统在平衡点(δs,0,ΔV′qs)处局部渐进稳定。且在平衡点处控制器的输出为稳态励磁控制bV′qs-ccosδs=Vset
S60结合虚拟同步机转子运动方程得到的虚拟功角δ、变流器调制波电压幅值E得到用于变流器控制的PWM调制波信号,如图4所示。该调制波与载波作比较生成的脉冲信号用于控制风电机组并网的三相变流器,实现风电机组并网的虚拟同步发电机控制。
附图说明
图1是风电机组并网单元系统结构图。
图2是虚拟同步控制下变流器对应传统同步机示意图。
图3是本发明提出的控制方法流程图。
图4整体控制流程示意图。
具体实施方式
下面结合实施例和附图对本发明作进一步说明,但不应以此限制本发明的保护范围。
本发明提供一种提高并网友好性的风电机组虚拟控制方法,该方法实施流程图如图3所示,其步骤包括:
S10测量变流器输入输出端口的电压电流频率特征量,计算得到虚拟同步发电机输入及输出功率;图1为适用于该方法的风电机组并网系统结构图。图中风电机组输出的频率幅值变化的交流电通过三相整流器得到幅值恒定的直流电,类似地,该三相整流器也可以采用虚拟同步电动机的控制策略;接着再经DC/AC变流器输出,通过LCL滤波环节后输送到电网。控制变流器内电力电子器件导通和关断的PWM脉冲由本发明控制方法得到的调制波信号E∠δ产生,以实现风电机组并网的VSG控制。
S20结合最大功率点跟踪建立背靠背风电机组并网系统的变流器的虚拟同步发电机的转子运动方程,如图4的上半部分所示,其中MPPT控制算法是通过结合有功/频率下垂控制器的下垂系数的调节信号,通过对下垂系数的改变,使得虚拟拖动转矩的指令值等于MPPT指令值,然后根据虚拟转轴的机械运动特性实现调速器功能。再由虚拟同步发电机的转子运动方程得到变流器PWM调制波相角;
S30建立虚拟同步机的虚拟励磁绕组电磁动态方程;
S40建立具有励磁控制的虚拟同步发电机单机无穷大系统数学模型并定义其平衡点:
具体地,该步骤包括:参照图2,结合步骤二、步骤三得到的转子运动方程、电磁动态方程,推导出虚拟同步发电机控制下的变流器系统数学模型表达式。根据虚拟同步发电机控制下的变流器系统数学模型表达式,定义具有励磁控制的虚拟同步发电机单机无穷大系统平衡点,并通过解非线性代数方程得到系统平衡点,用于虚拟励磁控制器的设计。
S50结合步骤四所得模型及其平衡点,如图4所示,结合前面的有功频率控制模型得出的ωr、Pset、δ等参数,通过构造李雅普诺夫函数,递归设计直接补偿阻尼系数的虚拟励磁控制器,即通过图4中递归计算得到用于控制风电机组并网变流器的PWM调制波信号幅值。
S60根据S50所设计的控制律,得到用于变流器控制的调制波信号幅值E,如图4所示。结合S20得到的虚拟功角δ,即可生成调制波信号E∠δ。该信号与载波作比较,生成的脉冲信号用于控制风电机组变流器开关管的导通和关断。

Claims (8)

1.一种提高并网友好性的风电机组虚拟控制方法,其特征在于,包括以下步骤:
(1)获取变流器输入输出端口的电压电流频率特征量,计算得到虚拟同步发电机输入及输出功率;
(2)结合最大功率点跟踪建立背靠背风电机组并网系统的变流器的虚拟同步发电机的转子运动方程,得到变流器PWM调制波相角;
(3)建立虚拟同步机的虚拟励磁绕组电磁动态方程;
(4)建立具有励磁控制的虚拟同步发电机单机无穷大系统数学模型并定义其平衡点;
(5)递归设计直接补偿阻尼系数的虚拟励磁控制器;
(6)结合虚拟同步机转子运动方程得到的虚拟功角δ、变流器调制波电压幅值E得出变流器调制信号,实现风电机组并网的虚拟同步发电机控制。
2.根据权利要求1所述的方法,其特征在于,所述获取所述虚拟同步发电机输出到电网端的电压电流以及电网频率,根据所测电压电流计算虚拟同步机输出功率,包括:
(1)获取虚拟同步机输出电压电流及电网频率:
通过电压电流互感器测量风电机组并网变流器在三相abc坐标系下的输出侧电压Uabc及输出电流Iabc;通过锁相环(PLL)模块测量电网频率。
(2)计算虚拟同步机输出功率:
利用前述所测输出电压电流,可直接计算虚拟同步机在abc坐标系输出功率,也可以将电压电流转换到dq坐标系再进行计算:
Pout=UaIa+UbIb+UcIc=UdId+UqIq (1)
式中Ud、Uq,Id、Iq分别为母线电压和母线电流在d轴和q轴分量。
3.根据权利要求1所述的方法,其特征在于,所述根据最大功率点跟踪(MPPT)建立背靠背风电机组并网系统的变流器的虚拟同步发电机的转子运动方程,得到变流器PWM调制波相角,包括:
(1)根据当前风速和风电机组的转子转速计算获得的MPPT指令值,结合有功频率下垂控制器的下垂系数的调节信号,通过对下垂系数的改变,使得虚拟拖动转矩的指令值等于MPPT指令值,然后根据虚拟转轴的机械运动特性实现调速器功能:
对于有功-频率下垂控制环节,即
Figure FDA0002826743640000021
式中,Pset为虚拟同步控制下变流器的参考有功输入(可为上层调度指令);Dp为下垂控制系数(由MPPT指令值决定);ωref为角速度参考值。
(2)建立虚拟同步机的转子运动特性方程为:
Figure FDA0002826743640000022
式中,δ[rad]为变流器虚拟功角,表示无穷大母线电压与变流器出口电压之间的夹角;Kd为阻尼系数;ωr[p.u.]为虚拟转子角速度与同步角速度的偏差;ωs=2πf为同步角速度;Pset[p.u.]为虚拟同步控制下变流器的参考功率指令(可为上层调度指令),对应于传统同步机的机械功率Pm;Pe[p.u.]为变流器输出电磁功率。
(3)基于上式(2)(3),建立有功频率控制模型,得到虚拟转子角速度及调制波信号相角为:
Figure FDA0002826743640000023
Figure FDA0002826743640000024
式中δ[rad]为变流器虚拟功角,也即控制变流器的PWM调制波信号相角。
4.根据权利要求1所述的方法,其特征在于,所述虚拟同步机的虚拟励磁绕组电磁动态方程,包括:
(1)采用与同步机励磁绕组电磁动态方程中参数相对应的方式,建立虚拟同步发电机励磁绕组电磁动态方程为:
Figure FDA0002826743640000025
式中,K′d是励磁绕组的时间常数(s);Vset[p.u.]为对应于系统稳态运行时的给定励磁电压;uf为对应于励磁电压的调节量;Vq是变流器出口电压,对应于发电机空载感应电动势;
Figure FDA0002826743640000031
是暂态电势,根据定义Vq与V′q存在关系为:
Figure FDA0002826743640000032
其中,xvir表示虚拟定子电抗,可用虚拟阻抗方法进行模拟;x′vir表示虚拟暂态同步电抗,同样可用虚拟阻抗方法模拟;Id为母线电流d轴分量,可表示为:
Figure FDA0002826743640000033
其中,U为无限大系统的母线电压,可看作常数,x′d∑=x′vir+xl表示为虚拟定子暂态电抗与线路电抗之和。
5.根据权利要求1所述的方法,其特征在于,建立具有励磁控制的虚拟同步发电机单机无穷大系统数学模型并定义其平衡点,包括:
(1)根据虚拟同步机虚拟励磁绕组电磁动态方程,将(7)代入(8)可得到Vq与V′q之间有如下关系成立:
Figure FDA0002826743640000034
(2)根据虚拟同步机电磁功率方程,可得到虚拟同步机有功功率表达式。
已知电动势、电压和电流的关系为
Figure FDA0002826743640000035
将(10)代入(1)得到虚拟同步机送入无穷大系统的有功功率表达式为:
Figure FDA0002826743640000036
(3)根据虚拟同步机转子运动方程(1)和虚拟励磁绕组电磁动态方程(2),最终建立具有励磁控制的虚拟同步发电机单机无穷大系统数学模型为:
Figure FDA0002826743640000041
(4)权利要求2中所述具有励磁控制的虚拟同步发电机单机无穷大系统数学模型简写为:
Figure FDA0002826743640000042
其中,
Figure FDA0002826743640000043
u=Vset+uf
(5)根据控制理论中对系统平衡点的定义,即系统
Figure FDA0002826743640000044
的平衡点是指满足f(x)≡0的点。定义系统(13)的状态变量为δ,ωr,V′q,并且其平衡点用(δsrs,V′qs)表示。根据该定义,系统(13)的平衡点应满足如下条件:
Figure FDA0002826743640000045
当给定变流器预设功率Pset和虚拟励磁电压Vset时,可通过求解(14)非线性代数方程得到系统平衡点。
6.根据权利要求1所述的方法,其特征在于,所述直接补偿阻尼系数的虚拟励磁控制器的递归设计方案,具体包括:
(1)考虑化简后的系统模型(13)中的子系统(δ,ωr),并将暂态电势V′q表示为V′q=V′qs+ΔV′q,则定义ΔV′q=V′q-V′qs,子系统(δ,ωr)可写为
Figure FDA0002826743640000051
其中Pe0=a1V′qssinδ-a2sin2δ、M=2H。
(2)将ΔV′q看作该子系统的虚拟控制律,希望通过调节变流器虚拟励磁使q轴暂态电势的变化ΔV′q满足a1ΔV′qsinδ=M·kDωr,由于存在调节误差,定义调节误差为:
e=M·kDωr-a1ΔV′qsinδ (16)
其中,kD为自定义变量且kD>0,MkD理想补偿变流器虚拟阻尼系数。由此可得:
Figure FDA0002826743640000052
定义
Figure FDA0002826743640000053
构造李雅普诺夫函数为
Figure FDA0002826743640000054
对上式求导得到:
Figure FDA0002826743640000055
选取控制律r为:
r=-ωsωr-kee (20)
其中ke为自定义变量且大于零,则
Figure FDA0002826743640000056
满足
Figure FDA0002826743640000057
由此,可得到直接改善虚拟同步机阻尼系数的励磁控制律u为:
Figure FDA0002826743640000061
7.根据权利要求3中得到的控制风电机组并网变流器的PWM调制波信号相角δ;再将控制律u代入式(9)(10),通过递归计算可得变流器出口电压参考值Vq,变流器调制波电压幅值E可由(23)求出:
Vq=Vdc·E (23)
从而得到变流器调制波电压幅值E。即可得到PWM调制波信号E∠δ,作用于风电机组并网变流器,实现风电机组并网的虚拟同步发电机控制。
8.根据权利要求3,分析转动惯量和阻尼系数对系统稳定性的影响,建立虚拟同步闭环传递函数,优化风电变流器的励磁控制参数,求解得到系统平衡点,提高风电机组并网自适应虚拟同步控制性能。
CN202011431695.0A 2020-12-10 2020-12-10 一种提高并网友好性的风电机组虚拟控制方法 Pending CN113328457A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011431695.0A CN113328457A (zh) 2020-12-10 2020-12-10 一种提高并网友好性的风电机组虚拟控制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011431695.0A CN113328457A (zh) 2020-12-10 2020-12-10 一种提高并网友好性的风电机组虚拟控制方法

Publications (1)

Publication Number Publication Date
CN113328457A true CN113328457A (zh) 2021-08-31

Family

ID=77413293

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011431695.0A Pending CN113328457A (zh) 2020-12-10 2020-12-10 一种提高并网友好性的风电机组虚拟控制方法

Country Status (1)

Country Link
CN (1) CN113328457A (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115882514A (zh) * 2023-02-16 2023-03-31 中国科学院电工研究所 新能源电力系统跟网构网一体化变流器集群聚合控制方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108199396A (zh) * 2018-02-06 2018-06-22 上海交通大学 储能逆变器虚拟励磁闭环控制系统及其设计方法
CN108270241A (zh) * 2018-02-06 2018-07-10 国网四川省电力公司电力科学研究院 风机并网逆变器虚拟同步发电机的控制方法
CN108493984A (zh) * 2018-02-06 2018-09-04 国网四川省电力公司电力科学研究院 适用于光伏并网系统的虚拟同步发电机控制方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108199396A (zh) * 2018-02-06 2018-06-22 上海交通大学 储能逆变器虚拟励磁闭环控制系统及其设计方法
CN108270241A (zh) * 2018-02-06 2018-07-10 国网四川省电力公司电力科学研究院 风机并网逆变器虚拟同步发电机的控制方法
CN108493984A (zh) * 2018-02-06 2018-09-04 国网四川省电力公司电力科学研究院 适用于光伏并网系统的虚拟同步发电机控制方法

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115882514A (zh) * 2023-02-16 2023-03-31 中国科学院电工研究所 新能源电力系统跟网构网一体化变流器集群聚合控制方法

Similar Documents

Publication Publication Date Title
CN108493984B (zh) 适用于光伏并网系统的虚拟同步发电机控制方法
CN108270241B (zh) 风机并网逆变器虚拟同步发电机的控制方法
Peresada et al. Power control of a doubly fed induction machine via output feedback
JP3435474B2 (ja) 変速風車
WO2022077847A1 (zh) 一种用于混合微电网mmc互联变换器的虚拟同步机控制方法
CN107968591B (zh) 基于预测公共点电压的并网逆变器虚拟惯性功率解耦控制方法
US7057908B2 (en) Method and arrangement in connection with network inverter
CN106877710B (zh) 基于虚拟同步电动机的三相pwm整流器多环路控制电路及控制方法
US8415818B2 (en) Method of and apparatus for compensation of oscillation effects in the event of mains asymmetry in a double-fed asynchronous machine
CN110994680B (zh) 一种自适应惯性系数的光伏虚拟同步机稳定控制方法
CN110429655B (zh) 基于同步机三阶模型的储能单元主动支撑控制方法及系统
US10819262B2 (en) Power generation system, system for suppressing sub-synchronous oscillation and method for controlling operation of power system
CN108199396B (zh) 储能逆变器虚拟励磁闭环控制系统及其设计方法
CN108988391B (zh) 基于转速控制的双馈风机转子侧变换器的稳定性分析方法
CN109217335A (zh) 一种海上风电vsc-hvdc输出系统的低频振荡阻尼控制方法
CN104638679A (zh) 一种采用自适应调节的虚拟同步发电机频率控制方法
CN112994111B (zh) 考虑功率控制的永磁同步风机接入弱电网稳定性分析方法
CN113489027B (zh) 抑制虚拟同步发电机频率振荡的虚拟电感控制方法和系统
WO2021110171A1 (zh) 一种基于p-u下垂特性的虚拟直流电机控制方法
Erazo-Damian et al. Stand-alone DFIG FOC sensitivity and stability under mismatched inductances
CN106533289B (zh) 一种非线性电压控制方法及系统
CN116316848A (zh) 一种基于虚拟同步发电机、虚拟同步电动机的微电网中直驱风机协同控制系统及方法
CN111917126A (zh) 基于无锁相环自同步控制的dfig不平衡电网电压补偿方法
Baesmat et al. Pole placement control for doubly-fed induction generators using compact representations in complex variables
EP3930180A1 (en) System and method for reducing voltage distortion from an inverter-based resource

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
RJ01 Rejection of invention patent application after publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20210831