[go: up one dir, main page]

CN114911300A - Photovoltaic system multi-pole tracking control method and device based on boundary positioning conductance increment control - Google Patents

Photovoltaic system multi-pole tracking control method and device based on boundary positioning conductance increment control Download PDF

Info

Publication number
CN114911300A
CN114911300A CN202210291054.2A CN202210291054A CN114911300A CN 114911300 A CN114911300 A CN 114911300A CN 202210291054 A CN202210291054 A CN 202210291054A CN 114911300 A CN114911300 A CN 114911300A
Authority
CN
China
Prior art keywords
irradiance
maximum power
photovoltaic
photovoltaic cell
current
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
CN202210291054.2A
Other languages
Chinese (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.)
Chongqing University
Electric Power Research Institute of State Grid Henan Electric Power Co Ltd
Original Assignee
Chongqing University
Electric Power Research Institute of State Grid Henan 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 Chongqing University, Electric Power Research Institute of State Grid Henan Electric Power Co Ltd filed Critical Chongqing University
Priority to CN202210291054.2A priority Critical patent/CN114911300A/en
Publication of CN114911300A publication Critical patent/CN114911300A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05FSYSTEMS FOR REGULATING ELECTRIC OR MAGNETIC VARIABLES
    • G05F1/00Automatic systems in which deviations of an electric quantity from one or more predetermined values are detected at the output of the system and fed back to a device within the system to restore the detected quantity to its predetermined value or values, i.e. retroactive systems
    • G05F1/66Regulating electric power
    • G05F1/67Regulating electric power to the maximum power available from a generator, e.g. from solar cell
    • 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/50Photovoltaic [PV] energy
    • Y02E10/56Power conversion systems, e.g. maximum power point trackers

