金夢(mèng)海灣及其鄰近海域水環(huán)境因子時(shí)空分布特征

匡翠萍，王 丹，趙 釩，劉會(huì)欣，朱 磊

（1.同濟(jì)大學(xué)土木工程學(xué)院，上海200092；2.河北省地礦局第八地質(zhì)大隊(duì)，河北秦皇島066001）

秦皇島是我國渤海海濱著名的旅游度假勝地，以美麗的海灘、怡人的夏季氣候聞名于國內(nèi)外［1］。金夢(mèng)海灣是秦皇島市重要的濱海浴場，為了開發(fā)旅游資源，自2011年起先后修建了蓮花島和海螺島2座人工島。隨著秦皇島城市化進(jìn)程的加快，陸域污染、沿海排污、外海養(yǎng)殖等對(duì)近岸海域生態(tài)環(huán)境產(chǎn)生較大負(fù)面影響，同時(shí)受近岸人工島的影響，周邊海域水動(dòng)力環(huán)境發(fā)生顯著變化［2-4］，致使綠潮等生態(tài)災(zāi)害不斷爆發(fā)，對(duì)當(dāng)?shù)睾Ｑ笊鷳B(tài)環(huán)境和社會(huì)經(jīng)濟(jì)等造成嚴(yán)重影響。

研究近岸海域水質(zhì)時(shí)空分布特征可以為有效削減污染負(fù)荷、改善水環(huán)境等提供重要的理論參考和數(shù)據(jù)支持［5］。不少學(xué)者對(duì)海灣水質(zhì)時(shí)空分布特征展開了諸多研究。陳曉玲等［6］根據(jù)實(shí)測資料，研究了香港海域葉綠素a濃度及其相關(guān)水質(zhì)因子的時(shí)空分布規(guī)律，結(jié)果表明香港海域東部的半封閉海灣，赤潮發(fā)生頻率較高，西部海域水體交換能力強(qiáng)，赤潮發(fā)生頻率較低。周暢浩等［7］研究了大亞灣海域顆粒態(tài)氮磷的組成分布及其關(guān)鍵控制因素，結(jié)果表明浮游植物生長和徑流輸入是顆粒態(tài)氮磷分布的主要影響因素；梁志宏等［8］基于水質(zhì)實(shí)測數(shù)據(jù)，系統(tǒng)分析了深圳灣水質(zhì)時(shí)空分布特征及其污染物主要來源；Zhou等［9］利用主成分分析和GIS等方法對(duì)香港東部海域水環(huán)境時(shí)空變化進(jìn)行研究，發(fā)現(xiàn)徑流污染等對(duì)水環(huán)境影響較大。Yao等［10］探討了膠州灣浮游植物色素及營養(yǎng)鹽的時(shí)空分布，結(jié)果表明春季和夏季膠州灣北部和東部營養(yǎng)鹽及浮游色素濃度較高，與人類活動(dòng)有密切關(guān)系。

數(shù)學(xué)模型能夠?qū)ξ廴疚锏臅r(shí)空分布及輸運(yùn)規(guī)律進(jìn)行數(shù)字化、可視化處理，在模擬預(yù)測水動(dòng)力、水環(huán)境等方面具有較大的優(yōu)勢(shì)［11-12］。自1925年Streeter和Phelps建立了BOD-DO耦合模型以來，水質(zhì)模型發(fā)展已有近百年歷史，陸莎莎等［13］將水質(zhì)模型的發(fā)展歷程分為了3個(gè)階段，第1階段為簡單的氧平衡模型，第2階段重點(diǎn)研究各水質(zhì)組分之間的相互作用，第3階段全面考慮了水沙邊界的影響。目前國內(nèi)水環(huán)境模型研究主要關(guān)注水動(dòng)力-水質(zhì)-水生態(tài)耦合模型的應(yīng)用［14］。鑒于此，本文基于MIKE 21軟件建立金夢(mèng)海灣及其鄰近海域水動(dòng)力與水質(zhì)耦合模型，模擬分析金夢(mèng)海灣主要污染物時(shí)空分布特性，在此基礎(chǔ)上采用主成分綜合得分對(duì)金夢(mèng)海灣及其鄰近海域的水質(zhì)污染狀況進(jìn)行評(píng)價(jià)，以期為該海域環(huán)境管理及生態(tài)災(zāi)害防治提供科學(xué)支撐。

