不同风机接入电力系统的潮流计算方法比较.docx

上传人:b****5 文档编号:28910641 上传时间:2023-07-20 格式:DOCX 页数:29 大小:206.64KB
下载 相关 举报
不同风机接入电力系统的潮流计算方法比较.docx_第1页
第1页 / 共29页
不同风机接入电力系统的潮流计算方法比较.docx_第2页
第2页 / 共29页
不同风机接入电力系统的潮流计算方法比较.docx_第3页
第3页 / 共29页
不同风机接入电力系统的潮流计算方法比较.docx_第4页
第4页 / 共29页
不同风机接入电力系统的潮流计算方法比较.docx_第5页
第5页 / 共29页
点击查看更多>>
下载资源
资源描述

不同风机接入电力系统的潮流计算方法比较.docx

《不同风机接入电力系统的潮流计算方法比较.docx》由会员分享,可在线阅读,更多相关《不同风机接入电力系统的潮流计算方法比较.docx(29页珍藏版)》请在冰豆网上搜索。

不同风机接入电力系统的潮流计算方法比较.docx

不同风机接入电力系统的潮流计算方法比较

不同风机接入电力系统的潮流计算方法比较中文摘要:

各类风机由于其结构特点的不同决定了其接入电力系统的潮流计算方法的差异,通过对普通异步风机、永磁直驱同步风机、双馈异步风机接入电力系统的潮流计算方法讨论和比较,可以更加清楚的掌握对各类风机的使用。

关键词:

普通异步风机永磁直驱同步风机双馈异步风机电力系统潮流计算正文:

近年来"随着人们对可持续发展战略的认识"风力发电在全球获

得了迅猛发展,截止2005年,全球风电装机容量已达,59332MW,排名前几位分别为德国18428MW,西班牙10027MW,美国9149MW,中国以1246MW排第8位/。

风电场建设包括风资源测量评估、风场选址、发电机选型、风场内部微观选址、风电场升压站建设、接入电网线路建设等一系列问题、其中"风电场选址与接入系统设计都是规划阶段的重要问题"其重要性不言而喻、对含风电场电力系统进行潮流计算"是规划阶段的基础工作"是风电场升压站、接入系统线路方案确定的考虑因素之一,另外"在风电系统动态仿真中"状态变量的初值也由潮流计算获得,因此"对含风电场电力系统进行潮流计算是非常有意义的。

风电场的潮流计算主要是对风电并网后根据给定的网络结构和运行条件确定整个网络的电气状态,主要是对风电并网后的电网中各个节点的电压幅值和相角、网络中功率的分布及功率损耗等,并进行越界检查,以了解和评价风电场并网后的运行情况。

常用于评估风电机组并网后对电网稳态运行的影响,同时也为分析风电场并网后分析对电网影响等其他理论研究工作的基础,具有重要的意义。

国内外对风电潮流的研究有着几十年的历史,风电场的潮流计算主要包括

含普通异步电机的风电场潮流计算、含永磁直驱风机的潮流计算和含双馈异步电机的风电场潮流计算。

风电场的潮流计算关键是正确建立风电机组的数学模型。

电力系统中节点分为PQ节点、PV节点和平衡节点。

一般异步发电机本身没有励磁调节装置,不具有调整节点电压的能力,因此不能像常规的同步发电机一样将它视为电压幅值恒定的PV节点,它只能依靠无功补偿装置才能保持风电场出口电压恒定。

同样风电场中的异步发电机向系统注入有功功率的同时还要从系统吸收一定的无功功率,吸收的无功大小与发电机的机端电压、发出的有功功率以及滑差有着密切相关,因此不能把它处理为恒功率的PQ节点。

一、异步电机计算潮流的方法首先给出了风电场的简化PQ模型、简化RX模型的潮流计算

模型,根据两种模型的缺点,提出了一种新的风电场潮流计算

模型即为扩展潮流模型。

扩展潮流模型将异步发电机直接纳入

电力网络,同时将异步发电机等值电路模型的中间激磁支路的端点作为一个虚拟PQ节点进行计算,建立约束方程;建立异

步发电机转子侧和定子侧之间的转矩平衡方程。

模型建立详细

潮流约束方程,推导雅可比矩阵元素,建立扩展风电潮流算法,

最后通过算例分析验证算法的有效性异步风力发电机潮流建

