CN110490354B - 一种基于roms模拟结果计算潮汐锋位置的优化方法 - Google Patents

一种基于roms模拟结果计算潮汐锋位置的优化方法 Download PDF

Info

Publication number
CN110490354B
CN110490354B CN201910579125.7A CN201910579125A CN110490354B CN 110490354 B CN110490354 B CN 110490354B CN 201910579125 A CN201910579125 A CN 201910579125A CN 110490354 B CN110490354 B CN 110490354B
Authority
CN
China
Prior art keywords
calculating
tide
term
tidal
roms
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
CN201910579125.7A
Other languages
English (en)
Other versions
CN110490354A (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 CN201910579125.7A priority Critical patent/CN110490354B/zh
Publication of CN110490354A publication Critical patent/CN110490354A/zh
Application granted granted Critical
Publication of CN110490354B publication Critical patent/CN110490354B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00—Administration; Management
    • G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism

Landscapes

  • Business, Economics & Management (AREA)
  • Engineering & Computer Science (AREA)
  • Strategic Management (AREA)
  • Economics (AREA)
  • Human Resources & Organizations (AREA)
  • General Business, Economics & Management (AREA)
  • General Physics & Mathematics (AREA)
  • Marketing (AREA)
  • Tourism & Hospitality (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Primary Health Care (AREA)
  • General Health & Medical Sciences (AREA)
  • Development Economics (AREA)
  • Game Theory and Decision Science (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Operations Research (AREA)
  • Quality & Reliability (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Length Measuring Devices With Unspecified Measuring Means (AREA)

Abstract

本发明涉及一种基于ROMS模拟结果计算潮汐锋位置的优化方法，包括下列步骤：考虑表层风摩擦作用、底层潮流摩擦作用、净热通量作用和浮力平流作用的情况下，得到计算潮汐锋的公式；计算得到底层潮流对锋的作用项A；计算得到表层风场对锋的作用项B；计算得到净热通量对锋的作用项C；计算得到浮力平流对锋的作用项M，为使计算过程简便同时减少M项误差，对M项的计算进行反算处理；计算得到潮汐锋的位置Hfront；将潮汐锋所处的位置深度Hfront与ROMS模式输出的地形水深Hdepth相减，得到差值ΔH，ΔH为零的地方即为海洋中潮汐锋真正存在的位置。

Description

一种基于ROMS模拟结果计算潮汐锋位置的优化方法
技术领域
本发明属于海洋科学领域，具体涉及一种基于ROMS模拟结果计算潮汐锋位置的优化方法。
背景技术
潮汐锋是海洋锋的一种，指在增温季节由潮混合与层化共同作用形成的近岸混合区与远岸层化区之间的过渡海域，对动力海洋学和海洋生物学研究以及探讨低氧的形成及变化机制都具有重要意义。表层风摩擦作用、底层潮混合作用、净热通量、浮力平流等都会对潮汐锋的生成和位置变化产生影响。大多数学者根据Simpson-Hunter层化参数KSH＝log(H/Umax 3)对潮汐锋的位置进行过计算，但Simpson-Hunter层化参数仅考虑了地形和潮流对潮汐锋位置的影响，并未将影响潮汐锋位置的其他因素考虑在内，不能够较为准确的计算出潮汐锋的位置。海洋数值模式是定量描述海洋物理现象及其变化的数值模型，也是海洋与气候研究、预测的核心工具，是对原位观测资料的有力补充和拓展，长远来看能够降低海洋研究成本、节省外出作业人力资源。ROMS(Regional Ocean Modeling System,区域海洋模式系统)，是一个采用自由表面、静力假定、sigma坐标、曲线正交网格的三维原始方程模式。目前该海洋模式已被广泛应用于中国近海水动力模拟研究中。本发明考虑了风摩擦作用、底层潮流摩擦作用、净热通量作用和浮力平流作用等对潮汐锋的影响因素，得到一种基于ROMS模拟结果计算潮汐锋位置的优化方法。该方法能够更为准确的刻画出潮汐锋的位置及变化情况，对分析潮汐锋的形成和变化机制具有重要意义。
发明内容
有鉴于此，本发明提供一种基于ROMS模拟结果计算潮汐锋位置的优化方法。本发明的技术方案如下：
一种基于ROMS模拟结果计算潮汐锋位置的优化方法，包括下列步骤：
1)考虑表层风摩擦作用、底层潮流摩擦作用、净热通量作用和浮力平流作用的情况下，计算潮汐锋的公式为：
Figure GDA0004100247540000021
其中Hfront为潮汐锋所处的水深，A为底层潮流对锋的作用项，B为表层风场对锋的作用项，C代表净热通量对锋的作用项，M代表浮力平流对锋的作用项；
2)通过ROMS模式输出底层潮流流速同时选择合适的常数参量，计算得到底层潮流对锋的作用项A；
3)利用关于风场的ERA数据资料同时选择合适的常数参量，计算得到表层风场对锋的作用项B；
4)利用关于热通量的ERA-Interim数据资料同时选择合适的常数参量，计算得到净热通量对锋的作用项C；
5)根据计算得到的A项、B项、C项数值，结合ROMS模式输出的密度和流场数据并选择合适的常数参量，计算得到浮力平流对锋的作用项M，为使计算过程简便同时减少M项误差，对M项的计算进行反算处理；
6)将计算得到的A项、B项、C项和M项值代入式(1)，可以计算得到潮汐锋的位置Hfront；
7)将潮汐锋所处的位置深度Hfront与ROMS模式输出的地形水深Hdepth相减，得到差值ΔH，ΔH为零的地方即为海洋中潮汐锋真正存在的位置。
本发明实施例提供的技术方案的有益效果是：
本发明的一种基于ROMS模拟结果计算潮汐锋位置的生态参数优化方法，进一步提高了潮汐锋位置计算准确度，弥补了Simpson-Hunter参数不能准确反映潮汐锋短时间尺度内位置变化的不足。同时通过反算的处理方式，大大简化了浮力平流作用项M的计算过程，减少了M项的噪音，提高了计算效率。充分利用模式结果和再分析数据对潮汐锋位置进行更为准确的计算，体现了模式和再分析数据对海洋科学研究的重要性，为更加准确地进一步探究潮汐锋的形成和变化机制提供了帮助。这是一次模式结果和再分析数据与计算方法很好结合的案例。
附图说明
图1本发明的基于ROMS模拟结果计算潮汐锋位置的优化方法流程。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚，下面对本发明实施方式作进一步地详细描述。
本发明实施例介绍了一种基于ROMS模拟结果计算潮汐锋的优化方法，目的在于通过一种改进的、高效的计算方法，更加准确地计算出潮汐锋的位置，为进一步研究潮汐锋的形成和变化机制提供帮助。具体实施步骤如下：
1)考虑表层风摩擦作用、底层潮流摩擦作用、净热通量作用和浮力平流作用的情况下，本发明计算潮汐锋的公式为：
Figure GDA0004100247540000031
其中Hfront为潮汐锋所处的水深，A为底层潮流对锋的作用项，B为表层风场对锋的作用项，C代表净热通量对锋的作用项，M代表浮力平流对锋的作用项。
2)ROMS模式输出底层潮流流速，通过计算得到底层潮流对锋的作用项A，A的计算公式如下：
A＝2mobU*b 3    (2)
式(2)中mob为下层挟卷参量，取为1.12，而U*b 2＝CdbUb 2，Ub为底层潮流流速，海底阻力系数Cdb＝1.5×10-3。
3)登陆ECMWF(欧洲中期天气预报中心,https://www.ecmwf.int/)网站，下载关于风场的ERA-Interim再分析数据资料，通过计算得到表层风场对锋的作用项B，B的计算公式如下：
B＝2mowU*w 3    (3)
式(3)中mow上层挟卷参量，取为1.25，
Figure GDA0004100247540000032
Uw为表层风速，海水平均密度
Figure GDA0004100247540000033
空气密度ρa＝1.18kg m-3，海面阻力系数Cdw＝1.15×10-3。
4)登陆ECMWF(欧洲中期天气预报中心,https://www.ecmwf.int/)网站，下载关于热通量的ERA-Interim再分析数据资料，通过计算得到净热通量对锋的作用项C，C的计算公式如下：
Figure GDA0004100247540000041
式(4)中，Q为净热通量，g＝9.8m s-2，热膨胀系数α＝3×10-4k-1，系统比热Cp＝3.93×103J kg-1k-1。
5)根据上述计算公式得到的A项、B项、C项数值，结合ROMS模式输出的密度和流场数据，计算得到浮力平流对锋的作用项M，M的计算公式如下：
Figure GDA0004100247540000042
式(5)中，ρ1、ρ2分别为上层海水密度，h1、h2分别为上、下层海水的厚度(h1、h2,根据ROMS模式输出的垂向密度差得到)，u1、v1分别为上层海水东分量流速和北分量流速，u2、v2分别为下层海水东分量流速和北分量流速。为了减小噪音和误差，将式(5)做如下转换：
Figure GDA0004100247540000043
由于近岸水体混合均匀(h2＝0)，此时有
Figure GDA0004100247540000044
6)将式(2)、(3)、(4)、(5)代入式(1)可以计算得到改进算法后的潮汐锋的位置Hfront。
7)海洋中潮汐锋位置的确定，将根据上述公式计算得到的潮汐锋所处的位置深度Hfront与ROMS模式输出的地形水深Hdepth相减，得到差值ΔH，即ΔH＝Hfron-Hdepth，ΔH为零的地方即为海洋中潮汐锋真正存在的位置。
8)绘图分析，以渤海辽东湾为例，利用绘图软件(如Matlab-2014a)对ΔH绘制平面图，将ΔH＝0的地方绘制成红线，同时将ROMS模式输出的作为潮汐锋判据的0.2℃表底温差等值线绘制为黑线，对计算结果进行检验，两者符合较好，认为该计算方法对潮汐锋的位置计算较为准确，能够用来研究潮汐锋的形成和变化机制。

