CN107204731A - 基于几何平均数的风光互补发电系统mppt控制方法 - Google Patents
基于几何平均数的风光互补发电系统mppt控制方法 Download PDFInfo
- Publication number
- CN107204731A CN107204731A CN201710426487.3A CN201710426487A CN107204731A CN 107204731 A CN107204731 A CN 107204731A CN 201710426487 A CN201710426487 A CN 201710426487A CN 107204731 A CN107204731 A CN 107204731A
- Authority
- CN
- China
- Prior art keywords
- wind
- mrow
- output
- generating system
- power generation
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 21
- 230000005611 electricity Effects 0.000 claims 4
- 238000001514 detection method Methods 0.000 claims 1
- 238000010248 power generation Methods 0.000 abstract description 74
- 238000010276 construction Methods 0.000 abstract 1
- 238000006243 chemical reaction Methods 0.000 description 29
- 238000005516 engineering process Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02S—GENERATION OF ELECTRIC POWER BY CONVERSION OF INFRARED RADIATION, VISIBLE LIGHT OR ULTRAVIOLET LIGHT, e.g. USING PHOTOVOLTAIC [PV] MODULES
- H02S10/00—PV power plants; Combinations of PV energy systems with other systems for the generation of electric power
- H02S10/10—PV power plants; Combinations of PV energy systems with other systems for the generation of electric power including a supplementary source of electric power, e.g. hybrid diesel-PV energy systems
- H02S10/12—Hybrid wind-PV energy systems
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05F—SYSTEMS FOR REGULATING ELECTRIC OR MAGNETIC VARIABLES
- G05F1/00—Automatic 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/66—Regulating electric power
- G05F1/67—Regulating electric power to the maximum power available from a generator, e.g. from solar cell
-
- H02J3/385—
-
- H02J3/386—
-
- Y—GENERAL 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
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E10/00—Energy generation through renewable energy sources
- Y02E10/50—Photovoltaic [PV] energy
- Y02E10/56—Power conversion systems, e.g. maximum power point trackers
-
- Y—GENERAL 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
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E10/00—Energy generation through renewable energy sources
- Y02E10/70—Wind energy
- Y02E10/76—Power conversion electric or electronic aspects
Landscapes
- Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Sustainable Development (AREA)
- Sustainable Energy (AREA)
- Electromagnetism (AREA)
- Physics & Mathematics (AREA)
- Power Engineering (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Automation & Control Theory (AREA)
- Control Of Eletrric Generators (AREA)
- Control Of Electrical Variables (AREA)
- Wind Motors (AREA)
Abstract
基于几何平均数的风光互补发电系统MPPT控制方法,首先根据MPPT控制方法分别计算风力、光伏发电系统的占空比。然后用几何平均数求出这两个占空比的平均值,用此平均值调整风光互补发电系统的输出电压,进而可以调整其输出功率。从而实现用一个DC/DC直流变换电路同时跟踪风力和光伏发电系统输出功率的目的,结构简单,既可以降低成本又能实现MPPT功能。保证系统稳定输出,达到资源优化配置。
Description
技术领域
本发明涉及风光互补发电系统MPPT(最大功率跟踪)控制方法,特别是实现用一个DC/DC直流变换电路同时跟踪风力和光伏发电系统最大功率的方法。
背景技术
随着科技的高速发展,人类对于能源的需求量不断增加。因此,加速清洁能源的普及与推广,已成为世界各国走可持续发展道路的必然选择。新能源发电供能系统得到了快速发展。太阳能和风能在时间上的互补性使得风光互补发电系统在资源分布上具有很好的匹配性。
在初始阶段,风光互补发电系统只是简单的将风力发电系统和光伏发电系统进行了组合。在此结构中风力发电系统经过AC/DC整流电路和DC/DC直流变换电路后连接到直流母线上;光伏发电系统经DC/DC直流变换后连接到直流母线上。但是其成本较高。随着电力电子技术、智能控制的高速发展以及诸多学者对风光互补发电技术的深入研究,很多先进的MPPT控制方法已日臻完善。
由于风能、太阳能会受到风速、光照强度等自然条件变化的影响。在使用一个DC/DC直流变换电路的风光互补发电系统中,适宜的MPPT控制方法至关重要,既可以降低成本又能实现MPPT功能。保证系统稳定输出,达到资源优化配置。
发明内容
本发明的目的是用一个DC/DC直流变换电路同时跟踪风力和光伏发电系统最大功率的目的。
本发明是基于几何平均数的风光互补发电系统MPPT控制方法,风力发电与光伏发电两个子系统共用一个DC/DC直流变换电路,其步骤为:
步骤1:在系统正常运行的情况下,采集风力发电系统的输出电流iw、光伏发电系统的输出电流is、二者并联后的输出电压E、系统中DC/DC直流变换电路的输出电压U0。检测风机的转速ω。根据采集到风力和光伏发电系统的输出电压和输出电流计算输出功率变化量△Pw=Pk-Pk-1、△Ps=Pk-Pk-1,其中Pk、Pk-1表示第k秒和第k-1秒时的输出功率。计算风机的转速变化量△ω=ωk-ωk-1,其中ωk、ωk-1表示第k秒和第k-1秒时的风机转速。计算输出电压的变化量△U0=Uk-Uk-1,其中Uk、Uk-1表示第k秒和第k-1秒时系统中DC/DC直流变换电路的输出电压;
步骤2:由风力发电MPPT控制方法得到风力发电系统中DC/DC直流变换电路第k秒时的电压占空比Dw[k];同理,能够得到光伏发电系统中DC/DC直流变换电路第k秒时的电压占空比Ds[k];
步骤3:将步骤2计算所得值Dw[k]与Ds[k]输入占空比计算器中,用几何平均数计算得出Dw[k]与Ds[k]的几何平均值Dsw[k]:
步骤4:用步骤3计算出的Dsw[k]对风光互补发电系统的输出电压进行调整。其中输出电压为:
其中:U0表示系统中DC/DC直流变换电路的输出电压;E表示风力发电与光伏发电并联后的输出电压;Dsw[k]表示风力发电系统中DC/DC直流变换电路第k秒时的电压占空比Dw[k]与光伏发电系统中DC/DC直流变换电路第k秒时的电压占空比Ds[k]的几何平均值。
步骤5:检测并计算dPw/dω值的变化,若dPw/dω=0,即风力发电系统达到最大功率;至此就能够完成对风力发电系统最大功率的跟踪;检测并计算dPs/dU0值的变化,若dPs/dU0=0,即光伏发电系统达到最大功率,从而实现了用一个DC/DC直流变换电路同时跟踪风光互补发电系统的最大功率。
本发明的有益之处是:由于风光互补发电系统会受到风速、光照强度等自然条件的制约和影响,系统的输出功率会不稳定。适宜的MPPT控制方法可以使得风光互补发电系统在自然条件变化时,可以快速跟踪系统的输出功率。首先利用风力和光伏发电MPPT控制方法得到风力和光伏发电系统中的电压占空比,用几何平均数求得这两个占空比的平均值,然后用该平均值调整风光互补发电系统的输出电压,从而同时跟踪风光互补发电系统的输出功率。因此,实现用一个DC/DC直流变换电路同时跟踪风光互补发电系统的输出功率,结构简单,降低成本。
附图说明
图1是风光互补发电系统结构图,图2是基于几何平均数的风光互补发电系统MPPT控制方法原理框图。
具体实施方式
如图1所示,本发明所述的风光互补发电系统由风力发电子系统1、AC/DC整流电路2、光伏发电子系统4、DC/DC直流变换电路5、DC/AC逆变电路6、控制器7、负载8等组成。光伏发电子系统4由光伏电池板3组成。此系统的特点是风力发电系统1经过AC/DC整流电路2与光伏子系统4的输出端并联,然后由DC/DC变换器5升压后再经过DC/AC逆变电路6接入负载8。
如图2所示,本发明是基于几何平均数的风光互补发电系统MPPT控制方法,具体包括以下步骤:
步骤1:在系统正常运行的情况下,采集风力发电系统的输出电流iw、风机的转速ω、光伏发电系统的输出电流is、二者并联后的输出电压E、DC/DC直流变换电路的输出电压U0。根据采集到的输出电压和输出电流计算风力和光伏发电系统输出功率变化量△Pw、△Ps。依据采集到的风机的转速和输出电压得到转速的变化量△ω及输出电压的变化值△U0;
步骤2:由风力发电MPPT控制方法得到风力发电系统中DC/DC直流变换电路第k秒时的电压占空比Dw[k]。同理,由光伏发电MPPT控制方法得到Ds[k];
步骤3:用几何平均数求出这两个占空比的平均值Dsw[k]。其中
步骤4:用几何平均值Dsw[k]对风光互补发电系统的输出电压进行调整。其输出电压为:
其中:U0表示系统中DC/DC直流变换电路的输出电压;E表示风力发电与光伏发电并联后的输出电压;Dsw[k]表示风力发电系统中DC/DC直流变换电路第k秒时的电压占空比Dw[k]与光伏发电系统中DC/DC直流变换电路第k秒时的电压占空比Ds[k]的几何平均值。
由此可以达到用一个DC/DC直流变换电路同时跟踪风力和光伏发电系统最大功率的目的。
如图2所示,通过采集风光互补发电系统的输出电压和电流等参数,计算系统输出功率的变化,再利用MPPT控制方法得出风力和光伏发电系统的占空比,然后用几何平均数计算这两个占空比的平均值,用该平均值调整风光互补发电系统的输出电压,从而实现用一个DC/DC直流变换电路同时跟踪风光互补发电系统的最大功率。
具体实施步骤如下:
步骤1:在系统正常运行的情况下,采集风力发电系统的输出电流iw、光伏发电系统的输出电流is、二者并联后的输出电压E、系统中DC/DC直流变换电路的输出电压U0。检测风机的转速ω。根据采集到风力和光伏发电系统的输出电压和输出电流计算输出功率变化量△Pw=Pk-Pk-1、△Ps=Pk-Pk-1,其中Pk、Pk-1表示第k秒和第k-1秒时的输出功率。计算风机的转速变化量△ω=ωk-ωk-1,其中ωk、ωk-1表示第k秒和第k-1秒时的风机转速。计算输出电压的变化量△U0=Uk-Uk-1,其中Uk、Uk-1表示第k秒和第k-1秒时系统中DC/DC直流变换电路的输出电压;
步骤2:由风力发电MPPT控制方法得到风力发电系统中DC/DC直流变换电路第k秒时的电压占空比Dw[k];同理,能够得到光伏发电系统中DC/DC直流变换电路第k秒时的电压占空比Ds[k];
步骤3:用几何平均数求出这两个占空比的平均值Dsw[k]。
步骤4:用占空比平均值Dsw[k]对风光互补发电系统的输出电压进行调整。其输出电压为:
步骤5:检测并计算dPw/dω值的变化,若dPw/dω=0,即风力发电系统达到最大功率。检测并计算dPs/dU0值的变化,若dPs/dU0=0,即光伏发电系统达到最大功率。进而可以实现用一个DC/DC直流变换电路同时跟踪风力、光伏发电系统的最大功率。
由以上步骤可以实现用一个DC/DC直流变换电路分时跟踪风力、光伏发电系统的输出功率,使风光互补发电系统工作在最大功率点。其结构简单,既降低了成本又能实现最大功率跟踪(MPPT)功能,保证系统稳定输出,达到资源最优化配置。
Claims (1)
1.基于几何平均数的风光互补发电系统MPPT控制方法,风力发电与光伏发电两个子系统共用一个DC/DC直流变换电路,其特征在于,其步骤为:
步骤1:在系统正常运行的情况下,采集风力发电系统的输出电流iw、光伏发电系统的输出电流is、二者并联后的输出电压E、系统中DC/DC直流变换电路的输出电压U0。检测风机的转速ω。根据采集到的风力和光伏发电系统的输出电压和输出电流计算输出功率变化量△Pw=Pk-Pk-1、△Ps=Pk-Pk-1,其中Pk、Pk-1表示第k秒和第k-1秒时的输出功率。计算风机的转速变化量△ω=ωk-ωk-1,其中ωk、ωk-1表示第k秒和第k-1秒时的风机转速。计算输出电压的变化量△U0=Uk-Uk-1,其中Uk、Uk-1表示第k秒和第k-1秒时系统中DC/DC直流变换电路的输出电压;
步骤2:由风力发电MPPT控制方法得到风力发电系统中DC/DC直流变换电路的电压占空比Dw[k];同理,由光伏发电MPPT控制方法得到光伏发电系统中DC/DC直流变换电路的电压占空比Ds[k];
步骤3:将步骤2计算所得值Dw[k]与Ds[k]输入占空比计算器中,用几何平均数计算得出Dw[k]与Ds[k]的几何平均值Dsw[k]:
<mrow>
<msub>
<mi>D</mi>
<mrow>
<mi>s</mi>
<mi>w</mi>
</mrow>
</msub>
<mo>&lsqb;</mo>
<mi>k</mi>
<mo>&rsqb;</mo>
<mo>=</mo>
<msqrt>
<mrow>
<msub>
<mi>D</mi>
<mi>s</mi>
</msub>
<mo>&lsqb;</mo>
<mi>k</mi>
<mo>&rsqb;</mo>
<msub>
<mi>D</mi>
<mi>w</mi>
</msub>
<mo>&lsqb;</mo>
<mi>k</mi>
<mo>&rsqb;</mo>
</mrow>
</msqrt>
</mrow>
步骤4:用步骤3计算出的Dsw[k]对风光互补发电系统的输出电压进行调整。其中输出电压为:
<mrow>
<msub>
<mi>U</mi>
<mn>0</mn>
</msub>
<mo>=</mo>
<mfrac>
<mi>E</mi>
<mrow>
<mn>1</mn>
<mo>-</mo>
<msub>
<mi>D</mi>
<mrow>
<mi>s</mi>
<mi>k</mi>
</mrow>
</msub>
<mo>&lsqb;</mo>
<mi>k</mi>
<mo>&rsqb;</mo>
</mrow>
</mfrac>
</mrow>
其中:U0表示系统中DC/DC直流变换电路的输出电压;E表示风力发电与光伏发电并联后的输出电压;Dsw[k]表示风力发电系统中DC/DC直流变换电路第k秒时的电压占空比Dw[k]与光伏发电系统中DC/DC直流变换电路第k秒时的电压占空比Ds[k]的几何平均值。
步骤5:检测并计算dPw/dω值的变化,若dPw/dω=0,即风力发电系统达到最大功率;至此就能够完成对风力发电系统最大功率的跟踪;检测并计算dPs/dU0值的变化,若dPs/dU0=0,即光伏发电系统达到最大功率,从而实现了用一个DC/DC直流变换电路同时跟踪风光互补发电系统的最大功率。
由以上步骤可以实现用一个DC/DC直流变换电路同时跟踪风力、光伏发电系统的输出功率,使风光互补发电系统工作在最大功率点。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710426487.3A CN107204731B (zh) | 2017-06-08 | 2017-06-08 | 基于几何平均数的风光互补发电系统mppt控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710426487.3A CN107204731B (zh) | 2017-06-08 | 2017-06-08 | 基于几何平均数的风光互补发电系统mppt控制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107204731A true CN107204731A (zh) | 2017-09-26 |
CN107204731B CN107204731B (zh) | 2018-11-23 |
Family
ID=59907336
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710426487.3A Active CN107204731B (zh) | 2017-06-08 | 2017-06-08 | 基于几何平均数的风光互补发电系统mppt控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107204731B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109245092A (zh) * | 2018-09-28 | 2019-01-18 | 东北大学 | 基于优化模型预测算法结合变步长扰动控制的mppt方法 |
CN111969582A (zh) * | 2020-08-17 | 2020-11-20 | 博阳能源科技有限公司 | 基于下垂并机单功率输入接口兼容光伏与风电的供电系统 |
CN111969581A (zh) * | 2020-08-17 | 2020-11-20 | 博阳能源科技有限公司 | 一种单功率输入接口兼容光伏与风电的供电系统 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2012018155A1 (ko) * | 2010-08-03 | 2012-02-09 | (주)하이레벤 | 태양광 발전 시스템의 엠피피티 인버터의 제어 장치 |
CN106125817A (zh) * | 2016-08-30 | 2016-11-16 | 三门峡速达节能新能源科技研究院 | 一种基于光伏电池四参数模型的光伏mppt方法 |
US20170018933A1 (en) * | 2013-12-30 | 2017-01-19 | King Fahd University Of Petroleum And Minerals | Parking lot shade for generating electricity having a photovoltaic system that tracks a maximum power point |
CN106712094A (zh) * | 2017-02-05 | 2017-05-24 | 兰州理工大学技术工程学院 | 一种光伏并网逆变器 |
CN106873709A (zh) * | 2017-03-17 | 2017-06-20 | 浙江大学 | 一种快速跟踪光伏发电最大功率点的mppt方法 |
-
2017
- 2017-06-08 CN CN201710426487.3A patent/CN107204731B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2012018155A1 (ko) * | 2010-08-03 | 2012-02-09 | (주)하이레벤 | 태양광 발전 시스템의 엠피피티 인버터의 제어 장치 |
US20170018933A1 (en) * | 2013-12-30 | 2017-01-19 | King Fahd University Of Petroleum And Minerals | Parking lot shade for generating electricity having a photovoltaic system that tracks a maximum power point |
CN106125817A (zh) * | 2016-08-30 | 2016-11-16 | 三门峡速达节能新能源科技研究院 | 一种基于光伏电池四参数模型的光伏mppt方法 |
CN106712094A (zh) * | 2017-02-05 | 2017-05-24 | 兰州理工大学技术工程学院 | 一种光伏并网逆变器 |
CN106873709A (zh) * | 2017-03-17 | 2017-06-20 | 浙江大学 | 一种快速跟踪光伏发电最大功率点的mppt方法 |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109245092A (zh) * | 2018-09-28 | 2019-01-18 | 东北大学 | 基于优化模型预测算法结合变步长扰动控制的mppt方法 |
CN111969582A (zh) * | 2020-08-17 | 2020-11-20 | 博阳能源科技有限公司 | 基于下垂并机单功率输入接口兼容光伏与风电的供电系统 |
CN111969581A (zh) * | 2020-08-17 | 2020-11-20 | 博阳能源科技有限公司 | 一种单功率输入接口兼容光伏与风电的供电系统 |
CN111969581B (zh) * | 2020-08-17 | 2021-09-24 | 博阳能源科技有限公司 | 一种单功率输入接口兼容光伏与风电的供电系统 |
CN111969582B (zh) * | 2020-08-17 | 2021-09-24 | 博阳能源科技有限公司 | 基于下垂并机单功率输入接口兼容光伏与风电的供电系统 |
Also Published As
Publication number | Publication date |
---|---|
CN107204731B (zh) | 2018-11-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104579120A (zh) | 一种光伏发电离网运行时的控制方法 | |
CN109787298A (zh) | 一种变流器及其控制方法 | |
CN104135030B (zh) | 一种用于智能电网的柔性孤岛‑并网控制装置及方法 | |
CN102403737A (zh) | 一种风光互补发电储能控制方法及系统 | |
CN104270022A (zh) | 一种光伏并网逆变电路、开关控制电路及控制方法 | |
CN107612019A (zh) | 一种组串式光伏逆变器有功功率控制方法及系统 | |
CN105552952A (zh) | 一种光储联合发电系统及其能量管理方法 | |
CN102638059B (zh) | 单相光伏并网发电系统功率解耦电路及其控制方法 | |
CN107204731B (zh) | 基于几何平均数的风光互补发电系统mppt控制方法 | |
CN113193755A (zh) | 一种基于拓扑集成的多端口变换器、控制方法及系统 | |
CN103746548B (zh) | 一种用于分布式电源并网检测的多功能模拟电源 | |
CN103904692A (zh) | 风光互补离网、并网双模式系统 | |
CN104156029B (zh) | 一种基于扰动自适应的mppt控制方法 | |
CN103248106A (zh) | 基于单周期控制策略的新型电动汽车充电机 | |
Werulkar et al. | Design of a constant current solar charge controller with microcontroller based soft switching buck converter for solar home lighting system | |
CN107017664A (zh) | 一种风光互补发电系统最大功率分时跟踪方法 | |
CN110336302A (zh) | 具有虚拟同步特性的光储联合并网系统及其控制方法 | |
CN202405798U (zh) | 一种风光互补发电储能控制器 | |
CN108832817A (zh) | 一种减少燃料电池低频电流纹波的功率变换控制器及方法 | |
CN107069820B (zh) | 分布式可再生能源发电并网功率波动控制系统及控制方法 | |
CN110676883A (zh) | 一种新型数模结合控制的逆变并网装置 | |
CN101924376A (zh) | 太阳能发电并网控制系统 | |
CN106611962B (zh) | 一种链式级联自同步太阳能光伏发电控制装置及控制方法 | |
CN103683312A (zh) | 一种集成智能光伏发电组件及其方法 | |
CN204103793U (zh) | 一种光伏并网逆变电路及开关控制电路 |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
TR01 | Transfer of patent right | ||
TR01 | Transfer of patent right |
Effective date of registration: 20230802 Address after: 730500 Muchang village, Taoyang Town, Lintao County, Dingxi City, Gansu Province Patentee after: GANSU XINNONG ECOLOGY ENERGY ENVIRONMENTAL PROTECTION SCIENCE & TECHNOLOGY Co.,Ltd. Address before: 730050 langong Ping Road, Lanzhou, Gansu Province, No. 287 Patentee before: LANZHOU University OF TECHNOLOGY |