(12)发明专利申请
(10)申请公布号 CN 111106618 A(43)申请公布日 2020.05.05
(21)申请号 201811252315.X(22)申请日 2018.10.25
(71)申请人 中国电力科学研究院有限公司
地址 100192 北京市海淀区清河小营东路
15号
申请人 国家电网有限公司
国网山东省电力公司(72)发明人 程鹏 李庆 张金平 李建立
贺敬 朱琼锋 (74)专利代理机构 北京安博达知识产权代理有
限公司 11271
代理人 徐国文(51)Int.Cl.
H02J 3/38(2006.01)H02J 3/01(2006.01)
权利要求书3页 说明书7页 附图1页
(54)发明名称
一种新能源发电装备接入电力系统的谐波分析方法及装置
(57)摘要
本发明涉及一种新能源发电装备接入电力系统的谐波分析方法及装置,包括:分别获取新能源发电装备接入电力系统的接口变流器等效模型的电量参数,包括独立电压源的瞬时电压值、受控电压源的瞬时电压值和阻抗的阻抗值;根据所述电量参数确定新能源发电装备接入电力系统的谐波电流。本发明提供的技术方案利用新能源发电装备接入电力系统的接口变流器等效模型中获取的电量参数,有效地分析了新能源发电装备输出电流的谐波特性,解决了新能源发电单元并网谐波电流难以评估的实际问题。CN 111106618 ACN 111106618 A
权 利 要 求 书
1/3页
1.一种新能源发电装备接入电力系统的谐波分析方法,其特征在于,所述方法包括:分别获取新能源发电装备接入电力系统的接口变流器等效模型的电量参数;根据所述电量参数确定新能源发电装备接入电力系统的谐波电流。2.如权利要求1所述的方法,其特征在于,所述新能源发电装备接入电力系统的接口变流器等效模型为串联的独立电压源、受控电压源和阻抗;
所述电量参数包括:独立电压源的瞬时电压值、受控电压源的瞬时电压值和阻抗的阻抗值。
3.如权利要求2所述的方法,其特征在于,所述独立电压源的瞬时电压值的确定过程包括:利用独立电压源频率和各次谐波电压幅值uA的对应关系,结合三相正弦波发生原理获取所述独立电压源的瞬时电压值;
按下式计算所述独立电压源在频率f∈(0,+∞)范围内的各次谐波电压幅值uA:上式中,s为拉普拉斯算子,j为复数因子,ω0为工频角频率,icref为接口变流器输出电流指令值,ug为电网相电压幅值;
其中,
kp、ki分别为电流控制闭环参数的比
例系数和积分系数,Td=(0.5~1.5)Ts,为控制系统延迟时间,Ts为控制系统采样时间。
4.如权利要求2所述的方法,其特征在于,所述受控电压源的瞬时电压值的确定过程包括:利用受控电压源的频率及其对应的受控电压源电压幅值uB的对应关系以及受控电压源的角频率及其对应的受控电压源电压幅值uC的对应关系,结合三相正弦波发生原理获取所述受控电压源的瞬时电压值;
当所述受控电压源的频率为(2k-1)ω0时,按下式计算所述受控电压源的频率对应的电压幅值uB:
当所述受控电压源的角频率为kωc±nωr,且k=6m,n=2m-1,m为正整数时,按下式计算所述受控电压源电压幅值uC:
当所述受控电压源的角频率为kωc±nωr,且k=6m,n=2m-1,m为正整数时,按下式计算所述受控电压源电压幅值uC:
上式中,k为正整数,ω0为工频频率,Tdt为调制死区时间,fsw为调制开关频率,Udc为直流电压,为输出电流初始相位角,t为时间,ωc为载波角频率,ωr为调制波角频率,Ji为第一类贝塞尔函数,下标i为阶数,m为调制度。
5.如权利要求2所述的方法,其特征在于,所述接口变流器模型中阻抗的阻抗值Zg的获取过程包括:
于新能源装备与电网连接点的切口处施加的1-1000Hz范围内的频率,测量获得小值电
2
CN 111106618 A
权 利 要 求 书
2/3页
压扰动信号和新能源装备输出的同频率电流反馈;
根据所述小值电压扰动信号和新能源装备输出的同频率电流反馈,按照欧姆定理获取所述接口变流器模型中阻抗的阻抗值。
6.如权利要求1所述的方法,其特征在于,所述根据所述电量参数确定新能源发电装备接入电力系统的谐波电流,包括:
按下式计算所述新能源发电装备接入电力系统的谐波电流id:id=(u-uabc1-uabc2)/Zg上式中,u为电力系统电压,uabc1为独立电压源的瞬时电压值,uabc2为受控电压源的瞬时电压值。
7.一种新能源发电装备接入电力系统的谐波分析装置,其特征在于,所述装置包括:获取模块,用于分别获取新能源发电装备接入电力系统的接口变流器等效模型的电量参数;
确定模块,用于根据所述电量参数确定新能源发电装备接入电力系统的谐波电流。8.如权利要求7所述的装置,其特征在于,所述获取模块,包括:等效单元,用于将所述新能源发电装备接入电力系统的接口变流器等效模型为串联的独立电压源、受控电压源和阻抗;
第一确定单元,用于利用独立电压源频率和各次谐波电压幅值uA的对应关系,结合三相正弦波发生原理获取所述独立电压源的瞬时电压值;
按下式计算所述独立电压源在频率f∈(0,+∞)范围内的各次谐波电压幅值uA:上式中,s为拉普拉斯算子,j为复数因子,ω0为工频角频率,icref为接口变流器输出电流指令值,ug为电网相电压幅值;
其中,
kp、ki分别为电流控制闭环参数的比
例系数和积分系数,Td=(0.5~1.5)Ts,为控制系统延迟时间,Ts为控制系统采样时间。
第二确定单元,用于利用受控电压源的频率及其对应的受控电压源电压幅值uB的对应关系以及受控电压源的角频率及其对应的受控电压源电压幅值uC的对应关系,结合三相正弦波发生原理获取所述受控电压源的瞬时电压值;
当所述受控电压源的频率为(2k-1)ω0时,按下式计算所述受控电压源的频率对应的电压幅值uB:
当所述受控电压源的角频率为kωc±nωr,且k=6m,n=2m-1,m为正整数时,按下式计算所述受控电压源电压幅值uC:
当所述受控电压源的角频率为kωc±nωr,且k=6m,n=2m-1,m为正整数时,按下式计算所述受控电压源电压幅值uC:
3
CN 111106618 A
权 利 要 求 书
3/3页
上式中,k为正整数,ω0为工频频率,Tdt为调制死区时间,fsw为调制开关频率,Udc为直流电压,为输出电流初始相位角,t为时间,ωc为载波角频率,ωr为调制波角频率,Ji为第一类贝塞尔函数,下标i为阶数,m为调制度。
第三确定单元,用于根据新能源装备与电网连接点的切口处施加的1-1000Hz范围内的频率测量得到的小值电压扰动信号和测量得到的新能源装备输出的同频率电流反馈,按照欧姆定理获取所述接口变流器模型中阻抗的阻抗值。
9.如权利要求7所述的装置,其特征在于,所述确定模块,用于:按下式计算所述新能源发电装备接入电力系统的谐波电流id:id=(u-uabc1-uabc2)/Zg上式中,u为电力系统电压,uabc1为独立电压源的瞬时电压值,uabc2为受控电压源的瞬时电压值。
4
CN 111106618 A
说 明 书
1/7页
一种新能源发电装备接入电力系统的谐波分析方法及装置
技术领域
[0001]本发明涉及新能源发电领域,具体涉及一种新能源发电装备接入电力系统的谐波分析方法及装置。
背景技术
[0002]随着电力电子装备的大量采用,电力系统的谐波问题日益严重,特别是在电力电子装备大量集中应用的新能源发电领域,发电装备的谐波问题愈发严重。因此,如何准确分析双馈发电机定子侧的谐波特性就显得尤为重要。[0003]针对新能源发电装备的谐波建模与分析,由于其多采用最大功率追踪控制而只向电力系统输入,最初将其等效为多频理想电流源叠加的形式,即根据外部测量的各次谐波含量在对应频率点处设置相应幅值的理想电流源,多个频率电流源输出相互叠加共同构成新能源发电装备输出电流,从而模拟其谐波电流输出特性。然而,随着电力电子装备应用规模的扩大以及研究的深入,有学者发现电力电子装备由于其闭环控制结构,不再是单纯电流源表现形式,而根据戴维南等效电路或诺顿等效电路,呈现电压源与阻抗串联、电流源与阻抗并联的复合形式。美国弗吉尼亚理工大学在两相同步旋转dq坐标系中提出了基于dq轴阻抗模型的变流器谐波等效模型。然而,在两相同步旋转dq坐标系中,d轴阻抗和q轴阻抗存在耦合,并且这两个变量物理意义模糊,难以与实际物理量相对应。美国伦斯理工学院利用谐波线性化方法,根据对称分量理论,推导了三相并网逆变器在正序、负序下的等效阻抗表达式,并从其时域模型出发,在复频域内得出一个纯等效阻抗模型,并且其中同频正负序的谐波阻抗完全相等,因此得到的二端口网络只含受控源,进一步等效成了更为简单的纯阻抗模型。这些研究成果多应用于新能源发电装备接入的振荡风险分析以及其控制策略的自优化,很少关注新能源发电装备的谐波行为。[0004]因此,现有研究成果在刻画新能源发电装备输出电流谐波特性方面存在不足,无法准确描述新能源发电装备接入电力系统的谐波电流。
发明内容
[0005]为了解决新能源发电装备接入电力系统的谐波特性分析不足的问题,本发明提供一种新能源发电装备接入电力系统的谐波分析方法及装置。通过本方法,确定了新能源发电装备接入电力系统的接口变流器等效模型的电量参数,并根据这些电量参数确定了新能源发电装备接入电力系统的谐波电流,本方法有效实现了对新能源发电装备谐波电流的模拟评估,有利于保障新能源发电装备接入电力系统的安全可靠性。
[0006]本发明提供的一种新能源发电装备接入电力系统的谐波分析方法,其改进之处在于,包括:
[0007]分别获取新能源发电装备接入电力系统的接口变流器等效模型的电量参数;[0008]根据所述电量参数确定新能源发电装备接入电力系统的谐波电流。[0009]优选的,所述新能源发电装备接入电力系统的接口变流器等效模型为串联的独立
5
CN 111106618 A
说 明 书
2/7页
电压源、受控电压源和阻抗;[0010]所述电量参数包括:独立电压源的瞬时电压值、受控电压源的瞬时电压值和阻抗的阻抗值。
[0011]进一步的,所述独立电压源的瞬时电压值的确定过程包括:利用独立电压源频率和各次谐波电压幅值uA的对应关系,结合三相正弦波发生原理获取所述独立电压源的瞬时电压值;
[0012]按下式计算所述独立电压源在频率f∈(0,+∞)范围内的各次谐波电压幅值uA:
[0013][0014]
上式中,s为拉普拉斯算子,j为复数因子,ω0为工频角频率,icref为接口变流器输出电流指令值,ug为电网相电压幅值;
其中,
kp、ki分别为电流控制闭环参数
[0015]
的比例系数和积分系数,Td=(0.5~1.5)Ts,为控制系统延迟时间,Ts为控制系统采样时间。[0016]进一步的,所述受控电压源的瞬时电压值的确定过程包括:利用受控电压源的频率及其对应的受控电压源电压幅值uB的对应关系以及受控电压源的角频率及其对应的受控电压源电压幅值uC的对应关系,结合三相正弦波发生原理获取所述受控电压源的瞬时电压值;
[0017]当所述受控电压源的频率为(2k-1)ω0时,按下式计算所述受控电压源的频率对应的电压幅值uB:
[0018]
当所述受控电压源的角频率为kωc±nωr,且k=6m,n=2m-1,m为正整数时,按下
式计算所述受控电压源电压幅值uC:
[0020]
[0019]
当所述受控电压源的角频率为kωc±nωr,且k=6m,n=2m-1,m为正整数时,按下
式计算所述受控电压源电压幅值uC:
[0022]
[0021]
上式中,k为正整数,ω0为工频频率,Tdt为调制死区时间,fsw为调制开关频率,Udc
为直流电压,为输出电流初始相位角,t为时间,ωc为载波角频率,ωr为调制波角频率,Ji为第一类贝塞尔函数,下标i为阶数,m为调制度。[0024]进一步的,所述接口变流器模型中阻抗的阻抗值Zg的获取过程包括:[0025]于新能源装备与电网连接点的切口处施加的1-1000Hz范围内的频率,测量获得小值电压扰动信号和新能源装备输出的同频率电流反馈;
[0026]根据所述小值电压扰动信号和新能源装备输出的同频率电流反馈,按照欧姆定理获取所述接口变流器模型中阻抗的阻抗值。
6
[0023]
CN 111106618 A[0027]
说 明 书
3/7页
优选的,所述根据所述电量参数确定新能源发电装备接入电力系统的谐波电流,
包括:
按下式计算所述新能源发电装备接入电力系统的谐波电流id:
[0029]id=(u-uabc1-uabc2)/Zg[0030]上式中,u为电力系统电压,uabc1为独立电压源的瞬时电压值,uabc2为受控电压源的瞬时电压值。
[0031]一种新能源发电装备接入电力系统的谐波分析装置,其改进之处在于,包括:[0032]获取模块,用于分别获取新能源发电装备接入电力系统的接口变流器等效模型的电量参数;[0033]确定模块,用于根据所述电量参数确定新能源发电装备接入电力系统的谐波电流。
[0034]优选的,所述获取模块,包括:[0035]等效单元,用于将所述新能源发电装备接入电力系统的接口变流器等效模型为串联的独立电压源、受控电压源和阻抗;[0036]第一确定单元,用于利用独立电压源频率和各次谐波电压幅值uA的对应关系,结合三相正弦波发生原理获取所述独立电压源的瞬时电压值;[0037]按下式计算所述独立电压源在频率f∈(0,+∞)范围内的各次谐波电压幅值uA:
[0038][0039][0028]
上式中,s为拉普拉斯算子,j为复数因子,ω0为工频角频率,icref为接口变流器输出电流指令值,ug为电网相电压幅值;
其中,
kp、ki分别为电流控制闭环参数
[0040]
的比例系数和积分系数,Td=(0.5~1.5)Ts,为控制系统延迟时间,Ts为控制系统采样时间。[0041]第二确定单元,用于利用受控电压源的频率及其对应的受控电压源电压幅值uB的对应关系以及受控电压源的角频率及其对应的受控电压源电压幅值uC的对应关系,结合三相正弦波发生原理获取所述受控电压源的瞬时电压值;[0042]当所述受控电压源的频率为(2k-1)ω0时,按下式计算所述受控电压源的频率对应的电压幅值uB:
[0043]
当所述受控电压源的角频率为kωc±nωr,且k=6m,n=2m-1,m为正整数时,按下
式计算所述受控电压源电压幅值uC:
[0045]
[0044]
当所述受控电压源的角频率为kωc±nωr,且k=6m,n=2m-1,m为正整数时,按下
式计算所述受控电压源电压幅值uC:
[0046]
7
CN 111106618 A[0047]
说 明 书
4/7页
上式中,k为正整数,ω0为工频频率,Tdt为调制死区时间,fsw为调制开关频率,Udc
为直流电压,为输出电流初始相位角,t为时间,ωc为载波角频率,ωr为调制波角频率,Ji为第一类贝塞尔函数,下标i为阶数,m为调制度。[0049]第三确定单元,用于根据新能源装备与电网连接点的切口处施加的1-1000Hz范围内的频率测量得到的小值电压扰动信号和测量得到的新能源装备输出的同频率电流反馈,按照欧姆定理获取所述接口变流器模型中阻抗的阻抗值。[0050]优选的,所述确定模块,用于:
[0051]按下式计算所述新能源发电装备接入电力系统的谐波电流id:[0052]id=(u-uabc1-uabc2)/Zg[0053]上式中,u为电力系统电压,uabc1为独立电压源的瞬时电压值,uabc2为受控电压源的瞬时电压值。
[0054]与最接近的现有技术比,本发明提供的技术方案具有以下优异效果:
[0055]本发明提供了一种新能源发电装备接入电力系统的谐波分析方法及装置,从建立新能源发电装备接入电力系统的接口变流器等效模型的角度,探索根据接口变流器等效模型的电量参数进行新能源发电装备接入电力系统的谐波电流的分析计算方法,解决了新能源发电装备并网谐波电流难以评估的实际问题,有利于保证新能源发电装备接入电力系统的安全可靠性。
[0048]
附图说明
[0056]图1是本发明提供的新能源发电装备接入电力系统的谐波分析方法的流程图;[0057]图2是本发明实施例中新能源发电装备接入电力系统的接口变流器等效模型示意图。
具体实施方式
[0058]下面结合附图对本发明的具体实施方式作详细说明。[0059]实施例一
[0060]本发明实施例提供一种新能源发电装备接入电力系统的谐波分析方法,流程图如图1所示,包括以下步骤:
[0061]分别获取新能源发电装备接入电力系统的接口变流器等效模型的电量参数;[0062]根据所述电量参数确定新能源发电装备接入电力系统的谐波电流。[0063]具体的,新能源发电装备接入电力系统的接口变流器等效模型,如图2所示,为串联的独立电压源、受控电压源和阻抗;接口变流器等效模型的电量参数包括:独立电压源的瞬时电压值、受控电压源的瞬时电压值和阻抗的阻抗值。[0064]独立电压源的瞬时电压值的确定过程包括:利用独立电压源频率和各次谐波电压幅值uA的对应关系,结合三相正弦波发生原理获取所述独立电压源的瞬时电压值;[0065]按下式计算所述独立电压源在频率f∈(0,+∞)范围内的各次谐波电压幅值uA:
8
CN 111106618 A[0066][0067]
说 明 书
5/7页
上式中,s为拉普拉斯算子,j为复数因子,ω0为工频角频率,icref为接口变流器输
出电流指令值,ug为电网相电压幅值;
其中,
kp、ki分别为电流控制闭环参数
[0068]
的比例系数和积分系数,Td=(0.5~1.5)Ts,为控制系统延迟时间,Ts为控制系统采样时间。[0069]受控电压源的瞬时电压值的确定过程包括:利用受控电压源的频率及其对应的受控电压源电压幅值uB的对应关系以及受控电压源的角频率及其对应的受控电压源电压幅值uC的对应关系,结合三相正弦波发生原理获取所述受控电压源的瞬时电压值;[0070]当所述受控电压源的频率为(2k-1)ω0时,按下式计算所述受控电压源的频率对应的电压幅值uB:
[0071]
当所述受控电压源的角频率为kωc±nωr,且k=6m,n=2m-1,m为正整数时,按下
式计算所述受控电压源电压幅值uC:
[0073]
[0072]
当所述受控电压源的角频率为kωc±nωr,且k=6m,n=2m-1,m为正整数时,按下
式计算所述受控电压源电压幅值uC:
[0075]
[0074]
上式中,k为正整数,ω0为工频频率,Tdt为调制死区时间,fsw为调制开关频率,Udc
为直流电压,为输出电流初始相位角,t为时间,ωc为载波角频率,ωr为调制波角频率,Ji为第一类贝塞尔函数,下标i为阶数,m为调制度。
[0077]接口变流器模型中阻抗的阻抗值Zg的获取过程包括:
[0078]于新能源装备与电网连接点的切口处施加的1-1000Hz范围内的频率,测量获得小值电压扰动信号和新能源装备输出的同频率电流反馈;
[0079]根据所述小值电压扰动信号和新能源装备输出的同频率电流反馈,按照欧姆定理获取所述接口变流器模型中阻抗的阻抗值。[0080]具体的,根据独立电压源的瞬时电压值、受控电压源的瞬时电压值和阻抗的阻抗值确定新能源发电装备接入电力系统的谐波电流,包括:
[0081]按下式计算所述新能源发电装备接入电力系统的谐波电流id:[0082]id=(u-uabc1-uabc2)/Zg[0083]上式中,u为电力系统电压,uabc1为独立电压源的瞬时电压值,uabc2为受控电压源的瞬时电压值。[0084]实施例二
[0085]本发明实施例还提供一种新能源发电装备接入电力系统的谐波分析装置,包括:
9
[0076]
CN 111106618 A[0086]
说 明 书
6/7页
获取模块,用于分别获取新能源发电装备接入电力系统的接口变流器等效模型的
电量参数;[0087]确定模块,用于根据所述电量参数确定新能源发电装备接入电力系统的谐波电流。
[0088]具体的,所述获取模块,包括:[0089]等效单元,用于将所述新能源发电装备接入电力系统的接口变流器等效模型为串联的独立电压源、受控电压源和阻抗;[0090]第一确定单元,用于利用独立电压源频率和各次谐波电压幅值uA的对应关系,结合三相正弦波发生原理获取所述独立电压源的瞬时电压值;[0091]按下式计算所述独立电压源在频率f∈(0,+∞)范围内的各次谐波电压幅值uA:
[0092][0093]
上式中,s为拉普拉斯算子,j为复数因子,ω0为工频角频率,icref为接口变流器输出电流指令值,ug为电网相电压幅值;
其中,
kp、ki分别为电流控制闭环参数
[0094]
的比例系数和积分系数,Td=(0.5~1.5)Ts,为控制系统延迟时间,Ts为控制系统采样时间。[0095]第二确定单元,用于利用受控电压源的频率及其对应的受控电压源电压幅值uB的对应关系以及受控电压源的角频率及其对应的受控电压源电压幅值uC的对应关系,结合三相正弦波发生原理获取所述受控电压源的瞬时电压值;[0096]当所述受控电压源的频率为(2k-1)ω0时,按下式计算所述受控电压源的频率对应的电压幅值uB:
[0097]
当所述受控电压源的角频率为kωc±nωr,且k=6m,n=2m-1,m为正整数时,按下
式计算所述受控电压源电压幅值uC:
[0099]
[0098]
当所述受控电压源的角频率为kωc±nωr,且k=6m,n=2m-1,m为正整数时,按下
式计算所述受控电压源电压幅值uC:
[0100][0101]
上式中,k为正整数,ω0为工频频率,Tdt为调制死区时间,fsw为调制开关频率,Udc
为直流电压,为输出电流初始相位角,t为时间,ωc为载波角频率,ωr为调制波角频率,Ji为第一类贝塞尔函数,下标i为阶数,m为调制度。[0103]第三确定单元,用于根据新能源装备与电网连接点的切口处施加的1-1000Hz范围内的频率测量得到的小值电压扰动信号和测量得到的新能源装备输出的同频率电流反馈,按照欧姆定理获取所述接口变流器模型中阻抗的阻抗值。
[0102]
10
CN 111106618 A[0104]
说 明 书
7/7页
具体的,所述确定模块,用于:
[0105]按下式计算所述新能源发电装备接入电力系统的谐波电流id:[0106]id=(u-uabc1-uabc2)/Zg[0107]上式中,u为电力系统电压,uabc1为独立电压源的瞬时电压值,uabc2为受控电压源的瞬时电压值。
[0108]本领域内的技术人员应明白,本申请的实施例可提供为方法、系统、或计算机程序产品。因此,本申请可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且,本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
[0109]本申请是参照根据本申请实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。[0110]这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
[0111]这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。[0112]以上仅为本发明的实施例而已,并不用于限制本发明,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均包含在申请待批的本发明的权利要求范围之内。
11
CN 111106618 A
说 明 书 附 图
1/1页
图1
图2
12
因篇幅问题不能全部显示,请点此查看更多更全内容