CN104864560A - Air conditioner electricity consumption pre-estimating method in office building - Google Patents

Air conditioner electricity consumption pre-estimating method in office building Download PDF

Info

Publication number
CN104864560A
CN104864560A CN201510225624.8A CN201510225624A CN104864560A CN 104864560 A CN104864560 A CN 104864560A CN 201510225624 A CN201510225624 A CN 201510225624A CN 104864560 A CN104864560 A CN 104864560A
Authority
CN
China
Prior art keywords
building
condition
full curve
temperature
curve sections
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
Application number
CN201510225624.8A
Other languages
Chinese (zh)
Other versions
CN104864560B (en
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.)
Suzhou Chint Enterprise Development Co.,Ltd.
Original Assignee
SHANGHAI SUNRISE POWER TECHNOLOGY 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 SHANGHAI SUNRISE POWER TECHNOLOGY Co Ltd filed Critical SHANGHAI SUNRISE POWER TECHNOLOGY Co Ltd
Priority to CN201510225624.8A priority Critical patent/CN104864560B/en
Publication of CN104864560A publication Critical patent/CN104864560A/en
Application granted granted Critical
Publication of CN104864560B publication Critical patent/CN104864560B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00Administration; Management
    • G06Q10/04Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"

Abstract

An air conditioner electricity consumption pre-estimating method in an office building relates to the energy-saving technical field and aims to solve the technical problem of reducing energy waste and saving electricity consumption cost. The method is that an air conditioner load curve, an indoor temperature curve and an outdoor temperature curve of the most three recent years in the office building are firstly acquired; secondly, according to the active power, the air conditioner rated power, the indoor temperature, an air conditioner preference temperature in the building and the number of people in the building, all continuous curve sections that meet three conditions are found out from the air conditioner load curve; thirdly, the fixed load value, the temperature coefficient, the load value of per capital of the whole minute point, the load coefficient of per capital are calculated according to those curve sections; and finally, the air conditioner electricity consumption load pre-estimating value is calculated. Thus, the method provided herein is suitable for air conditioners in office buildings.

Description

