CN108006863B - 一种混凝土辐射供冷系统负荷计算优化方法 - Google Patents

一种混凝土辐射供冷系统负荷计算优化方法 Download PDF

Info

Publication number
CN108006863B
CN108006863B CN201711153028.9A CN201711153028A CN108006863B CN 108006863 B CN108006863 B CN 108006863B CN 201711153028 A CN201711153028 A CN 201711153028A CN 108006863 B CN108006863 B CN 108006863B
Authority
CN
China
Prior art keywords
unit
concrete
heat storage
coefficient
thermal
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.)
Active
Application number
CN201711153028.9A
Other languages
English (en)
Other versions
CN108006863A (zh
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.)
Tianjin University
Original Assignee
Tianjin University
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 Tianjin University filed Critical Tianjin University
Priority to CN201711153028.9A priority Critical patent/CN108006863B/zh
Publication of CN108006863A publication Critical patent/CN108006863A/zh
Application granted granted Critical
Publication of CN108006863B publication Critical patent/CN108006863B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F24HEATING; RANGES; VENTILATING
    • F24FAIR-CONDITIONING; AIR-HUMIDIFICATION; VENTILATION; USE OF AIR CURRENTS FOR SCREENING
    • F24F5/00Air-conditioning systems or apparatus not covered by F24F1/00 or F24F3/00, e.g. using solar heat or combined with household units such as an oven or water heater
    • F24F5/0089Systems using radiation from walls or panels