1 研究區(qū)域概述

研究區(qū)域涵蓋金夢(mèng)海灣及其周邊海域（見圖1）。金夢(mèng)海灣是秦皇島市重要的濱海浴場，也是連接北戴河區(qū)和海港區(qū)的樞紐過度地帶。研究區(qū)域入海河流自北向南依次為新開河、湯河及新河，3條河流徑流量年際變幅較大，并且年內(nèi)分配很不均勻。新開河流經(jīng)秦皇島市區(qū)東部，是市區(qū)排水的主要河道，于東山東側(cè)注入渤海。河道全長11 km，其中主河道長4.05 km，流域面積43.9 km2，多年平均徑流量為740×104m3。湯河上游有兩大支流，東支發(fā)源于撫寧縣柳觀峪村西北，西支發(fā)源于撫寧縣溫泉堡西南的方家河村，河道全長28.5 km，流域面積184 km2，多年平均徑流量為3.68×107m3。新河發(fā)源于撫寧縣棲云寺山東麓，橫貫于北戴河森林濕地，河道全長15 km，河道常年自然淤積，常發(fā)生斷流現(xiàn)象［15-16］。

圖1 研究區(qū)域及模型網(wǎng)格Fig.1 Area studied and mesh of model

近年來，金夢(mèng)海灣的人工構(gòu)筑物建設(shè)密度較大，主要的人工建筑物有蓮花島、海螺島及3座潛堤。其中蓮花島于2011年審批建設(shè)，島嶼呈橢圓形布置，長軸約1 650 m，短軸約1 100 m，離岸約1 km；蓮花島北側(cè)近岸海域3座潛堤于2012年建設(shè)完工，離岸距離約為380 m，長約360 m；海螺島2013年底開始建設(shè)，離岸約600 m。

秦皇島海域經(jīng)過近幾年的環(huán)境綜合治理，水環(huán)境得到一定改善，但仍然存在一些突出問題［17］，生態(tài)災(zāi)害仍未得到有效控制，綠潮、赤潮、水母等海洋生態(tài)災(zāi)害頻發(fā)，研究金夢(mèng)海灣水環(huán)境時(shí)空分布特征對(duì)于海洋環(huán)境治理具有重要意義。本文選取的水質(zhì)指標(biāo)有化學(xué)需氧量COD、溶解氧DO、氨氮NH3-N、硝酸鹽氮NO-3-N、亞硝酸鹽氮NO-2-N和磷酸鹽PO3-4。

2 研究方法

2.1 數(shù)學(xué)模型建立

采用丹麥水力學(xué)研究所（Danish Hydraulic Institute，DHI）研發(fā)的MIKE 21 Flow Model軟件中的水動(dòng)力模塊（HD）和水質(zhì)模塊（ECO Lab）建立金夢(mèng)海灣海域潮流和水質(zhì)數(shù)學(xué)模型［18］。模型在空間上通過有限體積法進(jìn)行離散，顯示歐拉法求解，具備計(jì)算速度快、精度高、易收斂等特點(diǎn)［19-20］。

2.1.1 基本控制方程

潮流模型的控制方程是在笛卡爾坐標(biāo)系下，通過對(duì)三維連續(xù)方程和動(dòng)量方程沿深度進(jìn)行積分，得到如下二維淺水方程：

式中：t為時(shí)間；x、y為笛卡爾坐標(biāo)；h為總水深；d為靜止水深；η為水面高程；u、v分別為x、y方向的水深平均流速；S為源匯項(xiàng)；uS、vS分別為源/匯在x、y方向上的流入/流出速度；f為科氏力參數(shù)；g為重力加速度；Pa為大氣壓力；ρ為水的密度；ρ0為水的參考密度；τsx、τsy分別為x、y方向水面風(fēng)剪應(yīng)力分量；τbx、τby分別為x、y方向河床切應(yīng)力分量；Fu、Fv分別為x、y方向水平渦黏應(yīng)力量。

ECO Lab模塊基本方程如下：

式中：c為狀態(tài)變量質(zhì)量濃度；Fc為水平擴(kuò)散項(xiàng)；Cs為源質(zhì)量濃度；Pc為狀態(tài)變量的生化反應(yīng)過程。

2.1.2 計(jì)算網(wǎng)格及模型參數(shù)

