CN116044676A - Fan clearance detection method - Google Patents
Fan clearance detection method Download PDFInfo
- Publication number
- CN116044676A CN116044676A CN202211096211.0A CN202211096211A CN116044676A CN 116044676 A CN116044676 A CN 116044676A CN 202211096211 A CN202211096211 A CN 202211096211A CN 116044676 A CN116044676 A CN 116044676A
- Authority
- CN
- China
- Prior art keywords
- blade
- impeller
- headroom
- distance
- tower
- 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
Links
Images
Classifications
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F03—MACHINES OR ENGINES FOR LIQUIDS; WIND, SPRING, OR WEIGHT MOTORS; PRODUCING MECHANICAL POWER OR A REACTIVE PROPULSIVE THRUST, NOT OTHERWISE PROVIDED FOR
- F03D—WIND MOTORS
- F03D17/00—Monitoring or testing of wind motors, e.g. diagnostics
-
- 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/72—Wind turbines with rotation axis in wind direction
Landscapes
- Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Sustainable Development (AREA)
- Sustainable Energy (AREA)
- Chemical & Material Sciences (AREA)
- Combustion & Propulsion (AREA)
- Mechanical Engineering (AREA)
- General Engineering & Computer Science (AREA)
- Length Measuring Devices By Optical Means (AREA)
Abstract
本发明公开了一种风机净空检测方法,其步骤包括:(1)将至少一个拍摄设备对着风力发电机的叶轮的侧面;(2)转动所述叶轮,随着所述叶轮旋转形成的旋转平面,拍摄所述叶轮的各叶片转到最下端时的叶尖和塔架的图像;(3)测量各拍摄的图像中的所述叶片的叶尖到塔架的像素距离,计算各叶片的净空距离;还包括在步骤(3)之后,根据测量结果,比较各叶片的净空距离是否相等,从而确定所述叶轮是否平衡或检测各所述叶片的净空距离是否处于危险状态。本发明风机净空检测方法的检测过程简易、数据易计算,且检测用设备安装简单成本低。
The invention discloses a method for detecting the headroom of a wind turbine. The steps include: (1) pointing at least one photographing device at the side of the impeller of the wind power generator; plane, take the images of the blade tip and the tower when each blade of the impeller turns to the lowermost end; (3) measure the pixel distance from the blade tip of the blade in each captured image to the tower, and calculate the Clearance distance; also includes after step (3), according to the measurement results, comparing whether the clearance distances of each blade are equal, thereby determining whether the impeller is balanced or detecting whether the clearance distance of each blade is in a dangerous state. The detection process of the fan headroom detection method of the invention is simple, the data is easy to calculate, and the detection equipment is simple to install and low in cost.
Description
技术领域technical field
本发明涉及风机净空检测技术领域,尤指一种风机净空检测方法。The invention relates to the technical field of fan headroom detection, in particular to a fan headroom detection method.
背景技术Background technique
随着风力发电机单台功率的不断提高,风轮叶片也在不断往大型化方向发展。但是叶片越长越增加叶片扫塔的风险。叶片一旦发生扫塔,可能会造成叶片损坏,更甚则导致机组倒塔等恶劣事故发生。With the continuous improvement of the single power of wind turbines, the wind rotor blades are also constantly developing in the direction of large-scale. But the longer the blade, the more risk of blade sweeping. Once the blade sweeps the tower, it may cause damage to the blade, and even worse, it will cause serious accidents such as the unit falling down.
近年来通过视频或激光技术对叶片净空检测成为主流技术,叶片净空指的是叶片运行过程中,叶片距离塔筒的最小距离。但是这种主流技术需要将摄像机或激光发射器安装在机组的机舱或塔筒上进行测量,其安装工艺复杂、净空算法复杂,且成本高。In recent years, the detection of blade clearance through video or laser technology has become a mainstream technology. Blade clearance refers to the minimum distance between the blade and the tower during the operation of the blade. However, this mainstream technology requires cameras or laser transmitters to be installed on the cabin or tower of the unit for measurement. The installation process is complicated, the clearance algorithm is complicated, and the cost is high.
因此亟需一种检测过程简易、数据易计算,且检测用设备安装简单成本低的风机净空检测方法。Therefore, there is an urgent need for a fan headroom detection method with simple detection process, easy calculation of data, simple installation of detection equipment and low cost.
发明内容Contents of the invention
本发明的目的在于提供一种风机净空检测方法,检测过程简易、数据易计算,且检测用设备安装简单成本低。The object of the present invention is to provide a method for detecting the headroom of a blower fan, the detection process is simple, the data is easy to calculate, and the detection equipment is easy to install and low in cost.
为了实现上述目的,本发明的技术解决方案为:一种风机净空检测方法,其包括如下步骤:In order to achieve the above object, the technical solution of the present invention is: a fan headroom detection method, which includes the following steps:
(1)将至少一个拍摄设备对着风力发电机的叶轮的侧面;(1) Point at least one photographing device to the side of the wind turbine's impeller;
(2)转动所述叶轮,随着所述叶轮旋转形成的旋转平面,拍摄所述叶轮的各叶片转到最下端时的叶尖和塔架的图像;(2) Rotate the impeller, along with the plane of rotation formed by the rotation of the impeller, take images of the blade tip and the tower when each blade of the impeller turns to the lowermost end;
(3)测量各拍摄的图像中的所述叶片的叶尖到塔架的像素距离,计算各叶片的净空距离。(3) Measure the pixel distance from the tip of the blade to the tower in each captured image, and calculate the headroom distance of each blade.
优选地,还包括在步骤(3)之后,根据测量结果,比较各叶片的净空距离是否相等,从而确定所述叶轮是否平衡或检测各所述叶片的净空距离是否处于危险状态。Preferably, after step (3), comparing whether the clearance distances of the blades are equal according to the measurement results, so as to determine whether the impeller is balanced or detect whether the clearance distances of the blades are in a dangerous state.
优选地,所述步骤(1)中的相机位置设置为拍摄方向垂直于所述叶轮的旋转轴心线。Preferably, the position of the camera in the step (1) is set such that the shooting direction is perpendicular to the axis of rotation of the impeller.
优选地,所述步骤(2)中的旋转平面是由三个所述叶片与轮毂组成的所述叶轮绕所述轮毂的中心线旋转形成的平面。Preferably, the rotation plane in the step (2) is a plane formed by the rotation of the impeller composed of three blades and the hub around the center line of the hub.
优选地,所述步骤(3)中计算得到的净空距离C=c×L/l,其中c为测量出的所述叶尖至所述塔架的像素距离,L为与所述叶尖在同一水平面的所述塔架的横截面直径,l为所述塔架对应段的像素距离。Preferably, the clearance distance calculated in the step (3) is C=c×L/l, wherein c is the measured pixel distance from the blade tip to the tower, and L is the distance between the blade tip and the blade tip. The cross-sectional diameter of the tower at the same horizontal plane, l is the pixel distance of the corresponding section of the tower.
优选地,在拍摄的所述图像中,所述叶片处于最下端时是所述拍摄设备的拍摄线与所述叶片运动轨迹相切的时刻,因此需要修正计算得出的所述净空距离时所述叶片处于净空最小时刻的位置,修正距离为:△r1=R×(1-cosα),R是叶轮半径,α是所述相机拍摄方向与水平面之间的夹角。Preferably, in the captured image, when the blade is at the lowest end is the moment when the shooting line of the shooting device is tangent to the motion trajectory of the blade, so it is necessary to correct the calculated clearance distance. The position where the blade is at the moment of minimum headroom, the corrected distance is: Δr1=R×(1-cosα), R is the radius of the impeller, and α is the angle between the shooting direction of the camera and the horizontal plane.
优选地,当所述叶片的运行轨迹在图像中显示为椭圆形,则说明所述叶轮发生了偏航,所述叶轮偏航的角度β=arcsin(b/R),其中所述叶轮的半径为R,所述椭圆形轨迹的短轴长度为2b。Preferably, when the running track of the blade is shown as an ellipse in the image, it indicates that the impeller has yaw, and the yaw angle of the impeller is β=arcsin(b/R), wherein the radius of the impeller is R, and the minor axis length of the elliptical trajectory is 2b.
优选地,所述叶轮发生偏航时的真实净空距离C=C0/COSβ+L/2(1/COSβ-1),其中,L/2为与所述叶尖在同一水平面的所述塔架的横截面半径,C0为步骤(3)计算得到的所述净空距离。Preferably, the true clearance distance C=C0/COSβ+L/2(1/COSβ-1) when the impeller yaws, wherein, L/2 is the tower at the same level as the blade tip The radius of the cross section, C0 is the clearance distance calculated in step (3).
优选地,所述拍摄设备距离所述塔架的距离≤1000m。Preferably, the distance between the shooting device and the tower is ≤1000m.
优选地,所述拍摄设备距离所述塔架的距离为100m-300m。Preferably, the distance between the shooting device and the tower is 100m-300m.
优选地,所述拍摄设备采用相机或激光设备。Preferably, the photographing device adopts a camera or a laser device.
采用上述方案后,本发明风机净空检测方法具有以下有益效果:After adopting the above scheme, the air blower headroom detection method of the present invention has the following beneficial effects:
1、本发明通过在塔架下或底面设置相机或激光设备等的拍摄设备,从叶轮侧面拍摄各叶片的叶尖到塔架的净空距离,由各叶片的净空距离可判断是否处于危险状态,也可以通过比较三个叶片的净空距离是否相等,判断叶轮是否平衡,其检测过程简易、数据易计算,且检测用设备安装简单成本低;1, the present invention is by setting the shooting equipment such as camera or laser device under the tower frame or bottom surface, photographs the tip of each blade from the side of the impeller to the headroom distance of the tower frame, can judge whether to be in a dangerous state by the headroom distance of each blade, It is also possible to judge whether the impeller is balanced by comparing whether the clearance distances of the three blades are equal. The detection process is simple, the data is easy to calculate, and the detection equipment is simple to install and low in cost;
2、通过修正所计算的净空距离时叶片处于净空最小时刻的位置,可以更准确的得到真实的净空距离,从而更好的判断其是否处于危险状态,以及叶轮是否平衡;2. By correcting the position of the blade at the moment of the minimum clearance when the calculated clearance distance is calculated, the real clearance distance can be obtained more accurately, so as to better judge whether it is in a dangerous state and whether the impeller is balanced;
3、通过拍摄的叶片运行轨迹可以确定叶轮是否偏航,通过测得的偏航角采用简单的算法实现机组偏航状态下的真实的净空距离的检测。3. Whether the impeller is yaw can be determined through the captured blade running track, and the real headroom distance detection under the yaw state of the unit can be realized by using a simple algorithm through the measured yaw angle.
附图说明Description of drawings
图1是本发明计算叶片净空距离的拍摄设备拍摄方向的前视示意图;Fig. 1 is the front view schematic diagram of the shooting direction of the shooting device for calculating the blade clearance distance in the present invention;
图2是本发明计算叶片净空距离的拍摄设备拍摄方向的左视示意图;Fig. 2 is a schematic left view of the shooting direction of the shooting device for calculating the blade clearance distance in the present invention;
图3是本发明实际叶片净空距离的拍摄设备拍摄方向的前视示意图;Fig. 3 is the front view schematic diagram of the shooting direction of the shooting equipment of the actual blade clearance distance of the present invention;
图4是本发明实际叶片净空距离的拍摄设备拍摄方向的左视示意图;Fig. 4 is a schematic left view of the shooting direction of the shooting device of the actual blade clearance distance of the present invention;
图5是图4所示发生偏航时的偏航角算法示意图;Fig. 5 is a schematic diagram of the yaw angle algorithm when yaw occurs as shown in Fig. 4;
图6是本发明在叶轮偏航时的修正算法示意图。Fig. 6 is a schematic diagram of the correction algorithm of the present invention when the impeller is yaw.
下面结合附图,通过实施例对本发明做进一步的说明;Below in conjunction with accompanying drawing, the present invention will be further described by embodiment;
具体实施方式Detailed ways
下面根据附图所示实施方式阐述本发明。此次公开的实施方式可以认为在所有方面均为例示,不具限制性。本发明的范围不受以下实施方式的说明所限,仅由权利要求书的范围所示,而且包括与权利要求范围具有同样意思及权利要求范围内的所有变形。The present invention will be described below according to the embodiments shown in the accompanying drawings. It can be thought that embodiment disclosed this time is an illustration in every point, and is not restrictive. The scope of the present invention is not limited by the description of the following embodiments but only by the scope of the claims, and includes the same meaning as the scope of the claims and all modifications within the scope of the claims.
参考图1和图2所示,本发明风机净空检测方法包括如下步骤:With reference to Fig. 1 and shown in Fig. 2, fan headroom detection method of the present invention comprises the following steps:
(1)将至少一个拍摄设备对着风力发电机的叶轮的侧面进行拍摄,该拍摄设备可以采用相机或激光设备。如图2所示,本实施例采用相机1,将相机1对着叶轮2的侧面放置,本实施例的相机1拍摄方向为朝着叶轮2的一个叶片3在最下端时的方向,该叶轮2是由三个叶片3与轮毂4连接组成,并一体安装于塔架5的上端,三个叶片3安装于轮毂4的轴心方向,且三个叶片3沿圆周方向均布,相机1设置为拍摄方向垂直于叶轮2的旋转轴心线,相机1距离塔架5的距离选择≤1000m,优选距离为100m-300m;(1) Shooting at least one shooting device against the side of the impeller of the wind turbine, the shooting device may be a camera or a laser device. As shown in Figure 2, the present embodiment adopts the
(2)转动叶轮2,三个叶片3以轮毂4的中心线旋转形成旋转平面A,由相机1拍摄叶轮2的各叶片3转到最下端时的叶尖6和塔架5的图像,即沿拍摄线7拍摄叶尖6和塔架5;(2) Turn the
(3)测量各拍摄的图像中的叶片的叶尖6位于最下端时到塔架5的像素距离,计算各叶片3的净空距离。净空距离计算公式为:C=c×L/l,其中C为净空距离,c为测量出的叶尖6至塔架5的像素距离,L为与叶尖6在同一水平面的塔架5的横截面直径,l为塔架5对应段的像素距离。(3) Measure the pixel distance from the
通过步骤(3)比较各叶片3的净空距离C,可确定叶轮2是否平衡或检测各叶片3的净空距离C是否处于危险状态,当三个叶片3的净空距离C相等,则说明叶轮2平衡,否则为不平衡。By comparing the clearance distance C of each
参考图3所示,在拍摄的图像中,由于叶片3处于最下端时是相机1的拍摄线与叶片3运动轨迹相切的时刻,而并非叶片3在转动过程中处于净空的最小时。因此需要修正计算净空距离时叶片处于净空最小时的位置,修改后的净空距离为:As shown in FIG. 3 , in the photographed image, when the
△r1=R×(1-cosα),其中R是叶轮2的半径,α是相机1的拍摄方向与水平面之间的夹角。Δr1=R×(1-cosα), where R is the radius of the
参考图4所示,当叶片3的运行轨迹在图像中显示为椭圆形B时,则说明叶轮2发生了偏航,参考图5所示,叶轮2偏航的角度β=arcsin(b/R),其中R是叶轮2的半径,b是椭圆形轨迹的短半轴的长度。参考图6所示,在叶轮2发生偏航时的各叶片3真实的净空距离C=C0/COSβ+L/2(1/COSβ-1),其中C0为由步骤(3)计算得到的净空距离,L/2为与叶尖6在同一水平面的塔架5的横截面半径。该公式由以下计算步骤得出:Referring to Figure 4, when the running track of the
C=C1+C2(1)C=C1+C2(1)
C1=C0/COSβ(2)C1=C0/COSβ(2)
C2=L/2COSβ-L/2(3)C2=L/2COSβ-L/2(3)
由(1)、(2)、(3)可以得出C=C0/COSβ+L/2(1/COSβ-1)。From (1), (2) and (3), it can be obtained that C=C0/COSβ+L/2(1/COSβ-1).
本发明通过在塔架5下或底面设置相机1,也可以设置激光设备等拍摄设备,从叶轮2的左侧面拍摄各叶片3的叶尖6到塔架5的净空距离,由各叶片3的净空距离可判断是否处于危险状态,也可以通过比较三个叶片3的净空距离是否相等,判断叶轮2是否平衡,其检测过程简易、数据易计算,且检测用设备安装简单成本低;通过修正所计算的净空距离时叶片处于净空最小时刻的位置,可以更准确的得到真实的净空距离,从而更好的判断其是否处于危险状态,以及叶轮是否平衡;通过拍摄的叶片3运行轨迹可以确定叶轮是否偏航,通过测得的偏航角采用简单的算法实现机组偏航状态下的真实的净空距离的检测。The present invention can also be provided with photographing equipment such as laser equipment by setting
以上所述实施例仅仅是对本发明的优选实施方式进行描述,并非对本发明的范围进行限定,在不脱离本发明设计精神的前提下,本领域普通工程技术人员对本发明的技术方案作出的各种变形和改进,均应落入本发明的权利要求书确定的保护范围内。The above-mentioned embodiments are only descriptions of the preferred implementation modes of the present invention, and are not intended to limit the scope of the present invention. Variations and improvements should fall within the scope of protection defined by the claims of the present invention.
Claims (11)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211096211.0A CN116044676A (en) | 2022-09-08 | 2022-09-08 | Fan clearance detection method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211096211.0A CN116044676A (en) | 2022-09-08 | 2022-09-08 | Fan clearance detection method |
Publications (1)
Publication Number | Publication Date |
---|---|
CN116044676A true CN116044676A (en) | 2023-05-02 |
Family
ID=86117171
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202211096211.0A Pending CN116044676A (en) | 2022-09-08 | 2022-09-08 | Fan clearance detection method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116044676A (en) |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109958583A (en) * | 2017-12-22 | 2019-07-02 | 北京金风科创风电设备有限公司 | A method and device for measuring the tower clearance of a wind turbine |
US20210003114A1 (en) * | 2018-11-30 | 2021-01-07 | Beijing Goldwind Science & Creation Windpower Equipment Co., Ltd. | Method and device for determining tower clearance for wind turbine |
KR102322693B1 (en) * | 2020-05-14 | 2021-11-08 | 한국에너지기술연구원 | Apparatus for monitoring wind turbine blade and method thereof |
CN113864132A (en) * | 2020-06-30 | 2021-12-31 | 乌鲁木齐金风天翼风电有限公司 | Tower clearance monitoring method, device and equipment |
CN114623049A (en) * | 2022-01-29 | 2022-06-14 | 西安翔迅科技有限责任公司 | Wind turbine generator tower clearance monitoring method and computer program product |
CN114810514A (en) * | 2022-06-27 | 2022-07-29 | 东方电气风电股份有限公司 | Wind turbine generator blade tip clearance distance measuring method, device, equipment and storage medium |
-
2022
- 2022-09-08 CN CN202211096211.0A patent/CN116044676A/en active Pending
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109958583A (en) * | 2017-12-22 | 2019-07-02 | 北京金风科创风电设备有限公司 | A method and device for measuring the tower clearance of a wind turbine |
US20210003114A1 (en) * | 2018-11-30 | 2021-01-07 | Beijing Goldwind Science & Creation Windpower Equipment Co., Ltd. | Method and device for determining tower clearance for wind turbine |
KR102322693B1 (en) * | 2020-05-14 | 2021-11-08 | 한국에너지기술연구원 | Apparatus for monitoring wind turbine blade and method thereof |
CN113864132A (en) * | 2020-06-30 | 2021-12-31 | 乌鲁木齐金风天翼风电有限公司 | Tower clearance monitoring method, device and equipment |
CN114623049A (en) * | 2022-01-29 | 2022-06-14 | 西安翔迅科技有限责任公司 | Wind turbine generator tower clearance monitoring method and computer program product |
CN114810514A (en) * | 2022-06-27 | 2022-07-29 | 东方电气风电股份有限公司 | Wind turbine generator blade tip clearance distance measuring method, device, equipment and storage medium |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US20110206511A1 (en) | Wind turbine and method for measuring the pitch angle of a wind turbine rotor blade | |
US7870783B2 (en) | Method for measuring the turbulence intensity of a horizontal axis wind turbine | |
EP3645866B1 (en) | Computer-implemented method for re-calibrating at least one yaw-angle of a wind turbine, respective system, computer-implemented method for wind park optimization, and respective wind park | |
US10378517B2 (en) | Method for replacing the blades of a wind turbine to maintain safe operation | |
US20150252789A1 (en) | Method for Detecting Deflection of the Blades of a Wind Turbine | |
CN110520620B (en) | Method for reducing loads on a wind turbine blade pitch bearing | |
CN108869197B (en) | Method and system for accurately measuring height of fan through unmanned aerial vehicle | |
JP2010229824A (en) | Rotating blade diagnostic device | |
CN205135913U (en) | String of a musical instrument locating template and setting angle of blade measurement system of wind turbine generator system blade | |
WO2020108598A1 (en) | Method, device and system for correcting installation errors of wind vane | |
WO2011051778A1 (en) | System for determining wind turbine blade pitch settings | |
WO2022117616A1 (en) | Unmanned airborne visual diagnosis of an operating wind turbine generator | |
CN109973330A (en) | A method for detecting the influence of the wake of an upstream fan on a downstream fan | |
CN111102940B (en) | Blade pitch angle deviation detection method, device, storage medium and system | |
JP6147335B2 (en) | Method for controlling the pitch angle of at least one wind turbine blade | |
CN110857685B (en) | Method for detecting deviation of installation angle of blades of wind turbine generator | |
CN116044676A (en) | Fan clearance detection method | |
KR102322693B1 (en) | Apparatus for monitoring wind turbine blade and method thereof | |
WO2022174566A1 (en) | Method and system for identifying variable pitch out-of-synchronization failure of wind turbine generator system on basis of image processing | |
CN114127413B (en) | Method and apparatus for computer-implemented monitoring of one or more wind turbines in a wind farm | |
CN117889786B (en) | Vision-based automatic inspection method and system for wind power station | |
CN113565697A (en) | Impeller pneumatic unbalance optimization system and method based on laser and video measurement | |
CN115713558B (en) | A method for measuring wind turbine clearance based on ground vision | |
CN113864132B (en) | Tower clearance monitoring method, device and equipment | |
CN106704099A (en) | Method and device for controlling wind turbine |
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 | ||
CB02 | Change of applicant information | ||
CB02 | Change of applicant information |
Country or region after: China Address after: Room 181, 1st Floor, Building 12, Nanhaijiayuan Erli, Yinghai Town, Beijing Economic and Technological Development Zone, Daxing District, Beijing 100176 Applicant after: BEIJING ONE.WIND NEW ENERGY DEVELOPMENT CO.,LTD. Address before: 100176 room 2505, 20th floor, building 1, courtyard 7, Ronghua South Road, economic and Technological Development Zone, Daxing District, Beijing Applicant before: BEIJING ONE.WIND NEW ENERGY DEVELOPMENT CO.,LTD. Country or region before: China |