CN105977967B - 一种电力***中负电阻的消去方法 - Google Patents

一种电力***中负电阻的消去方法 Download PDF

Info

Publication number
CN105977967B
CN105977967B CN201610383224.4A CN201610383224A CN105977967B CN 105977967 B CN105977967 B CN 105977967B CN 201610383224 A CN201610383224 A CN 201610383224A CN 105977967 B CN105977967 B CN 105977967B
Authority
CN
China
Prior art keywords
negative resistance
resistance
branch
reactance
ground
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.)
Expired - Fee Related
Application number
CN201610383224.4A
Other languages
English (en)
Other versions
CN105977967A (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.)
South China University of Technology SCUT
Original Assignee
South China University of Technology SCUT
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 South China University of Technology SCUT filed Critical South China University of Technology SCUT
Priority to CN201610383224.4A priority Critical patent/CN105977967B/zh
Publication of CN105977967A publication Critical patent/CN105977967A/zh
Application granted granted Critical
Publication of CN105977967B publication Critical patent/CN105977967B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

本发明公开一种电力***中负电阻的消去方法,包括如下步骤:步骤S10,获取简化后电力网络***的潮流结果及新生成含负电阻串联支路的电阻值和电抗值,所述潮流结果是电力网络***化简之前的潮流计算结果;步骤S20,对含负电阻的串联支路进行等电流网络变换,变换为含有负电阻和正电阻的对地导纳支路;步骤S30,将对地导纳支路中的负电阻进行忽略处理,对其中的正电阻进行保留处理,对原含有负电阻的支路参数进行更新,形成新的支路。本发明方法基于潮流计算结果的消去负电阻方法,能够经过少量步骤将负电阻消去而不较大改变***的动态特性,从而满足了电磁暂态仿真软件对含负电阻网络的建模需求及含负电阻网络的仿真需求。

Description