模型平面坐標(biāo)系選用北京54坐標(biāo)系，采用非結(jié)構(gòu)三角網(wǎng)格進(jìn)行計(jì)算（見圖1）。模型共有13 431個(gè)節(jié)點(diǎn)和25 585個(gè)單元，權(quán)衡計(jì)算精度及計(jì)算時(shí)間，對(duì)河口海岸區(qū)進(jìn)行網(wǎng)格加密，網(wǎng)格空間步長最小為10 m，對(duì)外海區(qū)域采用稀疏網(wǎng)格，網(wǎng)格空間步長最大為2 500 m。

潮流模型中渦黏系數(shù)采用Smagorinsky方程來計(jì)算，其中常數(shù)取0.28。模型采用動(dòng)邊界處理技術(shù)，干水深、淹沒水深和濕水深分別取0.005 m、0.05 m和0.1 m。時(shí)間步長由模型自動(dòng)調(diào)節(jié)，為0.01～30 s。曼寧數(shù)根據(jù)底部泥沙粒徑和水深分布取值74 m1/3·s-1。水質(zhì)模型率定參數(shù)較多，由于篇幅所限，這里僅羅列部分參數(shù)取值，COD一階降解速率為0.02d-1，硝化反應(yīng)速率0.01d-1。

2.1.3 邊界條件

數(shù)學(xué)模型包含3條外海開邊界、3個(gè)入海河流開邊界以及陸域岸線閉邊界。水動(dòng)力模型的3條外海開邊界由Flather條件控制，其潮位和流速過程由渤海潮流模型（邊界為大連至煙臺(tái)）提供，入海河流上游開邊界由實(shí)測月平均流量進(jìn)行控制。水質(zhì)模型的外海開邊界為本底污染物質(zhì)量濃度，河流開邊界均為實(shí)測污染物質(zhì)量濃度值，文中COD質(zhì)量濃度為CODMn。初始水位和流速均為零，初始污染物質(zhì)量濃度設(shè)定為2018年3月實(shí)測污染物質(zhì)量濃度均值。陸域閉邊界采用不可滑移條件，即流速為零。

2.2 模型驗(yàn)證

2.2.1 水動(dòng)力模型驗(yàn)證

對(duì)金夢(mèng)海灣潮流模型分別進(jìn)行潮位和潮流驗(yàn)證，驗(yàn)證點(diǎn)位置如圖2所示。潮位驗(yàn)證（QHD點(diǎn)）采用2017年9月5—7日實(shí)測潮位過程，潮流驗(yàn)證（S1-S3點(diǎn)）采用2013年5月11日8：00至5月12日8：00實(shí)測潮流過程，潮流驗(yàn)證（S4點(diǎn)）采用2017年9月6日0：00至9月7日0：00實(shí)測潮流過程。驗(yàn)證結(jié)果如圖3所示，從驗(yàn)證結(jié)果可以看出，計(jì)算的潮位、流速和流向與實(shí)測值吻合較好，部分計(jì)算值與實(shí)測值略有偏差，這主要是由于秦皇島近岸海域地形復(fù)雜，且位于無潮點(diǎn)附近，潮波變化復(fù)雜，但總體上模型能較好地模擬出金夢(mèng)海灣近岸海域水動(dòng)力變化特征。

圖2 模型驗(yàn)證點(diǎn)位及T1分析點(diǎn)位置Fig.2 Location of verified stations and T1

圖3 潮位和流速、流向驗(yàn)證Fig.3 Verification of tidal level,current velocity,and current direction

2.2.2 水質(zhì)模型驗(yàn)證

水質(zhì)模型驗(yàn)證資料采用2018年3月13日在金夢(mèng)海灣（Q1-Q5點(diǎn)）5個(gè)站點(diǎn)的實(shí)測DO質(zhì)量濃度和2018年8月21日在金夢(mèng)海灣（Q1-Q3點(diǎn)）3個(gè)站點(diǎn)的實(shí)測水質(zhì)因子質(zhì) 量 濃 度（COD、NH3-N、NO-3-N、PO3-4）。驗(yàn)證點(diǎn)位置及水質(zhì)驗(yàn)證結(jié)果分別見圖2和圖4，其中站點(diǎn)Q2和Q3磷酸鹽PO3-4質(zhì)量濃度為一類海水水質(zhì)標(biāo)準(zhǔn)，實(shí)測值均小于0.015 mg·L-1。驗(yàn)證結(jié)果表明，計(jì)算值與實(shí)測值相對(duì)誤差較小，造成誤差的主要原因可能是站點(diǎn)監(jiān)測時(shí)間不同步，未考慮其他海上臨時(shí)性污染源等?？傮w上，模型能較好地模擬金夢(mèng)海灣及其鄰近海域的水質(zhì)變化，可以用于金 夢(mèng)海灣及其鄰近海域水動(dòng)力水質(zhì)模擬。

