人工岛工程对河口行洪冲淤的影响分析

盛天航，孙冬梅，张 杨

（天津大学 水利工程仿真与安全国家重点实验室，天津300072）

人工岛工程对河口行洪冲淤的影响分析

盛天航，孙冬梅，张 杨

（天津大学 水利工程仿真与安全国家重点实验室，天津300072）

以秦皇岛汤河河口人工岛在建工程为背景，采用大、小两重模型嵌套的方式，建立二维平面水动力-泥沙数学模型。在验证模型可行性的基础上，分析了50 a一遇设计洪水条件下工程区的水动力及泥沙冲淤特性，并进一步阐述了人工岛工程建设对工程区流场变化规律、河口冲淤规律以及工程区附近泥沙分布规律的影响。研究表明，在实施河道清淤的情况下，人工岛的建立并不会增加河口工程区附近的泥沙淤积，相反在清淤工程和梭形岛分流作用的影响下，泥沙淤积量有所减少。

人工岛；河口；数学模型；泥沙淤积；冲淤

近年来，随着经济建设的发展以及人口的快速增长，陆地资源逐渐枯竭，人均耕地大大减少，土地的供需矛盾日益尖锐［1］，各类用地已无法满足人们的需求，人类对海洋空间的探索不断增加，人工岛建设工程也日益增多［2］。建造人工岛不仅可以有效地开发利用海域资源来缓解陆地资源紧缺所带来的压力，还可以拓展人类生存空间以维持社会的全面、协调、可持续发展。人工岛在世界各沿海国家都极为普遍，诸如美国、日本、加拿大等，在人工岛的建设方面都积累了较为丰富的经验。然而人工岛的建造必然会改变河口工程区及其周边的水动力特性、泥沙输运及海岸演变的状况［3］，进而引发河口海岸、岸滩等冲淤特性的变化。因此，在人工岛建设之前，应充分评估其对周边河口海岸的影响，针对人工岛附近海域的水动力和泥沙特性以及海岸河口泥沙输运规律进行相应的分析研究，做出合理的规划和平面布设方案，最大程度减少人工岛工程对河口海岸环境的负面影响。

利用数值模拟技术对工程造成的损失进行评估和预测，对于减少工程在环境方面的负面影响和工程失误具有重要意义［4］，本文以汤河河口人工岛工程为背景，构建二维水动力-泥沙数学模型，模拟了在50 a一遇设计洪水条件下人工岛工程前后汤河河口附近的水动力和泥沙运动过程，并基于模拟结果分析人工岛对汤河河口次洪冲淤的影响。

1 研究区概况

1.1 人工岛概况

人工岛工程拟建在秦皇岛北戴河区东北、海港区西南海域，位于大、小两条汤河入海口。工程采用取围填海的方式进行建岛，总体呈海螺状，面积约为594 618.14 m3，包括防波堤、围堤、吹填造陆区以及游艇水域，北侧距离帆船基地约350 m，东侧距离西港区防波堤约150 m。人工岛工程区附近海域包括2个潮位观测站（H1、H2）和8个潮流观测站（V1～V8），位置见图1所示。

1.2 自然条件概况

大、小汤河均属于冀东沿海诸河水系，是典型的山溪性河流，源短流急，峰高量大，历时短，枯水季节流量很小。大汤河是秦皇岛海港区西部最大的河流，流域总面积约237.7 km3，长约30 km，流域内地势北高南低，平均坡度3.7%。大汤河主河道长约14 km，平均纵坡度1/500，上段为砂砾石河床，纵坡较陡，下段为沙质河床，纵坡较缓；小汤河是大汤河的分支，自北向南汇入大汤河并最终注入渤海。小汤河河堤全长12.2 km，流域面积53.75 km2，流域内地势北高南低，处于海滨冲击平原地带，河床岩性为细砂土和亚粘土。

研究区域位于大、小汤河入海口，属暖温带半湿润大陆性季风气候，受海洋影响较大，气温较温和，年平均气温10.2℃，7月份最高平均气温24.7℃，1月份最低平均气温-6.3℃。研究区降水丰富，多年平均降雨量为658.0 mm，主要集中在6～8月，占全年降水的70%以上降雨。此外，研究区依山傍海，附近海域受到潮汐影响，以潮力弱、潮差小为其显著特点。研究区汛期（6～9月）多年平均高潮位1.121 m，平均潮差0.74 m，最大潮差2.43。