一种电力***中负电阻的消去方法
技术领域
本发明涉及电力***仿真领域,具体涉及一种电力***中负电阻的消去方法。
背景技术
随着现代电力工业的迅速发展,电力***迅速向多机(上千台机)、大电网(几千条线路、几千条母线)、交直流联合输电及大区域联网运行发展,当进行电力***实时仿真或动模实验时,由于仿真规模和硬件设备的限制,必须对原有***做合理简化才能进行。实际上,对一个大电力***的动态性能研究一般只对某一个区域最感兴趣,该区域称之为研究***,而与此区域较远的区域,研究中只要计及其对研究区域的影响,其内部不必详细描写,往往可做降阶及简化,以节省研究的人力和物力,该拟作简化的区域称为外部***。通过***简化可以突出主要矛盾,对于掌握研究***的主要特征,是十分必要的。对于网络化简,CSR法是最为常用的方法,其基本思路是对非线性负荷中的定阻抗负荷、定电流负荷和定阻抗负荷进行分别处理,对于定阻抗负荷,是将其转化为等值阻抗,对于定电流负荷和定功率负荷是将其转化为等值电流移置于保留母线上,然后通过列写节点导纳矩阵并通过高斯消去法,形成新网络的各线路参数及各节点的注入功率。但由于定阻抗负荷的等值阻抗的影响,导致新生成线路的电阻值可能为负值,即出现负电阻。
对于线路负电阻,计算该线路的潮流时,该线路的有功损耗为负值,而对于真实的电力***,线路的有功损耗为正值,即含负电阻的线路的潮流计算与实际情况相违背,从而降低了仿真结果的可信度;对于当前主流的电磁暂态仿真软件无法对电阻值设置为负值,无法描述带负电阻的线路,进而无法对含负电阻的电力网络进行大量的电磁暂态仿真计算,无法满足当前对大电网的仿真模拟工作。目前对于如何处理含负电阻的电力网络的相关文献较少,对于网络化简产生负电阻原因的研究阐述也较少。
目前,对于线路负电阻的研究仅仅停留于分析产生原因的层面上,而且对于产生原因的分析也未深入探索。相关文献提出了一种基于潮流计算结果的逐点线性化的分析方法。该方法的思路是将电网化简过程中的非线性负荷进行逐点线性化处理,对线路等值参数进行跟踪,分别对PQ节点和PV节点进行分析,通过公式推导及简化分析出PQ节点和PV节点产生负电阻的机理。但该方法未对消去节点的非线性负荷的移置问题进行深入研究,导致分析出的产生负电阻的原因是相对片面的;也未联系实际的电力***,如线路电阻与电抗的大小关系,节点注入功率的有功部分与无功部分的大小关系等;也未提出明确的解决负电阻的方案。因此,急需一种处理含负电阻电力网络的方法,能够满足仿真软件的建模要求和仿真需求。
发明内容
本发明的目的在于针对上述现有技术中存在的问题,提出一种在不影响电力网络潮流分布的前提下,对带负电阻额电力网络进行简单的网络变换,从而有效提高网络化简精确度的电力***网络简化方法。
为达到上述发明的目的,本发明通过以下技术方案实现:
一种电力***中负电阻的消去方法,包括如下步骤:
步骤S10,获取简化后电力网络***的潮流结果及新生成含负电阻串联支路的电阻值和电抗值;
步骤S20,对含负电阻的串联支路进行等电流网络变换,变换为含有负电阻和正电阻的对地导纳支路;
步骤S30,将对地导纳支路中的负电阻进行忽略处理,对其中的正电阻进行保留处理,对原含有负电阻的支路参数进行更新,形成新的支路。
进一步,所述等电流变换包括如下步骤:
步骤S21,将含负电阻的串联支路等效变换为并联支路,计算出并联支路的电阻和电抗;
步骤S22,对并联支路的电阻分支通过等电流原则,转化为对地导纳支路。
进一步,所述步骤S22所述的对地导纳支路为不含负电阻的π型电路。
进一步,所述简化后电力网路***是通过CSR法化简的;所述潮流结果包括有节点电压幅值、电压相角、线路有功功率和无功功率。
进一步,所述新生成的含负电阻的串联电路上还串联线路间电抗。
进一步,所述线路间电抗为感性电抗。
进一步,所述串联支路变换为并联支路的等效变化式如下:
其中R为电阻,X为电抗,R1为变换后的电阻,X1为变换后的电抗,上述是将电阻与电抗串联支路变换为电阻与电抗并联支路。
所述π型对地导纳支路的两接地线路的阻抗,即第一接地线路阻抗Z2和第二接地线路阻抗Z3满足下式为:
其中,新生成线路两端的以相角和电压幅值表示的电压分别为V1∠θ1、V2∠θ2,其中V1、V2分别为第一端口、第二端口的电压幅值、∠θ1、∠θ2分别为第一端口和第二端口的电压相角。
本发明的一种电力***中负电阻的消去方法,能够有效的消去电力网络中的负电阻,并且易于实现,只需经过两次网络变换就将电力网络中的负电阻消去,而且不会明显的改变电网的潮流分布及电力网的动态特性。该消去方法是首先统计含负电阻支路的节点信息和线路信息,即节点名称和线路阻抗等参数;然后读取潮流结果文件中对应节点的电压幅值和电压相角,根据等电流原则,对含负电阻的串联支路变换成含负电阻的并联支路,再对负电阻分支进行变换,转化为线路两端的对地导纳支路,对对地导纳支路中的负电阻进行忽略处理,对正电阻进行保留;最后对电力***网络支路参数进行更新。这样,基于潮流计算结果的消去负电阻方法,能够经过少量步骤将负电阻消去而不较大改变***的动态特性。这样,满足了电磁暂态仿真软件对含负电阻网络的建模需求及含负电阻网络的仿真需求。也为研究和设计人员提供一种易于实现且高效的消去负电阻的方法。该发明的提出,对含负电阻网络的仿真计算及推动含负电阻网络的研究具有重大意义。
附图说明
图1为本发明的电力***中负电阻的消去方法的步骤图。
图2a为电力***网络优化后新生成的串联支路结构图。
图2b为与图2a等效的并联支路结构图。
图2c为与图2b等效的对地导纳支路结构图。
具体实施方式
下面结合附图和实施例对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部实施例。
从产生负电阻的原因分析,在形成节点导纳矩阵过程中,在电力***网络中,设节点a为待消去节点,与之关联的节点为b1,b2,…br。则待消去节点a与节点b1,b2,…br的节点导纳矩阵为:
其中y11为待消去节点a的自导纳,y22,y33,…y(r+1)(r+1)为保留节点自导纳,yij为节点i与节点j互导纳。对于待消去节点a的自导纳y11,还需要将节点a处的定阻抗负荷Sza转化成等值导纳yz1并计入y11中,即:
Sza=Pa×fza+jQa×gza (2)
y11′=G11‘+jB11‘=yz1+y11 (4)
其中Ua、fza、gza分别为待消去节点a的节点电压,定阻抗有功、无功负荷系数,Pa、Qa为节点a注入的非线性负荷,G11‘为消去节点a修正后的自电导,B11‘为消去节点a修正后的自电纳。
对于真实电力***,线路电抗为电感,则节点自导纳y11的虚部一般为负值,而对于yz1的虚部则取决于定阻抗无功负荷的性质,当无功负荷为感性时,其值为负;当无功负荷为容性时,其值为正。
令Y11=y11′,Y12=[y12…y1(r+1)],利用高斯法消去节点a后,设新的节点导纳矩阵为:
其中i,j=2,3…r+1。
对于新生成线路节点间的电导,由式(5)取实部可得:
其中G1j、Gi1、Gij分别为消去节点a与节点i、j间互电导,节点i与j间互电导;B1j、Bi1分别为消去节点a与节点i、j间互电纳。对于真实的电力***,线路电阻可认为远小于线路电抗,则式(6)中有:
G1jGi1-B1jBi1≈-B1jBi1≤0 (7)
ABS|G1jBi1+Gi1B1j|<<ABS|-B1jBi1| (8)
式(6)中G11‘和B11‘需将消去节点处的定阻抗负荷转化成的等值阻抗计入其中,真实电力***中,节点处注入的非线性负荷是与受端电网的需求相对应的。对于受端***,有功负荷的需求往往远大于无功负荷的需求,计算出的定阻抗负荷转化成的等值电导远大于等值电纳;且无功负荷的需求往往显感性,即B11‘≤0,有:
ABS|G11‘|>>ABS|B11‘| (9)
由式(7)、(8)、(9),可得:
-B1jBi1G11‘+(G1jBi1+Gi1B1j)B11‘<<0 (10)
当G11‘,B11‘满足一定条件时,由式(6)计算出的Gij‘≥0。新生成线路的电阻值为Gij‘的负倒数,其值为负数,即出现负电阻。
对于定阻抗负荷,可将其分为定阻抗有功负荷和定阻抗无功负荷,定阻抗负荷的大小及负荷的性质直接影响负电阻的产生。
对于定阻抗无功负荷的性质,当无功负荷为感性时,计算出的等值电纳为负值,则使得B11‘更小,由式(6)计算出的新生成线路的电阻值越容易为负值;当无功负荷为容性时,计算出的等值电纳为正值,使得B11‘增大,由式(6)计算出的新生成线路的电阻值能够有效的避免为负值。
非线性负荷是由定阻抗负荷、定电流负荷和定功率负荷共同构成的。对于负荷系数fZa、fIa、fpa、gZa、gIa、gpa,有fZa+fIa+fpa=1,gza+gpa+gIa=1。因此,当待消去节点处的注入非线性负荷固定时,1)定阻抗无功负荷为感性,负荷系数fZa、gZa的值越大,越容易产生负电阻;2)定阻抗无功负荷为容性,负荷系数fZa的值越大,gZa的值越小,越容易产生负电阻。
由上分析出了定阻抗无功负荷的性质及负荷系数的大小与线路电阻值大小的关系。如何消去网络中的负电阻是关键,若直接将串联支路的负电阻消去,则相当于将线路有功损耗直接忽略,对消去负电阻后的网络进行潮流计算必定导致线路潮流分布出现较大误差。
参看图1,本发明的电力***中负电阻的消去方法,包括有如下步骤:
一种电力***中负电阻的消去方法,包括如下步骤:
步骤S10,获取简化后电力网络***的潮流结果及新生成含负电阻串联支路的电阻值和电抗值;在本实施例中,采取从BPA软件中的BPA潮流结果文件中,读取含负电阻支路的潮流结果。
步骤S20,对含负电阻的串联支路进行等电流网络变换,变换为含有负电阻和正电阻的对地导纳支路;
步骤S30,将对地导纳支路中的负电阻进行忽略处理,对其中的正电阻进行保留处理,对原含有负电阻的支路参数进行更新,形成新的支路。
在对电力网络***进行化简之前,先对其做潮流计算,得到包含有节点电压幅值、电压相角、线路有功功率和无功功率的潮流结果,再对所述电力网络***进行化简,在本实施例中采用CSR法化简。
作为具体实施例,所述潮流结果以BPA软件中的潮流结果文件的形式保存,所述新生成含负电阻串联支路的电学参数以BPA软件数据文件的形式保存,所述电学参数包括线路参数、线路上的电阻值和电抗值。
可参看图2a,为电路***网络化简后新增线路的电路结构,在本实施例中,线路参数是线路两端的电压及相角,图中线路两端的以相角和电压幅值表示的电压分别为V1∠θ1、V2∠θ2,其中V1、V2分别为第一端口、第二端口的电压幅值、∠θ1、∠θ2分别为第一端口和第二端口的电压相角;线路中的电阻值和电抗值表示为串联在新生成线路上的负电阻R和线路间电抗X。
参看图2a至图2c,作为具体实施例,所述等电流变换包括如下步骤:
步骤S21,将含负电阻的串联支路等效变换为并联支路,计算出并联支路的电阻和电抗;所述串联支路变换为并联支路的等效变化式如下:
其中R为电阻,X为电抗,R1为变换后的电阻,X1为变换后的电抗,上述是将电阻与电抗串联支路变换为电阻与电抗并联支路,此时变换后的R1依然为负值。
步骤S22,对并联支路的电阻分支通过等电流原则,转化为对地导纳支路。即将R1变换为对地导纳支路,即在上述步骤S22形成的对地导纳支路为不含负电阻的π型电路,如图2c所示,变换后的π型对地导纳支路的第一接地线路阻抗Z2和第二接地线路阻抗Z3满足下式:
对于真实电力***,X>>R,则由式(11)、(12),可得|X1|<<|R1|,即流经R1的电流远小于流经X1的电流。由式(13)、(14)可知,图2c中流入大地的电流等于图2b中流经R1的电流,并且Z2、Z3中的电阻部分不全为负,则可以将计算出的Z2和Z3中的负电阻进行忽略处理,对于计算出的正电阻进行保留,这样可以使潮流分布影响达到最小而且能够有效的消除负电阻,该做法明显优于直接将串联支路的负电阻进行消去处理。
通过对直接忽略负电阻和经过两次电路
本发明的电力***网络化简方法将含负电阻的支路经过两次变换后,形成不含负电阻的π型支路。为了验证变换后支路是否对整个电力***的暂态特性产生影响,首先需要对整个***的潮流计算结果进行比对。选用IEEE 12机162母线电力***进行网络化简。为了消除同调发电机母线的化简及发电机的聚合引起的误差,待消去的节点均选择PQ节点及不对同调发电机进行聚合。网络化简的方案为消去52个网络节点。
对原始***按上述方案进行网络化简,得到化简后的电力***。统计网络化简后产生含负电阻的串联支路,对这些负电阻串联支路进行本专利的方案处理:串联支路变换为并联支路,然后电阻分支变换为对地分支,形成π型电路。变换前后相应支路参数如表1所示,其中支路参数均用标幺值表示。由表1可以看出,并联支路的电抗X1接近于串联支路的电抗X,电抗X2及电抗X3相对于电抗X1而言非常小,电抗X2的绝对值接近于电抗X3的绝对值。
表1变换前后新生成支路相关参数
表2原始***及变换***
对原始***,变换***进行潮流计算,比较与消去节点相邻节点的电压幅值及电压相角,对比结果如表2所示。由表2可以看出变换***与原始***的潮流计算结果有较小误差,偏差角度最大为0.04°,原因在于变换网络过程中,将计算出的对地导纳支路中的负电阻进行消去处理,从而引起传输线路的有功损耗出现偏差,最终导致节点电压相角有了稍许偏差,但这些变化都是在可以接受的范围内。
通过比较变换前后电力***的潮流计算结果,本专利提出的处理方式能够使等值化简***更接近于原始***,为如何处理电力网络中的负电阻提供一种有效的,易于实现的方法。
上述实施例仅用以说明本发明而并非限制本发明所描述的技术方案;因此,尽管本说明书参照上述的各个实施例对本发明已进行了详细的说明,但是,本领域的普通技术人员应当理解,仍然可以对本发明进行修改或者等同替换;而一切不脱离本发明的精神和范围的技术方案及其改进,其均应涵盖在本发明的权利要求范围当中。

