Most downloaded articles
2018, 46(2):128-133.
DOI: 10.3876/j.issn.1000-1980.2018.02.005
Abstract:
To improve the operational safety of earth rockfill dam, a longitudinal reinforcement is proposed. This rigid structure has triple functions as the seepage controlling, stress resistance and deformation resistance. Based on the classical mechanics, its anti-seepage effect, anti-deformation ability and the capacity to withstand external force is quantitatively analyzed. In addition, the quantitative relationships between the reinforcement thickness and the downstream water level, as well as the material properties of the rockfill dam are derived. For the theoretical derivation, the reinforced body is regarded as a fixed-end beam, loaded by water pressures from upstream and downstream, the active earth pressure from rockfill and the vertical drag force due to the differential settlement between core walls and rockfill. The top displacement together with the rotation angle and the bottom stress of the fixed end are derived from the deformation differential equation and validated through designing reviews. Subsequently, the construction method of buried steel pipe, which is adopted to perform the base grouting and also to meet the strength requirements, is put forward. The proposed method plays an important guiding role in rehabilitating and reinforcing the dangerous earth rockfill dams, as well as the design and construction of new dams.
To improve the operational safety of earth rockfill dam, a longitudinal reinforcement is proposed. This rigid structure has triple functions as the seepage controlling, stress resistance and deformation resistance. Based on the classical mechanics, its anti-seepage effect, anti-deformation ability and the capacity to withstand external force is quantitatively analyzed. In addition, the quantitative relationships between the reinforcement thickness and the downstream water level, as well as the material properties of the rockfill dam are derived. For the theoretical derivation, the reinforced body is regarded as a fixed-end beam, loaded by water pressures from upstream and downstream, the active earth pressure from rockfill and the vertical drag force due to the differential settlement between core walls and rockfill. The top displacement together with the rotation angle and the bottom stress of the fixed end are derived from the deformation differential equation and validated through designing reviews. Subsequently, the construction method of buried steel pipe, which is adopted to perform the base grouting and also to meet the strength requirements, is put forward. The proposed method plays an important guiding role in rehabilitating and reinforcing the dangerous earth rockfill dams, as well as the design and construction of new dams.
2019, 47(1):1-6.
DOI: 10.3876/j.issn.1000-1980.2019.01.001
Abstract:
To improve the accuracy of flood forecasting, 13 flood events between 2000 and 2017 were selected from the sections of Lutaizi and Lubeng, and the Xin’anjiang three-source model(XAJ model)was used to simulate the flood events. The Lutaizi basin can be divided into 9 natural sub-basins and the Lubeng basin can be divided into 4 natural sub-basins. Then, the XAJ model was calibrated for these basins. The node generalizability was carried out in the mainstream section of the Huaihe River with flood diversion and flood retarding areas. Combined with the diversion ratio method and the Muskingen method, the simulation method of flood diversion was established, and fairly good forecasting results were obtained. The research results show that: for the relative error of flood volume, the qualified rates of Lutaizi and Lubeng region are both up to 100%, and the Lutaizi reached a qualified rate of 100% in the flood peak, while the Lubeng region reached 92. 3%. In two floods in 2003 and 2007, the relative error and the certainty coefficient of flood peak have been significantly improved.
To improve the accuracy of flood forecasting, 13 flood events between 2000 and 2017 were selected from the sections of Lutaizi and Lubeng, and the Xin’anjiang three-source model(XAJ model)was used to simulate the flood events. The Lutaizi basin can be divided into 9 natural sub-basins and the Lubeng basin can be divided into 4 natural sub-basins. Then, the XAJ model was calibrated for these basins. The node generalizability was carried out in the mainstream section of the Huaihe River with flood diversion and flood retarding areas. Combined with the diversion ratio method and the Muskingen method, the simulation method of flood diversion was established, and fairly good forecasting results were obtained. The research results show that: for the relative error of flood volume, the qualified rates of Lutaizi and Lubeng region are both up to 100%, and the Lutaizi reached a qualified rate of 100% in the flood peak, while the Lubeng region reached 92. 3%. In two floods in 2003 and 2007, the relative error and the certainty coefficient of flood peak have been significantly improved.
2020, 48(1):1-7.
DOI: 10.3876/j.issn.1000-1980.2020.01.001
Abstract:
Based on the development of hydrology science and the basic framework of physical and virtual water coupled flow theory, this study indicated that the process hydrology is an interdiscipline that can reveal the coupled flow law and associated effects of physical and virtual water in complicated system. In addition, this study makes a detailed analysis of the scientific connotation, main characteristics, framework, discipline basis, method system and core tasks of process hydrology. Given the theoretical value and practical significance of process hydrology, it can be concluded that the development of process hydrology can provide measures and new views to cope with new challenges in hydrology that the sustainable developments of nature, economy and society faces.
Based on the development of hydrology science and the basic framework of physical and virtual water coupled flow theory, this study indicated that the process hydrology is an interdiscipline that can reveal the coupled flow law and associated effects of physical and virtual water in complicated system. In addition, this study makes a detailed analysis of the scientific connotation, main characteristics, framework, discipline basis, method system and core tasks of process hydrology. Given the theoretical value and practical significance of process hydrology, it can be concluded that the development of process hydrology can provide measures and new views to cope with new challenges in hydrology that the sustainable developments of nature, economy and society faces.
2021, 49(1):1-6.
DOI: 10. 3876/ j. issn.1000-1980.2021.01.001
Abstract:
This study introduces the basic situation of China‘s small and medium-sized rivers and flood disasters, and then key difficulties as well as urgent scientific and technological problems in the flood prevention of small and medium-sized rivers were analyzed. Thereafter, the study puts forward the strategy and measures for the flood forecasting and early warning of small and medium- sized rivers. Besides, prospects for the future forecasting, risk early warning, and refined risk management concepts and practices are also proposed. It will realize the change from the hydrological element forecasting to the flood impact forecasting, from the flood threshold-based warning to the flood risk warning, from the " prediction-response" emergency disposal to the grid refinement risk management, from the monitoring data provision to the information product service, in which the development level of the hydrological modernization is continuously improved.
This study introduces the basic situation of China‘s small and medium-sized rivers and flood disasters, and then key difficulties as well as urgent scientific and technological problems in the flood prevention of small and medium-sized rivers were analyzed. Thereafter, the study puts forward the strategy and measures for the flood forecasting and early warning of small and medium- sized rivers. Besides, prospects for the future forecasting, risk early warning, and refined risk management concepts and practices are also proposed. It will realize the change from the hydrological element forecasting to the flood impact forecasting, from the flood threshold-based warning to the flood risk warning, from the " prediction-response" emergency disposal to the grid refinement risk management, from the monitoring data provision to the information product service, in which the development level of the hydrological modernization is continuously improved.
2020, 48(1):22-28.
DOI: 10.3876/j.issn.1000-1980.2020.01.004
Abstract:
This study reviewed the environmental behavior of microplastics(including the adsorption of organic pollutants and the release of organic additives after entering water), and the interaction between microplastics and organic pollutants after entering into water body, including its mechanism, influencing factors and biological toxicity. It also analyzed the deficiencies of current researches, and believed that, in order to provide theoretical basis for further study on the behavior and function of microplastics, it is necessary to expand the research on the interaction between microplastics and organic pollutants, such as strengthening the research in the freshwater environment, establishing systematic research system, and paying attention to the effect of biofilm
This study reviewed the environmental behavior of microplastics(including the adsorption of organic pollutants and the release of organic additives after entering water), and the interaction between microplastics and organic pollutants after entering into water body, including its mechanism, influencing factors and biological toxicity. It also analyzed the deficiencies of current researches, and believed that, in order to provide theoretical basis for further study on the behavior and function of microplastics, it is necessary to expand the research on the interaction between microplastics and organic pollutants, such as strengthening the research in the freshwater environment, establishing systematic research system, and paying attention to the effect of biofilm
2019, 47(1):7-12.
DOI: 10.3876/j.issn.1000-1980.2019.01.002
Abstract:
In order to make full use of the continues pumping test data to estimate the hydraulic conductivity(K), the transmissivity(T)and the storage coefficient(μ*)of different head depression cones in a confined aquifer, a dimensionless analytical solution of drawdown transients with variable pumping rates was derived based on the Theis formula and the principle of superposition. Furthermore, a new type curve method was proposed. The in-situ pumping test was conducted to obtain the drawdown record perturbed by stepwise rates. Then, the drawdown curve was used to match the corresponding type curve. After that, this method selects one match point and records its coordinate values on both the type curve and the real drawdown curve. By substituting these recorded values into related formulas, the hydraulic parameters of the confined aquifer can be determined. This method was applied to a field pumping test conducted at a southern highway in Wuxi, Jiangsu Province, China. The type curve matching process is implemented in stages. Field test results reveal that the hydraulic conductivity of the confined aquifer with the 2-3 silt layer are 5. 12×10-4 cm/s, 2. 54×10-3 cm/s and 2. 83×10-3 cm/s respectively, which are corresponding to three-stage stepwise pumping rates.
In order to make full use of the continues pumping test data to estimate the hydraulic conductivity(K), the transmissivity(T)and the storage coefficient(μ*)of different head depression cones in a confined aquifer, a dimensionless analytical solution of drawdown transients with variable pumping rates was derived based on the Theis formula and the principle of superposition. Furthermore, a new type curve method was proposed. The in-situ pumping test was conducted to obtain the drawdown record perturbed by stepwise rates. Then, the drawdown curve was used to match the corresponding type curve. After that, this method selects one match point and records its coordinate values on both the type curve and the real drawdown curve. By substituting these recorded values into related formulas, the hydraulic parameters of the confined aquifer can be determined. This method was applied to a field pumping test conducted at a southern highway in Wuxi, Jiangsu Province, China. The type curve matching process is implemented in stages. Field test results reveal that the hydraulic conductivity of the confined aquifer with the 2-3 silt layer are 5. 12×10-4 cm/s, 2. 54×10-3 cm/s and 2. 83×10-3 cm/s respectively, which are corresponding to three-stage stepwise pumping rates.
2020, 48(2):95-101.
DOI: 10.3876/j.issn.1000-1980.2020.02.001
Abstract:
Although the Piecewise Muskingum method for river flood calculation has been widely used in China, the major limitation that hinders the further application of the Piecewise Muskingum method within distributed models still needs to be addressed. This method with empirical parameters cannot consider the impact of changing river cross-section characteristics and hydraulic roughness on the flood computation. In this study, the Muskingum-Cunge-Todini variable parameter method and the non-linear reservoir method, which could consider the channel cross-section geometry and hydraulic roughness, were used in the river flood calculation in the middle part of Huaihe river, from Wujiadu to Xiaoliuxiang, and the middle part of Hutuo River, from Huangbizhuang to Beizhongshan. In addition, the river flood calculation with the Piecewise Muskingum method also has been conducted as a comparison. Initial tests indicate that the Muskingum-Cunge-Todini variable parameter method and the non-linear reservoir method have high simulation accuracy in the studied cases. The correct rate of flood peak is above 86%, and the coefficient of certainty is greater than 0. 8. This study can provide references for the flood calculation considering the change of the channel cross section and hydraulic roughness characteristics.
Although the Piecewise Muskingum method for river flood calculation has been widely used in China, the major limitation that hinders the further application of the Piecewise Muskingum method within distributed models still needs to be addressed. This method with empirical parameters cannot consider the impact of changing river cross-section characteristics and hydraulic roughness on the flood computation. In this study, the Muskingum-Cunge-Todini variable parameter method and the non-linear reservoir method, which could consider the channel cross-section geometry and hydraulic roughness, were used in the river flood calculation in the middle part of Huaihe river, from Wujiadu to Xiaoliuxiang, and the middle part of Hutuo River, from Huangbizhuang to Beizhongshan. In addition, the river flood calculation with the Piecewise Muskingum method also has been conducted as a comparison. Initial tests indicate that the Muskingum-Cunge-Todini variable parameter method and the non-linear reservoir method have high simulation accuracy in the studied cases. The correct rate of flood peak is above 86%, and the coefficient of certainty is greater than 0. 8. This study can provide references for the flood calculation considering the change of the channel cross section and hydraulic roughness characteristics.
2015, 43(5):384-394.
DOI: 10.3876/j.issn.1000-1980.2015.05.002
Abstract:
Water resources constraints constitute a significant bottleneck in China. The origins of water are worth further research, as they relate to the Earth’s formation, the origins of the water globe, and the issue of deep water. First, we attempt to divide the hydrosphere into several belts and focus on the complex interaction between spherical layers. Then, water resources development issues are discussed, the concept of the comprehensive effect and related evaluation of the evolution of an entire basin is developed based on the spherical layer movement, and the evolution of the middle reaches of the Yellow River, the middle and upper reaches of the Yangtze River, and the inland river basin is analyzed based on this concept. Water resources development should consider three properties: limits, relativity, and ecology, as well as the instructive observation that water can carry a boat and also overturn a boat. The comprehensive development and utilization of water from six different sources, including rainwater, river water, lake water, groundwater, reservoir water, and sea water, are examined. We also discuss prevention and mitigation of six types of disasters related to the six water sources: flooding, waterlogging, drought, storm surges, geological disasters, and water pollution. It is pointed out that segmented control of water quality in the reservoir area, comprehensive development of surface water and groundwater, and construction of underground reservoirs are important. Finally, some suggestions for water resources development and disaster prevention are proposed: cooperation between multiple departments and multiple disciplines, implementation of water-saving measures, construction of an enhanced water resources network, determination of demand according to provision, and improvement of an investigation system for evaluation and development.
Water resources constraints constitute a significant bottleneck in China. The origins of water are worth further research, as they relate to the Earth’s formation, the origins of the water globe, and the issue of deep water. First, we attempt to divide the hydrosphere into several belts and focus on the complex interaction between spherical layers. Then, water resources development issues are discussed, the concept of the comprehensive effect and related evaluation of the evolution of an entire basin is developed based on the spherical layer movement, and the evolution of the middle reaches of the Yellow River, the middle and upper reaches of the Yangtze River, and the inland river basin is analyzed based on this concept. Water resources development should consider three properties: limits, relativity, and ecology, as well as the instructive observation that water can carry a boat and also overturn a boat. The comprehensive development and utilization of water from six different sources, including rainwater, river water, lake water, groundwater, reservoir water, and sea water, are examined. We also discuss prevention and mitigation of six types of disasters related to the six water sources: flooding, waterlogging, drought, storm surges, geological disasters, and water pollution. It is pointed out that segmented control of water quality in the reservoir area, comprehensive development of surface water and groundwater, and construction of underground reservoirs are important. Finally, some suggestions for water resources development and disaster prevention are proposed: cooperation between multiple departments and multiple disciplines, implementation of water-saving measures, construction of an enhanced water resources network, determination of demand according to provision, and improvement of an investigation system for evaluation and development.
2017, 45(3):277-282.
DOI: 10.3876/j.issn.1000-1980.2017.03.014
Abstract:
It is necessary to monitor the state of ZnO arresters to guarantee their safe operation. In order to select a best way to monitor the state of arresters, the advantages, disadvantages, and application of each monitoring method were analyzed. The results show that the harmonic analysis method is the most effective method. The properties of the harmonic analysis method in practical application could be improved by modifying the resistive leakage current and selecting an appropriate window function when conducting fast discrete Fourier transform on the collected signals. Through review of monitoring technologies in China and abroad, it can be found that in the future, the overall development trend of online monitoring methods of ZnO arresters is to develop fully automatic and multi-functional integrated online insulation monitoring systems.
It is necessary to monitor the state of ZnO arresters to guarantee their safe operation. In order to select a best way to monitor the state of arresters, the advantages, disadvantages, and application of each monitoring method were analyzed. The results show that the harmonic analysis method is the most effective method. The properties of the harmonic analysis method in practical application could be improved by modifying the resistive leakage current and selecting an appropriate window function when conducting fast discrete Fourier transform on the collected signals. Through review of monitoring technologies in China and abroad, it can be found that in the future, the overall development trend of online monitoring methods of ZnO arresters is to develop fully automatic and multi-functional integrated online insulation monitoring systems.
2020, 48(3):189-194.
DOI: 10.3876/j.issn.1000-1980.2020.03.001
Abstract:
On the basis of the original Xin’anjiang model, this study solved the matching problem between sub-catchment and grid channel, and enabled the model to calculate channel routing according to actual river conditions. A hybrid rainfall-runoff model(named XAJ-DCH)was further proposed by coupling the XAJ model with diffusion wave method and Muskingum-Cunge-Todini(MCT)method by accounting for the lateral inflow. The flood prediction capabilities of the traditional XAJ model and the XAJ-DCH model were compared at hourly scale for the Chengcun catchment. The performance of both models is comparable and satisfactory, and thus, the rationality of the new model is verified. By improving the channel routing method in the Xin’anjiang model, the spatiotemporal variability can be further considered. Furthermore, compared with the Xin’anjiang model, the XAJ-DCH model can not only be used to simulate the discharge at the outlet station, but also to predict the water level and discharge of the interior station in the catchment. In addition, the diffusion wave method was introduced into the XAJ-DCH model, which can improve the prediction accuracy in hydraulically mild slope area by accounting for the backwater effects.
On the basis of the original Xin’anjiang model, this study solved the matching problem between sub-catchment and grid channel, and enabled the model to calculate channel routing according to actual river conditions. A hybrid rainfall-runoff model(named XAJ-DCH)was further proposed by coupling the XAJ model with diffusion wave method and Muskingum-Cunge-Todini(MCT)method by accounting for the lateral inflow. The flood prediction capabilities of the traditional XAJ model and the XAJ-DCH model were compared at hourly scale for the Chengcun catchment. The performance of both models is comparable and satisfactory, and thus, the rationality of the new model is verified. By improving the channel routing method in the Xin’anjiang model, the spatiotemporal variability can be further considered. Furthermore, compared with the Xin’anjiang model, the XAJ-DCH model can not only be used to simulate the discharge at the outlet station, but also to predict the water level and discharge of the interior station in the catchment. In addition, the diffusion wave method was introduced into the XAJ-DCH model, which can improve the prediction accuracy in hydraulically mild slope area by accounting for the backwater effects.
2020, 48(4):377-384.
DOI: 10.3876/j.issn.1000-1980.2020.04.013
Abstract:
Firstly, the necessity of power electronization is explained from the development history of the power electronics and the challenges in power generation, transmission, distribution, and consumption. Secondly, the driving factors of power electronization, such as politics, economy, and technology, are elaborated in three aspects, including energy security and political construction, the resource and environmental benefits, and the transformation of power generation, transmission, distribution and consumption. Finally, the development trends and main restrictive factors of power electronization in power system are analyzed, and related study directions are summarized to provide some references and suggestions for future study.
Firstly, the necessity of power electronization is explained from the development history of the power electronics and the challenges in power generation, transmission, distribution, and consumption. Secondly, the driving factors of power electronization, such as politics, economy, and technology, are elaborated in three aspects, including energy security and political construction, the resource and environmental benefits, and the transformation of power generation, transmission, distribution and consumption. Finally, the development trends and main restrictive factors of power electronization in power system are analyzed, and related study directions are summarized to provide some references and suggestions for future study.
2022, 50(2):1-10.
DOI: 10.3876/j.issn.1000-1980.2022.02.001
Abstract:
In this study, a hydrological model parameter optimization scheme was constructed in SWAT based on streamflow assimilation using the ensemble Kalman filter(EnKF)approach, in which the model simulation and streamflow observation errors were reasonably quantified and the model parameter evolution and over-fitting issues were effectively treated. Based on this scheme, the capacity of streamflow assimilation on SWAT model parameter optimization was evaluated in the upper Huai River above Huaibin hydrological station. The results showed that the updated parameter ensemble gradually converged and stabilized during the data assimilation process. The runoff process based on the stabilized parameter ensemble was close to the measured runoff. The deterministic coefficient at the outlet of the basin reached 0. 88. The results indicated that the EnKF based streamflow assimilation had the ability to optimize the parameters in SWAT. It was feasible to calibrate model parameters using data assimilation method.
In this study, a hydrological model parameter optimization scheme was constructed in SWAT based on streamflow assimilation using the ensemble Kalman filter(EnKF)approach, in which the model simulation and streamflow observation errors were reasonably quantified and the model parameter evolution and over-fitting issues were effectively treated. Based on this scheme, the capacity of streamflow assimilation on SWAT model parameter optimization was evaluated in the upper Huai River above Huaibin hydrological station. The results showed that the updated parameter ensemble gradually converged and stabilized during the data assimilation process. The runoff process based on the stabilized parameter ensemble was close to the measured runoff. The deterministic coefficient at the outlet of the basin reached 0. 88. The results indicated that the EnKF based streamflow assimilation had the ability to optimize the parameters in SWAT. It was feasible to calibrate model parameters using data assimilation method.
2020, 48(2):102-108.
DOI: 10.3876/j.issn.1000-1980.2020.02.002
Abstract:
In order to obtain the characteristics of wind field and flow field as well as the influence of different types of wind fields on the flow field of Taihu Lake, 6 acoustic doppler current profilers(ADCP)and 4 meteorological stations were set up in different areas of Taihu Lake, and the synchronous data of flow field and wind field were obtained. The results show the floolwing: (a)The spatial distribution of wind field in Taihu Lake was not uniform, the wind speed varied greatly among 4 stations, and the difference of wind direction was expanded along with the decrease of wind speed. (b)The flow field in Taihu Lake had obvious characteristics of wind-generated flow. When the average wind speed was greater than 6. 7 m/s, there was a counterclockwise circulation that appeared in the western lake area, and there was no obvious flow stratification during this period. When the average wind speed was less than 3. 8 m/s, the circulation in the west fade away. However, the flow stratification was obviously developed, the current direction of each layer was relatively stable but different, and the compensation flow developed. (c)According to the observed data, the percentage of current speed less than 10 cm/s was more than 90%, and the current speed can reach 20 cm/s under strong wind conditions.
In order to obtain the characteristics of wind field and flow field as well as the influence of different types of wind fields on the flow field of Taihu Lake, 6 acoustic doppler current profilers(ADCP)and 4 meteorological stations were set up in different areas of Taihu Lake, and the synchronous data of flow field and wind field were obtained. The results show the floolwing: (a)The spatial distribution of wind field in Taihu Lake was not uniform, the wind speed varied greatly among 4 stations, and the difference of wind direction was expanded along with the decrease of wind speed. (b)The flow field in Taihu Lake had obvious characteristics of wind-generated flow. When the average wind speed was greater than 6. 7 m/s, there was a counterclockwise circulation that appeared in the western lake area, and there was no obvious flow stratification during this period. When the average wind speed was less than 3. 8 m/s, the circulation in the west fade away. However, the flow stratification was obviously developed, the current direction of each layer was relatively stable but different, and the compensation flow developed. (c)According to the observed data, the percentage of current speed less than 10 cm/s was more than 90%, and the current speed can reach 20 cm/s under strong wind conditions.
2015, 43(4):288-293.
DOI: 10.3876/j.issn.1000-1980.2015.04.002
Abstract:
Considering the hydrological calculation and prediction of ungauged basins, a distributed hydrological model, the THREW model, was adopted to study the upper basin of the Brahmaputra River from the Nuxia Hydrological Station. Surface meteorological observation, remote sensing of vegetation cover, snow cover area, and section hydrological monitoring data were used to calibrate the model, and the latest CMIP5 data were used to estimate the runoff evolution. The results showed that, in the upper basin of the Brahmaputra River from the Nuxia Hydrological Station, the Nash-Sutcliffe efficiency coefficient of monthly runoff simulation of the calibration period from 1991 to 1995 was 0. 75, the coefficient of monthly runoff simulation of the validation period from 1996 to 2000 was 0. 76, and the runoff markedly increased when the CO2 discharge in the AR5 of IPCC was maximal.
Considering the hydrological calculation and prediction of ungauged basins, a distributed hydrological model, the THREW model, was adopted to study the upper basin of the Brahmaputra River from the Nuxia Hydrological Station. Surface meteorological observation, remote sensing of vegetation cover, snow cover area, and section hydrological monitoring data were used to calibrate the model, and the latest CMIP5 data were used to estimate the runoff evolution. The results showed that, in the upper basin of the Brahmaputra River from the Nuxia Hydrological Station, the Nash-Sutcliffe efficiency coefficient of monthly runoff simulation of the calibration period from 1991 to 1995 was 0. 75, the coefficient of monthly runoff simulation of the validation period from 1996 to 2000 was 0. 76, and the runoff markedly increased when the CO2 discharge in the AR5 of IPCC was maximal.
总访问量： 今日访问：
Journal of Hohai University (Natural Sciences) ® 2026 All Rights ReservedAddress：No.1 Xikang Road, Nanjing 210098, P.R. ChinaPostcode：210098Tel：025-83786642
Website：https://jour.hhu.edu.cn/hhdxxbzren/home E-mail：xb1957@vip.163.com;xb@hhu.edu.cn
Supported by:Beijing E-Tiller Technology Development Co., Ltd.