CN110555602A - Method and system for ten-day runoff distribution under condition of runoff data shortage - Google Patents

Method and system for ten-day runoff distribution under condition of runoff data shortage Download PDF

Info

Publication number
CN110555602A
CN110555602A CN201910759536.4A CN201910759536A CN110555602A CN 110555602 A CN110555602 A CN 110555602A CN 201910759536 A CN201910759536 A CN 201910759536A CN 110555602 A CN110555602 A CN 110555602A
Authority
CN
China
Prior art keywords
runoff
day
data
division
distribution
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
CN201910759536.4A
Other languages
Chinese (zh)
Other versions
CN110555602B (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.)
Kunming Electric Power Transaction Center Co ltd
Original Assignee
Kunming Electric Power Trading Center LLC
Huazhong University of Science and Technology
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 Kunming Electric Power Trading Center LLC, Huazhong University of Science and Technology filed Critical Kunming Electric Power Trading Center LLC
Priority to CN201910759536.4A priority Critical patent/CN110555602B/en
Publication of CN110555602A publication Critical patent/CN110555602A/en
Application granted granted Critical
Publication of CN110555602B publication Critical patent/CN110555602B/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/06Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063Operations research, analysis or management
    • G06Q10/0631Resource planning, allocation, distributing or scheduling for enterprises or organisations
    • G06Q10/06312Adjustment or analysis of established resource schedule, e.g. resource or task levelling, or dynamic rescheduling
    • 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
    • G06Q50/00Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/06Energy or water supply
    • 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
    • Y02DCLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
    • Y02D10/00Energy efficient computing, e.g. low power processors, power management or thermal management