Landscapes

  • Engineering & Computer Science (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Sustainable Development (AREA)
  • Sustainable Energy (AREA)
  • Power Engineering (AREA)
  • Physics & Mathematics (AREA)
  • Electromagnetism (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Automation & Control Theory (AREA)
  • Photovoltaic Devices (AREA)

Abstract

The invention discloses a photovoltaic system multi-pole tracking control method and device based on boundary positioning conductance increment control, which comprises the following steps: s1, collecting operation data of each photovoltaic cell array; s2, establishing a photovoltaic cell irradiance non-measured calculation model; s3, calculating the irradiance received by each photovoltaic cell, and constructing a positive irradiance vector and an irradiance subscript vector; s4, constructing a photovoltaic array multi-pole maximum value interval positioning model and a maximum power point current analysis model through photovoltaic array dynamic characteristic analysis, calculating a maximum power point reference operation current, and performing multi-pole tracking control based on a conductance incremental method. According to the method, the maximum power point interval can be obtained only by measuring the operating voltage and current of each photovoltaic cell panel, so that energy loss is avoided; and the device calculates the current irradiance as I according to the real-time collected running voltage and current pma And the data source of x calculation avoids the application of an irradiance sensor.

Description

一种基于边界定位电导增量控制的光伏系统多极点追踪控制 方法及装置A method and device for multi-pole tracking control of photovoltaic system based on boundary positioning conductance incremental control

技术领域technical field

本发明属于光伏新能源开发及应用技术领域,具体涉及一种基于 边界定位电导增量控制的光伏系统多极点追踪控制方法及装置。The invention belongs to the technical field of photovoltaic new energy development and application, and in particular relates to a photovoltaic system multi-pole tracking control method and device based on boundary positioning conductance incremental control.

背景技术Background technique

在“碳达峰、碳中和”目标背景下,国家实施了整县屋顶分布式光 伏推进工程。光伏系统在实际运行中,太阳能电池往往受到丛云、树 木、烟囱、设备或者其他结构的阴影的不利影响。此时,处于局部阴 影不均匀光照条件下的光伏系统,整体输出功率能力下降严重,光伏 系统中的每个太阳能电池无法同时达到最大功率点,只要有10%的阵 列面积被遮挡,光伏系统总发电量会下降50%。因此在不均匀辐照度 条件下对太阳能电池串联电路进行多极点追踪控制,有利于实现太阳 能电池串联电路在不均匀光照条件下的能量高效转换。In the context of the goal of "carbon peaking and carbon neutrality", the state has implemented a county-wide rooftop distributed photovoltaic propulsion project. In practical operation of photovoltaic systems, solar cells are often adversely affected by shadows from clouds, trees, chimneys, equipment or other structures. At this time, under the condition of partial shadow and uneven illumination, the overall output power capacity of the photovoltaic system is seriously reduced, and each solar cell in the photovoltaic system cannot reach the maximum power point at the same time. As long as 10% of the array area is blocked, the total output power of the photovoltaic system Power generation would drop by 50%. Therefore, the multi-pole tracking control of the solar cell series circuit under the condition of uneven irradiance is beneficial to realize the efficient energy conversion of the solar cell series circuit under the condition of uneven illumination.

目前,光伏系统的多极点追踪控制多采用粒子群算法,这种方法 以功率突变为重启条件,功率突变并不是遮阴条件的充分条件,扩大 了重启范围;粒子群算法运行重启后的初始化条件大多是将光伏阵列 运行于开路状态或者短路状态进而获得所需要的信息,在环境变化剧 烈情况下,粒子群算法程序频繁重启造成光伏系统产生大量的能量浪 费。因此,研究一种新型的光伏系统多极点追踪控制方法是目前急需 解决的问题。At present, the multi-pole tracking control of photovoltaic systems mostly adopts the particle swarm algorithm. This method takes the sudden change of power as the restart condition. The sudden change of power is not a sufficient condition of the shading condition, which expands the restart range; the initialization condition after the particle swarm algorithm is restarted Most of the photovoltaic arrays are operated in an open-circuit state or a short-circuit state to obtain the required information. In the case of severe environmental changes, the frequent restart of the particle swarm algorithm program causes a large amount of energy waste in the photovoltaic system. Therefore, researching a new type of multi-pole tracking control method for photovoltaic systems is an urgent problem to be solved.

发明内容SUMMARY OF THE INVENTION

为了解决现有技术中存在的不足,本发明的目的是提供一种一种 基于边界定位电导增量控制的光伏系统多极点追踪控制方法及装置。 该方法能够实时确定光伏系统全局最大功率点所在区间边界,避免光伏 系统开、短路运行状态的多极点最大功率点追踪控制,有利于实现光伏 系统能量高效转换。In order to solve the deficiencies in the prior art, the purpose of the present invention is to provide a multi-pole tracking control method and device for a photovoltaic system based on boundary positioning conductance incremental control. This method can determine the interval boundary of the global maximum power point of the photovoltaic system in real time, avoid the multi-pole maximum power point tracking control of the photovoltaic system in the open and short-circuit operation states, and is conducive to the realization of efficient energy conversion of the photovoltaic system.

为了达到上述目的,本发明采用如下技术方案:In order to achieve the above object, the present invention adopts the following technical solutions:

一种基于边界定位电导增量控制的光伏系统多极点追踪控制方 法,包括如下步骤:A multi-pole tracking control method for photovoltaic systems based on boundary positioning conductance incremental control, comprising the following steps:

步骤S1、通过数据采集模块采集各光伏电池阵列的任意时刻的 运行数据和环境数据;Step S1, collecting the operation data and environmental data of each photovoltaic cell array at any time through the data collection module;

步骤S2、建立光伏电池辐照度非量测计算模型;Step S2, establishing a photovoltaic cell irradiance non-measurement calculation model;

步骤S3、通过步骤S2中所述的光伏电池辐照度非量测计算模型 对所述步骤S1中的光伏电池阵列运行数据进行计算,得出各个光伏 电池所接收的辐照度,将所述辐照度按照从低到高正向排列与计算, 构建正向辐照度向量和辐照度下标向量;Step S3: Calculate the operating data of the photovoltaic cell array in the step S1 by using the photovoltaic cell irradiance non-measurement calculation model described in the step S2, and obtain the irradiance received by each photovoltaic cell. The irradiance is arranged and calculated according to the forward order from low to high, and the forward irradiance vector and the irradiance index vector are constructed;

步骤S4、通过光伏阵列动态特性分析,构建光伏阵列多极点最 值区间定位模型和最大功率点电流解析模型,将正向辐照度向量输入 最值区间定位模型进行计算,得到全局最大功率区间所对应的辐照度; 将全局最大功率区间所对应的辐照度输入最大功率点电流解析模型, 得到最大功率点参考运行电流,并基于电导增量法进行多极点追踪控 制。Step S4, through the analysis of the dynamic characteristics of the photovoltaic array, construct a multi-pole maximum value interval positioning model and a maximum power point current analysis model of the photovoltaic array, input the forward irradiance vector into the maximum value interval positioning model for calculation, and obtain the global maximum power interval. Corresponding irradiance; Input the irradiance corresponding to the global maximum power interval into the maximum power point current analytical model to obtain the reference operating current of the maximum power point, and perform multi-pole tracking control based on the conductance increment method.

进一步的,所述步骤S1中运行数据包括各光伏电池阵列任意时 刻的运行电压V和电流I;所述环境数据包括光伏电池阵列横向附近 10米内的空气温度数据TairFurther, the operating data in the step S1 includes the operating voltage V and current I of each photovoltaic cell array at any time; the environmental data includes the air temperature data T air within 10 meters in the lateral vicinity of the photovoltaic cell array.

进一步的,所述步骤S2中建立光伏电池辐照度非量测计算模型 包括以下步骤:Further, establishing a photovoltaic cell irradiance non-measurement calculation model in the step S2 includes the following steps:

步骤S21、建立光伏电池辐照度非量测计算模型f(S)=I,所述计 算模型具体如下式所示:Step S21, establishing a photovoltaic cell irradiance non-measurement calculation model f(S)=I, and the calculation model is specifically shown in the following formula:

Figure BDA0003560115540000031
Figure BDA0003560115540000031

其中,

Figure BDA0003560115540000032
Isc为太阳能电池板的短路电流、Voc为 开路电压、Im为最大功率点电流,Vm为最大功率点电压,Isc、Voc、Im和Vm在标准测试条件下由光伏电池供应商提供;S为光伏电池板接 受的辐照度,e为自然对数,V和I为光伏电池运行电压和电流,Tair为环境温度;in,
Figure BDA0003560115540000032
I sc is the short circuit current of the solar panel, V oc is the open circuit voltage, Im is the maximum power point current, V m is the maximum power point voltage, I sc , V oc , I m and V m are determined by photovoltaic Provided by the battery supplier; S is the irradiance received by the photovoltaic panel, e is the natural logarithm, V and I are the operating voltage and current of the photovoltaic cell, and T air is the ambient temperature;

步骤S22、根据所述步骤S1中采集的电流I、电压V以及环境温 度Tair,构建辐照度求解方程f(S)-I=0。Step S22, according to the current I, the voltage V and the ambient temperature T air collected in the step S1, construct an irradiance solving equation f(S)-I=0.

进一步的,所述步骤S3具体包括以下步骤:Further, the step S3 specifically includes the following steps:

步骤S31、利用所述的辐照度求解方程f(S)-I=0结合所述各个光 伏电池的运行电压V和电流I,计算得出各个光伏电池所接收的辐照 度SiStep S31, using the irradiance to solve the equation f(S)-I=0 and combining the operating voltage V and current I of each photovoltaic cell to calculate the irradiance Si received by each photovoltaic cell;

步骤S32、将步骤1)中得到的辐照度Si按照从小到大进行排列, 建立正向辐照度向量S;Step S32, arranging the irradiance S i obtained in step 1) from small to large to establish a forward irradiance vector S;

步骤S33、将步骤2)中得到的辐照度向量S中各元素按照如下 不等式进行筛选,然后记录满足所述不等式辐照度Si的下标号,构成 所述的辐照度下标向量

Figure BDA0003560115540000041
并且统计满足不等式计算的次数q,即是 极点个数。Step S33, each element in the irradiance vector S obtained in step 2) is screened according to the following inequality, and then record the subscript that satisfies the inequality irradiance S i to form the irradiance subscript vector
Figure BDA0003560115540000041
And count the number of times q that satisfy the inequality calculation, that is, the number of poles.