Claims (4)

1.一种电力***中负电阻的消去方法,其特征在于,包括如下步骤:
步骤S10,获取简化后电力网络***的潮流结果及新生成含负电阻串联支路的电阻值和电抗值;
步骤S20,对含负电阻的串联支路进行等电流网络变换,变换为含有负电阻和正电阻的对地导纳支路;
所述等电流网络变换包括如下步骤:
步骤S21,将含负电阻的串联支路等效变换为并联支路,计算出并联支路的电阻和电抗;
所述串联支路变换为并联支路的等效变化式如下:
其中R为电阻,X为电抗,R1为变换后的电阻,X1为变换后的电抗,上述是将电阻与电抗串联支路变换为电阻与电抗并联支路,
步骤S22,对并联支路的电阻分支通过等电流原则,转化为对地导纳支路,
所述的对地导纳支路为不含负电阻的π型电路,
所述对地导纳支路的两接地线路的阻抗,即第一接地线路阻抗Z2和第二接地线路阻抗Z3满足下式为:
其中,新生成线路两端的以相角和电压幅值表示的电压分别为V1∠θ1、V2∠θ2,其中V1、V2分别为第一端口、第二端口的电压幅值、∠θ1、∠θ2分别为第一端口和第二端口的电压相角,
步骤S30,将对地导纳支路中的负电阻进行转换,转换流程为:已知负电阻所在节点的电压,因此可将负电阻转化为相应的有功负荷,并将此负荷归并入该节点的连接负荷中,对其中的正电阻进行保留处理,对原含有负电阻的支路参数进行更新,形成新的支路。
2.根据权利要求1所述的电力***中负电阻的消去方法,其特征在于,所述简化后电力网路***是通过CSR法化简的;所述潮流结果包括有节点电压幅值、电压相角、线路有功功率和无功功率。
3.根据权利要求1所述的电力***中负电阻的消去方法,其特征在于,所述新生成的含负电阻的串联电路还串联有线路间电抗。
4.根据权利要求3所述的电力***中负电阻的消去方法,其特征在于,所述线路间电抗为感性电抗。
CN201610383224.4A 2016-05-31 2016-05-31 一种电力***中负电阻的消去方法 Expired - Fee Related CN105977967B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610383224.4A CN105977967B (zh) 2016-05-31 2016-05-31 一种电力***中负电阻的消去方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610383224.4A CN105977967B (zh) 2016-05-31 2016-05-31 一种电力***中负电阻的消去方法