Landscapes

  • Engineering & Computer Science (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Sustainable Development (AREA)
  • Chemical & Material Sciences (AREA)
  • Combustion & Propulsion (AREA)
  • Mechanical Engineering (AREA)
  • General Engineering & Computer Science (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Investigating Or Analyzing Materials Using Thermal Means (AREA)

Abstract

一种混凝土辐射供冷系统负荷计算优化方法,包括:设定混凝土辐射供冷末端的热工参数和几何参数,其中,热工参数包括混凝土的导热系数和混凝土的蓄热系数,几何参数包括埋管间距及埋管深度;分别计算:第一热阻、第二热阻、热惰性指标、第一外表面蓄热系数、第二外表面蓄热系数和内表面蓄热系数;分别计算基本衰减系数和基本延迟时间;分别计算在不同周期的外扰作用下的衰减系数和延迟时间;将一天24小时的逐时房间冷负荷写为傅里叶展开式;计算一天24小时之中混凝土辐射供冷系统的逐时冷负荷。本发明解决了传统冷负荷计算方法无法用于混凝土辐射供冷系统的缺陷,有利于减小系统峰值负荷,实现削峰填谷的节能效果。

Description

一种混凝土辐射供冷系统负荷计算优化方法
技术领域
本发明涉及一种混凝土辐射供冷系统负荷计算优化方法。特别是涉及一种当单独使用混凝土辐射供冷系统供冷时,基于混凝土辐射供冷末端的蓄热特性,对混凝土辐射供冷系统的逐时冷负荷进行计算时采用的混凝土辐射供冷系统负荷计算优化方法。
背景技术
在建筑节能的大环境下,辐射供冷空调系统以其节能性与舒适性的优点在工程和学术界受到越来越多的关注。辐射供冷空调系统的末端形式有很多种,按照供冷末端与室内换热的主要方式进行分类,主要包括两大类:对流式末端和辐射式末端。对流式末端包括对流式冷天花、被动式冷梁及主动式冷梁,辐射式末端包括辐射式冷天花、毛细管末端及混凝土辐射供冷末端。其中,混凝土辐射供冷末端作为新型节能的空调末端形式,其供冷能力与性能也逐渐得到了学者们的关注与研究。
与其它辐射供冷末端形式不同,混凝土辐射供冷末端将供冷盘管直接安装在建筑构件如楼板之中,由于建筑构件的热惰性较大,使得混凝土辐射供冷末端有着显著的蓄热能力。在外界扰量突变或剧烈变化时,混凝土辐射供冷末端不像其它形式末端一样在短时间内对外扰做出响应,存在时间上的延迟及变化幅度上的衰减。混凝土辐射供冷末端的这种蓄热特性使得在应用混凝土辐射供冷系统供冷时,对混凝土辐射供冷系统的冷负荷计算不应采用传统的冷负荷计算方法。当房间内单独使用混凝土辐射供冷系统供冷时,由于混凝土辐射供冷末端的蓄热特性,热量由室内空气向冷水机组内冷水的传递过程中存在衰减延迟效应,因此冷水机组所需提供的供冷量与房间冷负荷之间存在差异。即在对混凝土辐射供冷系统的逐时冷负荷进行计算时,应考虑混凝土辐射供冷末端的蓄热特性带来的影响,对根据传统计算方法计算得到的逐时房间冷负荷进行修正。这种优化的冷负荷计算方法有助于减小混凝土辐射供冷系统的峰值负荷,且可以实现削峰填谷的节能效果。
发明内容
本发明所要解决的技术问题是,提供一种有利于减小系统峰值负荷,实现削峰填谷节能效果的混凝土辐射供冷系统负荷计算优化方法。
本发明所采用的技术方案是:一种混凝土辐射供冷系统负荷计算优化方法,包括如下步骤:
1)设定混凝土辐射供冷末端的热工参数和几何参数,其中,热工参数包括混凝土的导热系数λ,单位W/m.K,及混凝土的蓄热系数S,单位W/m2.K;几何参数包括埋管间距L,单位m,及埋管深度L1,单位m;
2)根据混凝土辐射供冷末端的热工参数和几何参数分别计算:
第一热阻R1,单位m2.K/W:
第二热阻R2,单位m2.K/W:
Figure BDA0001473476130000021
热惰性指标D:D=S·(R1+R2)
第一外表面蓄热系数y1,单位W/m2.K:
Figure BDA0001473476130000022
第二外表面蓄热系数y2,单位W/m2.K:
Figure BDA0001473476130000023
内表面蓄热系数yn,单位W/m2.K:
Figure BDA0001473476130000024
3)根据混凝土的蓄热系数S、热惰性指标D、第一外表面蓄热系数y1、第二外表面蓄热系数y2和内表面蓄热系数yn分别计算:
基本衰减系数γ0
Figure BDA0001473476130000025
基本延迟时间ξ0,单位h:
Figure BDA0001473476130000026
4)根据基本衰减系数γ0和基本延迟时间ξ0,分别计算在不同周期的外扰作用下的衰减系数和延迟时间:
周期为24h外扰作用下的衰减系数γ1:γ1=0.3γ0+0.6
周期为24h外扰作用下的延迟时间ξ1,单位rad:ξ1=0.3ξ0-0.24
周期为12h外扰作用下的衰减系数γ2:γ2=0.84γ0-0.1
周期为12h外扰作用下的延迟时间ξ2,单位rad:ξ2=0.44ξ0-0.15
周期为8h外扰作用下的衰减系数γ3:γ3=1.43γ0-0.8
周期为8h外扰作用下的延迟时间ξ3,单位rad:ξ3=0.51ξ0
5)已知一天24小时的逐时房间冷负荷Q(单位kW),将Q写为傅里叶展开式:
Figure BDA0001473476130000027
其中τ=1,2,3...24,表示一天24小时;A0、A1、A2和A3分别表示不同的系数;
Figure BDA0001473476130000028
Figure BDA0001473476130000029
分别表示不同的相位角,单位为rad;
6)根据步骤4)中在不同周期的外扰作用下的衰减系数和延迟时间和步骤5)中的傅里叶展开式计算一天24小时之中混凝土辐射供冷系统的逐时冷负荷q,单位kW:
Figure BDA00014734761300000210
其中τ=1,2,3...24,表示一天24小时;A0、A1、A2和A3分别表示不同的系数;
Figure BDA00014734761300000211
Figure BDA00014734761300000212
分别表示不同的相位角。
本发明的一种混凝土辐射供冷系统负荷计算优化方法,解决了传统冷负荷计算方法无法用于混凝土辐射供冷系统的缺陷,有利于减小系统峰值负荷,实现削峰填谷的节能效果。本发明具有以下优点:
1、提出了利用公式计算混凝土辐射供冷系统的逐时冷负荷的方法。
2、该方法考虑了混凝土辐射供冷系统冷负荷与房间冷负荷之间的差异,使得混凝土辐射供冷系统的逐时冷负荷计算更加准确。
3、在应用混凝土辐射供冷系统时,使用该方法所得的逐时冷负荷数据进行运行指导,有利于减小混凝土辐射供冷系统的峰值负荷,实现削峰填谷的节能效果。
附图说明
图1是本发明中混凝土辐射供冷末端的热工参数和几何参数示意图。
图中:
L:是埋管间距单位m L1:是埋管深度,单位m
1:是隔热层 2:冷水
3:冷水管 4:是混凝土
具体实施方式
下面结合实施例和附图对本发明的一种混凝土辐射供冷系统负荷计算优化方法做出详细说明。
本发明的一种混凝土辐射供冷系统负荷计算优化方法,包括如下步骤:
1)设定混凝土辐射供冷末端的热工参数和几何参数,如图1所示,其中,热工参数包括混凝土的导热系数λ,单位W/m.K,及混凝土的蓄热系数S,单位W/m2.K;几何参数包括埋管间距L,单位m,及埋管深度L1,单位m;
2)根据混凝土辐射供冷末端的热工参数和几何参数分别计算:
第一热阻R1,单位m2.K/W:
Figure BDA0001473476130000031
第二热阻R2,单位m2.K/W:
Figure BDA0001473476130000032
热惰性指标D:D=S·(R1+R2)
第一外表面蓄热系数y1,单位W/m2.K:
Figure BDA0001473476130000033
第二外表面蓄热系数y2,单位W/m2.K:
Figure BDA0001473476130000034
内表面蓄热系数yn,单位W/m2.K:
Figure BDA0001473476130000035
3)根据混凝土的蓄热系数S、热惰性指标D、第一外表面蓄热系数y1、第二外表面蓄热系数y2和内表面蓄热系数yn分别计算:
基本衰减系数γ0
Figure BDA0001473476130000041
基本延迟时间ξ0,单位h:
Figure BDA0001473476130000042
4)根据基本衰减系数γ0和基本延迟时间ξ0,分别计算在不同周期的外扰作用下的衰减系数和延迟时间:
周期为24h外扰作用下的衰减系数γ1:γ1=0.3γ0+0.6
周期为24h外扰作用下的延迟时间ξ1,单位rad:ξ1=0.3ξ0-0.24
周期为12h外扰作用下的衰减系数γ2:γ2=0.84γ0-0.1
周期为12h外扰作用下的延迟时间ξ2,单位rad:ξ2=0.44ξ0-0.15
周期为8h外扰作用下的衰减系数γ3:γ3=1.43γ0-0.8
周期为8h外扰作用下的延迟时间ξ3,单位rad:ξ3=0.51ξ0
5)已知一天24小时的逐时房间冷负荷Q(单位kW),将Q写为傅里叶展开式:
Figure BDA0001473476130000043
其中τ=1,2,3...24,表示一天24小时;A0、A1、A2和A3分别表示不同的系数;
Figure BDA0001473476130000045
分别表示不同的相位角,单位为rad;因此,根据上述Q的傅里叶展开式可以得到A0、A1、A2和A3以及
Figure BDA0001473476130000046
Figure BDA0001473476130000047
6)根据步骤4)中在不同周期的外扰作用下的衰减系数和延迟时间和步骤5)中的傅里叶展开式计算一天24小时之中混凝土辐射供冷系统的逐时冷负荷q,单位kW:
Figure BDA0001473476130000048
其中τ=1,2,3...24,表示一天24小时;A0、A1、A2和A3分别表示不同的系数;
Figure BDA0001473476130000049
Figure BDA00014734761300000410
分别表示不同的相位角。
实例:
1)设定一种混凝土辐射供冷末端的热工参数和几何参数,其中,热工参数包括混凝土的导热系数λ=1.74,单位W/m.K,及混凝土的蓄热系数S=17.2,单位W/m2.K;几何参数包括埋管间距L=0.2,单位m,及埋管深度L1=0.06,单位m;
2)根据混凝土辐射供冷末端的热工参数和几何参数分别计算得到:
第一热阻R1,单位m2.K/W:
Figure BDA00014734761300000411
第二热阻R2,单位m2.K/W:
Figure BDA00014734761300000412
热惰性指标D:D=S·(R1+R2)=1.087
第一外表面蓄热系数y1,单位W/m2.K:
Figure BDA0001473476130000051
第二外表面蓄热系数y2,单位W/m2.K:
Figure BDA0001473476130000052
内表面蓄热系数yn,单位W/m2.K:
Figure BDA0001473476130000053
3)根据混凝土的蓄热系数S、热惰性指标D、第一外表面蓄热系数y1、第二外表面蓄热系数y2和内表面蓄热系数yn分别计算得到:
基本衰减系数γ0
Figure BDA0001473476130000054
基本延迟时间ξ0,单位h:
Figure BDA0001473476130000055
4)根据基本衰减系数γ0和基本延迟时间ξ0,分别计算得到在不同周期的外扰作用下的衰减系数和延迟时间:
周期为24h外扰作用下的衰减系数γ1:γ1=0.3γ0+0.6=1.04
周期为24h外扰作用下的延迟时间ξ1,单位rad:ξ1=0.3ξ0-0.24=0.444
周期为12h外扰作用下的衰减系数γ2:γ2=0.84γ0-0.1=1.133
周期为12h外扰作用下的延迟时间ξ2,单位rad:ξ2=0.44ξ0-0.15=0.853
周期为8h外扰作用下的衰减系数γ3:γ3=1.43γ0-0.8=1.299
周期为8h外扰作用下的延迟时间ξ3,单位rad:ξ3=0.51ξ0=1.162
5)已知一天24小时的逐时房间冷负荷Q(单位kW)如表1
表1
τ(h) Q(kW) τ(h) Q(kW) τ(h) Q(kW)
1 21.2 9 19.0 17 31.1
2 20.5 10 20.5 18 28.7
3 20.1 11 20.6 19 26.3
4 19.9 12 23.3 20 24.0
5 22.0 13 24.7 21 22.2
6 24.4 14 26.0 22 21.1
7 23.1 15 27.0 23 21.8
8 21.7 16 28.1 24 21.6
将Q写为傅里叶展开式:
6)根据步骤4)中在不同周期的外扰作用下的衰减系数和延迟时间和步骤5)中的傅里叶展开式得到一天24小时之中混凝土辐射供冷系统的逐时冷负荷q(单位kW)的表达式:
Figure BDA0001473476130000062
将τ代入表达式,求得一天24小时之中混凝土辐射供冷系统的逐时冷负荷q(单位kW),如表2。
表2
τ(h) q(kW) τ(h) q(kW) τ(h) q(kW)
1 21.10 9 21.16 17 28.53
2 20.63 10 20.24 18 28.85
3 20.52 11 19.62 19 28.24
4 20.78 12 19.75 20 27.01
5 21.32 13 20.85 21 25.52
6 21.87 14 22.80 22 24.09
7 22.13 15 25.11 23 22.85
8 21.89 16 27.20 24 21.84