Figure BDA0003560115540000042
Figure BDA0003560115540000042

进一步的,所述步骤S4具体包括如下步骤:Further, the step S4 specifically includes the following steps:

步骤S41、通过光伏阵列动态特性分析,构建光伏阵列多极点最 值区间定位模型如下:Step S41, through the analysis of the dynamic characteristics of the photovoltaic array, construct a multi-pole maximum interval positioning model of the photovoltaic array as follows:

Figure BDA0003560115540000043
Figure BDA0003560115540000043

其中,i和j是辐照度下标向量

Figure BDA0003560115540000044
的元素,p是光伏阵列串联的电 池板个数,
Figure BDA0003560115540000045
为局部最大功率值大小的分界系数;where i and j are the irradiance index vectors
Figure BDA0003560115540000044
The element of , p is the number of solar panels connected in series in the photovoltaic array,
Figure BDA0003560115540000045
is the demarcation factor of the local maximum power value;

步骤S42、重复判断所述步骤S33中的不等式,以最终满足条件 的辐照度S下标号记为“max”,则全局最大功率点所在区域为(Smax-1, Smax),称为光伏阵列多极点最值区间;In step S42, the inequality in step S33 is repeatedly judged, and the irradiance S that finally meets the condition is marked as "max", then the region where the global maximum power point is located is (S max-1 , S max ), which is called PV array multi-pole maximum value interval;

步骤S43、构建最大功率点电流Ipmax解析模型,所述解析模型如 下式所示:Step S43, constructing an analytical model of the maximum power point current I pmax , the analytical model is shown in the following formula:

Figure BDA0003560115540000046
Figure BDA0003560115540000046

步骤S44、将全局最大功率区间所对应的辐照度输入步骤S43中 所述的最大功率点电流解析模型,得到最大功率点参考运行电流,并 基于电导增量法进行多极点追踪控制。Step S44, input the irradiance corresponding to the global maximum power interval into the maximum power point current analytical model described in step S43, obtain the maximum power point reference operating current, and perform multi-pole tracking control based on the conductance increment method.

进一步的,一种基于边界定位电导增量控制的光伏系统多极点追 踪控制装置,所述装置用于1-5任一项权利要求中所述的方法,包括:Further, a photovoltaic system multi-pole tracking control device based on boundary positioning conductance incremental control, the device is used for the method described in any one of claims 1-5, comprising:

数据采集模块,用于收集步骤S1中所述的光伏电池阵列运行数 据和环境数据;a data acquisition module, used to collect the photovoltaic cell array operation data and environmental data described in step S1;

光伏电池辐照度非量测计算模块,用于建立光伏电池辐照度计算 函数,通过利用光伏电池运行电压、电流和环境温度计算得出各个光 伏电池所接收的辐照度;The photovoltaic cell irradiance non-measurement calculation module is used to establish a photovoltaic cell irradiance calculation function, and calculate the irradiance received by each photovoltaic cell by using the photovoltaic cell operating voltage, current and ambient temperature;

光伏阵列多极点最值区间定位模块,用于将辐照度按照从低到高 正向排列与计算,构建正向辐照度向量、辐照度下标向量,将正向辐 照度向量和辐照度下标向量输入最值区间定位模型进行计算得到全 局最大功率区间所对应的辐照度;The PV array multi-pole maximum interval positioning module is used to arrange and calculate the irradiance from low to high forward, construct the forward irradiance vector and the irradiance index vector, and combine the forward irradiance vector and the irradiance index vector. The irradiance index vector is input into the maximum value interval positioning model to calculate the irradiance corresponding to the global maximum power interval;

最大功率点电流解析计算模块,用于构建所述步骤S43中的最大 功率点电流Ipmax解析模型,将全局最大功率区间所对应的辐照度输入 所述最大功率点电流解析模型。得到最大功率点参考运行电流;The maximum power point current analytical calculation module is used to construct the maximum power point current I pmax analytical model in the step S43, and input the irradiance corresponding to the global maximum power interval into the maximum power point current analytical model. Obtain the reference operating current of the maximum power point;

PWM追踪控制模块,用于将最大功率点参考运行电流作为电导 增量法的初值,通过软件编程嵌入到PWM控制器中,进行最大极点 追踪控制。The PWM tracking control module is used to take the reference operating current of the maximum power point as the initial value of the conductance increment method, and embed it into the PWM controller through software programming to perform maximum pole tracking control.

进一步的,所述光伏电池辐照度非量测计算模块根据光伏电池运 行电压V、电流I和环境温度Tair数据,利用所述步骤S21中的公式, 计算得到辐照度S=f-1(I)。Further, the photovoltaic cell irradiance non-measurement calculation module calculates the irradiance S=f −1 by using the formula in the step S21 according to the photovoltaic cell operating voltage V, current I and ambient temperature T air data. (I).

进一步的,所述光伏阵列多极点最值区间定位模块包括:Further, the photovoltaic array multi-pole maximum interval positioning module includes:

辐照度向量正向排列子模块,用于将Si从小到大排列建立辐照度 向量S;The irradiance vector forward arrangement sub-module is used to arrange Si from small to large to establish the irradiance vector S;

辐照度下标向量子模块,用于将辐照度向量S按顺序进行如下不 等式筛选;The irradiance subscript vector sub-module is used to screen the irradiance vector S according to the following inequality in sequence;

Figure BDA0003560115540000061
Figure BDA0003560115540000061

光伏阵列多极点最值区间定位子模块,用于通过光伏阵列动态特 性分析,构建所述步骤S41中的光伏阵列多极点最值区间定位模型。The photovoltaic array multi-pole maximum value interval positioning sub-module is used for constructing the photovoltaic array multi-pole maximum value interval positioning model in the step S41 by analyzing the dynamic characteristics of the photovoltaic array.

与现有技术相比,本发明具备的有益效果在于:Compared with the prior art, the present invention has the following beneficial effects:

(1)本发明能够迅速确定全局最大功率峰值点所在区间及其上 下边界,将全局最大功率区间所对应的辐照度输入最大功率点电流解 析模型,得到最大功率点参考运行电流,并作为电导增量法的初值, 在最大功率区间内运用电导增量法可迅速追踪到全局最大功率点;本 发明的全局最大功率区间边界能够实时更新,保证最大功率点区间跟 随改变,运行电流只要在该区间,电导增量控制原则会自动寻找到最 大功率点,具有良好的鲁棒性和自适应性;而且本发明的方法在光伏 系统启动时,不需要将整个电路运行至短路或者开路状态,只需要测 量各个太阳能电池板的运行电压电流,即可得到最大功率点区间,可 以避免能量损失。(1) The present invention can quickly determine the interval where the global maximum power peak point is located and its upper and lower boundaries, input the irradiance corresponding to the global maximum power interval into the maximum power point current analytical model, obtain the maximum power point reference operating current, and use it as the conductance The initial value of the incremental method can be quickly tracked to the global maximum power point by using the conductance incremental method in the maximum power interval; the global maximum power interval boundary of the present invention can be updated in real time to ensure that the maximum power point interval changes accordingly, and the operating current only needs to be within In this interval, the conductance increment control principle will automatically find the maximum power point, which has good robustness and adaptability; and the method of the present invention does not need to operate the entire circuit to a short-circuit or open-circuit state when the photovoltaic system is started. It is only necessary to measure the operating voltage and current of each solar panel to obtain the maximum power point range, which can avoid energy loss.

(2)本发明当运行电流I在最大功率点区间外时,输出ΔI(系 统实际输出电流与参考输出电流之差)使系统快速运行到最大功率点 区间,具有良好的快速性。(2) In the present invention, when the operating current I is outside the maximum power point range, the output ΔI (the difference between the actual output current of the system and the reference output current) enables the system to quickly run to the maximum power point range, with good rapidity.

(3)本发明的控制装置根据实时采集的运行电压和电流计算出 当前的辐照度作为Ipmax计算的数据来源,避免了辐照度传感器的运 用。(3) The control device of the present invention calculates the current irradiance as the data source for I pmax calculation according to the operating voltage and current collected in real time, avoiding the use of an irradiance sensor.

附图说明Description of drawings

图1是本发明基于边界定位电导增量控制的光伏系统多极点追 踪控制方法的流程图;Fig. 1 is the flow chart of the photovoltaic system multi-pole tracking control method based on boundary positioning conductance increment control of the present invention;

图2是本发明辐照度变化曲线图;Fig. 2 is the irradiance variation curve diagram of the present invention;

图3是本发明光伏电池温度变化曲线图;Fig. 3 is the temperature change curve diagram of photovoltaic cell of the present invention;

图4是本发明方法和全局扫描法对光伏系统多极点追踪性能对 比图。Fig. 4 is a comparison diagram of the multi-pole tracking performance of the method of the present invention and the global scanning method for a photovoltaic system.

具体实施方式Detailed ways

下面结合附图和实施例对本发明作进一步的详细说明。可以理解 的是,此处所描述的具体实施例仅仅用于解释本发明,而非对本发明 的限定。另外还需要说明的是,为了便于描述,附图中仅示出了与本 发明相关的部分而非全部。The present invention will be further described in detail below in conjunction with the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are only used for explaining the present invention, rather than limiting the present invention. In addition, it should be noted that, for the convenience of description, the drawings only show some but not all of the parts related to the present invention.

实施例Example

一种基于边界定位电导增量控制的光伏系统多极点追踪控制方 法,其流程图如图1所示,包括如下步骤:A multi-pole tracking control method for photovoltaic systems based on boundary positioning conductance incremental control, the flowchart of which is shown in Figure 1, including the following steps:

步骤S1、通过数据采集模块采集各光伏电池阵列的任意时刻的 运行数据和环境数据;所述运行数据包括各光伏电池阵列任意时刻的 运行电压V和电流I,所述环境数据包括光伏电池阵列横向附近10 米内的空气温度数据Tair,如图1中的步骤2所示;Step S1: Collect the operation data and environmental data of each photovoltaic cell array at any time through the data acquisition module; the operation data includes the operation voltage V and the current I of each photovoltaic cell array at any time, and the environmental data includes the horizontal direction of the photovoltaic cell array. Air temperature data T air within 10 meters nearby, as shown in step 2 in Figure 1;

步骤S2、首先,建立光伏电池辐照度非量测计算模型f(S)=I,如 图1中步骤3,所述计算模型具体如下式所示:Step S2, first, establish a photovoltaic cell irradiance non-measurement calculation model f(S)=I, as shown in step 3 in FIG. 1, the calculation model is specifically shown in the following formula:

Figure BDA0003560115540000081
Figure BDA0003560115540000081

其中,

Figure BDA0003560115540000082
Isc为太阳能电池板的短路电流、Voc为 开路电压、Im为最大功率点电流,Vm为最大功率点电压,Isc、Voc、Im和Vm在标准测试条件下由太阳能电池供应商提供;S为光伏电池板 接受的辐照度,e为自然对数,V和I为光伏电池运行电压和电流,Tair为环境温度;in,
Figure BDA0003560115540000082
I sc is the short circuit current of the solar panel, V oc is the open circuit voltage, Im is the maximum power point current, V m is the maximum power point voltage, I sc , V oc , Im and V m are measured by the solar Provided by the battery supplier; S is the irradiance received by the photovoltaic panel, e is the natural logarithm, V and I are the operating voltage and current of the photovoltaic cell, and T air is the ambient temperature;

本实施例选取太阳能电池板参数Isc、Voc、Im、Vm分别为25.44A、 66V、23.25A和54.2V;光伏电池板串联数量为3,3块光伏电池板 分别命名为PV1、PV2和PV3串联组成系统运行,其中,PV1按照 图2、图3中的S1和T1的趋势变化,PV2和PV3按照照图2、图3 中S2和T2的趋势变化;In this embodiment, the parameters I sc , V oc , Im and V m of the solar panels are selected to be 25.44A, 66V, 23.25A and 54.2V respectively; the number of photovoltaic panels in series is 3, and the three photovoltaic panels are named as PV1, PV2 and PV3 are connected in series to form a system operation, wherein PV1 changes according to the trend of S1 and T1 in Figure 2 and Figure 3, and PV2 and PV3 changes according to the trend of S2 and T2 in Figure 2 and Figure 3;