Office building air conditioning electricity predictor method
Technical field
The present invention relates to power-saving technology, particularly relate to a kind of technology of office building air conditioning electricity predictor method.
Background technology
In order to make full use of the green energy resource such as wind energy, solar energy, the non-renewable energy resources such as few thermoelectricity as much as possible, the power supply of every day is divided into three periods according to electricity consumption situation by power system, be respectively electrical network peak period, electrical network section, the electric-net valley period at ordinary times, in the electrical network peak period (such as 9 o'clock to 17 o'clock) that power supply requirement is higher, green energy resource proportion is lower, in the electrical network section at ordinary times that power supply requirement is relatively slightly low, then green energy resource proportion is then corresponding slightly high, and in power supply requirement minimum electric-net valley period, green energy resource proportion also reaches the highest.In order to encourage enterprise to utilize green energy resource, the electricity price of different power supply period is also different, and peak valley ordinary telegram valency is as one of electric energy measures to regulate rush-hour traffic, effectively can suppress peak value, improve peak-valley difference, avoid power consumption equipment capacity investment waste, also effectively make use of green energy resource simultaneously.
Office building air-conditioning is energy consumption rich and influential family, takes airconditioning control measure in advance, just can participate in the flat measures to regulate rush-hour traffic of electrical network peak valley, can cut operating costs, reduce requirement expense expenditure, the method of operation of all right elaborate scheme air-conditioning, reduces energy waste, saves electric cost.But also there is no effective method for taking airconditioning control measure in advance.
Summary of the invention
For the defect existed in above-mentioned prior art, technical problem to be solved by this invention is to provide one can reduce energy waste, saves the office building air conditioning electricity predictor method of electric cost.
In order to solve the problems of the technologies described above, a kind of office building air conditioning electricity predictor method provided by the present invention, it is characterized in that, concrete steps are as follows:
1) office building is obtained at air conditioner load curve, indoor temperature curve, the outdoor temperature curve of nearest 3 years;
2) according to indoor temperature curve, the outdoor temperature curve of nearest 3 years, the air conditioner load curve of nearest 3 years is searched all full curve sections of the A that satisfies condition, condition B, condition C;
Condition A: active power is greater than the S0 of air-conditioning rated power, and the variable quantity of active power is less than the K0 of air-conditioning rated power, and the difference of indoor temperature and outdoor temperature is less than W0, and the difference of indoor temperature and building air conditioning preferred temperature is less than W0, and the duration of this full curve section is greater than T0, and in building, personnel's number is less than the K0 of building setting personnel number;
Condition B: active power is greater than the S1 of air-conditioning rated power, and the variable quantity of active power is less than the K1 of air-conditioning rated power, and the difference of indoor temperature and building air conditioning preferred temperature is less than W0, and indoor temperature is less than outdoor temperature W1 degree, and the duration of this full curve section is greater than T1, and in building, personnel's number is less than the K1 of building setting personnel number;
Condition C: active power is greater than the S2 of air-conditioning rated power, and the variable quantity of active power is less than the K1 of air-conditioning rated power, and the difference of indoor temperature and building air conditioning preferred temperature is less than W0, and indoor temperature is less than outdoor temperature W2 degree, and the duration of this full curve section is greater than T2, and in building, personnel's number is greater than the K2 of building setting personnel number;
Wherein, the representative value of S0 is 30%, the representative value of the representative value of S1 to be the representative value of 50%, S2 be 80%, W0 is 1 degree Celsius, the representative value of W1 and W2 is 4 degrees Celsius, the representative value of T0 is 4 hours, and the representative value of T1 and T2 is 4 hours, and the representative value of K0 and K1 is 5%, the representative value of K2 is 60%, and the representative value in summer of building air conditioning preferred temperature is 26 degrees Celsius;
3) according to all full curve sections of the A that satisfies condition, calculate firm demand value, specific formula for calculation is:
P0=E0/D0
Wherein, P0 is firm demand value, and E0 is the meritorious integration of all full curve sections of A of satisfying condition, and D0 is the total duration of all full curve sections of A of satisfying condition;
4) according to all full curve sections of the B that satisfies condition, calculate temperature coefficient, specific formula for calculation is:
a1=(E1-P0×D1)/D1/H1
H1=∑△Tb/60/D1
Wherein, a1 is temperature coefficient, E1 is the meritorious integration of all full curve sections of B of satisfying condition, D1 is the total duration of all full curve sections of B of satisfying condition, H1 is the indoor/outdoor temperature-difference hourly average value of all full curve sections of B of satisfying condition, and ∑ △ Tb is the summation of indoor and outdoor temperature difference per minute of all full curve sections of B of satisfying condition;
5) to all full curve sections of the C that satisfies condition, calculate each whole point of hour per capita environment loads value, specific formula for calculation is:
A2[i]=(Px[i]-a1×△Tc[i]-P0)/m[i]
Wherein, i-th whole point hour per capita environment loads value in all full curve sections that A2 [i] is the C that satisfies condition, the burden with power at i-th whole point of hour in all full curve sections that Px [i] is the C that satisfies condition, the indoor and outdoor temperature at i-th whole point of hour in all full curve sections that △ Tc [i] is the C that satisfies condition is poor, personnel's number in the building at i-th whole point of hour in all full curve sections that m [i] is the C that satisfies condition;
6) according to all full curve sections of the C that satisfies condition, calculate per capita environment loads coefficient, specific formula for calculation is:
a2=∑△A2/n
Wherein, a2 is load coefficient per capita, and ∑ △ A2 is the summation of the whole point of hour per capita environment loads value satisfied condition in all full curve sections of C, and n is the number of the whole point of hour per capita environment loads value satisfied condition in all full curve sections of C;
7) calculate the discreet value of air conditioning electricity load, specific formula for calculation is:
P=P0+a1×△T+a2×m
Wherein, P is the discreet value of air conditioning electricity load, and △ T is the current indoor/outdoor temperature-difference of office building, and m is personnel's number in current building.
Office building air conditioning electricity predictor method provided by the invention, according to air conditioner load curve, temperature curve, building air conditioning preferred temperature, the interior personnel's number of building of history, calculate the discreet value of air conditioning electricity load, can the method for operation of elaborate scheme air-conditioning by this discreet value, reduce energy waste, save electric cost.
Accompanying drawing explanation
Fig. 1 is the calculation flow chart of the office building air conditioning electricity predictor method of the embodiment of the present invention.
Detailed description of the invention
Illustrate below in conjunction with accompanying drawing and embodiments of the invention are described in further detail; but the present embodiment is not limited to the present invention; every employing analog structure of the present invention and similar change thereof, all should list protection scope of the present invention in, the pause mark in the present invention all represent and relation.
As shown in Figure 1, a kind of office building air conditioning electricity predictor method that the embodiment of the present invention provides, it is characterized in that, concrete steps are as follows:
1) office building is obtained at air conditioner load curve, indoor temperature curve, the outdoor temperature curve of nearest 3 years;
2) according to indoor temperature curve, the outdoor temperature curve of nearest 3 years, the air conditioner load curve of nearest 3 years is searched all full curve sections of the A that satisfies condition, condition B, condition C;
Condition A: active power is greater than the S0 of air-conditioning rated power, and the variable quantity of active power is less than the K0 of air-conditioning rated power, and the difference of indoor temperature and outdoor temperature is less than W0, and the difference of indoor temperature and building air conditioning preferred temperature is less than W0, and the duration of this full curve section is greater than T0, and in building, personnel's number is less than the K0 of building setting personnel number;
Condition B: active power is greater than the S1 of air-conditioning rated power, and the variable quantity of active power is less than the K1 of air-conditioning rated power, and the difference of indoor temperature and building air conditioning preferred temperature is less than W0, and indoor temperature is less than outdoor temperature W1 degree, and the duration of this full curve section is greater than T1, and in building, personnel's number is less than the K1 of building setting personnel number;
Condition C: active power is greater than the S2 of air-conditioning rated power, and the variable quantity of active power is less than the K1 of air-conditioning rated power, and the difference of indoor temperature and building air conditioning preferred temperature is less than W0, and indoor temperature is less than outdoor temperature W2 degree, and the duration of this full curve section is greater than T2, and in building, personnel's number is greater than the K2 of building setting personnel number;
Wherein, the representative value of S0 is 30%, the representative value of the representative value of S1 to be the representative value of 50%, S2 be 80%, W0 is 1 degree Celsius, the representative value of W1 and W2 is 4 degrees Celsius, the representative value of T0 is 4 hours, and the representative value of T1 and T2 is 4 hours, and the representative value of K0 and K1 is 5%, the representative value of K2 is 60%, and the representative value in summer of building air conditioning preferred temperature is 26 degrees Celsius;
Wherein, in air-conditioner temperature load curve, building setting personnel number, building air conditioning preferred temperature, building, personnel's number all can obtain from existing building Property Management System;
3) according to all full curve sections of the A that satisfies condition, calculate firm demand value, specific formula for calculation is:
P0=E0/D0
Wherein, P0 is firm demand value, and E0 is the meritorious integration of all full curve sections of A of satisfying condition, and D0 is the total duration of all full curve sections of A of satisfying condition;
4) according to all full curve sections of the B that satisfies condition, calculate temperature coefficient, specific formula for calculation is:
a1=(E1-P0×D1)/D1/H1
H1=∑△Tb/60/D1
Wherein, a1 is temperature coefficient, E1 is the meritorious integration of all full curve sections of B of satisfying condition, D1 is the total duration of all full curve sections of B of satisfying condition, H1 is the indoor/outdoor temperature-difference hourly average value of all full curve sections of B of satisfying condition, and ∑ △ Tb is the summation of indoor and outdoor temperature difference per minute of all full curve sections of B of satisfying condition;
5) to all full curve sections of the C that satisfies condition, calculate each whole point of hour per capita environment loads value, specific formula for calculation is:
A2[i]=(Px[i]-a1×△Tc[i]-P0)/m[i]
Wherein, i-th whole point hour per capita environment loads value in all full curve sections that A2 [i] is the C that satisfies condition, the burden with power at i-th whole point of hour in all full curve sections that Px [i] is the C that satisfies condition, the indoor and outdoor temperature at i-th whole point of hour in all full curve sections that △ Tc [i] is the C that satisfies condition is poor, personnel's number in the building at i-th whole point of hour in all full curve sections that m [i] is the C that satisfies condition;
6) according to all full curve sections of the C that satisfies condition, calculate per capita environment loads coefficient, specific formula for calculation is:
a2=∑△A2/n
Wherein, a2 is load coefficient per capita, ∑ △ A2 is the summation (i.e. the summation of A2 [i]) of the whole point of hour per capita environment loads value satisfied condition in all full curve sections of C, and n is the number (i.e. the quantity of A2 [i]) of the whole point of hour per capita environment loads value satisfied condition in all full curve sections of C;
7) calculate the discreet value of air conditioning electricity load, specific formula for calculation is:
P=P0+a1×△T+a2×m
Wherein, P is the discreet value of air conditioning electricity load, and △ T is the current indoor/outdoor temperature-difference of office building, and m is personnel's number in current building.

