【正文】 2 2 2 2 212 ( )s b q qq D q u D q v D q w qD P P K D Et x y D?? ? ???? ? ? ? ? ?? ? ? ? ? ? ? ???? ? ? ? ? ??? （ 2— 8） 2 2 2 2 21l s b q llq l D q l u D q l v D w q l w W q ll E D P P K D Ft x y E D?? ? ? ??? ??? ? ? ? ? ?? ? ? ? ? ? ? ??? ??? ? ? ? ? ????? （ 2— 9） 其中， u、 v 和 w 為 x、 y 和 z 三個方向上各自速度的分量； T 為位溫 ； S是鹽度； f 是科氏力參數(shù)； P 為壓力； ? 是密度； 邊界條件 海面的邊界條件是在 σ = 0 處： 淮海工學 院二〇一三屆 本科畢業(yè)設計（論文） 第 7 頁 共 19 頁 ? ?0,sx symu v DK ??? ? ????? ??????? ； EP? ????? （ 2— 10） ? ? ? ?, , , , 0 ,nPhD Q x y t S W x y tcK???? ??????? （ 2— 11） 2 2 2 / 3 21。 0 。 0 。在我們劃分的三角形單元內(nèi)，三角形的三個節(jié)點可以用整數(shù) Nj(j)來表示，其中， j 的值按規(guī)定是順 時針方向由 l 到 3。 3 模型的配置與 驗證 模型配置 計算區(qū)域與網(wǎng)格剖分 本次模型的開邊界范圍為 和 ，其中涵蓋了整個模擬的海岸線 及其相關島嶼。 。 close all。rb39。,39。39。 temp = fscanf(fid1,39。 淮海工學 院二〇一三屆 本科畢業(yè)設計（論文） 第 11 頁 共 19 頁 for i=1:113 fprintf(fid2,39。% % % % % % % %\n39。,temp(4,4*(i1)+2),temp(4,4*(i1)+1),0,temp(4,4*(i1)+3),0,temp(4,4*(i1)+4),0,0)。它的公式如下： ? ? 1c osNo i i iit??? ? ? ? ?????? （ 3— 1） 其中： ζ i ， ω i ， θ i 分別為第 i 分潮的振幅、頻率及相位； No 是分潮總數(shù)。 MATLAB 將 矩陣 算法 、 數(shù)值 研究 、科學 數(shù)值 可視化 分析和 非 線性 動態(tài)系統(tǒng)的 建模 等 許 多強大 的 功能集 合淮海工學 院二〇一三屆 本科畢業(yè)設計（論文） 第 12 頁 共 19 頁 在一個 方便 使用的視 圖窗口 環(huán)境中，為科學 可視化分析 、工程 建模和 必須 采用有效 的 數(shù)值 分析 的眾科學 領域 提供了一 套 全面的解決方 法 ，并 且 在 一定 程度上擺脫了傳統(tǒng) 的 如 C 語言 、 Fortran 語句 等的 非交互式 程序設計語言 的編輯 環(huán)境 ，因此 MATLAB 代表了當今 在 科學計算軟件 領域 的 綜合 先進水平。 。 (2)以連云港海域連島北驗潮站的潮位實際測量值來驗證 FVCOM 計算模擬潮位的正確性，如圖 2 所示： 210121 815 22 29 36 43 50 57 64 71 78 85 92 99106 113 120 127 134 141 148 155 162 169 176 183 190 197 204 211 218 225 232 239T M D 計算值 F V C O M 計算值 圖 2 潮位驗證曲線 由圖可知，連島北驗潮站潮位的誤差最大不超過 3 厘米，這表示實測值與計算值比較吻合。 黃海暖流沿著北黃海直接至渤海海峽，所以蘇北淺灘上一部分泥沙也可被輸運到渤海和北黃海。模擬計算的粒子示蹤軌跡見圖 4(右為 其 軌跡 放大 圖 )，模擬 6h 后，每隔 2h，輸出一次該示蹤粒子的空間位置。N 附近水下 m 處垂直方向上進行周期約為 M2的上下運動。在 26h 后，該粒子在 176。實驗結果 還需要大量的實測數(shù)據(jù)檢驗，在此可為以后研究工作的開展提供參考。 工作 展望 (1) 考慮溫鹽度、蒸發(fā)、降水徑流、海面的熱通量、風場等參數(shù)綜合對本模式的影響； (2) 在連云港海域近岸及主要的海灣和海峽處，劃分的三角形網(wǎng)格需要進一步的加密處理，從而能獲得更細致的潮流、潮汐分布特征； (3) 根據(jù)模擬結果對連云港海域潮汐余流、淺水分潮、泥沙運動軌跡等的特征進行更深入的探討與分析。 其次， 本論文是在導師呂海濱副教授的悉心指導下完成的。在此，謹向呂老師表示衷心的感謝和最崇高的敬意！ 本論文的成功完成，離不開各位關心我的老師、朋友和同學的關心和幫助。在未來的日子里，我會更加努力的學習和工作，不辜負父母對我的殷殷期望！ 淮海工學 院二〇一三屆 本科畢業(yè)設計（論文） 第 18 頁 共 19 頁 參 考 文 獻 [1] 周志德 .淺談泥沙研究的困難 [J].中國水利水電科學研究院學報， 2020（ 1） ： 33~36. [2] 竇國仁 .論泥沙起動流速 [J].水利學報， 1960(4)： 22~31. [3] 竇國仁 .再論泥沙起動流速 [J].泥沙研究， 1999(6)： 1~9. [4] 龔政 ,邰佳愛 ,張東生 .陸源污染物對海州灣環(huán)境影響研究 [J].河海大學學報 :自然科學版 ,2020， 30(4)： 5~8. [5] 趙洪波 ,王廣聚 .田灣核電站海域潮流場數(shù)值模擬 [J].水道港口 ,2020，28(1)： 5~9. [6] 謝飛 ,逄勇 ,宋志堯 .海州灣海域三維潮流數(shù)值模擬 [J].河海大學學報 :自然科學版 ,2020， 35(6)： 718~721. [7] 呂海濱 ,周立 ,湯均博 ,等 .海州灣西部水動力數(shù)值模擬 [J].淮海工學院學報 :自然科學版 ,2020， 17(1)： 72~5. [8] 朱玉榮 .全新世渤、黃、東海陸架泥沙輸運演變過程模擬研究 [J].黃渤海海洋： 海洋科學進展 ， 2020， 19（ 2）： 25~38. [9] 李孟國 ,李文丹 ,時鐘 ,等 .田灣核電站海域潮流泥沙數(shù)值模擬研究 [J].泥沙研究 ,2020(2)： 16~23. [10] 朱學明 .中國近海潮汐潮流的數(shù)值模擬與研究 [D]:中國海洋大學 ,2020. [11] 邱照宇 .基于 FVCOM 模擬河流輸入對膠州灣水質(zhì)的影響 [D]：中國海洋大學， 2020. [12] 朱建榮．海洋數(shù)值計算方法和數(shù)值模式 [M]．北京：海洋 I 出版社， 2020：131~132. [13] 邊昌偉 .中國近岸泥沙在渤海 _黃 海和東海的輸運 [D]：中國海洋大學，2020. [14] 張二駿 , 張東生 , 高飛 , 等 . 連云港水域流場數(shù)值模擬 [C].連云港建港指揮部 :連云港回淤研究論文集 . 南京 : 河海大學出版社 , 1990： 141~153. [15] 胡春宏 ， 王延貴 ， 張燕菁 .泥沙研究的進展與展望 [J].水文， 2020， 26（ 3） ：37~41. [16] 李孟國 .海岸河口泥沙數(shù)學模型研究進展 [J].海 洋工程， 2020， 24（ 1）：淮海工學 院二〇一三屆 本科畢業(yè)設計（論文） 第 19 頁 共 19 頁 139~154. [17] 吳倫宇 .基于 FVCOM 的浪、流、泥沙模型耦合及應用 [D]:中國海洋大學 ，2020. [18] 劉高峰 .長江口水沙 運動及三維泥沙模型研究 [D]:華東師范大學 ,2020. [19] 陳斌 .長江口附近海域三維懸浮泥沙的數(shù)值模擬研究 [D]:中國科學院研究生院（海洋研究所） ,2020. [20] 陳曉宏 ,陳永勤 ,賴國友 .珠江口懸浮泥沙遷移數(shù)值模擬 [J].海洋學報 ，2020,25(2)： 120~126. [21] 萬遠揚，金中武，黃仁勇 . 泥沙模型研究述評與 前景展望 [J].南水北調(diào)與水利科技， 2020， 4（ 1） ： 48~51. [22] Chen, C. H. Liu, R. C. Beardsley, 2020a. An unstructured, finitevolume, threedimensional, primitive equation ocean model: application to coastal ocean and estuaries. J. Atm. amp。在論文的進展過程中，我們互相幫助，相互探討，沒有大家的努力就沒有此論文的成功。不僅使我體會到嚴謹?shù)?