Claims (1)

1.一种基于ROMS模拟结果计算潮汐锋位置的优化方法，包括下列步骤：
1)考虑表层风摩擦作用、底层潮流摩擦作用、净热通量作用和浮力平流作用的情况下，计算潮汐锋的公式为：
Figure FDA0004100247520000011
其中Hfront为潮汐锋所处的水深，A为底层潮流对锋的作用项，B为表层风场对锋的作用项，C代表净热通量对锋的作用项，M代表浮力平流对锋的作用项；
2)通过ROMS模式输出底层潮流流速，计算得到底层潮流对锋的作用项A，A＝2mobU*b 3，mob为下层挟卷参量，取为1.12，而U*b 2＝CdbUb 2，Ub为底层潮流流速，海底阻力系数Cdb＝1.5×10-3；
3)利用关于风场的ERA数据资料，计算得到表层风场对锋的作用项B，B＝2mowU*w 3，mow为上层挟卷参量，取为1.25，
Figure FDA0004100247520000012
Uw为表层风速，海水平均密度
Figure FDA0004100247520000013
Figure FDA0004100247520000014
空气密度ρa＝1.18kg m-3，海面阻力系数Cdw＝1.15×10-3；
4)利用关于热通量的ERA-Interim数据资料，计算得到净热通量对锋的作用项C，
Figure FDA0004100247520000015
Q为净热通量，g＝9.8m s-2，热膨胀系数α＝3×10-4k-1，系统比热Cp＝3.93×103J kg-1k-1；
5)根据计算得到的A项、B项、C项数值，结合ROMS模式输出的密度和流场数据，计算得到浮力平流对锋的作用项M，M的计算公式如下：
Figure FDA0004100247520000016
式中，ρ1、ρ2分别为上层和下层海水密度，h1、h2分别为上、下层海水的厚度，h1、h2根据ROMS模式输出的垂向密度差得到，u1、v1分别为上层海水东分量流速和北分量流速，u2、v2分别为下层海水东分量流速和北分量流速；为减小噪音和误差，将上式做如下转换：
Figure FDA0004100247520000021
由于近岸水体混合均匀，h2＝0，此时有
Figure FDA0004100247520000022
6)将计算得到的A项、B项、C项和M项值代入式(1)，计算得到潮汐锋的位置Hfront；
7)将潮汐锋所处的位置深度Hfront与ROMS模式输出的地形水深Hdepth相减，得到差值ΔH，ΔH为零的地方即为海洋中潮汐锋存在的位置。
CN201910579125.7A 2019-06-28 2019-06-28 一种基于roms模拟结果计算潮汐锋位置的优化方法 Active CN110490354B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910579125.7A CN110490354B (zh) 2019-06-28 2019-06-28 一种基于roms模拟结果计算潮汐锋位置的优化方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910579125.7A CN110490354B (zh) 2019-06-28 2019-06-28 一种基于roms模拟结果计算潮汐锋位置的优化方法