Claims (1)

1. an office building air conditioning electricity predictor method, is characterized in that, concrete steps are as follows:
1) office building is obtained at air conditioner load curve, indoor temperature curve, the outdoor temperature curve of nearest 3 years;
2) according to indoor temperature curve, the outdoor temperature curve of nearest 3 years, the air conditioner load curve of nearest 3 years is searched all full curve sections of the A that satisfies condition, condition B, condition C;
Condition A: active power is greater than the S0 of air-conditioning rated power, and the variable quantity of active power is less than the K0 of air-conditioning rated power, and the difference of indoor temperature and outdoor temperature is less than W0, and the difference of indoor temperature and building air conditioning preferred temperature is less than W0, and the duration of this full curve section is greater than T0, and in building, personnel's number is less than the K0 of building setting personnel number;
Condition B: active power is greater than the S1 of air-conditioning rated power, and the variable quantity of active power is less than the K1 of air-conditioning rated power, and the difference of indoor temperature and building air conditioning preferred temperature is less than W0, and indoor temperature is less than outdoor temperature W1 degree, and the duration of this full curve section is greater than T1, and in building, personnel's number is less than the K1 of building setting personnel number;
Condition C: active power is greater than the S2 of air-conditioning rated power, and the variable quantity of active power is less than the K1 of air-conditioning rated power, and the difference of indoor temperature and building air conditioning preferred temperature is less than W0, and indoor temperature is less than outdoor temperature W2 degree, and the duration of this full curve section is greater than T2, and in building, personnel's number is greater than the K2 of building setting personnel number;
Wherein, the representative value of S0 is 30%, the representative value of the representative value of S1 to be the representative value of 50%, S2 be 80%, W0 is 1 degree Celsius, the representative value of W1 and W2 is 4 degrees Celsius, the representative value of T0 is 4 hours, and the representative value of T1 and T2 is 4 hours, and the representative value of K0 and K1 is 5%, the representative value of K2 is 60%, and the representative value in summer of building air conditioning preferred temperature is 26 degrees Celsius;
3) according to all full curve sections of the A that satisfies condition, calculate firm demand value, specific formula for calculation is:
P0=E0/D0
Wherein, P0 is firm demand value, and E0 is the meritorious integration of all full curve sections of A of satisfying condition, and D0 is the total duration of all full curve sections of A of satisfying condition;
4) according to all full curve sections of the B that satisfies condition, calculate temperature coefficient, specific formula for calculation is:
a1=(E1-P0×D1)/D1/H1
H1=∑△Tb/60/D1
Wherein, a1 is temperature coefficient, E1 is the meritorious integration of all full curve sections of B of satisfying condition, D1 is the total duration of all full curve sections of B of satisfying condition, H1 is the indoor/outdoor temperature-difference hourly average value of all full curve sections of B of satisfying condition, and ∑ △ Tb is the summation of indoor and outdoor temperature difference per minute of all full curve sections of B of satisfying condition;
5) to all full curve sections of the C that satisfies condition, calculate each whole point of hour per capita environment loads value, specific formula for calculation is:
A2[i]=(Px[i]-a1×△Tc[i]-P0)/m[i]
Wherein, i-th whole point hour per capita environment loads value in all full curve sections that A2 [i] is the C that satisfies condition, the burden with power at i-th whole point of hour in all full curve sections that Px [i] is the C that satisfies condition, the indoor and outdoor temperature at i-th whole point of hour in all full curve sections that △ Tc [i] is the C that satisfies condition is poor, personnel's number in the building at i-th whole point of hour in all full curve sections that m [i] is the C that satisfies condition;
6) according to all full curve sections of the C that satisfies condition, calculate per capita environment loads coefficient, specific formula for calculation is:
a2=∑△A2/n
Wherein, a2 is load coefficient per capita, and ∑ △ A2 is the summation of the whole point of hour per capita environment loads value satisfied condition in all full curve sections of C, and n is the number of the whole point of hour per capita environment loads value satisfied condition in all full curve sections of C;
7) calculate the discreet value of air conditioning electricity load, specific formula for calculation is:
P=P0+a1×△T+a2×m
Wherein, P is the discreet value of air conditioning electricity load, and △ T is the current indoor/outdoor temperature-difference of office building, and m is personnel's number in current building.
CN201510225624.8A 2015-05-06 2015-05-06 Office building air conditioning electricity predictor method Active CN104864560B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510225624.8A CN104864560B (en) 2015-05-06 2015-05-06 Office building air conditioning electricity predictor method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510225624.8A CN104864560B (en) 2015-05-06 2015-05-06 Office building air conditioning electricity predictor method

