Marine meteorological service benefit evaluation method

Abstract

The invention discloses a marine meteorological service benefit evaluation method, which comprises the following steps of: taking a gust forecast grade sequence of a port forecast manuscript and a port area forecast manuscript for 3 hours by adopting a method of evaluating benefits for 3 hours in a scoring time period; in the gust sequence, the visibility sequence and the significant wave height sequence of the port live index station, a gust live grade sequence which is maximum in 3 hours, a visibility value sequence which is minimum in 3 hours and a significant wave height sequence which is maximum in 3 hours are obtained every 3 hours; port forecast, port area forecast and port index station live data are correspondingly substituted into a benefit evaluation calculation method to obtain benefit hours, the benefit hours are spliced to form a time-continuous benefit duration sequence, the time-continuous benefit duration sequence is compared with a single operation duration, and a final effective benefit duration sequence is selected. And finally, the effective service benefit durations obtained by all the forecast manuscripts are accumulated, so that the meteorological service benefit durations are obtained, and the meteorological service benefit durations are more accordant with a real production scene.

Description

Technical Field

[0001] The invention belongs to the technical field of meteorological service evaluation, and in particular relates to a method for evaluating the benefits of marine meteorological services. Background Art

[0002] Meteorological services play an important role in marine disaster prevention, marine economy, marine rights protection, climate change response and marine ecological environment protection. Severe weather such as strong winds, heavy fog, severe convection and typhoons will directly affect port operations and maritime shipping safety, leading to port congestion, suspension of operations and suspension of routes, causing great loss of life and property to port companies and cargo ships. The port and shipping meteorological service that came into being provides scientific, accurate and professional meteorological services for port operations, port-side business, oil refueling, key offshore engineering construction and maritime transportation, and is an important guarantee for promoting the development of the marine economy.

[0003] Previous researchers have proposed benefit estimation methods for port meteorological services, but they still have the following shortcomings: lack of objective and detailed meteorological data support, only forming estimation results with annual granularity; or only calculating for the single operation area of ​​the port, without forming a benefit evaluation algorithm for various industries in the marine economy; or only considering the impact of wind on port operations during the calculation process, without fully considering many influencing factors such as visibility and wave height, as well as actual factors such as the guarantee time required for a single operation. Summary of the invention

[0004] In order to make up for the shortcomings of the prior art, the purpose of the present invention is to provide a method for evaluating the benefits of marine meteorological services. By integrating the thresholds of meteorological and oceanic conditions such as wind force, visibility, and effective wave height in the production processes of various industries, a meteorological benefit evaluation value method is proposed for different service objects such as port areas, offshore refueling anchorages, passenger and cargo routes, iron ore and petrochemical products and other loading and unloading ports. It is more in line with the actual production scenario, and can quantify the increase in production and income of related marine industries, providing an effective experimental reference for the practical evaluation of the development of marine economic benefits.

[0005] The technical problem to be solved by the present invention can be achieved through the following technical solutions:

[0006] A method for evaluating the benefits of marine meteorological services comprises the following steps:

[0007] S1. Obtain the weather forecast text for the port and the weather forecast text for the port area (or the entire sea area, taking the port area as an example), and parse the forecast information into a forecast time series array;

[0008] S2. According to the gust level forecast time series array in the port forecast, extract the 3-hour gust level forecast subsequence in the scoring time period, and record the sequence as df i0 , dfi1 ...df ij ;According to the gust level forecast time series array in the port forecast draft, the 3-hour gust level forecast subsequence in the scoring time period is intercepted, and the sequence is recorded as DF i0 , DF i1 ...DF ij ;

[0009] S3. Take the actual gust wind speed data of the index station, obtain the actual data sequence of gust wind speed, the actual data sequence of hourly minimum visibility, and the actual data sequence of hourly maximum effective wave height in the scoring time period, and take the maximum value of the actual data sequence of gust wind speed every 3 hours in the scoring time period, convert the wind speed into wind level through the wind force level conversion table, and obtain the actual maximum value sequence of gust level every 3 hours; take the minimum value of the actual data sequence of hourly minimum visibility every 3 hours in the scoring time period, and obtain the actual minimum value sequence of visibility every 3 hours; take the minimum value of the actual data sequence of hourly maximum effective wave height every 3 hours in the scoring time period, and obtain the actual maximum value sequence of effective wave height every 3 hours;