Publications (2)

Publication Number Publication Date
CN110490354A CN110490354A (zh) 2019-11-22
CN110490354B true CN110490354B (zh) 2023-05-09

Family

ID=68546508

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910579125.7A Active CN110490354B (zh) 2019-06-28 2019-06-28 一种基于roms模拟结果计算潮汐锋位置的优化方法

Country Status (1)

Country Link
CN (1) CN110490354B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113642142B (zh) * 2021-06-08 2023-11-17 天津大学 一种基于海面热通量计算水体层化开始时刻的方法
CN113705046B (zh) * 2021-09-29 2023-08-15 自然资源部第一海洋研究所 一种提高海洋模式中海水温度计算准确性的方法

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103390248A (zh) * 2013-08-08 2013-11-13 牟林 一种海洋模型数值模拟的潮流能资源评估方法

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0820506B2 (ja) * 1986-09-10 1996-03-04 海洋科学技術センター 海洋音響トモグラフィーデータ伝送装置
JP2005273464A (ja) * 2004-03-23 2005-10-06 Hiroyuki Higa 海底潮流水車を用いた発電設備及び深層水汲み上げ設備
CN105894106B (zh) * 2015-01-05 2020-06-16 国家电网公司 一种海洋模式和气象模式的一体化耦合方法
CN105205200A (zh) * 2015-04-27 2015-12-30 河海大学 一种淤泥质海岸波生流数值模拟的验证方法
CN106446539A (zh) * 2016-09-19 2017-02-22 青岛地球软件技术有限公司 基于卫星定位系统的潮汐预报方法
CN113642142B (zh) * 2021-06-08 2023-11-17 天津大学 一种基于海面热通量计算水体层化开始时刻的方法

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103390248A (zh) * 2013-08-08 2013-11-13 牟林 一种海洋模型数值模拟的潮流能资源评估方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
常海洁 ; 王继光 ; .不同潮汐条件下的黄渤海海温模拟及特征分析.舰船电子工程.2017,(10),全文. *