其次,根据所述步骤S1中采集的电流I、电压V以及环境温度 Tair,构建辐照度求解方程f(S)-I=0。Next, according to the current I, the voltage V and the ambient temperature T air collected in the step S1 , an irradiance solving equation f(S)-I=0 is constructed.

步骤S3、首先,利用所述步骤S2中构建的辐照度求解方程 f(S)-I=0结合步骤S1中各个光伏电池的运行电压V和电流I,对f(S)=I 进行快速求解,得到光伏电池PV1、PV2和PV3的辐照度S1、S2、 S3;其次,将得到的所述S1、S2、S3按照从小到大顺序进行排列,建 立正向辐照度向量S,如图1中步骤4所示;最后,将所述的正向辐 照度向量S中的各元素按照如下不等式进行筛选,如图1中步骤4 所示,记录满足所述不等式辐照度Si的下标号,构成辐照度下标向量

Figure BDA0003560115540000091
并且统计满足不等式计算的次数q,也就是极点个数;所述辐照度下 标向量
Figure BDA0003560115540000092
的物理意义是物理意义是该向量中元素对应的电池板的Im附近存在局部极值;Step S3: First, use the irradiance constructed in the step S2 to solve the equation f(S)-I=0 and combine the operating voltage V and current I of each photovoltaic cell in the step S1 to quickly perform a fast calculation on f(S)=I. Solve to obtain the irradiance S 1 , S 2 , S 3 of the photovoltaic cells PV1, PV2 and PV3; secondly, arrange the obtained S 1 , S 2 , S 3 in ascending order to establish the forward irradiance The illuminance vector S, as shown in step 4 in Figure 1; finally, the elements in the forward irradiance vector S are screened according to the following inequality, as shown in step 4 in Figure 1, the record satisfies the The subscript of the inequality irradiance Si, which constitutes the subscript vector of the irradiance
Figure BDA0003560115540000091
And count the times q that satisfy the inequality calculation, that is, the number of poles; the irradiance subscript vector
Figure BDA0003560115540000092
The physical meaning of is that there is a local extremum near the Im of the panel corresponding to the element in the vector;

Figure BDA0003560115540000093
Figure BDA0003560115540000093

步骤S4、首先,通过光伏阵列动态特性分析,构建光伏阵列多 极点最值区间定位模型,如图1中步骤5所示,所述模型具体如下式 所示:Step S4, first, through the analysis of the dynamic characteristics of the photovoltaic array, build a photovoltaic array multi-pole maximum value interval positioning model, as shown in step 5 in Figure 1, the model is specifically shown in the following formula:

Figure BDA0003560115540000094
Figure BDA0003560115540000094

其中,i和j是向量

Figure BDA0003560115540000095
的元素,p是光伏阵列串联的电池板个数,
Figure BDA0003560115540000096
为局部最大功率值大小的分界系数。where i and j are vectors
Figure BDA0003560115540000095
The element of , p is the number of solar panels connected in series in the photovoltaic array,
Figure BDA0003560115540000096
is the demarcation factor of the local maximum power value.

实际中,(0.0005Sj-0.5)间没有数量级的差异,可忽略 ln(e+0.0005Sj-0.5)对

Figure RE-GDA0003724262310000096
的影响;通常情况下峰值电压Vm约是开 路电压Voc的0.8倍,则
Figure RE-GDA0003724262310000097
的典型值如表1所示。In practice, there is no order of magnitude difference between (0.0005S j -0.5), and the ln(e+0.0005S j -0.5) pair can be ignored.
Figure RE-GDA0003724262310000096
Influence; usually the peak voltage V m is about 0.8 times the open circuit voltage V oc , then
Figure RE-GDA0003724262310000097
Typical values are shown in Table 1.

表1局部最大功率分界系数的典型值Table 1 Typical values of local maximum power demarcation factor

Figure BDA00035601155400000911
Figure BDA00035601155400000911

Figure BDA0003560115540000101
Figure BDA0003560115540000101

其次,重复判断所述步骤S3中的不等式,以最终满足条件的辐 照度S下标号记为“max”,则全局最大功率点所在区域为(Smax-1,Smax), 称为光伏阵列多极点最值区间(GMPPZ,Global Maximum Power Point Zone),如图1中的步骤5所示;Next, the inequality in step S3 is repeatedly judged, and the irradiance S that finally meets the condition is marked as "max", then the region where the global maximum power point is located is (S max-1 , S max ), which is called photovoltaic Array multi-pole maximum zone (GMPPZ, Global Maximum Power Point Zone), as shown in step 5 in Figure 1;

然后,构建最大功率点电流Ipmax解析模型,如图1中步骤6中所 示,所述解析模型具体如下式所示:Then, an analytical model of the maximum power point current I pmax is constructed, as shown in step 6 in FIG. 1 , and the analytical model is specifically shown in the following formula:

Figure BDA0003560115540000102
Figure BDA0003560115540000102

最后,将全局最大功率区间所对应的辐照度输入所述的最大功率 点电流解析模型,得到最大功率点参考运行电流,并基于电导增量法 进行多极点追踪控制,如图1中的步骤8所示。Finally, input the irradiance corresponding to the global maximum power interval into the maximum power point current analytical model to obtain the maximum power point reference operating current, and perform multi-pole tracking control based on the conductance increment method, as shown in the steps in Figure 1 8 shown.