[0010] S4. Through field visits and investigations, record the industry operation stop conditions of each port, and set the wind power stop operation condition of the port as Wind stop , visibility stop working condition is Vis stop , the wave height stop operation condition is Wave stop , the wind operation threshold of the port is Wind threshold , evaluate j data one by one, and obtain j benefit time series r i0 , r i1 ...r ij ;

[0011] S5. Find the splicing time interval in the benefit time sequence to form the benefit time sequence res, and accumulate the elements in the benefit time sequence res, that is,

[0012] S6. Calculate the arithmetic accumulation of the benefit hours to obtain the benefit duration of the meteorological service, that is:

[0013]

[0014] Where N is the difference in days between the start and end dates.

[0015] Furthermore, in step S1, the port forecast draft production time is set to tm i , the port forecast was produced at TM i, the start time of the scoring period is tr0, and the end time of the scoring period is tr1; the time text of the forecast has a clear time definition and is intervals of 3 hours, the forecast release time is used as the initial time, and the meteorological text forecast content is parsed into a 3-hour forecast time series array according to the meteorological time terms. On the basis of the real-time 1-hour data of the automatic weather station, in order to match the forecast time series array, it is also processed into a real-time time series array of elements such as the real-time maximum wind speed every 3 hours, the real-time minimum visibility every 3 hours, and the real-time maximum effective wave height every 3 hours.

[0016] Furthermore, the specific operation in step S2 is: the gust level forecast time series array in the port forecast draft is generated according to the port forecast draft production time tm i The number of 3-hour gust level forecast subsequences from the start time tr0 of the scoring period to the end time tr1 of the scoring period is j = (tr1-tr0) / 3, and the sequence is recorded as df i0 , df i1 ...df ij ; The gust level forecast time series array in the port forecast is converted into the TM according to the port forecast production time i The number of 3-hour gust level forecast subsequences from the start time tr0 of the scoring period to the end time tr1 of the scoring period is j = (tr1-tr0) / 3, and the sequence is recorded as DF i0 , DF i1 ...DF ij .

[0017] Further, the specific operation in step S3 is: take the actual gust wind speed data of the index station, obtain the hourly gust wind speed actual data sequence, the hourly minimum visibility actual data sequence, and the hourly maximum effective wave height data sequence in the time period from tr0 to tr1, the number of which is tr1-tr0, and for this hourly gust wind speed actual data sequence, take the maximum value every 3 hours in the time period from tr0 to tr1, and convert the wind speed into wind level through the wind force level conversion table to obtain the maximum value sequence of the actual level of gust every 3 hours, the number of which is j=(tr1-tr0) / 3, and the sequence is recorded as dr io , dr i1 ...dr ij ; For the hourly minimum visibility live data sequence, take the minimum value every 3 hours in the time period from tr0 to tr1, the number is j = (tr1-tr0) / 3, and the sequence is recorded as Vis io , Vis i1 ...Vis ij; For the hourly maximum effective wave height data sequence, the maximum value is taken every 3 hours in the time period from tr0 to tr1, the number is j = (tr1-tr0) / 3, and the sequence is recorded as Wave io , Wave i1 ...Wave ij .

[0018] Furthermore, in step S4, the j data are evaluated one by one. First, the job stop condition is determined for the data at each time node. Taking the data at time node j as an example, if df ij >Wind stop or dr ij >Wind stop or Vis ij <Vis stop or Wave ij >Wave stop , then the benefit value r ij If df is 0, no benefit will be generated; otherwise, the judgment under normal operating conditions is entered. ij ≤Wind threshold 、dr ij ≤Wind threshold , while DF ij >Wind threshold , then the benefit value r ij is 3, otherwise r ij is 0; thus obtaining j benefit time series r i0 , r i1 ...r ij .

[0019] Furthermore, in step S5, r i1 , r i2 ...r ij The element value in the benefit time sequence is 3 or 0. Connect the adjacent time intervals with a value of 3 to form a benefit time sequence re i0 ...re ik , 0≤k<j, re i0 ...re ik Compare them one by one with the single operation time mt. If re ik ≥mt, then it is retained, and finally the benefit duration sequence res is formed i0 ...res ix , 0≤x<(k / mt).