Claims (1)

1.一种混凝土辐射供冷系统负荷计算优化方法,其特征在于,包括如下步骤:
1)设定混凝土辐射供冷末端的热工参数和几何参数,其中,热工参数包括混凝土的导热系数λ,单位W/m·K,及混凝土的蓄热系数S,单位W/m2.K;几何参数包括埋管间距L,单位m,及埋管深度L1,单位m;
2)根据混凝土辐射供冷末端的热工参数和几何参数分别计算:
第一热阻R1,单位m2.K/W:
Figure FDA0002266212960000011
第二热阻R2,单位m2.K/W:
Figure FDA0002266212960000012
热惰性指标D:D=S·(R1+R2)
第一外表面蓄热系数y1,单位W/m2.K:
Figure FDA0002266212960000013
第二外表面蓄热系数y2,单位W/m2.K:
Figure FDA0002266212960000014
内表面蓄热系数yn,单位W/m2.K:
Figure FDA0002266212960000015
3)根据混凝土的蓄热系数S、热惰性指标D、第一外表面蓄热系数y1、第二外表面蓄热系数y2和内表面蓄热系数yn分别计算:
基本衰减系数γ0
Figure FDA0002266212960000016
基本延迟时间ξ0,单位h:
Figure FDA0002266212960000017
4)根据基本衰减系数γ0和基本延迟时间ξ0,分别计算在不同周期的外扰作用下的衰减系数和延迟时间:
周期为24h外扰作用下的衰减系数γ1:γ1=0.3γ0+0.6
周期为24h外扰作用下的延迟时间ξ1,单位h:ξ1=0.3ξ0-0.24
周期为12h外扰作用下的衰减系数γ2:γ2=0.84γ0-0.1
周期为12h外扰作用下的延迟时间ξ2,单位h:ξ2=0.44ξ0-0.15
周期为8h外扰作用下的衰减系数γ3:γ3=1.43γ0-0.8
周期为8h外扰作用下的延迟时间ξ3,单位h:ξ3=0.51ξ0
5)已知一天24小时的逐时房间冷负荷Q,单位kW,将Q写为傅里叶展开式:
Figure FDA0002266212960000021
其中τ=1,2,3...24,表示一天24小时;A0、A1、A2和A3分别表示不同的系数;
Figure FDA0002266212960000022
Figure FDA0002266212960000023
分别表示不同的相位角,单位为rad;
6)根据步骤4)中在不同周期的外扰作用下的衰减系数和延迟时间和步骤5)中的傅里叶展开式计算一天24小时之中混凝土辐射供冷系统的逐时冷负荷q,单位kW:
其中τ=1,2,3...24,表示一天24小时;A0、A1、A2和A3分别表示不同的系数;
分别表示不同的相位角。
CN201711153028.9A 2017-11-19 2017-11-19 一种混凝土辐射供冷系统负荷计算优化方法 Active CN108006863B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711153028.9A CN108006863B (zh) 2017-11-19 2017-11-19 一种混凝土辐射供冷系统负荷计算优化方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711153028.9A CN108006863B (zh) 2017-11-19 2017-11-19 一种混凝土辐射供冷系统负荷计算优化方法