图4中展示了本发明的方法与全局扫描法对光伏系统多极点追 踪控制效对比效果。如图4所示,在0-0.1s时间段,光伏系统快速运 行到全局最大功率点区间,并结合电导增量法稳定运行在最值点;而 全局扫描法,将功率曲线运行至底部,导致功率损失加大,同时运行 到最大功率点的时间较长。0.3s时环境发生剧烈变化,多极点最值区 间边界发生改变,当前运行点与更新后的多极点最值区间相距较近, 经过短暂的暂态过程迅速稳定在新的最值点,在0.3s时刻几乎无波动 的过渡到新的工况,而全局扫描法有波动过程。0.5s时以同样的暂态 过程迅速过渡到新的平衡状态,保持全局最大功率输出。Fig. 4 shows the comparative effect of the method of the present invention and the global scanning method on the multi-pole tracking control effect of the photovoltaic system. As shown in Figure 4, in the 0-0.1s time period, the photovoltaic system quickly runs to the global maximum power point range, and combines with the conductance incremental method to stably run at the maximum point; while the global scanning method runs the power curve to the bottom, This results in increased power loss and a longer time to run to the maximum power point. At 0.3s, the environment changes drastically, the boundary of the multi-pole maximum interval changes, the current operating point is relatively close to the updated multi-pole maximum interval, and after a short transient process, it quickly stabilizes at the new maximum point, at 0.3 At time s, there is almost no fluctuation in transition to the new working condition, while the global scanning method has a fluctuation process. At 0.5s, it quickly transitions to a new equilibrium state with the same transient process, maintaining the global maximum power output.

一种用于上述基于边界定位电导增量控制的光伏系统多极点追 踪控制方法的装置,包括:A device for the above-mentioned multi-pole tracking control method of photovoltaic system based on boundary positioning conductance incremental control, comprising:

数据采集模块,用于收集步骤S1中所述的光伏电池阵列运行数 据和环境数据;所述光伏电池运行数据包括光伏电池运行电压、电流 数据,环境数据包括光伏电池横向附近10米内的空气温度数据;A data acquisition module for collecting the photovoltaic cell array operating data and environmental data described in step S1; the photovoltaic cell operating data includes photovoltaic cell operating voltage and current data, and the environmental data includes air temperature data within 10 meters of the lateral vicinity of the photovoltaic cell ;

光伏电池辐照度非量测计算模块,用于建立光伏电池辐照度计算 函数,通过利用光伏电池运行电压、电流和环境温度,再根据步骤 S2中建立的光伏电池辐照度非量测计算模型,计算得出各个光伏电 池所接收的辐照度S,所述辐照度S的计算式如下:The photovoltaic cell irradiance non-measurement calculation module is used to establish a photovoltaic cell irradiance calculation function. By using the photovoltaic cell operating voltage, current and ambient temperature, and then according to the photovoltaic cell irradiance non-measurement calculation established in step S2 model, the irradiance S received by each photovoltaic cell is calculated, and the calculation formula of the irradiance S is as follows:

S=f-1(I)S = f -1 (I)

光伏阵列多极点最值区间定位模块,其包括辐照度向量正向排列 子模块、辐照度下标向量子模块和光伏阵列多极点最值区间定位子模 块;The photovoltaic array multi-pole maximum interval positioning module, which includes the irradiance vector forward arrangement sub-module, the irradiance subscript vector sub-module and the photovoltaic array multi-pole maximum value interval positioning sub-module;

所述辐照度向量正向排列子模块,用于将辐照度Si从小到大排列 建立正向辐照度向量S;The irradiance vector forward arrangement submodule is used for arranging the irradiance S i from small to large to establish a forward irradiance vector S;

所述辐照度下标向量子模块,用于将所述正向辐照度向量S按顺 序进行如下不等式筛选The irradiance subscript vector submodule is used to perform the following inequality screening on the forward irradiance vector S in order

Figure BDA0003560115540000111
Figure BDA0003560115540000111

记录满足不等式辐照度Si的下标号,构成辐照度下标向量

Figure BDA0003560115540000112
并 且统计满足不等式计算的次数q,也就是极点个数;Record the subscripts that satisfy the inequality irradiance S i to form the irradiance subscript vector
Figure BDA0003560115540000112
And count the number of times q that satisfy the inequality calculation, that is, the number of poles;

所述光伏阵列多极点最值区间定位子模块:用于通过光伏阵列动 态特性分析,构建光伏阵列多极点最值区间定位模型,所述模型具体 如下式所示:The photovoltaic array multi-pole maximum interval positioning submodule: used to construct a photovoltaic array multi-pole maximum interval positioning model by analyzing the dynamic characteristics of the photovoltaic array, and the model is specifically shown in the following formula:

Figure BDA0003560115540000121
Figure BDA0003560115540000121

其中,i和j是辐照度下标向量

Figure BDA0003560115540000122
的元素,p是光伏阵列串联的电 池板个数,
Figure BDA0003560115540000123
为局部最大功率值大小的分界系数;where i and j are the irradiance index vectors
Figure BDA0003560115540000122
The element of , p is the number of solar panels connected in series in the photovoltaic array,
Figure BDA0003560115540000123
is the demarcation factor of the local maximum power value;

重复判断所述步骤S3中的不等式,以最终满足条件的辐照度S 下标号记为“max”,则全局最大功率点所在区域为(Smax-1,Smax),称 为光伏阵列多极点最值区间(GMPPZ,GlobalMaximum Power Point Zone);Repeatedly judging the inequality in step S3, and denoting the irradiance S that finally meets the conditions as "max", then the area where the global maximum power point is located is (S max-1 , S max ), which is called the photovoltaic array multi- Pole maximum interval (GMPPZ, GlobalMaximum Power Point Zone);

最大功率点电流解析计算模块,用于构建所述步骤S4中的最大 功率点电流Ipmax解析模型,将全局最大功率区间所对应的辐照度输入 所述最大功率点电流解析模型,得到最大功率点参考运行电流;The maximum power point current analytical calculation module is used to construct the maximum power point current I pmax analytical model in the step S4, and input the irradiance corresponding to the global maximum power interval into the maximum power point current analytical model to obtain the maximum power point reference running current;

所述构建的最大功率点电流Ipmax,具体如下式所示:The constructed maximum power point current I pmax is specifically shown in the following formula:

Figure BDA0003560115540000124
Figure BDA0003560115540000124

PWM追踪控制模块,用于将最大功率点电流作为光伏电池参考 运行电流,光伏电池参考运行电流作为电导增量追踪控制的初值,基 于电导增量法进行多极点追踪控制。The PWM tracking control module is used to use the maximum power point current as the reference operating current of the photovoltaic cell, and the reference operating current of the photovoltaic cell as the initial value of the conductance incremental tracking control, and perform multi-pole tracking control based on the conductance incremental method.