[0020] Compared with the prior art, the present invention has the following advantages:

[0021] (1) The marine meteorological service benefit evaluation method proposed in the present invention comprehensively considers the operation stop conditions (wind, visibility, effective wave height) in the production process of each industry, as well as the industry operation wind threshold in the operation area, the time required for a single operation and other conditions. The benefit duration is obtained by comparing the professional weather forecast of the operation area with the weather forecast of the port area where it is located, or the forecast draft of a larger spatial area (such as the public forecast draft is for the sea area within the jurisdiction of the city). The operation stop conditions, industry operation wind threshold and the time required for a single operation are obtained through field visits and face-to-face surveys with offshore operation users. The designed method is more in line with the actual operation scenario. The objective calculation method is clear, and the time granularity is fine to the day. The marine meteorological service benefits within any date range can be calculated. Therefore, whether it is the service benefits of a certain major weather process or the overall service benefits within a long period of time, it can be scientifically statistically obtained by this method, which has good social and economic benefits.

[0022] (2) The marine meteorological service benefit evaluation method proposed in the present invention has a wide coverage, covering all industries of the marine economy, and is applicable to ports, offshore refueling anchorages, passenger and cargo routes and port areas. A business system has been developed and the benefit evaluation method has been integrated into the business platform for daily customized professional meteorological business services. BRIEF DESCRIPTION OF THE DRAWINGS

[0023] Figure 1 is a flow chart of the method of the present invention;

[0024] Figure 2 This is a flow chart of the wind force level processing of the daily forecast draft of the present invention;

[0025] Figure 3 This is a flowchart of processing the daily small area (port) real-time data (wind force, visibility, and effective wave height) of the present invention;

[0026] Figure 4 The figure is a flow chart of the calculation of daily marine meteorological benefit evaluation of the present invention. DETAILED DESCRIPTION

[0027] In order to enable those skilled in the art to more clearly understand the technical solution of the present application, a method for evaluating the benefits of marine meteorological services of the present invention is further described below in conjunction with the accompanying drawings and specific embodiments.

[0028] like Figure 1 As shown in FIG. 1 , a method for evaluating the benefits of marine meteorological services includes the following specific steps:

[0029] (1) Obtain the port weather forecast text draft (public weather forecast, port area weather forecast, anchorage weather forecast, route weather forecast, port weather forecast) and the port area weather forecast text draft information, and parse the forecast information into a forecast time series array.

[0030] Take the comparison between the port forecast (small area) and the extreme wind forecast (large area) of the port area as an example. Assuming that the port forecast was made on a certain day at tm i , the port forecast was produced at TM i , the start time of the scoring period is tr0, and the end time of the scoring period is tr1.

[0031] The time text in the forecast has a clear time definition and is intervals of 3 hours, such as "morning" is 05-08, "am" refers to 08-11... Therefore, for each specific meteorological text forecast, the forecast release time can be used as the initial time, and the content of the meteorological text forecast can be parsed into a forecast time series array according to the meteorological time terms. On the basis of the real-time 1-hour data of the automatic weather station, in order to match the forecast time series array, it is also processed into a real-time time series array of elements such as the real-time maximum wind speed every 3 hours, the real-time minimum visibility every 3 hours, and the real-time maximum effective wave height every 3 hours.

[0032] (ii) Based on the gust level forecast time series array in the port forecast, extract the 3-hour gust level forecast subsequence in the scoring period; Based on the gust level forecast time series array in the port area forecast, extract the 3-hour gust level forecast subsequence in the scoring period.

[0033] The gust level forecast sequence in the port forecast draft is converted according to tm i The number of 3-hour gust level forecast subsequences in the time period from tr0 to tr1 is j = (tr1-tr0) / 3, and the sequence is recorded as df i0 , df i1 ...df ij . The gust level forecast sequence in the port forecast draft is converted according to TM i The number of 3-hour gust level forecast subsequences in the time period from tr0 to tr1 is j = (tr1-tr0) / 3, and the sequence is recorded as DF i0 , DF i1 ...DF ij .