Also Published As

Publication number Publication date
CN110490354A (zh) 2019-11-22

Similar Documents

Publication Publication Date Title
Zheng et al. A modeling study of the Satilla River estuary, Georgia. I: Flooding-drying process and water exchange over the salt marsh-estuary-shelf complex
CN105894106B (zh) 一种海洋模式和气象模式的一体化耦合方法
CN110490354B (zh) 一种基于roms模拟结果计算潮汐锋位置的优化方法
Marinov et al. Application of COHERENS model for hydrodynamic investigation of Sacca di Goro coastal lagoon (Italian Adriatic Sea shore)
CN110096792B (zh) 一种计算非定常Langmuir环流的动态模拟方法
Ribbat et al. Transport variability over the Hawkesbury Shelf (31.5–34.5 S) driven by the East Australian Current
Milne An experimental investigation of turbulence and unsteady loading on tidal turbines
Koropitan et al. Three-dimensional modeling of tidal circulation and mixing over the Java Sea
Lü et al. Tidally induced upwelling off Yangtze River estuary and in Zhejiang coastal waters in summer
Liu et al. A utilizing efficiency estimation method for wave-driven unmanned surface vehicle
Nayak et al. Tidal and Residual Circulation in the Gulf of Khambhat and its Surrounding on the West Coast of India
Simizu et al. Barotropic response of the Sea of Okhotsk to wind forcing
Song et al. Effects of wave–current interactions on bay–shelf exchange
CN117668141B (zh) 地球系统耦合模型海气系统高频变化优化垂向参数化方法
Zhang et al. A new hybrid vertical coordinate ocean model and its application in the simulation of the Changjiang diluted water
Wang et al. The impact of surface wave on the sediment erosion and deposition near the wellow river mouth, China.
Wang et al. Long-term mean circulation in the Japan Sea as reproduced by multiple eddy-resolving ocean circulation models
Huang et al. Suspended sediment dynamics and influencing factors during typhoons in Hangzhou Bay, China
Jones et al. Storm surge computations in estuarine and near-coastal regions: the Mersey estuary and Irish Sea area
CN113742901B (zh) 一种分析潮汐对无机氮跨陆架运输影响机理的方法
Hermawan et al. High spatial grid resolution of hydrodynamic numerical modeling for sea current energy site selection in Indonesia
Cao et al. Seamount-induced mixing revealed through idealized experiments and its parameterization in an Oceanic General Circulation Model
CN117577167B (zh) 一种太湖蓝藻漂移数值模拟改进方法
Gelati Modeling the impact on sediment texture of large-scale tidal power in the Bay of Fundy
Moeini et al. Extreme Wave Prediction Based on the Numerically Simulated Data in a Semi-Enclosed Basin

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