图1 人工岛地理位置Fig.1 Location of artificial island

人工岛工程区海岸带岬湾相间，岸线总体平顺，环海寺地咀和金山咀岬角规模较大，其他岬角向海突出有限，形成了发育程度不同的弧形海湾［5］。汤河口至沙河口沿岸基本上为人工岸线，码头西防波堤堤根处存在明显的沙质堆积体。金山咀至汤河口岸段属弧形海岸，岸滩平缓，潮间带宽10～15 m，岸滩沉积物主要为中细沙。该段海岸中部高潮线下发育多条相互平行的、自东向西南的顺岸沙坝，宽度逐渐增大，坡度逐渐变缓，显示出泥沙从东向西运移的趋势。根据该工程区海域各水文泥沙测站（见图1）观测资料，实测期间（2011年2～3月）工程区附近海域整体含沙量较低，平均含沙量介于0.008～0.025 kg/m3之间。工程区附近海域各垂线涨、落潮含沙量，大小潮各站含沙量，在垂线分布上的表层、中层和底层的含沙量以及在平面分布上自近岸向外的含沙量差别均不大。此外，按照、大、小潮期采集的悬沙水样粒径进行分析，工程区悬沙中值粒径介于0.008 7～0.011 2 mm之间，其中大潮平均中值粒径为0.009 9 mm，小潮平均中值粒径为0.012 2 mm，悬沙物质成分为粘土质粉砂。

2 二维水动力-泥沙数学模型的建立

工程区位于河口沿岸浅海地区，研究仅考虑了河口行洪、潮汐对工程区附近泥沙冲淤的影响，工程区内水平尺度远大于垂向尺度，可忽略流速、水深等水力参数在垂直方向的变化，并假设沿水深方向动水压强符合静水压强分布，将三维流动的基本方程沿垂直方向积分得到平面二维流动的基本方程，采用平面二维浅水方程［6-7］作为水动力模型控制方程。此外，天然河流在运动过程中推移质、悬移质和床沙之间会引起河床冲淤变化，可采用上述平面二维水动力模型耦合泥沙输移模型进行泥沙冲淤的数值模拟和分析。本次采用泥沙输运方程［8］作为控制方程与二维平面浅水方程进行耦合，泥沙控制方程主要包括：泥沙输移控制方程、悬沙运动造成的海床冲淤以及底沙（推移质）造成的海床冲淤数学模型等。

2.1 基本控制方程

二维浅水方程守恒型通式

底沙造成海床冲淤数学方程

式中：U为守恒量向量，U={h，hu，hv}T；F=[Fx，Fy]为通量向量，且x方向向量Fx={hu，hu2+gh2/2,huv}T，y方向向量Fy={hv，huv，hv2+gh2/2}T；W为源汇项向量，表达式为

式中：h表示水深；t为时间；g为重力加速度常数；u和v分别为沿x和y方向垂线平均速度；f为柯氏力系数；Ex，Ey分别为x，y方向的涡粘系数；C为谢才系数；n为曼宁粗糙系数；式（2）～（4）中：S为铅直方向积分的水体含沙浓度；Dx，Dy分别为x，y方向泥沙扩散系数；Fs=αω（S-S*）为泥沙源汇函数，且S*为水体挟沙力，ω为泥沙沉降速度，α为泥沙沉降几率；γ0为床面泥沙干容重；γb为床面底沙干容重；ηs为海底床面悬沙引起的冲淤厚度；ηb为沙底引起的海床冲淤厚度；qx、qy分别为沿x、y方向单位时间内单宽底沙输移量分量。

2.2 数学方程求解

水动力数学模型采用有限体积法（FVM），以守恒型二维浅水方程式（1）为出发点，将方程在任意控制体积作体积分，并利用Gauss原理将体积分化成面积分，重写方程组（1）为

对控制体单元取平均后，得到FVM半离散化方程为

由于Fx、Fy具有旋转不变性，可将Fx（U）、Fy（U）在法向上的投影转换为先投影U到法向上，即将U投影到n得到，且，再将其代入F得到F(Uˉ)，并根据公式进一步求得Fn（U），从而将二维法向量计算问题转化为一维局部坐标下的黎曼问题进行求解。经旋转变换之后的FVM体积法半离散化方程式为