[0034] (3) Take the actual gust wind speed data of the index station, obtain the actual data sequence of gust wind speed, the actual data sequence of hourly minimum visibility, and the actual data sequence of hourly maximum effective wave height in the scoring time period, and take the maximum value of the actual data sequence of gust wind speed according to the time period. After converting the wind speed into wind level through the wind force level conversion table, obtain the actual maximum value sequence of hourly gust level, the actual data sequence of hourly minimum visibility, and the actual data sequence of hourly maximum effective wave height.

[0035] Take the actual gust wind speed data of the index station, and obtain the actual data sequence of gust wind speed in the time period from tr0 to tr1, the actual data sequence of hourly minimum visibility, and the actual data sequence of hourly maximum effective wave height, the number of which is tr1-tr0. For this actual data sequence of gust wind speed, take the maximum value every 3 hours, and convert the wind speed into wind level through the wind force level conversion table (such as Table 1), and obtain the maximum value sequence of the actual level of gust wind every 3 hours, the number of which is j=(tr1-tr0) / 3, and the sequence is recorded as dr io , dr i1 ...dr ij ; Take the minimum value of the minimum visibility real-time data sequence every 3 hours, the number is j = (tr1-tr0) / 3, and the sequence is recorded as Vis io , Vis i1 ...Vis ij ; Take the maximum value of the hourly maximum effective wave height data sequence every 3 hours, the number is j = (tr1-tr0) / 3, and the sequence is recorded as Wave io , Wave i1 ...Wave ij .

[0036] Table 1 Wind level and wind speed comparison table

[0037] Wind level Wind speed (m / s, meters / second) 0 0～0.2 1 0.3～1.5 2 1.6～3.3 3 3.4～5.4 4 5.5～7.9 5 8.0～10.7 6 10.8～13.8 7 13.9～17.1 8 17.2～20.7 9 20.8～24.4 10 24.5～28.4 11 28.5～32.6 12 32.7～36.9 13 37.0～41.4 14 41.5～46.1 15 46.2～50.9 16 51.0～56.0 17 56.1～61.2

[0038] (IV) Through field visits and investigations, record the industry operation stop conditions of each port. For example, the wind power stop operation condition of a port is Wind stop (such as 10 levels), visibility stop working condition is Vis stop (e.g. 2000m), the wave height stop operation condition is Wave stop (such as 3 meters), the wind operation threshold of the port is Wind threshold (such as level 7), evaluate j data one by one, and thus obtain j benefit time series r i0 , r i1 ...r ij .

[0039] Evaluate j data one by one. First, determine the job stop condition for the data at each time node. Take the data at time node j as an example. If df ij >Wind stop or dr ij >Wind stop or Vis ij <Vis stop or Wave ij >Wave stop , then the benefit value r ijIf df is 0, no benefit will be generated; otherwise, the judgment under normal operating conditions is entered. ij ≤Wind threshold 、dr ij ≤Wind threshold , while DF ij >Wind threshold , then the benefit value r ij is 3, otherwise r ij is 0; thus obtaining j benefit time series r i0 , r i1 ...r ij .

[0040] (V) Find the splicing time interval in the benefit time sequence to form a benefit time sequence, where r i1 , r i2 ...r ij The element value in the benefit time sequence is 3 or 0. Connect the adjacent moments with a value of 3 to form a benefit time sequence re i0 ...re ik , 0≤k<j, re i0 ...re ik Compare them one by one with the single operation time mt. If re ik ≥mt, then it is retained, and finally the benefit duration sequence res is formed i0 ...res ix , 0≤x<(k / mt). Accumulate the elements in the benefit duration sequence res, that is,

[0041] (VI) The arithmetical accumulation of the benefit hours is used to obtain the benefit duration of meteorological services, namely:

[0042]

[0043] Where N is the difference in days between the start and end dates.

