CN116151085B - 风暴潮与天文潮共同作用下的电厂取水卷吸率计算方法 - Google Patents
风暴潮与天文潮共同作用下的电厂取水卷吸率计算方法 Download PDFInfo
- Publication number
- CN116151085B CN116151085B CN202310148210.4A CN202310148210A CN116151085B CN 116151085 B CN116151085 B CN 116151085B CN 202310148210 A CN202310148210 A CN 202310148210A CN 116151085 B CN116151085 B CN 116151085B
- Authority
- CN
- China
- Prior art keywords
- sea area
- model
- power plant
- typhoon
- water
- 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
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/25—Design optimisation, verification or simulation using particle-based methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/04—Power grid distribution networks
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/06—Wind turbines or wind farms
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Data Mining & Analysis (AREA)
- Mathematical Analysis (AREA)
- General Engineering & Computer Science (AREA)
- Computational Mathematics (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Geometry (AREA)
- Operations Research (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明涉及风暴潮与天文潮共同作用下的电厂取水卷吸率计算方法,包括,建立大尺度海域水动力数学模型;台风的气压场及风场计算;大尺度海域台风风暴潮模型验证;建立嵌套的中尺度工程海区水动力数学模型;嵌套模型的潮流验证;建立中尺度工程海区粒子追踪数学模型;电厂取水卷吸率的计算。本发明针对受台风影响的滨海核电厂为研究对象,建立适用于风暴潮模拟的大尺度海域模型与用于工程近区流场模拟的中尺度海区模型嵌套的二维数学模型,采用示踪粒子方法研究台风引起的风暴潮与天文潮共同作用下漂浮物的运移轨迹,分析漂浮物在电厂取水明渠内的分布情况,以此基础计算电厂取水卷吸率。为核电厂取水安全预警及应急措施的规划和设计提供科学依据。
Description
技术领域
本发明涉及一种风暴潮与天文潮共同作用下的电厂取水卷吸率计算方法,是一种环境水力学计算方法,是一种适用于滨海电厂在台风引起的风暴潮与天文潮共同作用下电厂取水卷吸率的计算方法。
背景技术
在核电厂运行过程中,需要通过取水系统从河流、海洋等天然水域取水以冷却机组。如果取水系统发生堵塞,将可能导致机组被迫降功率或停机,严重时甚至引发核安全事故。取水明渠拦导工程措施在核电取水工程中陆续得到应用,它是核电厂冷源安全保障的重要部分。取水明渠内的流速、流向分布等水力特性对进入明渠的漂浮物和海生物的卷吸效应影响显著,是拦导工程设计中需要考虑的重要因素。对于沿海台风频发,不利的海洋水文、气象条件下漂浮物(海生物)对冷源安全造成的突发影响如何予以评估、如何计算台风引起的风暴潮与天文潮共同作用下滨海电厂取水卷吸率的计算,是一个需要解决的问题。
发明内容
为了克服现有技术的问题,本发明提出了一种风暴潮与天文潮共同作用下的电厂取水卷吸率计算方法。所述的方法针对滨海核电厂,利用大尺度海域模型与中尺度海区模型嵌套的二维数学模型结合台风模型及示踪粒子的方法,提出了一种台风引起的风暴潮与天文潮共同作用下电厂取水卷吸率的计算方法。为核电厂取水安全预警以及应急措施的规划和设计提供科学依据。
本发明是一种风暴潮与天文潮共同作用下的电厂取水卷吸率计算方法,所述方法的步骤如下:
步骤1,建立大尺度海域水动力数学模型:根据厂址所在海域的位置及台风路径确定模型范围,建立适用于台风风暴潮模拟计算的大范围海域水动力数学模型,平面直角坐标系下模型的控制方程如下:
连续方程:
运动方程:
x方向:
y方向:
上式中:t为时间,x,y为笛卡尔坐标;η为水位,h=d+η为总水深,d为静水深;为x方向的水深平均速度,/>为y方向的水深平均速度:/>u为x方向速度,v为y方向速度,z为笛卡尔坐标;us为源在x方向的排水速度;vs为源在y方向的排水速度;f为柯氏力参数,f=2Ωsinφ,其中:Ω为旋转角速度,φ为地理纬度;g为重力加速度;ρ为水的密度,ρ0为参考水密度;pa为大气压力;A为涡粘性系数;τb=(τbx,τby)为底部应力,cf为阻力系数,/>床底上部流速;τs=(τsx,τsy)为水面风应力,ρa为空气密度,cd为空气阻力系数,/>水面上空10m处风速;S为源汇项;
大尺度海域数学模型的开边界水位条件由全球潮汐预报模型结果提取得到,台风的气压场及风场作为输入条件;大尺度海域数学模型的计算结果将为后续嵌套的中尺度海区模型提供计算所需的边界条件;
步骤2,台风的气压场及风场计算:分析电厂所在工程海域的台风资料,应用台风模型计算台风的气压场及风场;用于计算台风期间气压场变化的台风气压模型如下:
式中:p(r)为距台风中心r距离处的气压值;pc为台风中心气压,pn为台风外围气压(正常气压),Rmw为台风最大风速半径;B为台风的形状参数,通常取值范围为1<B<2.5;
气压场确定之后,用于计算理论台风风场的公式如下:
式中:Vg(r)为距台风中心r距离处的风速;ρA为大气的密度;
步骤3,大尺度海域台风风暴潮模型验证:针对台风引起的风暴潮进行模拟计算,通过对台风期间风暴潮增减水过程实测值与计算值的比较,进行大尺度海域数学模型的验证,确保模型的可靠性;
步骤4,建立嵌套的中尺度工程海区水动力数学模型:根据电厂取排水工程所在海域的位置、岸线及水下地形等条件确定模型范围;建立适用于工程海区流场模拟计算的中尺度工程海区数学模型;模型的控制方程与步骤1相同;中尺度工程海区模型中水位和流量的边界条件由大尺度海域模型的计算结果提取得到;
步骤5,嵌套模型的潮流验证:应用大尺度海域数学模型,进行大范围海域内常规气象条件下天文潮为主的水动力模拟计算,从计算结果中提取中尺度工程海区模型所需的潮位和流量边界条件;在此基础上,应用中尺度工程海区模型模拟计算常规气象条件下厂址海域潮流场;通过对厂址近区潮位、流速、流向的实测值与计算值的比较,进行嵌套模型的验证,确保嵌套模型的可靠性;
步骤6,建立中尺度工程海区粒子追踪数学模型:以嵌套模型的中尺度工程海区水动力模型为基础,采用示踪粒子的方法,建立中尺度工程海区漂浮物运移轨迹模拟的粒子追踪数学模型;模拟计算中主要考虑水流的携带作用,不考虑漂浮物形状等的影响;粒子运动轨迹模拟控制方程如下:
dXt=a(t,Xt)dt+b(t,Xt)ξtdt
式中:Xt为t时刻粒子位置,a(t,Xt)为流、风影响产生的漂流移动项,b(t,Xt)为分子运动、湍流产生的扩散项,ξt为随机数;
步骤7,电厂取水卷吸率的计算:以台风的气压场和风场作为输入条件,模拟台风引起的风暴潮与天文潮共同作用下厂址工程海区内粒子的运移轨迹,分析粒子在电厂取水明渠内的分布情况,统计进入取水明渠内的粒子数,计算电厂取水卷吸率随时间的变化,所述的卷吸率ψ为进入取水明渠内的粒子数k与当前已投放粒子总数K之比:
本发明产生的有益效果是:本发明针对受台风影响的滨海核电厂为研究对象,建立适用于风暴潮模拟的大尺度海域模型与用于工程近区流场模拟的中尺度海区模型嵌套的二维数学模型,采用示踪粒子方法研究台风引起的风暴潮与天文潮共同作用下漂浮物的运移轨迹,分析漂浮物在电厂取水明渠内的分布情况,在此基础上计算电厂取水卷吸率。为核电厂取水安全预警以及应急措施的规划和设计提供科学依据。
附图说明
下面结合附图和实施例对本发明作进一步说明。
图1是本发明的实施例所述方法的流程示意图;
图2是本发明的实施例所述使用举例中的台风路径示意图;
图3是本发明的实施例所述使用举例中的大尺度海域模型与嵌套的中尺度工程海区模型的模拟范围示意图;
图4是本发明的实施例所述使用举例中台风期间瞬时气压场(图4-a)和风场(图4-b);
图5是本发明的实施例所述使用举例中厂址附近台风风暴潮增减水验证图;
图6是本发明的实施例所述使用举例中取排水工程建设前常规气象条件下潮位观测站(W1)、流速和流向观测站(V1、V2、V3)位置示意图;
图7是本发明的实施例所述使用举例中常规气象条件下潮位验证图;
图8是本发明的实施例所述使用举例中常规气象条件下流速、流向验证图;
图9是本发明的实施例所述使用举例中电厂取排水工程建设后取水明渠位置示意图
图10是本发明的实施例所述使用举例中台风风暴潮和天文潮共同作用下某时刻示踪粒子短时运移轨迹及取水明渠内粒子分布与无台风影响时结果的对比图;
图11是本发明的实施例所述使用举例中台风期间风暴潮和天文潮共同作用下进入取水明渠内的粒子数随时间变化过程线与无台风影响时结果的对比图;
图12是本发明的实施例所述使用举例中台风期间风暴潮和天文潮共同作用下电厂取水卷吸率随时间变化过程线与无台风影响时结果的对比图。
具体实施方式
实施例:
本实施例是一种风暴潮与天文潮共同作用下的电厂取水卷吸率计算方法,流程如图1所示。本实施例所述方法包括如下步骤:
步骤1:建立大尺度海域水动力数学模型:根据厂址所在海域的位置及台风路径确定模型范围,建立适用于台风风暴潮模拟计算的大范围海域水动力数学模型。平面直角坐标系下模型的控制方程如下:
连续方程:
运动方程:
x方向:
y方向:
上式中:t为时间,x,y为笛卡尔坐标;η为水位,h=d+η为总水深,d为静水深;为x方向的水深平均速度,/>为y方向的水深平均速度:/>u为x方向速度,v为y方向速度,z为笛卡尔坐标;us为源在x方向的排水速度;vs为源在y方向的排水速度;f为柯氏力参数,f=2Ωsinφ,其中:Ω为旋转角速度,φ为地理纬度;g为重力加速度;ρ为水的密度,ρ0为参考水密度;pa为大气压力;A为涡粘性系数;τb=(τbx,τby)为底部应力,cf为阻力系数,/>床底上部流速;τs=(τsx,τsy)为水面风应力,ρa为空气密度,cd为空气阻力系数,/>水面上空10m处风速;S为源汇项。
大尺度海域数学模型的开边界水位条件由全球潮汐预报模型结果提取得到,台风的气压场及风场作为输入条件;大尺度海域数学模型的计算结果将为后续嵌套的中尺度海区模型提供计算所需的边界条件。
应用实例:
某受台风影响的滨海核电厂,应用实例的厂址位置及台风路径如图2所示。利用大尺度海域模型与中尺度工程海区模型嵌套的二维数学模型,结合台风模型计算得到的气压场和风场,模拟计算台风风暴潮与天文潮共同作用下工程近区水动力特性。在此基础上采用示踪粒子的方法,应用粒子追踪数学模型模拟风暴潮与天文潮共同作用下厂址工程海区内漂浮物的运移轨迹,计算电厂取水卷吸率。本实施例所述的大尺度水域是指范围在100万平方千米以上的水域和气象范围。所述的中尺度水域是指0.5~5万平方千米的水域,相对而言的小尺度水域是指100平方千米以下的水域。
应用实例中建立大尺度海域水动力数学模型的步骤包括:根据厂址所在海域位置及台风路径确定大尺度海域模型范围,应用实例的模型范围如图3所示(大约230万平方千米)。采用世界上广泛应用的MIKE21建立适用于台风风暴潮模拟计算的大范围海域水动力数学模型。模型的开边界水位条件由用于天文潮预报的全球潮汐模型的结果提取得到。
步骤2,台风的气压场及风场计算:分析电厂所在工程海域的台风资料,应用台风模型计算台风的气压场及风场;用于计算台风期间气压场变化的台风气压模型如下:
式中:p(r)为距台风中心r距离处的气压值;pc为台风中心气压,pn为台风外围气压(正常气压),Rmw为台风最大风速半径;B为台风的形状参数,通常取值范围为1<B<2.5;
气压场确定之后,用于计算理论台风风场的公式如下:
式中:Vg(r)为距台风中心r距离处的风速;ρA为大气的密度。
应用实例的台风期间气压场及风场计算的步骤包括:分析电厂所在工程海域的台风资料,依据台风的路径、最大风速、中心气压、台风外围气压等数据应用Holland台风模型计算台风的气压场。在气压场确定之后,应用理论台风风场公式计算风场。应用实例的台风气压场及风场如图4所示。计算得到的台风期间气压场及风场为大尺度海域水动力数学模型提供输入条件。
步骤3,大尺度海域台风风暴潮模型验证:对台风引起的风暴潮进行模拟计算,通过对台风风暴潮增减水过程的实测值与计算值的比较,进行大尺度海域数学模型的验证,确保模型的可靠性。
应用实例的大尺度海域台风风暴潮模型验证的步骤包括:应用大尺度海域水动力数学模型对台风引起的风暴潮进行模拟计算,通过台风期间厂址附近观测站(位置见图4)的风暴潮增减水过程实测值与计算值的比较,对大尺度海域数学模型进行验证。应用实例的风暴潮增减水验证如图5所示,良好的验证结果证明了大尺度海域水动力数学模型的可靠性。
步骤4,建立嵌套的中尺度工程海区水动力数学模型:根据电厂取排水工程所在海域的位置、岸线及水下地形等条件确定模型范围(大约1.2万平方千米)。建立适用于工程海区流场模拟计算的中尺度工程海区数学模型。中尺度工程海区模型所需的水位和流量边界条件由大尺度海域模型的计算结果提取得到。
应用实例建立嵌套的中尺度工程海区水动力数学模型的步骤包括:根据厂址取排水工程所在海域的位置、岸线及水下地形等条件确定中尺度工程海区模型范围,应用实例的模型范围见图3所示。采用MIKE21建立适用于工程海区流场模拟计算的中尺度工程海区水动力数学模型。模型所需的水位和流量边界条件由大尺度海域模型的计算结果提取得到。
步骤5,嵌套模型的潮流验证:应用大尺度海域数学模型,进行大范围海域内常规气象条件下天文潮为主的水动力模拟计算,从计算结果中提取中尺度工程海区模型所需的潮位和流量边界条件。在此基础上,应用中尺度工程海区模型模拟计算常规气象条件下厂址海域潮流场。通过对厂址近区潮位、流速、流向的实测值与计算值的比较,进行嵌套模型的验证,确保嵌套模型的可靠性。
应用实例的嵌套模型潮流验证的步骤包括:首先应用大尺度海域数学模型进行大范围海域内常规气象条件下天文潮为主的水动力模拟计算,从计算结果中提取中尺度工程海区模型边界位置处的潮位和流量过程作为嵌套的中尺度工程海区模型的边界条件。在此基础上,应用中尺度工程海区模型模拟计算常规气象条件下厂址海域潮流场。通过厂址近区观测站的潮位、流速和流向的实测值与计算值比较,进行嵌套模型的验证。应用实例的工程建设前厂址附近观测站位置如图6所示,潮位验证结果如图7所示,流速、流向的验证结果如图8所示。良好的验证结果证明了嵌套模型的可靠性。
步骤6,建立中尺度工程海区粒子追踪数学模型:以嵌套模型的中尺度工程海区水动力模型为基础,采用示踪粒子的方法,建立用于漂浮物运移轨迹模拟的中尺度工程海区粒子追踪数学模型。
应用实例建立中尺度工程海区粒子追踪数学模型的步骤包括:针对取排水工程建设后的地形(取、排水明渠位置如图9所示),以嵌套模型的中尺度工程海区水动力模型为基础,采用示踪粒子的方法,应用MIKE21的粒子追踪模块,建立用于漂浮物运移轨迹模拟的中尺度工程海区粒子追踪数学模型。应用实例中粒子的投放位置为位于海湾顶部断面上200个均匀分布的点位,断面位置如图9所示。粒子的投放方式为每个计算时间步长(5分钟)在每个点位各投放1个粒子,连续投放约6天,共计投放粒子的总数为344600个。
步骤7,电厂取水卷吸率的计算:以台风的气压场和风场作为输入条件,模拟台风风暴潮与天文潮共同作用下厂址工程海区内粒子的运移轨迹,分析粒子在电厂取水明渠内的分布情况,统计进入取水明渠内的粒子数,计算电厂取水卷吸率随时间的变化,所述的卷吸率ψ为进入取水明渠内的粒子数k与当前已投放粒子总数K之比:
应用实例的电厂取水卷吸率计算分析的步骤包括:以台风的气压场和风场作为输入条件,模拟台风风暴潮与天文潮共同作用下厂址工程海区内粒子的运移轨迹。统计进入电厂取水明渠内的粒子数随时间的变化。定义电厂取水卷吸率为进入取水明渠内的粒子数与当前已投放粒子总数之比,以此计算电厂取水卷吸率随时间的变化。应用实例的台风期间风暴潮和天文潮共同作用下某时刻示踪粒子短时运移轨迹及取水明渠内粒子分布如图10所示,图中显示了与无台风影响时结果的对比。结果表明台风对粒子的运移及取水明渠内粒子聚集情况的影响显著。应用实例的台风风暴潮与天文潮共同作用下进入电厂取水明渠内的粒子数随时间变化过程线如图11所示,计算得到的电厂取水卷吸率随时间变化过程线如图12所示,图中显示了与无台风影响时结果的对比。风暴潮与天文潮共同作用下最大取水卷吸率为3.42%,无台风影响时天文潮作用下最大取水卷吸率为3.02%,台风影响的最大取水卷吸率增量为0.73%。
最后应说明的是,以上仅用以说明本发明的技术方案而非限制,尽管参照较佳布置方案对本发明进行了详细说明,本领域的普通技术人员应当理解,可以对本发明的技术方案(比如模型的选择、公式的运用、步骤的先后顺序等)进行修改或者等同替换,而不脱离本发明技术方案的精神和范围。
Claims (1)
1.一种风暴潮与天文潮共同作用下的电厂取水卷吸率计算方法,其特征在于,所述方法的步骤如下:
步骤1,建立大尺度海域水动力数学模型:根据厂址所在海域的位置及台风路径确定模型范围,建立适用于台风风暴潮模拟计算的大范围海域水动力数学模型,平面直角坐标系下模型的控制方程如下:
连续方程:
运动方程:
x方向:
y方向:
上式中:t为时间,x,y为笛卡尔坐标;η为水位,h=d+η为总水深,d为静水深;为x方向的水深平均速度,/>为y方向的水深平均速度:/>z为笛卡尔坐标;us为源在x方向的排水速度;vs为源在y方向的排水速度;f为柯氏力参数,f=2Ωsinφ,其中:Ω为旋转角速度,φ为地理纬度;g为重力加速度;ρ为水的密度,ρ0为参考水密度;pa为大气压力;/>A为涡粘性系数;τb=(τbx,τby)为底部应力,/>cf为阻力系数,/>床底上部流速;τs=(τsx,τsy)为水面风应力,ρa为空气密度,cd为空气阻力系数,/>水面上空10m处风速;S为源汇项;
大尺度海域数学模型的开边界水位条件由全球潮汐预报模型结果提取得到,台风的气压场及风场作为输入条件;大尺度海域数学模型的计算结果将为后续嵌套的中尺度海区模型提供计算所需的边界条件;
步骤2,台风的气压场及风场计算:分析电厂所在工程海域的台风资料,应用台风模型计算台风的气压场及风场;用于计算台风期间气压场变化的台风气压模型如下:
式中:p(r)为距台风中心r距离处的气压值;pc为台风中心气压,pn为台风外围气压(正常气压),Rmw为台风最大风速半径;B为台风的形状参数,通常取值范围为1<B<2.5;
气压场确定之后,用于计算理论台风风场的公式如下:
式中:Vg(r)为距台风中心r距离处的风速;ρA为大气的密度;
步骤3,大尺度海域台风风暴潮模型验证:针对台风引起的风暴潮进行模拟计算,通过对台风期间风暴潮增减水过程实测值与计算值的比较,进行大尺度海域数学模型的验证,确保模型的可靠性;
步骤4,建立嵌套的中尺度工程海区水动力数学模型:根据电厂取排水工程所在海域的位置、岸线及水下地形等条件确定模型范围;建立适用于工程海区流场模拟计算的中尺度工程海区数学模型;模型的控制方程与步骤1相同;中尺度工程海区模型中水位和流量的边界条件由大尺度海域模型的计算结果提取得到;
步骤5,嵌套模型的潮流验证:应用大尺度海域数学模型,进行大范围海域内常规气象条件下天文潮为主的水动力模拟计算,从计算结果中提取中尺度工程海区模型所需的潮位和流量边界条件;在此基础上,应用中尺度工程海区模型模拟计算常规气象条件下厂址海域潮流场;通过对厂址近区潮位、流速、流向的实测值与计算值的比较,进行嵌套模型的验证,确保嵌套模型的可靠性;
步骤6,建立中尺度工程海区粒子追踪数学模型:以嵌套模型的中尺度工程海区水动力模型为基础,采用示踪粒子的方法,建立中尺度工程海区漂浮物运移轨迹模拟的粒子追踪数学模型;模拟计算中主要考虑水流的携带作用,不考虑漂浮物形状等的影响;粒子运动轨迹模拟控制方程如下:
dXt=a(t,Xt)dt+b(t,Xt)ξtdt
式中:Xt为t时刻粒子位置,a(t,Xt)为流、风影响产生的漂流移动项,b(t,Xt)为分子运动、湍流产生的扩散项,ξt为随机数;
步骤7,电厂取水卷吸率的计算:以台风的气压场和风场作为输入条件,模拟台风引起的风暴潮与天文潮共同作用下厂址工程海区内粒子的运移轨迹,分析粒子在电厂取水明渠内的分布情况,统计进入取水明渠内的粒子数,计算电厂取水卷吸率随时间的变化,所述的卷吸率ψ为进入取水明渠内的粒子数k与当前已投放粒子总数K之比:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310148210.4A CN116151085B (zh) | 2023-02-22 | 2023-02-22 | 风暴潮与天文潮共同作用下的电厂取水卷吸率计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310148210.4A CN116151085B (zh) | 2023-02-22 | 2023-02-22 | 风暴潮与天文潮共同作用下的电厂取水卷吸率计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN116151085A CN116151085A (zh) | 2023-05-23 |
CN116151085B true CN116151085B (zh) | 2023-08-11 |
Family
ID=86350422
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310148210.4A Active CN116151085B (zh) | 2023-02-22 | 2023-02-22 | 风暴潮与天文潮共同作用下的电厂取水卷吸率计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116151085B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117892659B (zh) * | 2024-01-29 | 2024-06-25 | 中国水利水电科学研究院 | 一种波流共同作用下电厂取水堵塞物相对密集度计算方法 |
CN117993716A (zh) * | 2024-02-20 | 2024-05-07 | 中国水利水电科学研究院 | 基于反向流场溯源的电厂取水卷吸风险概率评估方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102288729A (zh) * | 2011-06-28 | 2011-12-21 | 河海大学 | 一种淤泥质海岸浮泥形成条件的定量计算方法 |
RU2652642C1 (ru) * | 2016-12-19 | 2018-04-28 | Василий Николаевич Храмушин | Способ построения гидродинамической модели источника штормовых нагонов и экстремальных течений под воздействием подвижного атмосферного тайфуна |
CN112001069A (zh) * | 2020-08-10 | 2020-11-27 | 广州大学 | 一种轴非对称台风风场模拟的方法 |
CN112085270A (zh) * | 2020-09-04 | 2020-12-15 | 东南大学 | 基于随机统计模型与水动力模型耦合的风暴潮极端水位预测方法 |
CN112257352A (zh) * | 2020-10-21 | 2021-01-22 | 黄河水利委员会黄河水利科学研究院 | 一维水动力模型和二维水动力模型的耦合方法及系统 |
CN113901557A (zh) * | 2021-10-20 | 2022-01-07 | 中国水利水电科学研究院 | 一种基于降低取水卷吸效应的电厂取水明渠布设分析方法 |
-
2023
- 2023-02-22 CN CN202310148210.4A patent/CN116151085B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102288729A (zh) * | 2011-06-28 | 2011-12-21 | 河海大学 | 一种淤泥质海岸浮泥形成条件的定量计算方法 |
RU2652642C1 (ru) * | 2016-12-19 | 2018-04-28 | Василий Николаевич Храмушин | Способ построения гидродинамической модели источника штормовых нагонов и экстремальных течений под воздействием подвижного атмосферного тайфуна |
CN112001069A (zh) * | 2020-08-10 | 2020-11-27 | 广州大学 | 一种轴非对称台风风场模拟的方法 |
CN112085270A (zh) * | 2020-09-04 | 2020-12-15 | 东南大学 | 基于随机统计模型与水动力模型耦合的风暴潮极端水位预测方法 |
CN112257352A (zh) * | 2020-10-21 | 2021-01-22 | 黄河水利委员会黄河水利科学研究院 | 一维水动力模型和二维水动力模型的耦合方法及系统 |
CN113901557A (zh) * | 2021-10-20 | 2022-01-07 | 中国水利水电科学研究院 | 一种基于降低取水卷吸效应的电厂取水明渠布设分析方法 |
Non-Patent Citations (1)
Title |
---|
江苏海域"达维"台风风暴潮数值模拟试验研究;刘文坚;李瑞杰;林祥;董啸天;李玉婷;;水运工程(第03期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN116151085A (zh) | 2023-05-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN116151085B (zh) | 风暴潮与天文潮共同作用下的电厂取水卷吸率计算方法 | |
Jeong et al. | Salinity intrusion characteristics analysis using EFDC model in the downstream of Geum River | |
Shaji et al. | Storm surge studies in the North Indian Ocean: A review | |
Abbs | Sea-breeze interactions along a concave coastline in southern Australia: Observations and numerical modeling study | |
CN110414041B (zh) | 一种基于gis技术建立风暴潮漫滩和淹没分析的方法及系统 | |
Steedman et al. | Wind-driven circulation of Cockburn Sound | |
CN112085270A (zh) | 基于随机统计模型与水动力模型耦合的风暴潮极端水位预测方法 | |
CN114510765B (zh) | 一种台风风暴潮作用下的海堤漫溃堤动态淹没预报方法 | |
Xu et al. | Effects of tide-surge interactions on storm surges along the coast of the Bohai Sea, Yellow Sea, and East China Sea | |
CN107589220B (zh) | 一种适用于狭长航道的船舶尾气监测站选址方法 | |
Li et al. | Observed sub-hectometer-scale low level jets in surface-layer velocity profiles of landfalling typhoons | |
Zhang et al. | Simulation and analysis of Back siltation in a navigation channel using MIKE 21 | |
Yue et al. | Effects of coastal reclamation on the storm surge in the Bohai Bay | |
CN111611641B (zh) | 一种考虑余流的港口建设方法 | |
CN117892659B (zh) | 一种波流共同作用下电厂取水堵塞物相对密集度计算方法 | |
Schoenbaechler et al. | TxBLEND model calibration and validation for the Laguna Madre Estuary | |
Schoenbaechler et al. | TxBLEND model calibration and validation for the Nueces estuary | |
Yuk et al. | Modelling of storm-induced seawater flooding in the Suyeong River area, South Korea: A case study due to the storm surge and waves during Typhoon Sanba | |
CN117993716A (zh) | 基于反向流场溯源的电厂取水卷吸风险概率评估方法 | |
Kwon et al. | Simulation of vertical water temperature structure variation in the Yellow Sea typhoon MEARI | |
CN117709131B (zh) | 近岸风暴潮高精度批量数值模拟方法及系统 | |
Jinshan et al. | Numerical simulation of storm surge in the submerged radial sandbars of the Southern Yellow Sea, China | |
Averkiev et al. | Calculation of extreme water levels in the eastern part of the Gulf of Finland | |
CN112434423B (zh) | 一种同心圆网格与新型台风场模式结合的风暴潮模拟方法 | |
Wang et al. | Development of an oil spill forecast system for offshore China |
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 |