式中：C0为无量纲谢才系数；γ、γs为水容重和泥沙颗粒容重；ωb为底沙颗粒沉降速度；为水流平均速度矢量；为底沙颗粒临界启动流速；当|时，当|时，m=0；泥沙数学模型初始条件如式

式中：S0（x，y，t0）为初始时刻t0的已知值；模型计算水域与陆地交界的固边界Γ1上有

计算水域与陆地交界的固边界Γ2上有

图2 模型计算范围Fig.2 Model computational domain

式中：S*（x，y，t）为实测或分析计算所得已知值；n为陆地边界单位法向量矢量，该式表述泥沙沿固边界的法向量通量为零。

2.3 模型范围及网格

为保证局部流场计算符合潮流场的整体物理特征，采用大、小两重模型嵌套进行计算（图2）。其中大模型包含整个渤海、开边界位于大连老虎滩和烟台两个潮位观测站连线上，小模型是以工程区为中心、沿岸宽33 km的扇形区域，小模型潮位开边界条件由大模型提供，开边界设置在大、小汤河入河口处。

大、小模型控制方程使用有限体积法进行离散，采用三角形网格，很好地拟合了复杂岸线和建筑物边界（图3）。大模型共计网格节点2 283个，计算时间步长10 s，小模型网格节点6 153个，最大空间步长1 000 m，最小空间步长10 m，计算时间步长为0.5 s。小模型计算时底部糙率取值0.012～0.022之间。

图3 计算网格Fig.3 Computational grid

2.4 模型的验证

采用2011年2月25日9：00～26日13：00大潮和3月5日8：00～6日12：00小潮观测资料对模型潮位、流速、流向、含沙量和冲淤进行了验证［13］。

图4 大潮潮位验证Fig.4 Spring tide level validation

图5 大潮流速流向验证Fig.5 Spring tide velocity and flow direction validation

2.4.1 潮流验证

为了验证潮流模型的合理性，分别对图1中H1、H2潮位点和V1～V8潮流点进行潮位及流速、流向的验证，限于篇幅，仅给出2个大潮点的验证情况，如图4至图5所示。

图6 大潮含沙量验证Fig.6 Spring tide sediment concentration validation

图7 D1断面水深变化图Fig.7 Water depth of D1 section

图8 泥沙冲淤验证Fig.8 Silt erosion and deposition validation

由实测与计算结果比较可见，各测点计算值与实测潮位、流速、流向在连续的变化过程中都比较接近，所建潮流模型可以比较全面地反应工程区附近海域的流动规律，可以进一步为分析工程后流场和泥沙运动模拟提供必要的水动力条件。

2.4.2 泥沙验证

含沙量观测资料与潮流同步，图6为大潮V1、V2测点实测含沙量与计算结果比较情况。

验证结果表明，各测站计算与实测含沙量在连续变化过程中都比较接近，位相基本相符，工程水域内绝大多数测点的含沙量结果符合交通运输部《海岸与河口潮流泥沙模拟技术规程》的要求。因此，所建立的泥沙运动模型能够合理反应工程区附近海域的含沙量分布，可以进一步用于冲淤数学模型的验证。

2.4.3 冲淤验证

因本区域缺少用于地形冲淤验证的试挖坑回淤资料，且对于短期施工造成的地形冲淤变化，因地形资料、水动力条件及工程进度情况的缺失，也很难进行模拟。但为了保证所建立的模型在冲淤性质及地形变化量级上的准确性，根据大范围海图等深线对比分析结果，选取人工岛所在海域附近的D1断面2003～2009年水深变化进行地形冲淤验证（见图7）。在此期间，该断面附近无大型工程建设，地形冲淤变化可以只考虑自然动力因素的影响，适于作为地形冲淤验证的资料，冲淤验证结果见图8所示。

验证结果表明，计算值能够真实反应D1断面地形淤积变化的整体趋势，量级一致。因此，所建立的地形冲淤数学模型是合理的，可以进一步用于分析工程后地形冲淤变化。

图9 人工岛嵌入位置示意图Fig.9 Sketch of artificial island location

3 人工岛对行洪冲淤的影响分析