[0044] In summary, the present invention adopts a method of benefit evaluation for 3 hours within the scoring period, which means that the wind forecast value in a small area is ≤ the operating threshold, while the wind forecast value in a large area is > the operating threshold, and the operating benefit hours generated therefrom. Taking the 2023 port forecast as an example, the conditions are the port forecast draft made at 05:00 in the morning of the same day and the port area forecast draft made at 05:00 in the morning of the same day, and the benefit evaluation is conducted for the single operation time of 6 hours in the time period from 08:00 on the same day to 08:00 on the next day under the operating condition of 7-force wind. For a certain day in 2023, take the gust forecast level sequence within the range of 08:00 on the same day and 08:00 on the next day in the port forecast draft and port area forecast draft made at 05:00 on the morning of the same day (the number of forecasts is 24 / 3=8), and then obtain the maximum gust actual level sequence every 3 hours, the minimum visibility value sequence every 3 hours and the maximum visibility value sequence every 3 hours from the gust sequence, visibility sequence and effective wave height sequence of the port actual index station within the range of 08:00 on the same day and 08:00 on the next day in a 3-hour section. The effective wave height sequence (there are 24 actual data, and the maximum value is taken in 3 hours to obtain 24 / 3=8), and the 8 port forecasts, port area forecasts and actual data (maximum gust, minimum visibility, maximum effective wave height) are correspondingly substituted into the benefit evaluation calculation method, and 8 benefit hours (3 hours or 0 hours) can be obtained. The benefit time is spliced ​​to form a time-continuous benefit duration sequence, and then compared with the single operation duration of 6 hours, and the final effective benefit duration sequence is selected. Columns are accumulated as the effective service benefit duration of the single-day forecast (for example, 08:00-11:00 is effective, 11:00-14:00 is effective, 14:00-17:00 is effective, 17:00-20:00 is invalid, 20:00-23:00 is effective, 23:00 is effective-02:00 the next day, 02:00 the next day-05:00 the next day is invalid, 05:00 the next day-08:00 the next day is valid), first splice to form 3 benefit duration sequences, namely 08:00-17:00 (duration 9 hours), 20:00-02:00 the next day (duration 6 hours), It is valid from 05:00 to 08:00 the next day (3 hours), so it is compared with the single operation duration of 6 hours. The effective benefit duration is greater than or equal to 6 hours, that is, the effective benefit duration sequence is 08:00-17:00 (9 hours) and 20:00-02:00 the next day (6 hours), and the effective service benefit duration of a single-day forecast is 9+6=15 hours). Finally, the effective service benefit duration obtained from all forecast drafts in 2023 is accumulated to obtain the meteorological service benefit duration in 2023.

[0045] Example 1

[0046] like Figure 2-Figure 4 As shown, taking the small area (port) weather forecast versus the large area (port area) weather forecast as an example, this is a method for calculating a single-day benefit evaluation.

[0047] like Figure 2The figure shows the daily forecast wind force level processing flow: the first part is the port area forecast wind force level processing, which includes 4 steps A, B, C, and D; the second part is the port forecast wind force level processing, which includes 4 steps E, F, G, and H;

[0048] A: Get the production time of the corresponding port forecast: the date is known (such as September 2, 2024). At the same time, since the forecast production time of the meteorological department is basically fixed, it is also known which port forecast is based on, such as the morning of the day. Then the production date of the port forecast is: September 2, 2024 + 0 (the day), and the production time is: 5:00 (morning), that is, 5:00 on September 2, 2024. The terms such as the morning of the day are matched one by one through the configuration file.

[0049] B: Search the port forecast database and obtain the text content of the wind forecast: Search the forecast database by the forecast production time and the port name, obtain the forecast content (including weather forecast and wind forecast, etc.), and use string clipping technology to obtain the wind forecast content.

[0050] C: Analyze the wind forecast text and process it into a maximum wind level forecast sequence with 3-hour time matching: Given the forecast draft production time, such as 5:00 on September 2, 2024, the time text of the meteorological text forecast draft has a clear time definition and is at intervals of 3 hours, such as "morning" is 05-08, "morning" refers to 08-11... Therefore, for each specific meteorological text forecast draft, the forecast draft release time can be used as the initial time, and the wind forecast text content can be parsed into a gust forecast time series array according to meteorological time terms. For example, "This morning, the northeast wind is 5-6, and it will turn to southeast wind 5-6 in the morning, with gusts of 7", this paragraph of text can be parsed into a time array of [5:00 on September 2, 2024, 8:00 on September 2, 2024], a maximum wind array of [5.5, 7], and a wind direction of [northeast wind, southeast wind].