模一牛顿拉夫逊法潮流计算,牛顿法是解非线性方程的有效的

方法。

它把非线性方程的求解变成反复的求解线性方程,逐步

接近非线性方程的解的过程,通常称为逐次线性化过程。

而且

牛顿一拉夫逊法求解潮流计算具有平方项收敛的速度,能够使

潮流快速收敛。

极坐标的牛顿拉夫逊法求解潮流问题的步骤如下:

(1)计算有功功率和无功功率的不平衡量假设系统中有n号节点,第1-m号节点为PQ节点,第m+1~n-1号为PV节点,n号为平衡节点。

PQ节点不平衡量为:

m

(弟1)

A/J=心-理=珞-片£Vj(G-j+Bysin^-)=0

m

&幺=Qis-Qj-益匕呵血切cos%)=0

PV节点的功率不平衡量为:

(3-2)

(2)计算雅可比矩阵

修正方程如下:

(3-3)

当i=j时:

%"sin%-B.cos吗)

%=-05®+血sin%)

H9=輕(G”cos0ti+Btjsin耳)

弘二-将5sin©-B卩co%)

(3-4)

(3-5)

nH

小啊-Q

(3)潮流方程的求解化成下面的方程

a

1

+J

■F

=0

(3-6)

L」

根据上述步骤进行潮流计算,直至潮流熟练,输出结果

风电机组容量等值计算,等值后的风电机组容量为:

PJV

P严字

*Q讦工Q*其中W'U(3-7)

卜如1

k・F

式中N表示风电场中风机的台数#毘、比为等值后的有功功率、无功功率、视在功率。

风力发电机的电气参数等值如下:

R犷RJNXs^Xs!

N辱二和N(3-8)

其中,&、为定子电阻、定子电抗.R八耳为转予电阻、转子电抗,g为激磁电抗”gXt£>尬、%、*嗣分别为等值后的定子电阻*定子电抗,转子电阻、转子电抗、激碓电抗・

简化结构异步电机潮流算法:

对于常规的含普通异步电机的风电场潮流计算时,一般只考虑

风力发电机和所连接的电网,忽略了风力机部分以及风力机与

发电机之间的联系。

同时,在进行潮流计算时,只考虑风机发出的功率全部输送到与风机连接电网中,不考虑风机由于浆矩

角的变化对输出功率的影响。

在进行风电场的潮流计算时,由

于风机的机械部分的变化跟不上电气部分的变化速度,一般将

风电场节点作为PQ节点进行处理,但PQ节点的无功功率与

风机自身的参数以及风场节点电压有关。

忽略空间和时间谐波、

忽略励磁饱和、忽略铁损,将转子电阻和转子电抗折算到定子

侧,图3.1(a)给出异步电机稳态H等值电路,其中「乩

是定子和转子电阻,」,是定子和转子电抗,'是激磁回路

电抗,S为异步发电机的滑差,儿是对地并联电容。

由于激磁

电抗远大于定子电抗,且定子电阻较小。

因此把励磁支路外移,

合并定子和转子支路。

采用简化等值电路时,门型简化

等值电路如图3.1(b)所示。

大型风电场中的普通异步发电机以超同步转速情况下发电运行,发电机将风力机提供的机

械功率转化为电磁功率输出,同时从电网或无功补偿装置中吸收无功功率来维持励磁电流的大小

X0X舟

■hh—<

—HI

——*

L1

!

(a)发电机n型穩态等值电路

(b)简比1■型等值电爵图3』异步风力境电机的等惬啦賂

RfXf

图3—l(b)所示简化等值电路中,定子电抗与转子电抗合

并」"式。

定义电纳容性为正,电抗感性为负,将激磁支路的电抗与对地电容并联等值电抗为■-I式。

■*—

Xk=x^xr(3-9)

Xp^XcflXm(3*10)

假定异步风机流向电网功率为正,则将异步风机视为发电机处

理,定义风机的输出有功功率:

尸⑺、发电机滑差以及无功功率a

如下:

(3-U)

G-12)

(3-13)

无功功率对节点电压求导得:

dV

(3-14)

 

将风电场节点作为PQ节点,根据风机出力曲线确定不同风速

下的风机的出力,根据潮流建模公式代入潮流计算,采用算例

分析,计算出潮流结果。