上述仅为本发明的较佳实施例及所运用技术原理。本领域技术人 员会理解,本发明不限于这里所述的特定实施例,对本领域技术人员 来说能够进行各种明显的变化、重新调整和替代而不会脱离本发明的 保护范围。因此,虽然通过以上实施例对本发明进行了较为详细的说 明,但是本发明不仅仅限于以上实施例,在不脱离本发明构思的情况 下,还可以包括更多其他等效实施例,而本发明的范围由所附的权利 要求范围决定。The above are only the preferred embodiments of the present invention and the applied technical principles. Those skilled in the art will understand that the present invention is not limited to the specific embodiments described herein, and various obvious changes, readjustments and substitutions can be made to those skilled in the art without departing from the protection scope of the present invention. Therefore, although the present invention has been described in detail through the above embodiments, the present invention is not limited to the above embodiments, and can also include more other equivalent embodiments without departing from the concept of the present invention. The scope is determined by the scope of the appended claims.

Claims (8)

1. A photovoltaic system multi-pole tracking control method based on boundary positioning conductance increment control is characterized by comprising the following steps:
step S1, collecting the operation data and the environment data of each photovoltaic cell array at any time through a data collection module;
step S2, establishing a photovoltaic cell irradiance non-measured calculation model;
step S3, calculating the photovoltaic cell array operation data in the step S1 through the photovoltaic cell irradiance non-measurement calculation model in the step S2 to obtain irradiance received by each photovoltaic cell, arranging and calculating the irradiance according to a low-to-high positive direction, and constructing a positive irradiance vector and an irradiance subscript vector;
step S4, constructing a photovoltaic array multi-pole maximum interval positioning model and a maximum power point current analysis model through photovoltaic array dynamic characteristic analysis, inputting a forward irradiance vector into the maximum interval positioning model for calculation, and obtaining irradiance corresponding to a global maximum power interval; and inputting irradiance corresponding to the global maximum power interval into a maximum power point current analysis model to obtain a maximum power point reference running current, and performing multi-pole tracking control based on a conductance increment method.
2. The photovoltaic system multi-pole tracking control method based on the boundary positioning conductance increment control of claim 1, wherein the operation data in the step S1 includes an operation voltage V and a current I of each photovoltaic cell array at any time; the environmental data includes air temperature data T within 10 meters of lateral proximity of the photovoltaic cell array air
3. The photovoltaic system multi-pole tracking control method based on the border-positioned conductance increment control as claimed in claim 2, wherein said step S2 of establishing a non-measured calculation model of irradiance of the photovoltaic cell comprises the steps of:
step S21, establishing a non-measured calculation model f (S) -I of irradiance of the photovoltaic cell, where the calculation model is specifically shown as follows:
Figure FDA0003560115530000021
wherein,
Figure FDA0003560115530000022
I sc short-circuit current, V, for solar panels oc Is an open circuit voltage, I m Is the maximum power point current, V m Is the maximum power point voltage, I sc 、V oc 、I m And V m Provided by the photovoltaic cell supplier under standard test conditions; s is the irradiance received by the photovoltaic cell panel, e is the natural logarithm, V and I are the operating voltage and current of the photovoltaic cell, T air Is ambient temperature;
step S22, collecting current I, voltage V and environment temperature T according to the step S1 air And constructing an irradiance solving equation f (S) -I ═ 0.
4. The photovoltaic system multi-pole tracking control method based on the boundary positioning conductance increment control of claim 3, wherein the step S3 specifically comprises the following steps:
step S31, solving equation f (S) -I ═ 0 by using the irradiance, and calculating irradiance S received by each photovoltaic cell by combining operating voltage V and current I of each photovoltaic cell i
Step S32, obtaining irradiance S in step 1) i Arranging according to the sequence from small to large, and establishing a positive irradiance vector S;
step S33, screening elements in the irradiance vector S obtained in the step 2) according to the following inequality, and then recording irradiance S meeting the inequality i Forming said irradiance subscript vector
Figure FDA0003560115530000031
And counting the times q of inequality calculation, namely the number of poles.
Figure FDA0003560115530000032
5. The photovoltaic system multi-pole tracking control method based on the boundary positioning conductance increment control according to claim 4, wherein the step S4 specifically comprises the following steps:
step S41, through analysis of dynamic characteristics of the photovoltaic array, a positioning model of the maximum and minimum multi-pole interval of the photovoltaic array is constructed as follows:
Figure FDA0003560115530000033
wherein i and j are vectors
Figure FDA0003560115530000034
P is the number of solar panels connected in series by the photovoltaic array,
Figure FDA0003560115530000035
a dividing coefficient is the size of the local maximum power value;
step S42, repeatedly judging the inequality in the step S33, and marking the irradiance S which finally meets the condition as the maximum power point as the maximum power (max), wherein the area where the global maximum power point is located is (S) max-1 ,S max ) The maximum value interval of multiple poles of the photovoltaic array is called;
step S43, constructing the maximum power point current I pmax An analytical model, the analytical model being represented by the following formula:
Figure FDA0003560115530000036
and step S44, inputting the irradiance corresponding to the global maximum power interval into the maximum power point current analysis model in the step S43 to obtain a maximum power point reference running current, and performing multi-pole tracking control based on a conductance increment method.
6. A photovoltaic system multi-pole tracking control device based on boundary positioning conductance increment control, which is used for the method of any one of claims 1-5, and comprises:
a data acquisition module for collecting the photovoltaic cell array operation data and the environmental data in step S1;
the photovoltaic cell irradiance non-measurement calculation module is used for establishing a photovoltaic cell irradiance calculation function and calculating irradiance received by each photovoltaic cell by utilizing the operating voltage, current and ambient temperature of the photovoltaic cell;
the photovoltaic array multi-pole maximum interval positioning module is used for arranging and calculating irradiance in a forward direction from low to high, constructing a forward irradiance vector and an irradiance subscript vector, and inputting the forward irradiance vector and the irradiance subscript vector into a maximum interval positioning model to calculate so as to obtain irradiance corresponding to a global maximum power interval;
a maximum power point current analysis and calculation module for constructing the maximum power point current I in the step S43 pmax And the maximum power point current analysis model is used for inputting the irradiance corresponding to the global maximum power interval into the maximum power point current analysis model. Obtaining a maximum power point reference operating current;
and the PWM tracking control module is used for taking the maximum power point reference running current as an initial value of a conductance increment method, embedding the initial value into the PWM controller through software programming and carrying out maximum pole tracking control.
7. The device as claimed in claim 6, wherein the photovoltaic system multi-pole tracking control device based on the border-positioned conductance increment control is characterized in that the photovoltaic cell irradiance non-measured calculation module is used for calculating the current I and the ambient temperature T according to the photovoltaic cell operating voltage V and the ambient temperature T air Using the formula in step S21, the irradiance S ═ f is calculated -1 (I)。
8. The photovoltaic system multipole point tracking control device based on boundary positioning conductance increment control according to claim 6, wherein the photovoltaic array multipole point maximum interval positioning module comprises:
irradiance vector positive direction arrangement sub-dieBlock for coupling S i Arranging the irradiance vectors S from small to large;
the irradiance subscript vector submodule is used for screening the irradiance vector S according to the following inequality in sequence;
Figure FDA0003560115530000051
and the photovoltaic array multipole point maximum interval positioning submodule is used for constructing the photovoltaic array multipole point maximum interval positioning model in the step S41 through photovoltaic array dynamic characteristic analysis.
CN202210291054.2A 2022-03-23 2022-03-23 Photovoltaic system multi-pole tracking control method and device based on boundary positioning conductance increment control Pending CN114911300A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210291054.2A CN114911300A (en) 2022-03-23 2022-03-23 Photovoltaic system multi-pole tracking control method and device based on boundary positioning conductance increment control

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210291054.2A CN114911300A (en) 2022-03-23 2022-03-23 Photovoltaic system multi-pole tracking control method and device based on boundary positioning conductance increment control