[0051] D: According to the scoring period, extract the 3-hour maximum wind level forecast subsequence within the time interval from the maximum wind level forecast sequence: the start time of the scoring period is tr0, and the end time of the scoring period is tr1. If tr0 is 08:00 on the current day, and if tr1 is 08:00 on the next day, and the date is known (such as September 2, 2024), then tr0 is 08:00 on September 2, 2024, and tr1 is 08:00 on September 3, 2024; and the known time array is [5:00 on September 2, 2024, 8:00 on September 2, 2024, 11:00 on September 2, 2024, 14:00 on September 2, 2024, 17:00 on September 2, 2024, 20:00, September 2, 2024 23:00, September 3, 2024 2:00, September 3, 2024 5:00, September 3, 2024 8:00...], the maximum wind array is [5.5, 7, 8, 8, 8, 7, 7, 8, 9, 7...]; then according to the interception condition [tr0, tr1), [8:00 on September 2, 2024, September 2, 11:00, September 2, 2024, 14:00, September 2, 2024, 17:00, September 2, 2024, September 2, 20:00, September 2, 2024 23:00, September 3, 2024 2:00, September 3, 2024, 5:00 on September 3, 2024] and [7, 8, 8, 8, 7, 7, 8, 9] can be intercepted.

[0052] E: Similar to step A, change “port area forecast” to “port forecast”.

[0053] F: Similar to step B, change “port area forecast” to “port forecast”.

[0054] G: Similar to step C, change “port area forecast” to “port forecast”.

[0055] H: Similar to step D, change “port area forecast” to “port forecast”.

[0056] like Figure 3 As shown in the figure, it is the daily small area (port) real-time data (wind force, visibility, and effective wave height) processing flow: the first part is the port real-time wind force level processing, which includes 3 steps I, J, and K; the second part is the port real-time visibility, which includes 2 steps L and M; the third part is the port real-time effective wave height processing, which includes 2 steps N and O;

[0057] I: Search the real-time database, and obtain the real-time sequence of maximum wind speed matched by 1 hour according to the port index station and the scoring period: The real-time database of the meteorological system stores the maximum wind speed data with a granularity of 1 hour. Each port is matched with a corresponding wind index station. According to the port wind index station and the start and end time of the scoring period are tr0 and tr1 (such as tr0 is 08:00 on September 2, 2024, and tr1 is 08:00 on September 3, 2024), 24 maximum wind speed data of 1 hour can be obtained according to the query condition (tr0, tr1] to form a sequence.

[0058] J: Group the hourly maximum wind speed sequence in the scoring period into 3-hour groups in chronological order, and take the maximum value of each 3-hour wind speed group data: the start time of the scoring period is tr0, and the end time of the scoring period is tr1. If tr0 is 08:00 on September 2, 2024, and tr1 is 08:00 on September 3, 2024, then the time grouping sequence is [11:00 on September 2, 2024, 14:00 on September 2, 2024, 17:00 on September 2, 2024, 20:00 on September 2, 2024, 23:00 on September 2, 2024] , September 3, 2024 2:00, September 3, 2024 5:00, September 3, 2024 5:00], 24 1-hour maximum wind speed data can also be cut into 3-hour groups of data accordingly, such as (September 2, 2024 8:00, September 2, 2024 11:00] and match the first time element in the time grouping sequence, and take the maximum value of the 3 maximum wind speed data in the time interval of (September 2, 2024 8:00, September 2, 2024 11:00], so 8 wind speed data will eventually be formed, and match them one by one with the time grouping sequence.

[0059] K: Use the meteorological wind speed to wind scale table to convert the maximum wind speed into wind scale, and finally form a 3-hour maximum wind scale sequence R_list within the scoring period, with count as the number: Refer to Table 1, convert 8 wind speed data into 8 wind scale data, and match them one by one with the time grouping sequence.

[0060] L: The internal database of the meteorological system has the 1-hour minimum visibility element. The search method is similar to step I, but change "port wind index station" to "port visibility index station". Find the hourly minimum visibility sequence.