人工岛工程必然迫使周边流场发生变化，进而影响泥沙的运动，导致海床地形的冲淤变化。利用验证合理的水动力模型，在工程区嵌入人工岛（见图9）进行工程前、后水动力过程及冲淤过程的模拟分析。模型计算条件根据汤河控制流域范围，按照《秦皇岛市水文手册》推荐的小流域设计洪水经验公式进行50 a一遇设计洪峰的推求，并采用洪峰放大倍比对典型洪水过程进行放大，得到大、小汤河50 a一遇设计洪水过程线（见图10）作为模型入流条件。工程区开边界为50 a一遇的设计高潮位1.601 m；河口平均含沙量0.015 kg/m3，多年平均来水含沙量0.16 kg/m3。工程区所在海域的悬沙中值粒径为0.01 mm，泥沙发生絮凝沉降［14-15］，相应海水的悬沙沉速取平均值0.04 cm/s，并作为本次预测过程中的泥沙沉速［13］。

图10 设计洪水过程Fig.10 Designed flood hydrograph

3.1 清淤工程

根据拟建人工岛项目规划，为了满足游艇泊位及航道水深要求，计划对老京山铁路桥以下段河道进行航道清淤工程，如图11所示。其中大、小汤河汇合口段以上段两支河道河底高程清至-3.0 m，汇合口以下至海螺岛周边河底高程清至-4.5 m以下。为保证京山铁路桥安全，桥梁下游150 m范围内不清淤。

3.2 雍水及流态分析

工程建成后，在考虑到汤河河口及人工岛附近清淤工程的基础上，对工程区内水位及水流流态分布变化情况进行分析。工程前后汤河河道水面线变化如图12所示。

图11 大、小汤河河底高程Fig.11 Bed elevation of the Tang river

图12 大、小汤河河道水面线Fig.12 Water surface profile of the Tang river

可见，工程建设后，大、小汤河河道及入海口附近均会产生一定的雍水，但由于清淤工程的实施，泄水能力较现状河道明显加强，大、小汤河河道附近的水位与工程前水位相比明显下降。

人工岛具有一定的阻水作用，致使工程区附近水位雍高，但由于清淤工程，雍高幅度并不大（见图13，最大雍高水位仅为0.03 m）。此外，人工岛具有分流作用，致使水流流态和流速发生变化（见图14），变化范围在-0.245～0.38 m/s，主要表现为入海洪水遇岛分成东、西两支水流，岛西与海堤之间为主要的水流通道，岛东由于受到灯塔阻碍作用，在灯塔北侧流速减小，南侧流速加强。

图13 50 a一遇洪水洪峰时刻水位分布图Fig.13 Peak flood water level distribution under the condition of 50 a frequency design flood

3.3 冲淤分析

大、小汤河河道段在工程前后均以冲刷为主。工程前，大汤河最大冲刷厚度为0.306 m，位于大、小汤河汇合口（起点距1 127 m），最大流速3.134 m/s，小汤河最大次洪冲刷厚度0.334 m，最大流速3.134 m/s。小汤河以下至汇合口的冲刷量逐渐增大，在汇合口700 m范围内仍以冲刷为主，之后进入海域，流速迅速下降，开始淤积，淤积厚度在0.12～0.14 m；工程后，大汤河最大冲刷厚度0.481 m，最大流速4.877 m/s，小汤河最大冲刷厚度0.518 m，最大流速5.276 m/s，且越靠近汇合口，冲刷量越小。从汇合口延至350 m范围内，大、小汤河水流交汇，流速仍然较大，越靠近海域，淤积量越大，人工岛北端淤积厚度0.067 m，最大流速0.5 m/s（见图15）。由工程前后大、小汤河河道沿程的冲淤量变化分析，工程前后汤河及河口附近冲淤变化比较明显，受到下游河道清淤的影响，河道上游部分工程后冲刷量较工程前有所增大，而在河道下游部分工程后冲刷量明显小于工程前。

图14 50 a一遇洪水流场分布图Fig.14 Peak flood flow distribution under the condition of 50 a frequency design flood