圖4 近岸海域污染物質(zhì)量濃度驗(yàn)證Fig.4 Verification of pollutant concentration in coastal waters

2.3 主成分分析

主成分分析方法是將多個(gè)具有相關(guān)性的變量轉(zhuǎn)化成少數(shù)幾個(gè)不相關(guān)的綜合指標(biāo)的分析與統(tǒng)計(jì)方法，其中主成分保留了原始變量的大多數(shù)信息，這樣能夠簡化復(fù)雜的研究，提高分析效率。主成分分析的一般步驟為：首先將原始數(shù)據(jù)標(biāo)準(zhǔn)化，建立變量的相關(guān)系數(shù)矩陣，再計(jì)算相關(guān)系數(shù)矩陣的特征值和特征向量以及各主成分的貢獻(xiàn)率和累積貢獻(xiàn)率，在此基礎(chǔ)上確定主成分個(gè)數(shù)，最后求出主成分綜合得分。本文以主成分綜合得分研究金夢(mèng)海灣及其鄰近海域水質(zhì)指標(biāo)的時(shí)空變化格局，并分析其影響因素［21-25］。

3 結(jié)果與討論

3.1 潮流特征

金夢(mèng)海灣及其鄰近海域典型時(shí)刻（漲急、漲憩、落急和落憩）潮流場如圖5所示。漲急、落急時(shí)刻，近岸海域潮流流速略小于外海流速，近岸流速約為0.22 m·s-1，外海流速約為0.40m·s-1，其中湯河口人工島群周邊海域流速較小，約為0.1～0.2 m·s-1。由于岬角效應(yīng)，金山嘴附近流速大且流向偏轉(zhuǎn)較大。漲憩和落憩時(shí)刻，潮流流速較小，從外海向河口海域遞減。總體而言，研究海域潮流流速較小，水動(dòng)力較弱，潮流呈現(xiàn)明顯的沿岸往復(fù)流運(yùn)動(dòng)，漲潮流方向?yàn)閺腘E（東北）向SW（西南），落潮流方向?yàn)閺腟W向NE，漲潮流略強(qiáng)于落潮流。

圖5 金夢(mèng)海灣及其鄰近海域典型時(shí)刻潮流場Fig.5 Typical current fields in JMB and its adjacent waters

3.2 污染物時(shí)空分布特征

基于驗(yàn)證良好的MIKE 21水動(dòng)力水質(zhì)耦合模型，對(duì)2018年3月至8月金夢(mèng)海灣及其鄰近海域污染物變化情況進(jìn)行模擬分析。

3.2.1 污染物季節(jié)性變化特征

根據(jù)模型的計(jì)算結(jié)果，在近岸海域選取分析點(diǎn)T1（位置見圖2），不同月份污染物濃度變化如圖6。

圖6 近岸海域不同污染物質(zhì)量濃度月際變化Fig.6 Monthly variation of different pollutant concentrations in coastal water