简化RX模型:

采用RX模型进行潮流计算的基本思想是两个迭代步骤:

(1)

通过常规潮流计算计算出发电机的端电压;

(2)通过异步发

电机的滑差迭代计算出发电机的滑差。

采用简化RX模型避免

了潮流计算的两步迭代,能够节省潮流计算的计算时间。

图3.1(a)为普通异步发电机的稳态等值电路图,计算出普通异步发电机的电磁功率。

由异步发电机原理知道,异步发电

机发出的电磁功率P*。

与转子电流■人、滑差s等有关,其表达式为

勺―岸岭(S)"(3-15)

由此可见异步发电机的输送到电网的电磁功率E随着滑差s的变化而变化,同时风力机的转速、叶尖速比、风能利用系数以及机械功率也与滑差s有关。

根据系统的功率平衡的关系,风力机的机械功率等于输送到电网的有功功率相等。

由于初始

的功率不相等,通过滑差s的迭代使两个功率最终达到平衡。

采用牛顿拉夫逊计算潮流,引入风力机的机械功率与送到电网的有功功率的差和滑差修正量,修正方程为:

胪&並&錨&叱v7_eF妝卩

色E3/03霜縮劈"

+

■I

 

式(3-17)式中Pw为机械功率,为机械功率与电磁功率差值•式(348),(3-19)中雅可比矩阵的新增元素,但新增元素中的电压了为激臓支路电压.而潮流计尊中使用的风电场的节点电压.

(3-18)

(3-19)

采用此用简化RX模型来计算風电场的潮流,由于新增雅可比矩阵中元素的

程序框图如下:

洁)

it}

o

二、双馈异步风机的潮流计算方法

本节在对双馈风机进行潮流建模时,将双馈风电机组的风电场

作为恒定功率因数的PQ节点进行处理,功率因数值设为0

98。

然后根据双馈异步发电机的稳态等值电路以及双馈风机

建立约束方程,推导雅可比矩阵元素,建立潮流计算模型

的发电系统示意图分析各个功率之间的关系,

及转矩平衡关系,

糧搦怕豊寿屋,计宜出蚯禅

峪节点电压,需定带点盼!

2的fl鏑和功事

••正方稼式Xfti节点也尸广—合金严电机蛆“叭也昨朴』砒

廿JT节点功草度比昭血皐、凤怎崭竿盪斗电血、功宰不出宅fflft奄呻・

H■図黑节点干平场■:

虫只・和计声輿力握电机朗氏MKwMQ忖』卩

图4—l(a)为双馈异步发电机的稳态等值电路,图4.1

(b)为双馈异步电机系统发电示意图。

在图4.1(a),叫所为定子电压,丹以为转子电压弘为转子端到中间节点的

复功率,氐为定子端到中间节点的复功率,其他参数同普通异

步电机等值电路的参数意义相同。

在图4—1(b)中,风电场的节点设为i节点,'为网侧变流器出口电压’分别为风机输入到电网的有功功率和无功功率,'为变压器的阻抗。

系统中的普通节点的潮流计算约束方程按照节点类型建立约束方程,推导雅可比矩阵元素。

风电场节点的潮流计算变量的

约束方程根据双馈风电机组的功率平衡关系以及风力机与发电机之间的转矩平衡条件来建立,并推导出对应的雅可比矩阵

元素。

令儿,以下定义类似。

节点i流向节点W的复功率:

瓦喘哉)

(4-5)

(4-6)

_念俘-祁呷对艸”冷即涸n触’&怦-也血触卜恥如血触

鬲方TR+X?

节点w流向定子侧的复功率为:

W1Wwi

-上怪少竺进sin^vj,COS6^J卜念唧曲瓯

aR+x?

T訴#节点w流向转子侧的复功率为:

“羊・匚!

£

阴护-爲心血卄曲血加.^XrVj山v时aj$脣+舟血阪$“》鬲可"必耐旷%P畀色%(")

转子端流向节点w的复功率为:

=搗+%逐驚沪^牆”窘严(他

网侧变流器端口流出的复功率为:

旷S<^Rr^JXrj

确喘弋閱a跖*肝堆函sin為Xj\晦-呛函嗨詁附御融尙不旷+;需一一(⑷

节点i流向变流器的复功率为:

"%

 

(4-10)

_夕墜咐缺吨十灯乃金咼%「乃■(疔-剤处cos%)■斯附価血血~昭X*7哈靜

(1)风机与外网之间的连接

根据基尔霍夫电流定律,流入节点电流之和为零。

■分别

注入i节点的有功功率和无功功率。

节点i对应于电网有功功率平衡约束方程为:

I丿卄)

I丹

件11)

对应变量的雅可比矩阵的元素为:

 

t_p-rrRySIU酩COS供w

JVt*tp

则务'-2帀-才$%抽触耘卩加cos龜亠眄卜肝I仗血%

“与怒%-窃鈕%

PVw

(4-12)

(4-13)

(4-U)

(4•⑸

*(4-16)

(447)

 

 

节点i对应于电网的无功功率平衡约束方程为:

亠盅■血

)賓呼%2%卜施%血备.帀(呼比險COS链卜砒检询纶

申X*

(4-18)

对应变量的雅可比矩阵元素为:

r-1/1/涉W二尽A®住伸

如一fwR^+Xj

爲血昨t■与cos卷

J=W

伽e他2+72"山夢=7孙7艸

&r碍%血细4'XyllFwCO骑*L2并)旳"临f対11%+才『{冷^003血-2片)“阳、

少=莎?

…+科吗(仇2)

JQVw=¥i

(4-20)

件21)

R件卅

「—Xfgs%+与咖血

^2772

JTLypiAqn

(4-23)

(4-24)

(2)双馈电机的等值电路

在双馈异步电机稳态等值电路中,

虚拟内节点w连接定子支路

和转子支路,以及激磁支路,图4.2为虚拟节点的电路结构。

相比较普通异步电机,转子电压值不为零。

Rs心人孚齐牛

f—EZpw-——«

图虚拟内节点结构

根据节点功率平衡,建立内节点

w有功功率平衡约束方程:

邺勺+坪006$厂予严备丸(4,25)斛胡wr

推导对应变量的雅可比矩阵元素为:

 

PVr

(4-28)

g2)

+i2^

bylS^Xr-RrVVV辞帥睐#+2遍并005跖■中工#由口fj^Vwr;-心g

囲+阳]

74

"叭"」P$7師

i—pR§cos&^j・X在血几孑

鸟巧沖岛-tv^-x^sinewi碍兀w心ws如7価&nr

+*ri7i—"■

j&&血為汁龙$cos為f

虚拟节点w的无功平衡约束方程为:

必"%-瓷-弘

比(⑴沖隔-岡+胪丹in硼培s2XrV23砌00如比血需)—,*「石示动"件列

对应变量的雅可比矩阵的元素为:

(3)双馈风机的外部电路

对于风电场出口节点来讲,不仅连接外电网和双馈异步发电机

的定子侧,还通过变压器连接着网侧变流器,变压器的阻抗为

rysioOwj~&cos&\^i

(4-37)

件38)

(4-40)

J_VV次F血“MT-心605^wr

破必—无叼—

JQdw=^JQffi~JQ^r

M&sin琳计屁cos%i

=艸―両—

t_r/尺八仙0时+对00«日册眇"—

「昭如附5必血&血-%)%2?

Xrrw和&心备+%COS%

p叶工J)_3_—r一-~_口+vr

2sJ^XrV^t站cos漏-IsRfXf鈕漏.2心0氐*浮小甘

氐贬・•为风电场出口处的结构,根据基尔霍夫电流原理,

对1节点注入功率之和为0。

建立1节点与发电机之间功率平

衡约束方程。

4-3风电场的岀口结构

 

对i节点与发电机部分建立有功功率平衡约束方程:

A/L,=-P-P-«P-

协rKref

砒2-%咖树+灿心in知研甲,陆亦时柵险盹竹—、

P呵°°(441)

对应的雅可比矩阵元素为:

『一w与鈕弟揺严%

丿刊厂7弹^—5

Am+

JP&i=~JP&w^JP0c

r念(2巧7务co砌■胡+心7泸in舟v血[2岭-喰严0£珈)+忑『临亡曲氐

JPVi工

(4-42)

(4-43)

(4-44)

件45)

PVc

T一吐叫卄兀血9和二氏与co%-与血血

(4-46)

輕+球

节点i对与发电机部分建立无功功率平衡约束方程:

(4-47)

昨+X*