对人工岛周边洪水过程冲淤变化情况进行分析，工程前、后次洪过程冲淤厚度等值线分布见图16所示。工程前，洪水进入海域后受到海水顶托作用，流速减少，泥沙淤积较多，主要沿西侧的海堤流动，并随距离增加，泥沙淤积量逐渐减小；工程后，由于人工岛的阻碍作用，河口附近流速减小，泥沙开始淤积。但由于清淤工程的影响，泥沙淤积较工程前有所减少。加之人工岛北端呈梭状，在河口处起到了一定的分流作用，致使岛西、岛东两侧流速增加，也在一定程度上减少了岛两侧泥沙的淤积量。此外，岛东侧西港区存在灯塔，对淤积产生一定影响，岛东侧水流受到灯塔阻碍作用，淤积量较多，而在水流通过灯塔后流速增大，淤积量减少。

图15 河道沿程的冲淤量Fig.15 Sediment erosion-accumulation along the river

对比工程前、后工程区附近冲淤量分析可知，工程后，泥沙淤积量并未因人工岛阻水作用而显著增加，相反在清淤工程及人工岛分流作用的影响下，岛两侧泥沙淤积量较工程前有所减少，最大淤积厚减少量为0.049 m，且工程前后工程区附近的泥沙淤积分布规律变化较大，工程前泥沙淤积主要发生在河口附近，工程后人工岛北端、东北侧及西南侧则淤积较多。这主要是由于局部流态的变化，改变了水流的挟沙能力，从而导致次洪过程工程区附近的淤积规律发生了一定的变化。

图16 冲淤厚度等值线分布图Fig.16 Sediment erosion-accumulation distribution

4 结语

本文以汤河河口拟建的人工岛工程为背景，通过构建二维水动力-泥沙数学模型，对50 a一遇设计洪水条件下工程区的水动力及泥沙冲淤特性进行了模拟，并将模型应用于人工岛工程，分析了工程建设对工程区水动力及泥沙冲淤变化的影响。结果表明构建的二维水动力-泥沙模型基本反映了工程区内潮流运动特点及含沙量的分布情况，模拟精度符合要求，可利用模型进行人工岛工程建设对周边海域影响的分析。从人工岛工程所引起的水位、流场变化及次洪冲淤的变化情况可知，工程后工程区附近雍水较高，流速和流态均发生变化。但在河道清淤工程的影响下，汤河及河口附近冲淤变化比较明显，工程后冲刷量小于工程前，冲刷主要发生在河道中泓部分；此外，工程实施后，在清淤工程及人工岛分流作用的共同影响下，人工岛附近泥沙淤积量较工程前有所减少，人工岛工程不会对河口冲淤造成影响，人工岛的修建并未引起泥沙淤积量的增加。人工岛北端、东北侧及西南侧淤积较多，人工岛东侧水流受到灯塔阻碍作用，也会导致泥沙淤积的增加。因此为了两岸堤防的稳定安全和完整，冲淤工程在实施过程中应尽量注意冲刷较大河段的堤防稳定，必要时增加适当防护工程。防波堤北端及西侧要充分考虑冲淤的影响，采取相应工程措施。灯塔对泥沙淤积具有一定影响，在工程建设中，也可考虑去掉灯塔，以有效减少人工岛东侧及西港区之间的淤积量。

［1］姜志德，姜爱林.中国土地资源现状及可持续利用对策研究［J］.华中理工大学学报（社会科学版），1998（3）：61-64. JIANG Z D，JIANG A L.Study on the present situation and countermeasures of land resources sustainable utilization in China［J］. Journal of HUST：Social Science Edition，1998（3）：61-64.

［2］林元军，吴家鸣.人工岛工程建设对海洋环境影响数值模拟分析方法探讨［J］.广东造船，2008（4）：35-37. LIN Y J，WU J M.Numerical methods for analyzing influence of artificial island project on marine environment［J］.Guangdong Shipbuilding magazine，2008（4）：35-37.

［3］Rozemeijer M J C，Kok J D，Boon J，et al.Towards an assessment of the effects of some artificial island-designs on the dutch coastal system（summary）［J］.Deutsche Hydrografische zeits chrift，1999，51（10）：41-42.

［4］魏龙，王义刚，黄慧明，等.冀东南堡人工岛工程潮流泥沙数值模拟［J］.水运工程，2012（6）：43-47. WEI L，WANG Y G，HUANG H M，et al.Flow and sediment numerical simulation on artificial island of south fort in eastern Hebei province［J］.Port&Waterway Engineering，2012（6）：43-47.