COD是水體污染的主要參數(shù)之一，反映了水體受還原性物質(zhì)污染的程度［26］。由圖6a可見，4月份COD質(zhì)量濃度較高，這可能也與浮游植物的生長有關(guān)，相關(guān)研究表明，渤海的浮游植物量在春秋兩季相對(duì)較高，浮游植物的死亡及排泄作用加強(qiáng)，導(dǎo)致COD質(zhì)量濃度升高［27-28］。COD質(zhì)量濃度季節(jié)變化高峰出現(xiàn)在7月，這與豐水期河流攜帶泥沙同時(shí)伴隨大量有機(jī)物進(jìn)入近岸海域有關(guān)。DO反映了生物生長和環(huán)境質(zhì)量狀況，是海水中各種物理、化學(xué)、生物過程的重要參數(shù)［29-30］。DO質(zhì)量濃度呈下降趨勢(shì)，這主要和溫度及有機(jī)物耗氧分解等有關(guān)。3月份水溫低，氧氣溶解度較大，同時(shí)有機(jī)物耗氧速率較低，因此3月份DO質(zhì)量濃度最高。海水中的溶解性無機(jī)氮（DIN）是浮游植物的主要營養(yǎng)鹽之一，由三部分組成，影響氮營養(yǎng)鹽的因素較多，其季節(jié)變化機(jī)制較復(fù)雜［31］。從圖6 c-6 e可以看出，金夢(mèng)海灣各種形態(tài)氮濃度存在明顯的季節(jié)變化，3～6月NH3-N為DIN的主要形態(tài)，7～8月NO-3-N含量較高，由于極不穩(wěn)定，水體中含量較低。磷酸鹽也是海域環(huán)境主要污染物之一，金夢(mèng)海灣PO3-4含量呈上升趨勢(shì)，其影響因素也較為復(fù)雜，陸源排放、底泥釋放等都有可能導(dǎo)致PO3-4增加［32-33］。

3.2.2 污染物空間分布特征

圖7～12為2018年8月12和13日漲落急時(shí)刻污染物濃度分布圖，從空間分布上可以看出，近岸海域污染物濃度場隨潮流呈周期性變化，河流污染物排入近岸海域后，高濃度區(qū)主要集中在河口近岸海域，隨漲潮流向SW方向運(yùn)移，隨落潮流向NE方向運(yùn)移。COD質(zhì)量濃度分布由近岸向外海遞減，并且近岸海域質(zhì)量濃度梯度較外海大。DO質(zhì)量濃度由近岸向外海呈現(xiàn)出先減后增的趨勢(shì)，其影響因素較為復(fù)雜，其中河口低氧區(qū)的形成主要與人類活動(dòng)有關(guān)。氮磷營養(yǎng)鹽的空間分布呈現(xiàn)出湯河口人工島群周邊海域營養(yǎng)鹽濃度明顯高于其他河口區(qū)域，這可能與湯河口海域接受的陸源污染負(fù)荷更大有關(guān)，另一方面工程改變了周邊海域水動(dòng)力環(huán)境，水體交換能力降低，導(dǎo)致對(duì)污染物的混合及稀釋能力降低。

圖7 金夢(mèng)海灣及其鄰近海域漲落急時(shí)刻COD質(zhì)量濃度分布Fig.7 Spatial distribution of COD concentration at maximum flood and ebb

圖8 金夢(mèng)海灣及其鄰近海域漲落急時(shí)刻DO質(zhì)量濃度分布Fig.8 Spatial distribution of DO concentration at maximum flood and ebb

圖10 金夢(mèng)海灣及其鄰近海域漲落急時(shí)刻-N質(zhì)量濃度分布Fig.10 Spatial distribution of -N concentration at maximum flood and ebb

3.2.3 主成分綜合得分分析

基于數(shù)學(xué)模型計(jì)算結(jié)果，選取35個(gè)站點(diǎn)（見圖13）3月、5月和8月各水質(zhì)指標(biāo)的月均值進(jìn)行主成分綜合得分分析，其中3月代表枯水期、5月代表平水期、8月代表豐水期。采用SPSS軟件進(jìn)行分析，結(jié)果顯示，KMO（Kaiser-Meyer-Olkin）檢驗(yàn)值為0.703（大于0.50），Bartlett球形檢驗(yàn)值為0（小于0.05），表明變量間不相互獨(dú)立，符合主成分分析要求。經(jīng)過計(jì)算，根據(jù)特征根大于1的原則提取了2個(gè)主成分，累計(jì)方差貢獻(xiàn)率為89.443%，說明這2個(gè)主成分能夠反應(yīng)全部數(shù)據(jù)的大部分信息。

圖13 分析點(diǎn)位置圖Fig.13 Location of analysis points

圖11 金夢(mèng)海灣及其鄰近海域漲落急時(shí)刻-N質(zhì)量濃度分布Fig.11 Spatial distribution of -N concentration at the maximum flood and ebb

圖12 金夢(mèng)海灣及其鄰近海域漲落急時(shí)刻PO3-4質(zhì)量濃度分布Fig.12 Spatial distribution of PO3-4 concentration at maximum flood and ebb