Publications (2)

Publication Number Publication Date
CN105977967A CN105977967A (zh) 2016-09-28
CN105977967B true CN105977967B (zh) 2018-09-25

Family

ID=57010681

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610383224.4A Expired - Fee Related CN105977967B (zh) 2016-05-31 2016-05-31 一种电力***中负电阻的消去方法

Country Status (1)

Country Link
CN (1) CN105977967B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113098409A (zh) * 2021-03-25 2021-07-09 中国矿业大学 一种多个不同形式负电阻的串联及并联方式
CN113054747B (zh) * 2021-04-13 2022-06-07 国网福建省电力有限公司 一种配电台区非线性负荷的非侵入式检测方法

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103336866A (zh) * 2013-06-27 2013-10-02 华南理工大学 一种电磁暂态仿真中含负电阻支路的处理方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5525349B2 (ja) * 2010-06-28 2014-06-18 株式会社東芝 電力系統監視制御装置

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103336866A (zh) * 2013-06-27 2013-10-02 华南理工大学 一种电磁暂态仿真中含负电阻支路的处理方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
交直流互联电网动态等值的实用化方法;张一荻等;《电力自动化设备》;20130228;第33卷(第2期);第120-125段 *

Also Published As

Publication number Publication date
CN105977967A (zh) 2016-09-28