[0061] M: Similar to step J, except that the hourly visibility sequence is grouped into 3-hour groups and the minimum value is taken. The sequence found is the 3-hour minimum visibility sequence.

[0062] N: The internal database of the meteorological system contains the element of the maximum effective wave height in one hour. The search method is similar to step I, but the "port wind index station" is changed to "port wave index station". Find the maximum effective wave height sequence by hour.

[0063] O: Similar to step J, and the maximum value of the hourly maximum effective wave height sequence is taken when it is grouped by 3 hours, and the sequence found is the 3-hour maximum effective wave height sequence.

[0064] like Figure 4 As shown in the figure, the specific steps of daily marine meteorological benefit evaluation calculation are as follows: the first part is to cyclically search for benefit hour sequence; the second part is to splice benefit hour sequence to form effective benefit hour sequence, i.e., step P; the third part is to perform arithmetic accumulation on effective benefit hour sequence. Method of splicing benefit hour sequence to form effective benefit hour sequence:

[0065] P: Find the splicing time interval in the r sequence, splice the time adjacent elements with a value of 3, and finally form a new re sequence: that is, if the wind forecasts at 11:00 on September 2, 2024, 14:00 on September 2, 2024, and 17:00 on September 2, 2024 are all effective, splice the time into 11:00 on September 2, 2024-20:00 on September 2, 2024, with a duration of 9 hours. The specific calculation steps are as follows:

[0066] Calculation step P1: Select a sequence with a duration of 3 from the r sequence, called the original sequence. The elements in the sequence are objects, which contain three attributes: start time, end time, and duration.

[0067] Calculation step P2: If there is only one object in the original sequence, return this object.

[0068] Calculation step P3: If the number of elements in the sequence is ≥ 2, first store the first object of the original sequence into the result sequence.

[0069] Calculation step P4: When iterating to the second or more items, it is necessary to compare the current object in the original sequence with each element in the result sequence. If there is an element in the result sequence whose end time is exactly equal to the start time of the current object, then extract this element from the result sequence, reassign the end time, recalculate the duration, and replace the old element in the result sequence; if there is no such situation in the result sequence, store the current object in the original sequence in the result sequence.

[0070] After the above steps P1 to P4, the temporal adjacent elements with a value of 3 can be temporally spliced.

[0071] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention, rather than to limit it. Although the present invention has been described in detail with reference to the aforementioned embodiments, those skilled in the art should understand that they can still modify the technical solutions described in the aforementioned embodiments, or replace some or all of the technical features therein with equivalents. However, these modifications or replacements do not cause the essence of the corresponding technical solutions to deviate from the scope of the technical solutions of the embodiments of the present invention.

Claims

1. A method for evaluating the benefits of marine meteorological services, characterized in that: The following steps are involved: S1. Obtain the port weather forecast text and the port area weather forecast text information, and parse the forecast information into a forecast time series array; S2. According to the gust level forecast time series array in the port forecast, extract the 3-hour gust level forecast subsequence in the scoring time period, and record the sequence as df i0 , df i1 ...df ij ;According to the gust level forecast time series array in the port forecast draft, the 3-hour gust level forecast subsequence in the scoring time period is intercepted, and the sequence is recorded as DF i0 , DF i1 ...DF ij ; S3. Take the actual gust wind speed data of the port index station, obtain the hourly gust wind speed actual data sequence, the hourly minimum visibility actual data sequence, and the hourly maximum effective wave height data sequence in the scoring time period, and take the maximum value of the gust wind speed actual data sequence every 3 hours in the scoring time period, convert the wind speed into wind level through the wind force level conversion table, and obtain the gust actual level maximum value sequence every 3 hours; take the minimum value of the hourly minimum visibility actual data sequence every 3 hours in the scoring time period, and obtain the 3-hour actual minimum visibility value sequence; take the minimum value of the hourly maximum effective wave height data sequence every 3 hours in the scoring time period, and obtain the 3-hour actual maximum value sequence of effective wave height; S4. Through field visits and investigations, record the industry operation stop conditions of each port, and set the wind power stop operation condition of the port as Wind stop , visibility stop working condition is Vis stop , the wave height stop operation condition is Wave stop , the wind operation threshold of the port is Wind threshold , evaluate j data one by one, and obtain j benefit time series r i0 , r i1 ...r ij ; S5. Find the splicing time period in the benefit time sequence to form the benefit time sequence res, and accumulate the elements in the benefit time sequence res, that is, S6. Calculate the arithmetic accumulation of the benefit hours to obtain the benefit duration of the meteorological service, that is: Where N is the difference in days between the start and end dates.