Publications (2)

Publication Number Publication Date
CN108006863A CN108006863A (zh) 2018-05-08
CN108006863B true CN108006863B (zh) 2020-01-31

Family

ID=62053111

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711153028.9A Active CN108006863B (zh) 2017-11-19 2017-11-19 一种混凝土辐射供冷系统负荷计算优化方法

Country Status (1)

Country Link
CN (1) CN108006863B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111780377B (zh) * 2020-06-16 2021-07-09 珠海格力电器股份有限公司 一种风机转速控制方法、装置及空调设备

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102980272A (zh) * 2012-12-08 2013-03-20 珠海派诺科技股份有限公司 一种基于负荷预测的空调系统节能优化方法
CN103104509A (zh) * 2013-02-25 2013-05-15 天津大学 变频水泵全工况运行状态获取方法
CN107062383A (zh) * 2017-04-27 2017-08-18 华电电力科学研究院 一种考虑建筑物热惰性及太阳辐射影响的实时负荷计算方法
CN206488398U (zh) * 2016-12-22 2017-09-12 长安大学 一种tabs辐射供冷装置及其防结露控制系统

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2010286204A (ja) * 2009-06-15 2010-12-24 Haishima Kogyo Kk 家屋における年間冷・暖房エネルギー消費量低減方法及び外皮内の安全性を確保する方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102980272A (zh) * 2012-12-08 2013-03-20 珠海派诺科技股份有限公司 一种基于负荷预测的空调系统节能优化方法
CN103104509A (zh) * 2013-02-25 2013-05-15 天津大学 变频水泵全工况运行状态获取方法
CN206488398U (zh) * 2016-12-22 2017-09-12 长安大学 一种tabs辐射供冷装置及其防结露控制系统
CN107062383A (zh) * 2017-04-27 2017-08-18 华电电力科学研究院 一种考虑建筑物热惰性及太阳辐射影响的实时负荷计算方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
基于频域回归法的混凝土辐射供冷楼板非稳态传热模型的建立及验证;田喆等;《天津大学学报(自然科学与工程技术版)》;20151124;第49卷(第8期);第848-854页 *
混凝土辐射供冷楼板热动态响应特性的研究;雒志明;《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》;20171115;C038-157 *