Publications (2)

Publication Number Publication Date
CN104864560A true CN104864560A (en) 2015-08-26
CN104864560B CN104864560B (en) 2017-08-25

Family

ID=53910574

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510225624.8A Active CN104864560B (en) 2015-05-06 2015-05-06 Office building air conditioning electricity predictor method

Country Status (1)

Country Link
CN (1) CN104864560B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117649062A (en) * 2024-01-30 2024-03-05 中国建筑西北设计研究院有限公司 Engineering configuration method and device for multi-heat source combined operation energy station system capacity

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101782258A (en) * 2009-01-19 2010-07-21 中华电信股份有限公司 Energy-saving method for air conditioner
JP2011185534A (en) * 2010-03-09 2011-09-22 Panasonic Corp Air conditioner
CN102779228A (en) * 2012-06-07 2012-11-14 华南理工大学 Method and system for online prediction on cooling load of central air conditioner in marketplace buildings
CN104112077A (en) * 2014-07-22 2014-10-22 上海申瑞继保电气有限公司 Method for calculating per capita air conditioner electricity consumption of office building
CN104181898A (en) * 2014-09-01 2014-12-03 东北电力大学 Intelligent control method and system for interactive home appliances on basis of time-of-use electricity price response

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101782258A (en) * 2009-01-19 2010-07-21 中华电信股份有限公司 Energy-saving method for air conditioner
JP2011185534A (en) * 2010-03-09 2011-09-22 Panasonic Corp Air conditioner
CN102779228A (en) * 2012-06-07 2012-11-14 华南理工大学 Method and system for online prediction on cooling load of central air conditioner in marketplace buildings
CN104112077A (en) * 2014-07-22 2014-10-22 上海申瑞继保电气有限公司 Method for calculating per capita air conditioner electricity consumption of office building
CN104181898A (en) * 2014-09-01 2014-12-03 东北电力大学 Intelligent control method and system for interactive home appliances on basis of time-of-use electricity price response

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117649062A (en) * 2024-01-30 2024-03-05 中国建筑西北设计研究院有限公司 Engineering configuration method and device for multi-heat source combined operation energy station system capacity

