CN107944608B - 基于卫星遥感的海面漂移物与溢油漂移扩散预报方法 - Google Patents
基于卫星遥感的海面漂移物与溢油漂移扩散预报方法 Download PDFInfo
- Publication number
- CN107944608B CN107944608B CN201711133263.XA CN201711133263A CN107944608B CN 107944608 B CN107944608 B CN 107944608B CN 201711133263 A CN201711133263 A CN 201711133263A CN 107944608 B CN107944608 B CN 107944608B
- Authority
- CN
- China
- Prior art keywords
- drift
- forecasting
- oil
- path
- diffusion
- 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
- 238000009792 diffusion process Methods 0.000 title claims abstract description 69
- 239000003305 oil spill Substances 0.000 title claims abstract description 64
- 238000013277 forecasting method Methods 0.000 title claims abstract description 6
- 238000000034 method Methods 0.000 claims abstract description 31
- 230000007246 mechanism Effects 0.000 claims abstract description 19
- 150000003839 salts Chemical class 0.000 claims abstract description 16
- 239000002245 particle Substances 0.000 claims description 38
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 15
- 230000002093 peripheral effect Effects 0.000 claims description 13
- 238000005516 engineering process Methods 0.000 claims description 10
- 238000004458 analytical method Methods 0.000 claims description 5
- 238000005295 random walk Methods 0.000 claims description 5
- 238000007689 inspection Methods 0.000 claims description 4
- 238000012800 visualization Methods 0.000 claims description 4
- 230000009182 swimming Effects 0.000 claims description 2
- 239000003921 oil Substances 0.000 description 42
- 238000004364 calculation method Methods 0.000 description 8
- 230000008569 process Effects 0.000 description 8
- 230000009471 action Effects 0.000 description 5
- 238000002474 experimental method Methods 0.000 description 5
- 238000011161 development Methods 0.000 description 4
- 230000004907 flux Effects 0.000 description 4
- 230000000007 visual effect Effects 0.000 description 4
- 230000008859 change Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000003672 processing method Methods 0.000 description 3
- 239000002344 surface layer Substances 0.000 description 3
- 238000012937 correction Methods 0.000 description 2
- 238000009826 distribution Methods 0.000 description 2
- 239000003550 marker Substances 0.000 description 2
- 239000000047 product Substances 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 230000002411 adverse Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000005094 computer simulation Methods 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 238000013480 data collection Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000003111 delayed effect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000004836 empirical method Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 230000008020 evaporation Effects 0.000 description 1
- 238000001704 evaporation Methods 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 239000000295 fuel oil Substances 0.000 description 1
- 230000004927 fusion Effects 0.000 description 1
- 238000007654 immersion Methods 0.000 description 1
- 239000013067 intermediate product Substances 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000009364 mariculture Methods 0.000 description 1
- 238000013508 migration Methods 0.000 description 1
- 230000005012 migration Effects 0.000 description 1
- 238000012821 model calculation Methods 0.000 description 1
- 238000003786 synthesis reaction Methods 0.000 description 1
- 238000012876 topography Methods 0.000 description 1
- 238000009827 uniform distribution Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
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
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
- G06F16/20—Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
- G06F16/29—Geographical information databases
-
- 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
- G06Q50/02—Agriculture; Fishing; Forestry; Mining
-
- 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
- G06Q50/10—Services
- G06Q50/26—Government or public services
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Theoretical Computer Science (AREA)
- Tourism & Hospitality (AREA)
- Strategic Management (AREA)
- Physics & Mathematics (AREA)
- Human Resources & Organizations (AREA)
- General Physics & Mathematics (AREA)
- Economics (AREA)
- General Business, Economics & Management (AREA)
- Marketing (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Health & Medical Sciences (AREA)
- Development Economics (AREA)
- Databases & Information Systems (AREA)
- Life Sciences & Earth Sciences (AREA)
- Marine Sciences & Fisheries (AREA)
- Remote Sensing (AREA)
- General Engineering & Computer Science (AREA)
- Educational Administration (AREA)
- Agronomy & Crop Science (AREA)
- Animal Husbandry (AREA)
- Data Mining & Analysis (AREA)
- Mining & Mineral Resources (AREA)
- Game Theory and Decision Science (AREA)
- Entrepreneurship & Innovation (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Testing Or Calibration Of Command Recording Devices (AREA)
- Indicating Or Recording The Presence, Absence, Or Direction Of Movement (AREA)
Abstract
本发明涉及一种基于卫星遥感的海面漂移物与溢油漂移扩散预报方法,包括:S1:收集资料;S2:建立三维温、盐、流数值预报模型,获取海表流场;S3:建立海面漂移物漂移动态预报模式,得到漂移物的漂移路径结果;S4:建立溢油漂移扩散动态预报模式,得到溢油的漂移路径与扩散结果;S5:获取不同预报机构的漂移物的漂移路径结果,结合步骤S3获得的漂移物的漂移路径结果,获取最优的漂移物的漂移路径结果;同时获取不同预报机构的漂移物的漂移路径结果,结合步骤S4获得的溢油的漂移路径与扩散结果,获取最优的溢油的漂移路径与扩散结果。与现有技术相比,本发明克服一家预报机构的局限性,综合多条漂移扩散预报结果,提高预报精度。
Description
技术领域
本发明涉及海上搜救技术领域、溢油漂移扩散预报技术领域,尤其是涉及一种基于卫星遥感的海面漂移物与溢油漂移扩散预报方法。
背景技术
伴随着海洋开发规模的扩大及我国海洋石油工业和海上石油运输业的迅猛发展,海上石油勘探、开发、海底管线铺设的规模不断扩大,海事活动的频繁,人员落水、船舶以及各种海上溢油等海上突发事故时有发生,海难事故更为人们所关注。海难事故在政治、经济、军事上都会给世界沿海各国带来巨大灾难,不仅是人员伤亡和财产的损失,还会给社会发展带来不良影响。一旦发生海难事故,会造成严重人员伤亡和巨大经济损失,海上溢油事故还会破坏海洋生态环境,影响海水养殖生产。海洋与渔业、交通、海事等部门对海上突发事故应急保障预报服务提出了很高的要求,既有海洋环境保障预报服务的要求,又有落水人员和溢油漂移轨迹精细化保障预报的要求,以实现应对海上突发事故“精准预报、高效救援”的目标,为此,迫切需要开展覆盖该区域海上突发事故应急辅助决策服务技术的研究。
由于针对海面漂移物和溢油漂移扩散的快速预报技术涉及高效的整合物理海洋环境相关的动力学计算以及全自动的数据下载以及实时的图形可视化系统,目前关于海面漂移物和溢油漂移扩散预报的相关技术多基于经验判断和观察追踪,成本较高且预报效率较低,并且由于没有耦合海洋动力学模型致使在物理海洋环境急剧变化的情况下无法准确、快速、高效的预报海面漂移物和溢油漂移路径以及扩散情况,这在很大程度上影响了指挥与协调工作的快速性和准确性,往往会使海上突发事故紧急处理行动贻误时机。
发明内容
本发明的目的就是为了克服上述现有技术存在的缺陷而提供一种基于卫星遥感的海面漂移物与溢油漂移扩散预报方法,克服一家预报机构的局限性,综合多条漂移扩散预报结果,提高预报精度,为海洋与渔业、交通、海事等部门海上突发事故紧急处理等相关决策提供更多种的漂移扩散路径参考。
本发明的目的可以通过以下技术方案来实现:
基于卫星遥感的海面漂移物与溢油漂移扩散预报方法包括以下步骤:
S1:收集海面漂移物及溢油的历史资料,并获取相应时刻的卫星遥感数据;
S2:建立近海及周边海域三维温、盐、流数值预报模型,获取海表流场;
S3:建立海面漂移物漂移动态预报模式,根据步骤S1的历史资料、卫星遥感数据和步骤S2的海表流场得到漂移物的漂移路径结果;
S4:建立溢油漂移扩散动态预报模式,根据步骤S1的历史资料、卫星遥感数据和步骤S2的海表流场得到溢油的漂移路径与扩散结果;
S5:获取不同预报机构的漂移物的漂移路径结果,结合步骤S3获得的漂移物的漂移路径结果,获取最优的漂移物的漂移路径结果;同时获取不同预报机构的漂移物的漂移路径结果,结合步骤S4获得的溢油的漂移路径与扩散结果,获取最优的溢油的漂移路径与扩散结果。
所述步骤S2中对近海及周边海域三维温、盐、流数值预报模型的海洋与风浪数据进行数据同化,为近海及周边海域三维温、盐、流数值预报模型提供初始场和检验场。
所述漂移物包括落水人员和遇险船只,所述步骤S3中海面漂移物漂移动态预报模式包括漂移物漂移运动公式,满足以下公式:
Xn+1=Xn+UΔt+Δα
U=k·W
式中,Xn为漂移物在t时刻的位置向量,Xn+1为漂移物在t+Δt时刻的位置向量,Δt为随机游动时间步长,U为在瞬时位置Xn时的风压漂流速度,Δα为α方向上的随机运动距离,α代表海面上x或y方向,k为设定的系数,W为海面风速,R(t)为(-1,1)之间的均匀分布随机数,Kα为α方向上的扰动系数。
所述步骤S4中溢油漂移扩散动态预报模式包括溢油中每个油粒子的水平二维空间位置公式,满足以下公式:
式中,Y为溢油中每个油粒子的位置向量,为油粒子在海面上x、y方向上的位置坐标,为油粒子在海面上x、y方向上的初始位置坐标,ux、uy为在海面上x、y方向上的流速分量,Wx、Wy为在海面上x、y方向上的海面10米风速,A为风拖曳因子,β为风拖曳角,Δt为随机游动时间步长,为在海面上x、y方向上的湍流扩散。
所述风拖曳因子A取值范围为2.5-3.5%,所述风拖曳角β取值范围为0°-25°。
所述步骤S5中最优的漂移物的漂移路径结果包括漂移物最优的位置向量(X)best,漂移物最优的位置向量(X)best满足以下公式:
式中,X为步骤S3获得的漂移物的漂移路径结果中的漂移物位置向量,μ0为对应X的漂移物置信因子,N为预报机构的总数,Xi′为对应第i个预报机构的漂移物的漂移路径结果中的漂移物位置向量,μi为对应第i个预报机构的漂移物置信因子,μ0,μi∈(0,1);
最优的溢油的漂移路径与扩散结果包括溢油中每个油粒子最优的位置向量(Y)best,溢油中每个油粒子最优的位置向量(Y)best满足以下公式:
式中,Y为步骤S4获得的溢油的漂移路径与扩散结果中的溢油中每个油粒子的位置向量,为对应Y的溢油置信因子,Yi′为对应第i个预报机构的溢油的漂移路径与扩散结果中的溢油中每个油粒子的位置向量,为对应第i个预报机构的溢油置信因子,
还包括步骤S6:通过可视化技术将最优的漂移物的漂移路径结果和最优的溢油的漂移路径与扩散结果直观展示出来,并生成包含发生地海域风场、浪场和流场的海洋环境分析报告。
所述步骤S6中,在地图上显示单线条路径和概率预报路径,所述单线条路径包括步骤S3获得的与不同预报机构的漂移物的漂移路径结果以及步骤S4获得的与不同预报机构的溢油的漂移路径与扩散结果,所述概率预报路径包括最优的漂移物的漂移路径结果和最优的溢油的漂移路径与扩散结果,同时不同时刻的单线条路径和概率预报路径存储为GIF动画格式的显示文件。
与现有技术相比,本发明具有以下优点:
(1)步骤S2中采用观测的地波雷达海表流场、卫星遥感融合海表温度和大浮标剖面速等数据,针对温盐流模型开展基于4D-Var同化技术的研究、应用、检验,可以较好地解决观测仪器的不稳定性对对同化数据的影响,同时优化并行算方法,解决4D-Var同化技术计算时间长,不利于业务化的问题,提高模型计算效率和模型预报准确率。
(2)步骤S3在海流数值模型提供表面流的基础上,通过拉格朗日法进行水质点的跟踪,同时在模型和卫星遥感观测风场的支持下,对各类型的定位物(如不同类型及吨位的船只、人员)进行风漂移纠正。对于不同类型的漂移物进行轨迹预报,需要在海上进行追踪实验,以便验证系统的准确性,根据实验获得的数据,对该系统进行进一步的完善优化。
(3)对于潮汐及环流动力所引起的海流,步骤S2的三维海流模型可以在事故发生前一个月甚至更早就给出计算结果。但是对于由风引起的风海流,由于风场的时效性短,只能是利用现场数据或气象模型预测数据进行实时计算。这种情况下,通过模型在较短时间内给出预报报表有一定难度,本发明步骤S4将已计算结果作为中间产品,通过现场的实时情况,如风速大小、漂移物形态等,引入半经验的风导系数和随机漂移数等方法,重新热启动计算,可短时间内给出有时效性的结果。
(4)鉴于一家预报机构的局限性,本发明步骤S5基于多家预报机构不同海洋数值模式及数值风浪模式来实现集合预报,能够给出多条漂移扩散预报结果,在此基础上计算多家漂移路径置信因子,进而也给出加权平均后的漂移扩散预报结果,为海洋与渔业、交通、海事等部门海上突发事故紧急处理等相关决策提供更多种的漂移扩散路径参考。
(5)通过步骤S6开发可视化GIS的海上突发事故应急辅助决策系统,实现海上失踪船只、落水人员及溢油事故漂移方位及轨迹信息等动态预报的可视化操作。
附图说明
图1为本发明方法流程图;
图2为本发明方法中数据流向示意图;
图3为高精度的近海及周边海域三维温、盐、流数值预报模型的示意图;
图4为高精度的近海及周边海域三维温、盐、流数值预报模型的模型参数调整示意图。
具体实施方式
下面结合附图和具体实施例对本发明进行详细说明。本实施例以本发明技术方案为前提进行实施,给出了详细的实施方式和具体的操作过程,但本发明的保护范围不限于下述的实施例。
如图1和图2所示,基于卫星遥感的海面漂移物与溢油漂移扩散预报方法包括以下步骤:
S1:资料收集
收集近十年来海上船只事故、人员失踪及海上溢油资料,并获取相应时刻的卫星遥感数据,建立海面漂移物和溢油漂移扩散预报平台的数据库,进行案例分析。
S2:建立高精度的近海及周边海域三维温、盐、流数值预报模型,获取对应预报机构的海表流场。
如图3所示,基于最新水深、岸线、地形等海洋基础信息,以全球模型提供外边界条件,借助高性能的计算机条件,以大小网格嵌套技术,建立高分辨率、高精度的近海及周边海域三维温盐流业务化数值预报模型。
本模型用大小网格嵌套方法引入非均一网格技术,建立具有港湾分辨率的三维海流数值模型,并采用欧共体Mercator全球模式或Hycom全球模式提供边界条件,以及海表的净热通量、水气通量和动量通量则采用国家海洋预报中心的MM5或WRF的产品进行计算。
步骤S2中对海洋与风浪数据进行数据同化,为近海及周边海域三维温、盐、流数值预报模型提供初始场和检验场。实际应用时,建立数值同化模块同化各类观测数据,提供近海及周边海域内丰富的浮标、卫星遥感及地波雷达观测数据(海面风、温和盐流观测值),采用4D-Var同化技术,进行数值模型的数据同化、调式及验证。针对不同的数据源,采用不同的方法同化这些数据,最终提高数值模型的温、盐、流预报精度及减小漂移物位置预报误差。
如图4所示,为了提高模型计算精度,本模型全面考虑潮汐、风场、温度、盐度、径流、热通量、降雨、蒸发等动力过程的作用和影响,优化数值模型参数条件,在垂直方向采用S坐标系,对水平扩散项进行纠正,以减少模型中海底地形剧烈变化所引起的虚假垂向混合。优化改进Yamada二阶半紊流扩散模型,解决引入内潮所导致的湍动能量变化,提供更为合理的垂向扩散系数,使模型能准确模拟近海区域。优化改进Yamada二阶半紊流扩散模型是采用数值实验的方式,采用不同的对流项处理方法,模拟了海湾潮流场,并将数模结果与物模模拟结果进行对比,分析不同对流项处理方式对海湾区域流态模拟结果的影响,以此为依据提出适合港湾的对流项处理方法。
S3:建立海面漂移物漂移动态预报模式,根据步骤S1的历史资料、卫星遥感数据和步骤S2的海表流场得到漂移物的漂移路径结果。
在高精度三维温盐流数值模型的基础上,通过拉格朗日法进行水质点的跟踪,实现水质点轨迹预测预报。在不同强度的各类风场作用下,进行海上现场实验和计算机模拟实验,建立不同类型漂移物的轨迹相对水质点轨迹的订正经验式及相关系数,以适应对不同类型漂移物(如船只、人员)的轨迹的动态预测预报。
基于拉格朗日粒子追踪方法研发搜救漂移预测模块,将落水人员、海面艇筏当作标识粒子,计算粒子随流体机械运动的轨迹,跟踪所释放粒子的位置随时间的变化。即将落水人员作为海面风和表层流作用下的被动漂移物体,只考虑其在水平方向上的漂移运动,得到海面漂移物漂移动态预报模式的漂移物漂移运动公式,满足以下公式:
Xn+1=Xn+UΔt+Δα
式中,Xn为漂移物在t时刻的位置向量,Xn+1为漂移物在t+Δt时刻的位置向量,Δt为随机游动时间步长,U为质点在瞬时位置Xn时的风压漂流速度,Δα为α方向上的随机运动距离,α代表海面上x或y方向,x与y方向相垂直。
在海流数值模型提供表面流的基础上,通过拉格朗日法进行水质点的跟踪,同时在模型和卫星遥感观测风场的支持下,对各类型的定位物(如不同类型及吨位的船只、人员)进行风漂移纠正。其中,风导漂移指由拖拽或风动压力产生的漂移物漂移,风动压力是指处于一定运动状态下的漂移物水面上所受的空气压力。风压漂流速度U的计算公式为
U=k·W
式中,k为设定的系数,W为海面风速,随着漂移物的不同,排水状态以及漂移物的浸没比例等因素变化而不同。风导漂移根据经验关系建立模型计算所得。
海洋中的湍流是造成海流运动具有不确定性的重要原因,海洋中湍流是由各种不同尺度的涡漩叠合而成的流动,这些涡漩的大小及旋转轴的方向分布是随机的。漂移物在湍流涡漩作用下的随机运动距离Δα可以由以下公式计算
式中,R(t)为(-1,1)之间的均匀分布随机数,Kα为α方向上的扰动系数。
S4:建立溢油漂移扩散动态预报模式,根据步骤S1的历史资料、卫星遥感数据和步骤S2的海表流场得到溢油的漂移路径与扩散结果。
建立能够描述近海及周边海域溢油漂移扩散动力学过程、适用于对溢油轨迹、油膜厚度及其影响区域进行动态预测的溢油漂移扩散动态预报模式,可以实现不同油品(重油、轻油)产生连续点源(长时间连续溢油)和瞬时点源(短时间突发溢油)溢油的漂移路径和扩散范围。
溢油的漂移是指在风、表层及次表层流合成的环境动力作用下的拉格朗日漂移过程。油膜用大量的标识油粒子表示,每个标识油粒子代表一定的油量,油粒子在Δt时间内将溢油运动过程分成两个主要的部分,即平流运动过程和扩散运动过程。用确定性方法模拟平流过程,用随机走动法模拟扩散过程。通过追踪粒子的坐标和轨迹,统计出油膜的漂移轨迹和面积。
则溢油漂移扩散动态预报模式的溢油中每个油粒子的水平二维空间位置公式,满足以下公式:
式中,Y为溢油中每个油粒子的位置向量,为油粒子在海面上x、y方向上的位置坐标,为油粒子在海面上x、y方向上的初始位置坐标,ux、uy为在海面上x、y方向上的流速分量,Wx、Wy为在海面上x、y方向上的海面10米风速,A为风拖曳因子,取值范围为2.5-3.5%,例如3%,β为风拖曳角,取值范围为0°-25°,例如:15°,Δt为随机游动时间步长,为在海面上x、y方向上的湍流扩散。
S5:获取不同预报机构的漂移物的漂移路径结果,结合步骤S3获得的漂移物的漂移路径结果,获取最优的漂移物的漂移路径结果;同时获取不同预报机构的漂移物的漂移路径结果,结合步骤S4获得的溢油的漂移路径与扩散结果,获取最优的溢油的漂移路径与扩散结果。
基于多家预报机构不同海洋数值模式及数值风浪模式来实现集合预报,能够给出多条漂移扩散预报结果,在此基础上结合历史海上船只事故、人员失踪及海上溢油资料,综合分析多家漂移路径对应的置信因子以及本发明步骤S3、S4获取的漂移路径对应的置信因子,进而也给出加权平均后的最优漂移扩散预报结果。
本次采用的加权平均方法具体公式如下:
步骤S5中最优的漂移物的漂移路径结果包括漂移物最优的位置向量(X)best,漂移物最优的位置向量(X)best满足以下公式:
式中,X为步骤S3获得的漂移物的漂移路径结果中的漂移物位置向量,μ0为对应X的漂移物置信因子,N为预报机构的总数,Xi′为对应第i个预报机构的漂移物的漂移路径结果中的漂移物位置向量,μi为对应第i个预报机构的漂移物置信因子,μ0,μi∈(0,1);
最优的溢油的漂移路径与扩散结果包括溢油中每个油粒子最优的位置向量(Y)best,溢油中每个油粒子最优的位置向量(Y)best满足以下公式:
式中,Y为步骤S4获得的溢油的漂移路径与扩散结果中的溢油中每个油粒子的位置向量,为对应Y的溢油置信因子,Yi′为对应第i个预报机构的溢油的漂移路径与扩散结果中的溢油中每个油粒子的位置向量,为对应第i个预报机构的溢油置信因子,
S6:通过可视化技术将最优的漂移物的漂移路径结果和最优的溢油的漂移路径与扩散结果直观展示出来,并生成包含发生地海域风场、浪场和流场的海洋环境分析报告。具体有:在地图上显示单线条路径、概率预报路径、浓度场分布三种,以及支持对位置点进行拖拽修改,避免漂移路径经过陆地现象,单线条路径包括步骤S3获得的与不同预报机构的漂移物的漂移路径结果以及步骤S4获得的与不同预报机构的溢油的漂移路径与扩散结果,概率预报路径包括最优的漂移物的漂移路径结果和最优的溢油的漂移路径与扩散结果,同时不同时刻的单线条路径和概率预报路径存储为GIF动画格式的显示文件。
建立基于可视化GIS的集成风场、海流场、波浪场及漂移物迁移的海上突发事故应急辅助决策系统,通过风场、海流边界条件、事故发生的时间和地点等信息的输入,实现海上失踪船只、落水人员及溢油事故漂移方位及轨迹信息等动态预报的便捷可视化操作,能够动态显示各时刻漂移物和溢油所在地点及漂移扩散等,同时形成发生地海域风场、浪场、流场等海洋环境分析报告,并实现电子邮件、传真、短信等多种方式的一键快速推送,为有关部门提供海上突发事故决策依据。
Claims (6)
1.基于卫星遥感的海面漂移物与溢油漂移扩散预报方法,其特征在于,包括以下步骤:
S1:收集海面漂移物及溢油的历史资料,并获取相应时刻的卫星遥感数据;
S2:建立近海及周边海域三维温、盐、流数值预报模型,获取海表流场;
S3:建立海面漂移物漂移动态预报模式,根据步骤S1的历史资料、卫星遥感数据和步骤S2的海表流场得到漂移物的漂移路径结果;
S4:建立溢油漂移扩散动态预报模式,根据步骤S1的历史资料、卫星遥感数据和步骤S2的海表流场得到溢油的漂移路径与扩散结果;
S5:获取不同预报机构的漂移物的漂移路径结果,结合步骤S3获得的漂移物的漂移路径结果,获取最优的漂移物的漂移路径结果;同时获取不同预报机构的漂移物的漂移路径结果,结合步骤S4获得的溢油的漂移路径与扩散结果,获取最优的溢油的漂移路径与扩散结果;
所述漂移物包括落水人员和遇险船只,所述步骤S3中海面漂移物漂移动态预报模式包括漂移物漂移运动公式,满足以下公式:
Xn+1=Xn+UΔt+Δα
U=k·W
式中,Xn为漂移物在t时刻的位置向量,Xn+1为漂移物在t+Δt时刻的位置向量,Δt为随机游动时间步长,U为在瞬时位置Xn时的风压漂流速度,Δα为α方向上的随机运动距离,α代表海面上x或y方向,k为设定的系数,W为海面风速,R(t)为(-1,1)之间的均匀分布随机数,Kα为α方向上的扰动系数;
所述步骤S4中溢油漂移扩散动态预报模式包括溢油中每个油粒子的水平二维空间位置公式,满足以下公式:
2.根据权利要求1所述的基于卫星遥感的海面漂移物与溢油漂移扩散预报方法,其特征在于,所述步骤S2中对近海及周边海域三维温、盐、流数值预报模型的海洋与风浪数据进行数据同化,为近海及周边海域三维温、盐、流数值预报模型提供初始场和检验场。
3.根据权利要求1所述的基于卫星遥感的海面漂移物与溢油漂移扩散预报方法,其特征在于,所述风拖曳因子A取值范围为2.5-3.5%,所述风拖曳角β取值范围为0°-25°。
4.根据权利要求1所述的基于卫星遥感的海面漂移物与溢油漂移扩散预报方法,其特征在于,所述步骤S5中最优的漂移物的漂移路径结果包括漂移物最优的位置向量(X)best,漂移物最优的位置向量(X)best满足以下公式:
式中,X为步骤S3获得的漂移物的漂移路径结果中的漂移物位置向量,μ0为对应X的漂移物置信因子,N为预报机构的总数,Xi′为对应第i个预报机构的漂移物的漂移路径结果中的漂移物位置向量,μi为对应第i个预报机构的漂移物置信因子,μ0,μi∈(0,1);
最优的溢油的漂移路径与扩散结果包括溢油中每个油粒子最优的位置向量(Y)best,溢油中每个油粒子最优的位置向量(Y)best满足以下公式:
5.根据权利要求1所述的基于卫星遥感的海面漂移物与溢油漂移扩散预报方法,其特征在于,还包括步骤S6:通过可视化技术将最优的漂移物的漂移路径结果和最优的溢油的漂移路径与扩散结果直观展示出来,并生成包含发生地海域风场、浪场和流场的海洋环境分析报告。
6.根据权利要求5所述的基于卫星遥感的海面漂移物与溢油漂移扩散预报方法,其特征在于,所述步骤S6中,在地图上显示单线条路径和概率预报路径,所述单线条路径包括步骤S3获得的不同预报机构的漂移物的漂移路径结果以及步骤S4获得的不同预报机构的溢油的漂移路径与扩散结果,所述概率预报路径包括最优的漂移物的漂移路径结果和最优的溢油的漂移路径与扩散结果,同时不同时刻的单线条路径和概率预报路径存储为GIF动画格式的显示文件。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711133263.XA CN107944608B (zh) | 2017-11-16 | 2017-11-16 | 基于卫星遥感的海面漂移物与溢油漂移扩散预报方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711133263.XA CN107944608B (zh) | 2017-11-16 | 2017-11-16 | 基于卫星遥感的海面漂移物与溢油漂移扩散预报方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107944608A CN107944608A (zh) | 2018-04-20 |
CN107944608B true CN107944608B (zh) | 2021-12-17 |
Family
ID=61932395
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711133263.XA Active CN107944608B (zh) | 2017-11-16 | 2017-11-16 | 基于卫星遥感的海面漂移物与溢油漂移扩散预报方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107944608B (zh) |
Families Citing this family (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR101740720B1 (ko) * | 2015-11-17 | 2017-05-29 | 한국해양과학기술원 | 위험유해물질의 확산 예측 장치 및 방법 |
CN108875185A (zh) * | 2018-06-06 | 2018-11-23 | 牟林 | 一种基于多源数据建立海上遇险目标风漂模型的方法 |
CN110533217A (zh) * | 2019-07-22 | 2019-12-03 | 北京金交信息通信导航设计院 | 一种基于遥感的海上溢油污染漂移预测系统 |
CN110399676B (zh) * | 2019-07-24 | 2021-12-28 | 国家海洋环境预报中心 | 西北太平洋三维溢油业务化应急预报及评估系统 |
CN111368414B (zh) * | 2020-02-28 | 2022-04-08 | 天津大学 | 一种基于厚度信息的海上溢油轨迹模拟方法 |
CN111239003B (zh) * | 2020-03-26 | 2022-03-22 | 上海勘测设计研究院有限公司 | 油膜漂移模拟装置 |
CN111855942A (zh) * | 2020-07-25 | 2020-10-30 | 山东交通学院 | 一种基于3s技术的船舶与海洋溢油污染的监测系统 |
CN112232705B (zh) * | 2020-11-05 | 2023-09-01 | 交通运输部规划研究院 | 一种基于不确定性分析的区域溢油风险分析方法和装置 |
CN112445856B (zh) * | 2020-12-01 | 2023-04-21 | 海南长光卫星信息技术有限公司 | 一种海面高度影响关联分析方法和装置 |
CN115186569B (zh) * | 2022-07-07 | 2023-01-31 | 中科三清科技有限公司 | 漂浮物的漂移模拟方法、装置、存储介质及电子设备 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101706972A (zh) * | 2009-11-16 | 2010-05-12 | 大连海事大学 | 海上溢油的三维可视化算法 |
CN102024200A (zh) * | 2010-09-03 | 2011-04-20 | 焦作大学 | 河流溢油扩散动态预测技术 |
CN102156817A (zh) * | 2011-04-11 | 2011-08-17 | 国家海洋信息中心 | 海上溢油漂移路径及扩散快速预报方法 |
CN105160681A (zh) * | 2015-09-09 | 2015-12-16 | 浙江海洋学院 | 一种海面溢油监测方法及装置 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20150130460A1 (en) * | 2013-11-13 | 2015-05-14 | Schlumberger Technology Corporation | Methods for separating oil and water on multidimensional nuclear magnetic resonance maps |
-
2017
- 2017-11-16 CN CN201711133263.XA patent/CN107944608B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101706972A (zh) * | 2009-11-16 | 2010-05-12 | 大连海事大学 | 海上溢油的三维可视化算法 |
CN102024200A (zh) * | 2010-09-03 | 2011-04-20 | 焦作大学 | 河流溢油扩散动态预测技术 |
CN102156817A (zh) * | 2011-04-11 | 2011-08-17 | 国家海洋信息中心 | 海上溢油漂移路径及扩散快速预报方法 |
CN105160681A (zh) * | 2015-09-09 | 2015-12-16 | 浙江海洋学院 | 一种海面溢油监测方法及装置 |
Non-Patent Citations (1)
Title |
---|
珠江口溢油预测预警与应急决策技术研究;牛志刚;《中国优秀硕士学位论文全文数据库 工程科技I辑》;20130215(第2期);第16-79,103-126页 * |
Also Published As
Publication number | Publication date |
---|---|
CN107944608A (zh) | 2018-04-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107944608B (zh) | 基于卫星遥感的海面漂移物与溢油漂移扩散预报方法 | |
CN105760575B (zh) | 渤海海上溢油输移、扩展数值预报系统的建立方法 | |
CN110399676B (zh) | 西北太平洋三维溢油业务化应急预报及评估系统 | |
CN105159954A (zh) | 基于gis与虚拟现实技术的三维电子海图生成方法 | |
Rao et al. | Estimation of extreme water levels due to cyclonic storms: a case study for Kalpakkam coast | |
Zhurbas et al. | Rotation of floating particles in submesoscale cyclonic and anticyclonic eddies: a model study for the southeastern Baltic Sea | |
CN110910448A (zh) | 一种海上溢油漂移轨迹溯源方法 | |
Ponomarev et al. | Simulation of mesoscale circulation in the Tatar Strait of the Japan Sea | |
CN111626501B (zh) | 一种评估互花米草适生区域的方法 | |
CN112633545A (zh) | 一种基于闽江下游突发水污染事故风险预警方法及系统 | |
García-Sánchez et al. | Very high resolution tools for the monitoring and assessment of environmental hazards in coastal areas | |
Zhurbas et al. | Hindcast of the mesoscale eddy field in the Southeastern Baltic Sea: Model output vs satellite imagery | |
Simecek-Beatty | Oil spill trajectory forecasting uncertainty and emergency response | |
National Research Council et al. | Coastal meteorology: a review of the state of the science | |
Prants et al. | Lagrangian Oil Spill Simulation in Peter the Great Bay (Sea of Japan) with a High-Resolution ROMS Model | |
Elliott et al. | The need for operational forecasting during oil spill response | |
González et al. | Operational Oceanography System applied to the Prestige oil-spillage event | |
Brushett et al. | Evaluation of met-ocean forecast data effectiveness for tracking drifters deployed during operational oil spill response in Australian waters | |
Kamza et al. | Prediction of the flooding area of the northeastern Caspian Sea from satellite images | |
Borrione et al. | Understanding altimetry signals in the Northeastern Ligurian sea using a multi-platform approach | |
Helmi et al. | Intergrated approach of tsunami vulnerability assessment at coastal area of kalianda sub district, south lampung district, lampung Province, Indonesia | |
French-McCay et al. | Evaluation of field-collected data measuring fluorescein dye movements and dispersion for dispersed oil transport modeling | |
Yichuan et al. | Quantifying wind-induced impacts on particulate Cu footprint in the Yangtze Estuary | |
Pereira | Coupled Modelling of the Minho and Lima Estuaries: Hydrological Response to Climate Changes | |
Morozov et al. | Polar frontal zone of the Barents Sea Western Trough based on the direct measurements in 2007 |
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 | ||
CB02 | Change of applicant information |
Address after: 350001 unit 6, 4 building, 245 Xiu Shan Road, Xin Dian Town, Jinan District, Fuzhou, Fujian, 608 Applicant after: Four creation technology limited company Address before: 350012 unit 6, 4 building, 245 Xiu Shan Road, Xin Dian Town, Jinan District, Fuzhou, Fujian, 608 Applicant before: Fujian Strong Software Co., Ltd. |
|
CB02 | Change of applicant information | ||
GR01 | Patent grant | ||
GR01 | Patent grant |