列甲*诊8S鼠卜刷冷融务项甲比險006%;-丹%®址(4/8)

对应变量的雅可比矩阵的元素为:

QVi

心叫厂牍05%

j勺血%-心倔益

JQ6c-呼竽

斥+於

JQ0i=“巒"JQ9c

-兀(2片-诊cos曲胡-禺金srn妳wXT[2^j-VgCco56icyRrVgCsiTi^c

2r

X

+

JQVc

JQVyv

件49)

(4-50)

(<51)

(4・52)

(4-53)

(4-54)

在稳态运行时,忽略变换器的开关损耗。

根据能量守恒原理,

背靠背变流器输出有功功率之和等于零,即向两侧流出(或注

入)的有功功率不变。

网侧变流器输出无功为给定值一般取零,以防止PV控制方式下,电网需求无功过小,网侧变流器输出

无功倒流至异步电机。

对变流器建立约束方程,并推导雅可比矩阵元素。

4・4为双馈风机背靠背变流器稳态运行的功率关系。

背靠背变流器有功功率的约束方程如下:

野^=叫厂仏

灯(咲-呛函8班沪小卜应咻U1匍昭呼-$咐^期外耳+^齐伶%血乳,(455)

对应变量的雅可比矩阵的元素为:

乜Rfcos伦产AYsinffciPVf^gc

件56)

FVf

PVw

2

fzXrsin^

=*r

坷(2耳-sVwcos+sin6^

与(2金-彳5%)*X*jsin^i

JP0i_KgcKid2“27p&c

耳十X*

t_j/1/S^r^rw+^rw_r

JP9r-斥⑷ni?

2„2-7啄

R?

+s2x}

j_r$cosX"艸sin2jX#❻卩⑷cos

(4-57)

(4-58)

件59)

(4-60)

(4-61)

(4-62)

 

由于变流器之间不考虑无功功率的传输,网侧变流器采用单位

功率因数的控制方式运行,0网的值设为o。

对应的网侧变流

器的无功功率约束方程为:

^t[嚎丹風血詹应囁吩sin^j

%主QCref%=Qcref一―一j」二°件斶

对应变量的雅可比矩阵为:

(4-64)

j臥=齐片Xrsin^j-/?

rcos^

(4-65)

r_vXtcos^j+Jffsin^-QVi=rgc~

砾+輕

r勺(2金叫驱%•卜今乍in%Q%■f2.x2

rLp"Fyl

(4-66)

(4-67)

 

(4)转矩平衡简化模型

双馈异步发电机的电磁功率为:

cos瓯*卩厨sin

Pem=_S

(4・68)

转矩平衡方程为:

=-昱且=0(4-69)

©叫

双馈风机的机械功率采用最大功率跟踪计算,其表达式为:

仏以测(1-$)‘(4-70)

和普通异步电机相比,在双馈异步电机的转矩平衡方程时,考虑了转子电压对功率平衡的影响。

p^_业空空OS站+耳%的n轴

(4-72)

(4-73)

(4・74)

(4-75)

建立潮流计

对应变量的雅可比矩阵元素为:

Ji—Mik竺险+比『诳祈乐:

/如2叫匸并cos為)+Xj^/sin跖

JTVw£

Ts一血/x(I)

根据潮流计算约束方程以及推导雅可比矩阵元素,算模型:

双馈风电机组潮流计算的程序框图:

三、永磁直驱同步风机的潮流计算牛顿法是解非线性方程的有效的方法。

它把非线性方程的求解变成反复的求解线性方程,逐步接近非线性方程的解的过程,通常称为逐次线性化过程。

这是牛顿法的核心。

用牛顿法解三相潮流问题的步骤如下:

(1)、计算功率不平衡方程

AS习迪斗虻\纠\酣宀A/fA/仏硝"呵碱魅…边遂咤Y(3,)

式中表示了一个有n+1个母线的系统功率不平衡矩阵,其中有m个

PQ母线,n-m个PV母线,1个平衡节点。

PQ节点的功率不平衡量为该节点的功率给定值与当前电压计算出来

的实际功率差,可表示为:

冊認-)用岌拧|(q尸曲霁+矿刪。

1jHnfv

*

月tma

(3-2)

其中i=1,2,3….,p=a,b,c。

而对PV节点来说,节点电压幅值是给定的,不再作为变量。

同时,该点无法预先给定无功功率。

这样,该点的无功不平衡量也就失去了约束作用。

因此,在迭代过程中无须计算与PV节点有关的无功功率方程式。

砒認tfI££咛I(%心矿+矿血矿)g

户]tHN