Publications (1)

Publication Number Publication Date
CN114911300A true CN114911300A (en) 2022-08-16

Family

ID=82763203

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210291054.2A Pending CN114911300A (en) 2022-03-23 2022-03-23 Photovoltaic system multi-pole tracking control method and device based on boundary positioning conductance increment control

Country Status (1)

Country Link
CN (1) CN114911300A (en)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20100117623A1 (en) * 2008-11-11 2010-05-13 Fife John M System and method of determining maximum power point tracking for a solar power inverter
CN103475267A (en) * 2013-08-14 2013-12-25 浙江工业大学 Method for improving generating efficiency of photovoltaic array
CN103472884A (en) * 2013-08-14 2013-12-25 浙江工业大学 Global maximum power point tracking method for photovoltaic array
CN106774607A (en) * 2016-12-02 2017-05-31 东北电力大学 A kind of photovoltaic maximum power tracking of on-line amending open-circuit voltage

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20100117623A1 (en) * 2008-11-11 2010-05-13 Fife John M System and method of determining maximum power point tracking for a solar power inverter
CN103475267A (en) * 2013-08-14 2013-12-25 浙江工业大学 Method for improving generating efficiency of photovoltaic array
CN103472884A (en) * 2013-08-14 2013-12-25 浙江工业大学 Global maximum power point tracking method for photovoltaic array
CN106774607A (en) * 2016-12-02 2017-05-31 东北电力大学 A kind of photovoltaic maximum power tracking of on-line amending open-circuit voltage

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
乐全明等: "不均匀光照条件下太阳能电池串联电路特性分析及GMPPT控制", 《光子学报》 *

Similar Documents

Publication Publication Date Title
KR101326420B1 (en) System and method of determining maximum power point tracking for a solar power inverter
CN104680424B (en) The voltage power situation predictor method in large-sized photovoltaic power station
CN101783621B (en) Global maximum power point tracking method of photovoltaic generating system and system device
CN108614612B (en) Method and system for tracking maximum power of solar photovoltaic cell
CN103472884B (en) A kind of overall maximum power point of photovoltaic array tracking
CN103730910B (en) The Dynamic Equivalence that a kind of large-scale photovoltaic power station is grid-connected
CN103838291B (en) A kind of maximum power point tracing method of solar cell
CN105930918B (en) Global Distribution-Particle Swarm Optimization Algorithm Applied to Multimodal MPPT
CN104238622A (en) Output power acquisition method based on photovoltaic array current characteristic under partial shadow
TWI391807B (en) A maximum power tracking system and method for photovoltaic power generation systems
CN103475267A (en) Method for improving generating efficiency of photovoltaic array
CN116800172B (en) Photoelectric conversion method and energy storage system
CN115759391A (en) A high-precision photovoltaic power generation power prediction analysis method and system
CN106951024B (en) A kind of quick MPPT method of solar components
CN118591983A (en) Method and apparatus for tracking the maximum power point of a photovoltaic module and monitoring the degradation of the photovoltaic module
Yang et al. Techno-economic-environmental optimization of hybrid photovoltaic-thermoelectric generator systems based on data-driven approach
CN114911300A (en) Photovoltaic system multi-pole tracking control method and device based on boundary positioning conductance increment control
CN106970677A (en) A kind of quick MPPT methods of solar inverter
CN104898758A (en) Photovoltaic array maximum power tracing apparatus and method
CN115169087A (en) Household photovoltaic power station scheme generation method and corresponding intelligent design platform
CN114995578A (en) Tracking method for maximum power of photovoltaic array under local shadow
Salem et al. Design of a perturb and observe and neural network algorithms-based maximum power point tracking for a gridconnected photovoltaic system.
CN107340799B (en) A kind of maximum power point of photovoltaic array tracking based on shadow Detection
Belhadj Djilali et al. Improved Incremental Conductance MPPT Technique Designed to Addressing Drift Problem in a Photovoltaic System
Cao Photovoltaic applications in demanding situations: estimation and optimisation of solar ressources for autonomous power supplies

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

Application publication date: 20220816

RJ01 Rejection of invention patent application after publication