Landscapes

  • Business, Economics & Management (AREA)
  • Human Resources & Organizations (AREA)
  • Engineering & Computer Science (AREA)
  • Economics (AREA)
  • Strategic Management (AREA)
  • Theoretical Computer Science (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Health & Medical Sciences (AREA)
  • Marketing (AREA)
  • General Physics & Mathematics (AREA)
  • General Business, Economics & Management (AREA)
  • Tourism & Hospitality (AREA)
  • Physics & Mathematics (AREA)
  • Public Health (AREA)
  • Primary Health Care (AREA)
  • Water Supply & Treatment (AREA)
  • Development Economics (AREA)
  • Educational Administration (AREA)
  • General Health & Medical Sciences (AREA)
  • Game Theory and Decision Science (AREA)
  • Operations Research (AREA)
  • Quality & Reliability (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The invention discloses a ten-day runoff distribution method and a ten-day runoff distribution system under the condition of runoff data shortage, belonging to the field of hydrological research in runoff data shortage or data-free areas, wherein the method comprises the following steps: carrying out runoff division with ten days as a unit on the known runoff data according to the constructed quadratic programming model; judging whether the ten-day runoff obtained by the division reaches a set standard or not according to the evaluation parameters, if so, finishing the division; if not, the division is continued. The method can fully utilize the existing monthly runoff (annual runoff data) data to supplement the runoff data of smaller time units, effectively reduces the inherent error of the existing distribution method, does not need parameter limitation, enables the distribution result to be more fit with the actual situation, can be applied to the data-free hydrology research field to solve the problem of related runoff simulation distribution, is simple and efficient, is convenient to implement, and has wide application value.

Description

method and system for ten-day runoff distribution under condition of runoff data shortage
Technical Field
The invention belongs to the field of hydrological research in areas with or without runoff data shortage, and particularly relates to a ten-day runoff distribution method and system under the condition of runoff data shortage.
background
Hydrological research in areas without data or data shortage is one of the hot spots and difficult problems of modern international hydrological water resource research, and as far as China is concerned, many reservoirs and hydropower stations have no systematic runoff data record at the beginning of construction, and most runoff data take months as the minimum unit. However, in practical engineering application or scientific research, more detailed data support is often required, such as ten days and day data, and at this time, the monthly runoff data needs to be scientifically and reasonably divided.
However, the current dividing method is relatively rough, the common method is to divide the monthly runoff into every ten days or every day (equal divisions), although the method is simple, the method has obvious defects, the formation of the surface runoff is influenced by various factors such as climate, land utilization type and the like, so the time sequence of the runoff and the relation of rainfall-runoff generally have obvious non-stationarity, and the equal division causes the runoff value jump between ten days and ten days (day and day) to be often overlarge or zero, which is contrary to the reality.
of course, many random numbers generated by random simulation methods are used to replace actual runoff data, but the simulation of statistical characteristics and the establishment of random models (including order determination, parameter estimation and the like) are involved, a series of factors such as trend terms, skip terms, period terms, random terms and the like need to be considered, the steps are relatively complex, the calculation burden is increased, and the existing runoff data is not fully utilized.
Generally speaking, the existing runoff distribution method has the problem that the simulation data is inconsistent with the actual situation due to obvious errors in distributing the ten-day runoff.
Disclosure of Invention
aiming at the defects of the prior art, the invention aims to provide a method and a system for allocating the ten-day runoff under the condition of runoff data shortage, and aims to solve the problem that the simulation data is inconsistent with the actual situation due to obvious errors when the ten-day runoff is allocated by adopting the conventional method.
in order to achieve the above object, the present invention provides a ten-day runoff allocation method under a condition of runoff data shortage, including:
(1) Carrying out runoff division with ten days as a unit on the known runoff data according to the constructed quadratic programming model;
(2) judging whether the ten-day runoff obtained by the division reaches a set standard or not according to the evaluation parameters, if so, finishing the division; if not, returning to execute the step (1).
Further, the step (1) is specifically to perform the runoff division with ten days as a unit on the known runoff data by adopting a large-scale linear programming solver according to the constructed quadratic programming model.
Further, the quadratic programming model is as follows:
wherein i represents the number of runoff in ten days, N represents the number of divided ten days in the runoff sequence, and QiIs the ten-day runoff of the time period i.
Further, the quadratic programming model satisfies the following constraint conditions:
and (3) total balance constraint: q3j+Q3j+1+Q3j+2=3qjwherein j represents a monthly runoff number, qjThe monthly runoff for time period j; non-negative constraints: qi≥0。
Further, the evaluation parameters are:
Wherein Cvm is a monthly runoff sequence deviation coefficient, Cvx is a runoff sequence deviation coefficient after the optimized allocation, and K is an evaluation parameter, which reflects the variation of the runoff sequence and the monthly runoff sequence deviation parameter.
in another aspect, the present invention provides a ten-day runoff distribution system under the condition of runoff data shortage, including:
The ten-day runoff allocation module is used for carrying out runoff division with ten-day as a unit on the known runoff data according to the constructed quadratic programming model;
The distribution result evaluation module is used for judging whether the ten-day runoff obtained by the division reaches a set standard or not according to the evaluation parameters, and if so, the division is finished; if not, continuing to adopt the ten-day runoff allocation module to carry out ten-day runoff allocation.
Further, the ten-day runoff allocation module performs runoff division on known runoff data by taking ten days as a unit by adopting a large-scale linear programming solver according to the constructed quadratic programming model.
further, the quadratic programming model is as follows:
wherein i represents the number of runoff in ten days, N represents the number of divided ten days in the runoff sequence, and Qiis the ten-day runoff of the time period i.
Further, the quadratic programming model satisfies the following constraint conditions:
And (3) total balance constraint: q3j+Q3j+1+Q3j+2=3qjwherein j represents a monthly runoff number, qjThe monthly runoff for time period j; non-negative constraints: qi≥0。
Further, the evaluation parameters are:
wherein Cvm is a monthly runoff sequence deviation coefficient, Cvx is a runoff sequence deviation coefficient after the optimized allocation, and K is an evaluation parameter, which reflects the variation of the runoff sequence and the monthly runoff sequence deviation parameter.
Through the technical scheme, compared with the prior art, the invention has the following beneficial effects:
(1) Compared with the conventional annual ten-day runoff distribution method without data, the method disclosed by the invention can fully utilize the conventional monthly runoff data (annual runoff data) data to supplement the runoff data of smaller time units, so that the inherent error of the conventional distribution method is effectively reduced, the parameter limitation is not required, and the distribution result is more suitable for the actual situation.
(2) according to the method, the large-scale linear programming solver is called, so that the global optimal solution of the quadratic programming model can be quickly and accurately solved, the optimal allocation of the annual flow every ten days is realized, and the method is simple, efficient, convenient to implement and wide in application value.
drawings
Fig. 1 is a flow chart of a ten-day runoff allocation method under a condition of runoff data shortage according to the present invention;
FIG. 2(a) is a schematic diagram of a ten-day traffic sequence allocation obtained by applying the method of the present invention;
Fig. 2(b) is a schematic diagram of a ten-day traffic sequence allocation obtained by applying the conventional allocation method.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, the present invention is described in further detail below with reference to the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are merely illustrative of the invention and are not intended to limit the invention.
Referring to fig. 1, in an aspect, an embodiment of the present invention provides a ten-day runoff allocation method in a runoff data shortage situation, including:
(1) and carrying out runoff division with ten days as a unit on the known runoff data according to the constructed quadratic programming model.
specifically, according to known monthly runoff data or part of known ten-day runoff data, a mathematical programming model is established by taking the minimum sum of squares of fluctuation values of adjacent ten-day runoff in a ten-day runoff sequence as a target, and because a linear programming solver gurobi adopts an advanced optimization technology, the advantages of a multi-core processor are fully utilized, efficient and stable solution is realized, and a convenient and light interface is provided for direct calling, therefore, the invention adopts a large-scale linear programming solver gurobi to perform runoff division on the known runoff data by taking ten-day as a unit, so that the mathematical programming model is adapted to the large-scale linear programming solver to obtain an optimal solution, and is reconstructed into a strict quadratic programming model:Wherein i represents the number of runoff in ten days, N represents the number of divided ten days in the runoff sequence, and QiIs the ten-day runoff of the time period i. The quadratic programming model in the invention meets the following constraint conditions: and (3) total balance constraint: q3j+Q3j+1+Q3j+2=3qjwherein j represents a monthly runoff number, qjThe monthly runoff for time period j; non-negative constraints: qi≥0。
(2) Judging whether the ten-day runoff obtained by the division reaches a set standard or not according to the evaluation parameters, if so, finishing the division; if not, returning to execute the step (1).
Specifically, the ten-day runoff evaluation parameters obtained by the allocation method are as follows:
Wherein Cvm is a month runoff sequence deviation coefficient, Cvx is an optimized post-allocation ten-day runoff sequence deviation coefficient, K is an evaluation parameter, which reflects the variation of the ten-day runoff sequence and the month runoff sequence deviation parameter, and the smaller the value of K, the more similar the post-allocation ten-day runoff sequence and the previous month runoff sequence, and according to the previous allocation experience, it is generally considered that the similarity between the ten-day runoff sequence allocation and the month runoff meets the requirement when K is less than 0.01.
In order to verify the effectiveness of the method, a runoff sequence of a certain year of a hydrological station in a certain data shortage area is selected for distribution, fig. 2(a) and fig. 2(b) are distribution sequences generated by the distribution method and the existing distribution method respectively, and it can be seen that the ten-day runoff sequence of fig. 2(a) has more stable runoff fluctuation and more practical fluctuation trend compared with fig. 2(b), wherein the K value of fig. 2(a) is 0.06%, the K value of fig. 2(b) is 2.90%, and the similarity of the ten-day runoff sequence distribution is improved by nearly two orders of magnitude compared with that of the ten-day runoff sequence distribution, and meets the requirement of an evaluation standard, thereby further proving that the ten-day runoff sequence generated by the distribution method of the invention has more representative and practical significance.
Another aspect of the embodiments of the present invention provides a ten-day runoff distribution system under a condition of runoff data shortage, including: the ten-day runoff allocation module is used for carrying out runoff division with ten-day as a unit on the known runoff data according to the constructed quadratic programming model; the distribution result evaluation module is used for judging whether the ten-day runoff obtained by the division reaches a set standard or not according to the evaluation parameters, and if so, the division is finished; if not, continuing to adopt the ten-day runoff allocation module to carry out ten-day runoff allocation. The ten-day runoff allocation module of the invention adopts a large-scale linear programming solver to divide the runoff with the ten-day units for the known runoff data according to the constructed quadratic programming model. The quadratic programming model adopted is as follows:wherein i represents the number of runoff in ten days, N represents the number of divided ten days in the runoff sequence, and QiIs the ten-day runoff of the time period i. The quadratic programming model satisfies the following constraint conditions: and (3) total balance constraint: q3j+Q3j+1+Q3j+2=3qjwherein j represents a monthly runoff number, qjThe monthly runoff for time period j; non-negative constraints: qiIs more than or equal to 0. The evaluation parameters used were:wherein Cvm is a monthly runoff sequence deviation coefficient, Cvx is a runoff sequence deviation coefficient after the optimized allocation, and K is an evaluation parameter, which reflects the variation of the runoff sequence and the monthly runoff sequence deviation parameter.
It will be understood by those skilled in the art that the foregoing is only a preferred embodiment of the present invention, and is not intended to limit the invention, and that any modification, equivalent replacement, or improvement made within the spirit and principle of the present invention should be included in the scope of the present invention.

Claims (10)

1. A ten-day runoff allocation method under the condition of runoff data shortage is characterized by comprising the following steps:
(1) Carrying out runoff division with ten days as a unit on the known runoff data according to the constructed quadratic programming model;
(2) Judging whether the ten-day runoff obtained by the division reaches a set standard or not according to the evaluation parameters, if so, finishing the division; if not, returning to execute the step (1).
2. The method for ten-day runoff distribution under the condition of runoff data shortage according to claim 1, wherein the step (1) is specifically to perform runoff division in ten-day units on the known runoff data by adopting a large-scale linear programming solver according to the constructed quadratic programming model.
3. the method for ten-day runoff distribution under the condition of runoff data shortage according to claim 1 or 2, wherein the quadratic programming model is as follows:
Wherein i represents the number of runoff in ten days, N represents the number of divided ten days in the runoff sequence, and Qiis the ten-day runoff of the time period i.
4. The method for ten-day runoff distribution under the condition of runoff data shortage according to claim 3, wherein the quadratic programming model meets the following constraint conditions:
And (3) total balance constraint: q3j+Q3j+1+Q3j+2=3qjWherein j represents a monthly runoff number, qjThe monthly runoff for time period j; non-negative constraints: qi≥0。
5. The method for ten-day runoff distribution in the case of runoff data shortage according to any one of claims 1 to 4, wherein the evaluation parameters are:
Wherein Cvm is a monthly runoff sequence deviation coefficient, Cvx is a runoff sequence deviation coefficient after the optimized allocation, and K is an evaluation parameter, which reflects the variation of the runoff sequence and the monthly runoff sequence deviation parameter.
6. A ten day runoff distribution system under the condition of runoff data shortage, which is characterized by comprising:
The ten-day runoff allocation module is used for carrying out runoff division with ten-day as a unit on the known runoff data according to the constructed quadratic programming model;
The distribution result evaluation module is used for judging whether the ten-day runoff obtained by the division reaches a set standard or not according to the evaluation parameters, and if so, the division is finished; if not, continuing to adopt the ten-day runoff allocation module to carry out ten-day runoff allocation.
7. The system according to claim 6, wherein the ten-day runoff distribution module performs runoff division in ten-day units on the known runoff data by using a large-scale linear programming solver according to the constructed quadratic programming model.
8. The system according to claim 6 or 7, wherein the quadratic programming model is:
Wherein i represents the number of runoff in ten days, N represents the number of divided ten days in the runoff sequence, and QiIs the ten-day runoff of the time period i.
9. the system according to claim 8, wherein the quadratic programming model satisfies the following constraints:
And (3) total balance constraint: q3j+Q3j+1+Q3j+2=3qjwherein j represents a monthly runoff number, qjThe monthly runoff for time period j; non-negative constraints: qi≥0。
10. A ten day runoff distribution system according to any one of claims 1 to 9 wherein the evaluation parameters are:
Wherein Cvm is a monthly runoff sequence deviation coefficient, Cvx is a runoff sequence deviation coefficient after the optimized allocation, and K is an evaluation parameter, which reflects the variation of the runoff sequence and the monthly runoff sequence deviation parameter.
CN201910759536.4A 2019-08-16 2019-08-16 Method and system for ten-day runoff distribution under condition of runoff data shortage Active CN110555602B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910759536.4A CN110555602B (en) 2019-08-16 2019-08-16 Method and system for ten-day runoff distribution under condition of runoff data shortage

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910759536.4A CN110555602B (en) 2019-08-16 2019-08-16 Method and system for ten-day runoff distribution under condition of runoff data shortage

Publications (2)

Publication Number Publication Date
CN110555602A true CN110555602A (en) 2019-12-10
CN110555602B CN110555602B (en) 2023-04-07

Family

ID=68737741

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910759536.4A Active CN110555602B (en) 2019-08-16 2019-08-16 Method and system for ten-day runoff distribution under condition of runoff data shortage

Country Status (1)

Country Link
CN (1) CN110555602B (en)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103577895A (en) * 2013-11-08 2014-02-12 河海大学 Method for forecasting monthly runoff through secondary coupling under condition of data shortage
US20180128939A1 (en) * 2016-11-09 2018-05-10 Guizhou Normal University Adjustable karst underground water and soil loss simulation apparatus
CN108053083A (en) * 2018-01-16 2018-05-18 河南创辉水利水电工程有限公司 A kind of hydro plant with reservoir non-flood period combined optimization power generation dispatching method

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103577895A (en) * 2013-11-08 2014-02-12 河海大学 Method for forecasting monthly runoff through secondary coupling under condition of data shortage
US20180128939A1 (en) * 2016-11-09 2018-05-10 Guizhou Normal University Adjustable karst underground water and soil loss simulation apparatus
CN108053083A (en) * 2018-01-16 2018-05-18 河南创辉水利水电工程有限公司 A kind of hydro plant with reservoir non-flood period combined optimization power generation dispatching method

Non-Patent Citations (6)

* Cited by examiner, † Cited by third party
Title
周婷等: "基于小波分解的优化支持向量机模型在水库年径流预测中的应用", 《华北水利水电大学学报(自然科学版)》 *
徐长江等: "资料短缺地区径流模拟研究" *
洪家宁等: "反褶积模型在日径流模拟中的应用" *
罗吉忠等: "基于SWAT模型的缺资料流域径流模拟研究", 《西南民族大学学报(自然科学版)》 *
胡彩虹等: "汾河水库上游流域径流年内分配时程变化规律研究", 《干旱区资源与环境》 *
靳晟等: "设计年径流及其年内分配计算软件的开发研究", 《人民黄河》 *

Also Published As

Publication number Publication date
CN110555602B (en) 2023-04-07

Similar Documents

Publication Publication Date Title
CN110501733B (en) Self-adaptive grid VRS generation and service method
CN104036327A (en) Fast bus load forecasting method for smart distribution network
CN109766349A (en) The anti-weighing method of task, device, computer equipment and storage medium
CN103810538A (en) Electricity sale quantity distribution method and device
Veselov et al. Methodological approach for harmonization of the investment and pricing policy options in the electric power industry
CN103853939A (en) Combined forecasting method for monthly load of power system based on social economic factor influence
CN108133005A (en) A kind of environmental model analogy method, terminal device and storage medium based on memory database
CN104951903A (en) Custom asset evaluation method based on financial asset management system architecture
CN110555602B (en) Method and system for ten-day runoff distribution under condition of runoff data shortage
CN113902160A (en) Water regulation method applied to hydro-junction engineering
CN103246539B (en) Based on the generalized information system module loading method of user's use habit
Zagona et al. RiverWare
CN111090710B (en) Distribution network multi-time scale graph-model abnormal automatic maintenance method
CN109872250B (en) Method and system for settlement processing of initial cardinality deviation electric quantity of electric power market
CN109989747B (en) Oil field injection and production well group development non-uniformity quantitative characterization method
CN112036604A (en) Medium runoff forecasting method considering multiple time sequence process factors
CN112953846A (en) Domain name bandwidth cost optimization type delivery method and device and computer equipment
CN108599976B (en) Network quality evaluation model construction method and device and computer equipment
CN111080115A (en) Water supply hydraulic model-based water consumption distribution method and device
CN115953251B (en) Insurance commission accounting method and system
CN110162839A (en) A kind of discrimination method and system of Sediment variation influence factor contribution rate
CN115965056A (en) Medium-and-long-term electric quantity variable weight combination prediction method and system
CN115269564B (en) Centralized creative migration method for large-scale system
CN103559266A (en) Multi-mode matching method and device
Chunye Research on Measuring Method of Function Size of Information System

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: 20230506

Address after: No.73 Tuodong Road, Guandu District, Kunming, Yunnan 650011

Patentee after: KUNMING ELECTRIC POWER TRANSACTION CENTER CO.,LTD.

Address before: No. 73 Tuodong Road, Panlong District, Kunming City, Yunnan Province, 650011

Patentee before: KUNMING ELECTRIC POWER TRANSACTION CENTER CO.,LTD.

Patentee before: HUAZHONG University OF SCIENCE AND TECHNOLOGY