Similar Documents

Publication Publication Date Title
Balamurugan et al. Review of power flow studies on distribution network with distributed generation
JP5815382B2 (ja) 平衡配電系統のための電力潮流解析
CN106532711B (zh) 随迭代和节点类型改变雅可比矩阵的牛顿法潮流计算方法
CN103969552B (zh) 一种分布式发电***的谐波源定位分析方法
CN102609598A (zh) 一种对大型电力***进行电磁暂态仿真的方法
CN105160058B (zh) 一种用于电网电磁暂态计算的网络边界等值方法
CN109617080B (zh) 基于改进的雅可比矩阵的直角坐标牛顿法潮流计算方法
CN106484962B (zh) 一种谐振式无线输电***稳态特性的符号分析方法
CN105977967B (zh) 一种电力***中负电阻的消去方法
CN104636829A (zh) 一种提高电力***温度最优潮流计算效率的解耦算法
CN103956735A (zh) 一种分布式发电***的谐波潮流分析方法
CN105977962B (zh) 一种基于改进节点法的输电网-配电网联合故障分析方法
CN112966364B (zh) 一种用于特征值计算的光伏电站等值模型建模方法及装置
CN106410811B (zh) 首次迭代小阻抗支路端点改变雅可比矩阵的潮流计算方法
Hua et al. Practical method to determine the harmonic contribution of a specific harmonic load
Su et al. The harmonic calculation model of non-linear load
CN113162002A (zh) 一种计及宽频测量环节的直流行波保护方法及***
CN110569558B (zh) 适用于微电网实时仿真的混合电磁暂态仿真方法
CN103838962A (zh) 一种计及pmu量测的分步线性状态估计方法
Campello et al. Representation of multiport rational models in an electromagnetic transients program: Networks with lumped and distributed parameters
CN106712029B (zh) 小阻抗支路pq端点变雅可比矩阵的牛顿法潮流计算方法
Su et al. Advancements on the integration of electromagnetic transients simulator and transient stability simulator
CN105811413B (zh) 一种避免负电阻出现的电力***网络化简方法
Ramzan et al. Accurate harmonic analysis of distribution systems
CN104269872B (zh) 一种变压器三相节点导纳矩阵奇异的处理方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20180925