為反映金夢(mèng)海灣及其鄰近海域的水質(zhì)狀況，通過計(jì)算不同季節(jié)主成分綜合得分并構(gòu)建時(shí)空分布圖（見圖14）來評(píng)估污染水平，主成分綜合得分越大，水質(zhì)污染狀況越嚴(yán)重。經(jīng)過計(jì)算，3月、5月和8月主成分綜合得分中位數(shù)分別為-1.4、-0.2和1.5，金夢(mèng)海灣及其鄰近海域8月水質(zhì)污染最嚴(yán)重，3月污染較輕，這與陸源排污和降水有直接關(guān)系，秦皇島每年7、8月進(jìn)入主汛期，大量農(nóng)田化肥、生活垃圾等隨雨水匯入河流，入海河流污染物濃度升高，造成近岸海域水體污染嚴(yán)重。另一方面，也與海域內(nèi)浮游植物的種類、密度和分布等有一定的關(guān)系，豐水期比枯水期水溫高，近岸河口區(qū)光照充足，浮游植物生長代謝旺盛，消耗大量氧氣，同時(shí)浮游植物殘骸及代謝產(chǎn)物等被分解為營養(yǎng)鹽加劇海域水質(zhì)污染。

不同時(shí)期主成分綜合得分空間分布圖可以直觀展示金夢(mèng)海灣及其鄰近海域水質(zhì)狀況（見圖14），受海水溫度、陸源排放及海洋潮汐等影響，近岸海域水質(zhì)綜合得分有明顯的空間變化特征。從區(qū)域上看，平水期新河近岸海域水質(zhì)污染狀況較嚴(yán)重，枯水期和豐水期湯河口海域水質(zhì)污染狀況最嚴(yán)重。從總體上看，河口水域水質(zhì)污染較嚴(yán)重，主成分綜合得分呈現(xiàn)出由近岸向遠(yuǎn)海逐漸減小的趨勢(shì)，進(jìn)一步表明金夢(mèng)海灣及其鄰近海域污染物主要來自陸源徑流的輸入，因此改善近岸海域水環(huán)境必須從源頭控制，減少河流向海洋的污染輸入，治海先治河、治河先治水、治水先治污。另一方面，湯河口人工島群的建設(shè)對(duì)周圍海域的生態(tài)環(huán)境產(chǎn)生一定的影響，原有的水動(dòng)力條件發(fā)生改變，因此，海岸工程在規(guī)劃設(shè)計(jì)時(shí)要充分考慮多個(gè)工程開發(fā)的綜合影響以及對(duì)水動(dòng)力環(huán)境影響的疊加效應(yīng)。

圖14 不同季節(jié)主成分綜合得分空間分布Fig.14 Spatial distributions of comprehensive scores of principal components in different seasons

4 結(jié)論

（1）金夢(mèng)海灣近岸海域潮流流速較小，水動(dòng)力較弱，潮流總體特征為沿岸往復(fù)流，漲潮西南向，落潮東北向。污染物濃度隨潮流呈周期性變化，漲潮向西南向輸運(yùn)，落潮向東北向輸運(yùn)。

（2）應(yīng)用主成分綜合得分評(píng)價(jià)了金夢(mèng)海灣及其鄰近海域的水質(zhì)污染狀況，在時(shí)間尺度上，豐水期水質(zhì)劣于平水期和枯水期；在空間尺度上，水質(zhì)污染程度整體表現(xiàn)為從外海向河口逐漸加重的趨勢(shì)，空間分布特征反映了金夢(mèng)海灣水質(zhì)污染的主要驅(qū)動(dòng)因素來自陸源污染，近岸海域工程、潮流動(dòng)力等控制了污染物在海域的遷移擴(kuò)散。

（3）湯河口近岸海域潮流動(dòng)力較弱且工程較密集，建議進(jìn)一步加強(qiáng)湯河等河流的水質(zhì)監(jiān)測和污染防治工作，加大河流污染治理力度，防止金夢(mèng)海灣及其鄰近海域水質(zhì)進(jìn)一步惡化，減少生態(tài)災(zāi)害。

作者貢獻(xiàn)說明：

匡翠萍：項(xiàng)目負(fù)責(zé)人，論文撰寫及修改。

王丹：數(shù)據(jù)分析，論文撰寫及修改。

趙釩：數(shù)值模擬。

劉會(huì)欣：現(xiàn)場實(shí)測。

朱磊：現(xiàn)場實(shí)測。