2. A method for evaluating the benefits of marine meteorological services according to claim 1, characterized in that: In step S1, the port forecast draft production time is tm i , the port forecast was produced at TM i , the start time of the scoring period is tr0, and the end time of the scoring period is tr1; the time text of the forecast has a clear time definition and is intervals of 3 hours, the forecast release time is used as the initial time, and the meteorological text forecast content is parsed into a 3-hour forecast level time series array according to the meteorological time terms; based on the real-time 1-hour data of the automatic weather station, in order to match the forecast level time series array, it is also processed into a real-time time series array of elements such as the real-time maximum wind speed every 3 hours, the real-time minimum visibility every 3 hours, and the real-time maximum effective wave height every 3 hours.

3. A method for evaluating the benefits of marine meteorological services according to claim 2, characterized in that: The specific operation in step S2 is: the gust level forecast time series array in the port forecast draft is converted into a time series array according to the port forecast draft production time tm i The number of 3-hour gust level forecast subsequences from the start time tr0 of the scoring period to the end time tr1 of the scoring period is j = (tr1-tr0) / 3, and the sequence is recorded as df i0 , df i1 ...df ij ; The gust level forecast time series array in the port forecast is converted into the TM according to the port forecast production time i The number of 3-hour gust level forecast subsequences from the start time tr0 of the scoring period to the end time tr1 of the scoring period is j = (tr1-tr0) / 3, and the sequence is recorded as DF i0 , DF i1 ...DF ij .

4. A method for evaluating the benefits of marine meteorological services according to claim 3, characterized in that: The specific operation in step S3 is: take the actual gust wind speed data of the index station, obtain the actual gust wind speed data sequence, the hourly minimum visibility actual data sequence, and the hourly maximum effective wave height data sequence in the time period from tr0 to tr1, the number of which is tr1-tr0, and take the maximum value of the actual gust wind speed data sequence every 3 hours in the time period from tr0 to tr1, and convert the wind speed into wind level through the wind force level conversion table to obtain the actual gust level maximum value sequence every 3 hours, the number of which is j=(tr1-tr0) / 3, and the sequence is recorded as dr io , dr i1 ...dr ij ; For the minimum visibility live data sequence, take the minimum value every 3 hours in the time period from tr0 to tr1, the number is j = (tr1-tr0) / 3, and the sequence is recorded as Vis io , Vis i1 ...Vis ij ; For the hourly maximum effective wave height data sequence, the maximum value is taken every 3 hours in the time period from tr0 to tr1, the number is j = (tr1-tr0) / 3, and the sequence is recorded as Wave io , Wave i1 ...Wave ij .

5. A method for evaluating the benefits of marine meteorological services according to claim 4, characterized in that: In step S4, the j data are evaluated one by one. First, the job stop condition is determined for the data at each time node. Taking the data at time node j as an example, if df ij >Wind stop or dr ij >Wind stop or Vis ij <Vis stop or Wave ij >Wave stop , then the benefit value r ij If df is 0, no benefit will be generated; otherwise, the judgment under normal operating conditions is entered. ij ≤Wind threshold 、dr ij ≤Wind threshold , while DF ij >Wind threshold , then the benefit value r ij is 3, otherwise r ij is 0; thus obtaining j benefit time series r i0 , r i1 ...r ij .

6. A method for evaluating the benefits of marine meteorological services according to claim 5, characterized in that: In step S5, r i1 , r i2 ...r ij The element value in the benefit time sequence is 3 or 0. Connect the adjacent time intervals with a value of 3 to form a benefit time sequence re i0 ...re ik , 0≤k<j, re i0 ...re ik Compare them one by one with the single operation time mt. If re ik ≥mt, then it is retained, and finally the benefit duration sequence res is formed i0 ...res ix , 0≤x<(k / mt).

Images