［5］王中起，韩志远，严冰.秦皇岛海域水文泥沙特征分析［J］.水道港口，2010，31（4）：247-252. WANG Z Q，HAN Z Y，YAN B.Hydrographic and sediment characteristics in Qinhuangdao sea area［J］.Journal of Waterway and Harbor，2010，31（4）：247-252.

［6］赵棣华，戚晨，庾维德，等．平面二维水流水质有限体积法及黎曼近似解模型［J］.水科学进展，2000，11（4）：368-374．ZHAO D H，QI C，YU W D，et al.Finite volume method and Riemann solver for depth-averaged two-dimensional flow-pollutants coupled medel［J］.Advances in Water Science，2000，11（4）：368-374．

［7］赵棣华，姚琪，蒋艳，等．通量向量分裂格式的二维水流-水质模拟［J］.水科学进展，2002，13（6）：701-706. ZHAO D H，YAO Q，JIANG Y，et al.FVS scheme in 2-D depth-averaged flow-pollutants modeling［J］.Advances in Water Science，2002，13（6）：701-706．

［8］王光谦，胡春宏.泥沙研究进展［M］.北京：中国水利水电出版社，2006.

［9］窦国仁.潮汐水流中的悬沙运动及冲淤计算［J］.水利学报，1963（4）：13-24. DOU G R.Suspended sediment movement&Erosion and Sediment Calculation in Tidal Flow［J］.Journal of Hydraulic Engineering，1963（4）：13-24.

［10］窦国仁，董风舞，窦希萍.潮流和波浪的挟沙能力［J］.科学通报，1995，40（5）：443-446. DOU G R，DONG F W，DOU X P，et al.The sediment carrying capacity in combined actions of waves and tidal current［J］.Chinese Science Bulletin，1995，40（5）：443-446.

［11］窦国仁，董风舞，窦希萍，等.河口海岸泥沙数学模型研究［J］.中国科学（A辑），1995，25（9）：995-1 001. DOU G R，DONG F W，DOU X P.Numerical simulation research on the sediment of estuarial and coastal areas［J］.Science in China：Series A，1995，25（9）：995-1 001.

［12］窦国仁，赵士清，黄亦芬.河道二维全沙数学模型的研究［J］.水利水运科学研究，1987（2）：1-12. DOU G R，ZHAO S Q，HUANG Y F.Study on two-dimensional total sediment transport mathematic model［J］.Journal of Nanjing Hydraulic Research Institute，1987（2）：1-12.

［13］天津水运工程科学研究院海岸与海洋资源利用研究中心.秦皇国际游轮游艇港海螺岛项目波浪潮流泥沙数模试验研究报告［R］.天津：交通运输部天津水运工程科学研究所，2012．

［14］黄建维.海岸与河口黏性泥沙运动规律的研究和应用［M］.北京：海洋出版社，2008.

［15］Mehta A J，Parchure T M，Ariathurai R.Resuspension potential of deposited cohesive sediment beds［J］.V.S.Kennedy（Ed.），Estuarine Comparisons，1982，60：591-609.

Analysis of artificial island construction influences on individual flood erosion and deposition in estuary

SHENG Tian⁃hang,SUN Dong⁃mei,ZHANG Yang
（State Key Laboratory of Hydraulic Engineering Simulation and Safety，Tianjin University，Tianjin 300072，China）

Given that the artificial island project is under construction in Qinhuangdao,large and small double nested model was adopted to establish plane two-dimensional hydrodynamic-sediment coupled mathematical model in Tanghe estuary artificial island project.Based on the feasible results of the model validation,the erosion and deposition characteristics of hydrodynamic force and sediment in this project were analyzed under the condition of 50 a frequency design flood.The artificial island construction influences on the flow field change rules,the individual flood erosion and deposition rule in the estuary,as well as the sediment accumulation distribution law near the project areas were also investigated.This study indicates that sediment accumulation is not increasing due to the artificial island construction,rather,it is decreasing under the influence of the distributary and sediment dredging.

artificial island;estuary;mathematical model;sediment accumulation;flood erosion and deposition

TV 142；O 242.1

A

1005-8443（2016）01-0018-09

2015-04-27；

2015-07-29

盛天航(1989-)，男，内蒙古人，硕士研究生，主要从事水文、城市雨洪管理、河道水动力及水质环境数值模拟方面的研究。

Biography：SHENG Tian-hang（1989-），male，master student.