中,FT.Z

只有当迭代结束后,即各节点的电压向量求得后,再求PV节点应当

维持的无功功率。

(3斗)

‘冋炳罗0声心旷_坊尸006矿)祁J

目/J

NJ=*0「«产00^声+耳尸曲1耳尸)详j

njt£)7(qKaBqR+qRsinqk)_2^qR彳川曲

"JJ=lfl7(q严0050"+码严sinqa洋j

其中.

工匕%qR口實+耳尸血貂)月严

屮二干匕%(^鈕^尸一马円cosqv冷丿

厶网二十工叮©尸血矿*_马尸5陽“)十耳严矿"

(4)含风电场的电力系统三相潮流的求解最终能化成求解下列

方程

M=J*AF(3-5)

其中少为节点三相电压的修正列向量:

Ar=[A^Ar/r]r

昭M,△昭…0為,△殆,城屮

△叩號聊/珞幷/咕…仏岛/久,△临哄人唱化J

风电场在牛顿法中的处理

风电场等各种分布式电源可以建立成PQ节点,PV节点,PI节点和P-Q(V)节点这四种节点类型。

对PQ类型的分布式电源只需将它们简单处理成功率值是负的负荷即可。

本节主要分析其他三种类型的分布式电源在程序中的处理。

PV节点可以直接代入牛顿法中处理。

若PV母线与系统通过n(n=1,2,3…)相线路连接,则母线上各节点注入功率为母线总注入功率的n分之一。

在每次迭代后,可以求出节点的电压相角和无功功率。

若计算出的节点无功越限,则将其转换成对应的PQ节点,Q值等于该分布式电源能发出的最大无功值。

如果在后续迭代中,又出现该节点电压越界,重新将其转换成PV节点。

(2)P恒定,电流幅值I恒定的PI节点

PI节点不可以直接代入牛顿法中处理,所以在每次迭代前须做一定的处理。

若PI母线与系统通过n(n=l,2,3…)相线路连接,则母线上各节点注入功率为母线总注入功率的n分之一。

相应的无功

功率可以由上一次迭代得到的电压,给定的电流幅值和有功功率计算得出:

2+严』丽+尸両0-6)

其中,为第k+1次迭代的分布式电源的无功功率值;/斤五分别为第k次迭代得到的电压的实部和虚部5十必北);I为恒定的分布式电源的电流相量的幅值;尸为恒定的有功功率值。

因此在进行潮流计算时,第抖1次迭代前可以把PI节点的无功注入量求出,在第k+1迭代过程中便可将PI节点处理成有功和无功输出分别为P和•的PQ节点。

在每次迭代后,可以求出节点的电压相角和无功功率。

PI型的分布式电源也有无功输出的限制,但从式

(3—6)可以看出,’的标幺值一般在1.0附近,P和,是

两个必需维持的值,所以影响最后计算出来的值的只是该PI

节点的给定有功功率和电流幅值,即P和,若给定得合理,则计算得出的无功功率不会越限。

(3)P恒定,V不定,Q受P、V限定的P-Q(V)节点P-Q(V)节点不可以直接代入牛顿法中处理,所以在每次迭代前须做一

定的处理。

P-Q(V)节点给定的输出有功功率只为异步电机的输

出有功功率,节点电压U在每次迭代后都得到修正,节点的注入无功

r(U2-4帀卅)

s=:

(3-7)

功率Q计算公式如下:

2PeXtr

Q-

4-+a

■Pe

(38)

(3-9)

Qc

QN-Unit

(3-10)

Qn=Qh7咸*Tp

(3-H)

其中,s为异步电机的转差率;

(342)

狗为发电机定子电抗工•与转子电抗*」

之和:

几为励磁电抗;r为转子电阻;Q'为异步电机的吸收无功;

606熾

展开阅读全文
相关资源
猜你喜欢
相关搜索

当前位置:首页 > 医药卫生 > 基础医学

copyright@ 2008-2022 冰豆网网站版权所有

经营许可证编号:鄂ICP备2022015515号-1