Also Published As

Publication number Publication date
CN108006863A (zh) 2018-05-08

Similar Documents

Publication Publication Date Title
Shen et al. Dynamic thermal performance of pipe-embedded building envelope utilizing evaporative cooling water in the cooling season
Li et al. A study on pipe-embedded wall integrated with ground source-coupled heat exchanger for enhanced building energy efficiency in diverse climate regions
Shen et al. Energy saving potential of pipe-embedded building envelope utilizing low-temperature hot water in the heating season
Zhai et al. Experimental investigation and performance analysis of a fin tube phase change cold storage unit for high temperature cooling application
Xie et al. An active pipe-embedded building envelope for utilizing low-grade energy sources
CN103438611B (zh) 一种太阳能地源热泵系统优化设计方法
Zhang et al. Simulation on the performance and free cooling potential of the thermosyphon mode in an integrated system of mechanical refrigeration and thermosyphon
CN101588701B (zh) 机柜的温控方法和系统
Yan et al. Model validation and application of the coupled system of pipe-encapsulated PCM wall and nocturnal sky radiator
CN108006863B (zh) 一种混凝土辐射供冷系统负荷计算优化方法
CN111829059B (zh) 一种供热系统动态建模方法、模型及调控系统
Allerhand et al. Energy and thermal comfort performance evaluation of PCM ceiling panels for cooling a renovated office room
Yang et al. Experimental study on the thermal performances of a non-pump-driven thermo-activated building system based on flat-plate heat pipe array
Yan et al. Parametric analysis on performances of the pipe-encapsulated PCM (PenPCM) wall system coupled with gravity heat-pipe and nocturnal radiant cooler
CN202282922U (zh) 一种节能散热机柜
Gao et al. Modeling the airside dynamic behavior of a heat exchanger under frosting conditions
CN102901270A (zh) 一种冷媒直接蒸发式地埋管换热系统
Liu et al. 2D Axisymmetric Model Research of Helical Heat Exchanger inside Pile Foundations
Yu et al. Influence of intermittent operation on soil temperature and energy storage duration of ground-source heat pump system for residential building
Huang et al. Heat and mass transfer of evaporative wall and its cooling load calculation
CN202209819U (zh) 一种冷媒直接蒸发式地埋管换热系统
CN112984602B (zh) 大空间建筑在太阳辐射下地板辐射供暖系统及分区方法
CN203163546U (zh) 水源换热器及水源热泵空调
CN213334752U (zh) 一种建筑墙面节能系统
Zhang et al. Design and application of a seasonal solar soil heat storage system applied in greenhouse heating

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
CP02 Change in the address of a patent holder
CP02 Change in the address of a patent holder

Address after: 300452 Binhai Industrial Research Institute Campus of Tianjin University, No. 48 Jialingjiang Road, Binhai New Area, Tianjin

Patentee after: Tianjin University

Address before: 300192 92 Weijin Road, Nankai District, Tianjin

Patentee before: Tianjin University