Also Published As

Publication number Publication date
CN104864560B (en) 2017-08-25

Similar Documents

Publication Publication Date Title
CN108039710B (en) Step characteristic-based air conditioner load-participating power grid day-ahead scheduling method
Mathieu et al. Using residential electric loads for fast demand response: The potential resource and revenues, the costs, and policy recommendations
CN103574845B (en) A kind of ice-storage system optimal control method based on cooling load prediction
CN109059195B (en) Control method and control system for central air conditioner for reducing load peak value of power grid
CN101929721B (en) Predicting method of central air conditioner energy-conservation control autoregressive (AR) model load predicting system
CN104214912A (en) Aggregation air conditioning load scheduling method based on temperature set value adjustment
CN104534556A (en) Heat supply control method based on energy consumption monitoring
CN104456845A (en) Public building central air-conditioning preopen time calculating method
CN104864559A (en) Per capital air conditioner energy consumption calculating method in public building
CN104238531A (en) Railway station energy management system and energy-saving control method
CN101090335A (en) Remote regulating method and system for indoor temp. and load of domestic air conditioner
CN109974218A (en) A kind of multi-online air-conditioning system regulation method based on prediction
CN103499136A (en) Ice storage control system with next-day energy consumption simulating function
CN103323664B (en) Electric power subitem energy consumption timesharing requirement method for early warning
CN110260479A (en) Tail end of central air conditioner monitoring and energy consumption management system and management method based on APP
CN105546768A (en) Energy-saving method and system of central air conditioner
CN203501384U (en) Central air-conditioner monitoring system
CN105783210B (en) A kind of multi-online air-conditioning system and its control method applied to subway station
Xue et al. Interactive building load management for smart grid
CN205505299U (en) Central air conditioning intelligence temperature control system
CN103471178B (en) Heat supply energy-saving control system and energy-saving control method
CN104864560A (en) Air conditioner electricity consumption pre-estimating method in office building
CN105864963B (en) A kind of polymerization air conditioner load control method based on transformation time priority list
CN107563547B (en) Comprehensive energy management and control method for optimizing depth of energy consumption of user side
CN105549389A (en) Household energy management algorithm based on building thermodynamic model

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into 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: 20200521

Address after: 215200 south of Lianyang road and east of Chang'an Road, Wujiang Economic and Technological Development Zone, Suzhou City, Jiangsu Province (Science and technology entrepreneurship Park)

Patentee after: Wujiang science and Technology Pioneer Park Management Service Co., Ltd

Address before: 200233, building 14, building 470, No. 4, Guiping Road, Shanghai, Xuhui District

Patentee before: SHANGHAI SUNRISE POWER TECHNOLOGY Co.,Ltd.

TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20201102

Address after: No. 599, South Road, Zhenze Town, Wujiang District, Suzhou City, Jiangsu Province 215200

Patentee after: Suzhou Chint Enterprise Development Co.,Ltd.

Address before: 215200 south of Lianyang road and east of Chang'an Road, Wujiang Economic and Technological Development Zone, Suzhou City, Jiangsu Province (Science and technology entrepreneurship Park)

Patentee before: Wujiang science and Technology Pioneer Park Management Service Co.,Ltd.