Method and device for determining emergency plan area

Abstract

The invention discloses a method and a device for determining an emergency planning area, and relates to the technical field of nuclear industry. The method comprises the following steps: acquiring K * n target meteorological sequences; for the first type of accidents, determining K * n groups of first dose result sets; determining a second dose result set from the K * n groups of first dose result sets; aiming at each target distance in D target distances of the second dose result set, determining a third dose result set of C doses corresponding to the target distance to obtain C * D groups of third dose result sets; and according to the C * D groups of third dose result sets and a preset target dose result set, determining an emergency plan area suitable for the first type of accidents. According to the embodiment of the invention, the method can reduce the calculation amount, reduce the consumed calculation time and reduce the storage space while guaranteeing the accuracy of the determination of the emergency plan area.

Description

Technical Field

[0001] The present application belongs to the field of nuclear industry technology, and specifically relates to a method and device for determining an emergency planning area. Background Art

[0002] The emergency response plan of a nuclear power plant is the last link in the nuclear safety defense-in-depth strategy. The determination of the emergency plan area is an important part of the emergency response plan and is of great significance to the safety of the public around the nuclear power plant. "Emergency Preparedness and Emergency Response of Nuclear Power Plant Operating Units" (HAD002 / 01-2019) stipulates that the operating unit should clarify the scope of the emergency plan area to be established in its emergency plan based on the analysis of possible nuclear accidents and their consequences.

[0003] Since the dose results under different meteorological conditions are different, in the prior art, the dose results after an accident under different meteorological conditions are usually calculated using meteorological data of one or several years to determine the emergency planning area. However, the above method is computationally intensive when determining the emergency planning area, and consumes a large amount of computing time and storage space. Summary of the invention

[0004] The technical problem to be solved by the present application is to provide a method and device for determining an emergency plan area in view of the above-mentioned deficiencies in the prior art. The method can reduce the amount of calculation, the computing time and the storage space while ensuring the accuracy of the emergency plan area determination.

[0005] In a first aspect, an embodiment of the present application provides a method for determining an emergency planning area, comprising:

[0006] S11. Obtain n meteorological sequences corresponding to K target types respectively as target meteorological sequences, and obtain K×n target meteorological sequences, wherein the meteorological sequences include values ​​of V meteorological variables at T hours at the meteorological stations corresponding to the target nuclear power plant sites, and n, K, V and T are all positive integers;

[0007] S12. For the first type of accident, determine the first dose result set of C doses at D target distances and P target directions corresponding to the release of the source term of the first type of accident under each target meteorological sequence, so as to obtain K×n sets of first dose result sets, the first type of accident is any one in the accident set, and D, P, and C are all positive integers;

[0008] S13, determining a second dose result set of C doses at D target distances and P target positions from the K×n groups of first dose result sets;

[0009] S14. For each of the D target distances in the second dose result set, according to the dose result sets of C doses at P target orientations corresponding to the target distance in the second dose result set, determine a third dose result set of C doses corresponding to the target distance from the second dose result set to obtain a C×D group of third dose result sets;

[0010] S15. Determine an emergency plan area applicable to the first type of accident based on the third dose result set of group C×D and the preset target dose result set.

[0011] In some implementations of the first aspect, before obtaining n meteorological sequences corresponding to the K target types as target meteorological sequences and obtaining K×n target meteorological sequences, the method further includes:

[0012] Obtain Q meteorological sequences of the target nuclear power plant site, where Q is a positive integer;

[0013] The Q meteorological sequences are classified to obtain K target types, each target type includes multiple meteorological sequences, and K is less than Q; wherein the Q meteorological sequences are meteorological sequences for each hour throughout the year.

[0014] In some implementations of the first aspect, Q meteorological sequences are classified to obtain K target types, including:

[0015] Select K meteorological sequence features from the Q meteorological sequence features that correspond one to one to the Q meteorological sequences as K initial classification centers;

[0016] For the remaining QK meteorological sequence features among the Q meteorological sequence features, respectively determine the first similarity between the meteorological sequence features and each initial classification center;

[0017] According to the first similarity, the QK meteorological sequences are classified into K initial classification centers respectively to obtain K target types.

[0018] In some implementations of the first aspect, for the remaining QK meteorological sequence features among the Q meteorological sequence features, determining a first similarity between the meteorological sequence features and each initial classification center respectively includes:

[0019] For the remaining QK meteorological sequence features among the Q meteorological sequence features, the first similarity between the meteorological sequence features and each initial classification center is determined respectively according to the target similarity algorithm;

[0020] The target similarity algorithm includes any one of Euclidean distance, cosine similarity, edit distance and Manhattan distance.

[0021] In some embodiments of the first aspect, according to the first similarity, the Q - K meteorological sequences are respectively classified into K initial classification centers to obtain K target types, including:

[0022] For the Q - K meteorological sequences, the meteorological sequences are classified into the target types where the initial classification centers corresponding to their maximum first similarities are located to obtain K target types;

[0023] Alternatively, according to the first similarity, the Q - K meteorological sequences are respectively classified into K initial classification centers to obtain K target types, including:

[0024] S21. According to the first similarity, the Q - K meteorological sequences are respectively classified into K initial classification centers to obtain K initial types;

[0025] S22. For each initial type among the K initial types, the average value of the multiple meteorological sequence features corresponding to the initial type is used as the initial classification center to obtain K initial classification centers;

[0026] S23. On the basis of iteration, steps S21 and S22 are repeatedly executed, and the number of executions is counted until the number of executions is greater than or equal to the preset number of executions. The K initial types finally obtained after iteration are determined as K target types.

[0027] In some embodiments of the first aspect, determining a first dose result set of C doses at D target distances and P target azimuths after the first - type accident release source term under each target meteorological sequence, including:

[0028] Determining an initial dose result set of C doses at Y distances and W azimuths after the first - type accident release source term under each target meteorological sequence, where 0 < D ≤ Y, 0 < P ≤ W, and both Y and W are integers;

[0029] Determining D target distances related to the emergency planning area from the Y distances, and determining P target azimuths related to the emergency planning area from the W azimuths;

[0030] Determining a first dose result set of C doses at D target distances and P target azimuths from the initial dose result set.

[0031] In some embodiments of the first aspect, determining a second dose result set of C doses at D target distances and P target azimuths from the K×n groups of first dose result sets, including:

[0032] For each target meteorological sequence, obtaining the first weight value corresponding to the target meteorological sequence;

[0033] For each target meteorological sequence, multiplying the first weight value by the first dose result set corresponding to the target meteorological sequence as a fourth dose result set to obtain K×n sets of fourth dose result sets;

[0034] Sorting the fourth dose result set of K×n groups to obtain a first sorting;

[0035] The fourth dose result set located at a preset percentile in the first sorting is used as the second dose result set.

[0036] In some implementations of the first aspect, the first weight values ​​corresponding to the target meteorological sequences under the same target type are the same;

[0037] Obtaining a first weight value corresponding to a target meteorological sequence includes:

[0038] Acquire a first quantity and a total quantity of meteorological sequences of the target nuclear power plant site, wherein the first quantity is the quantity of meteorological sequences included in the target type corresponding to the target meteorological sequence;

[0039] Substitute the first quantity and the total quantity into formula (1) to calculate and obtain the first weight value corresponding to the target meteorological sequence;

[0040] Formula (1) includes:

[0041]

[0042] Wherein, M represents the first weight value; N k represents a first quantity; N k Indicates the total quantity.

[0043] In some embodiments of the first aspect, according to a dose result set of C doses at P target orientations corresponding to the target distances in the second dose result set, determining a third dose result set of C doses corresponding to the target distances from the second dose result, comprises:

[0044] Determine from the second dose result set a fifth dose result set of C doses corresponding to each target orientation at the target distance;

[0045] The fifth dose result set with the largest value is determined as the third dose result set of C doses corresponding to the target distance.

[0046] In some embodiments of the first aspect, determining an emergency plan area applicable to a first type of accident according to C×D third dose result sets and a preset target dose result set includes:

[0047] For each target distance, according to the third dose result set of C doses corresponding to the target distance and the preset target dose result set, C initial radius ranges corresponding to the C doses are determined to obtain C×D initial radius ranges;

[0048] For each target distance, determine a first radius range from C initial radius ranges corresponding to the target distance to obtain D first radius ranges;

[0049] For each target distance, the intersection of D first radius ranges is used as the target radius range of the emergency planning area for the first type of accident;

[0050] The area formed by taking the center of the reactor in the target nuclear power plant as the center of the circle and the maximum target distance within the target radius as the radius is determined as the emergency planning area suitable for the first type of accidents.

[0051] In some embodiments of the first aspect, determining C initial radius ranges corresponding to the C doses one by one according to a third dose result set of C doses corresponding to the target distance and a preset target dose result set includes:

[0052] For the first dose result corresponding to each dose in the third dose result set, when the first dose result is less than or equal to the target dose result corresponding to the dose in the preset target dose result set, the initial radius range is determined to be greater than 0 and less than or equal to the target distance corresponding to the first dose result, or, when the first dose result is greater than the target dose result, the initial radius range is determined to be greater than the target distance corresponding to the first dose result, so as to determine C initial radius ranges corresponding to the C doses one by one.

[0053] In some implementations of the first aspect, determining the first radius range from C initial radius ranges corresponding to the target distance includes:

[0054] When the C initial radius ranges corresponding to the target distances are all the same, the initial radius range is determined as the first radius range;

[0055] Alternatively, when there are at least two different initial radius ranges among the C initial radius ranges corresponding to the target distance, the second radius range among the C initial radius ranges is determined as the first radius range;

[0056] Among them, any value within the second radius range is greater than or equal to all values ​​within the other C-1 initial radius ranges except the second radius range.

[0057] In some embodiments of the first aspect, after determining the emergency plan area applicable to the first type of accident according to C×D third dose result sets and the preset target dose result set, the method further includes:

[0058] S16. All accidents in the accident set except the first category accidents are regarded as the first category accidents;

[0059] S17, repeating steps S12 to S16 until an emergency plan area corresponding to each accident in the accident set is obtained;

[0060] S18. Determine the largest emergency planning area among all emergency planning areas as the emergency planning area applicable to the accident collection.

[0061] In some embodiments of the first aspect, the V meteorological variables include at least one of precipitation, stability, wind speed in the U direction, and wind speed in the V direction;

[0062] And / or, type C dose includes at least one of a 2-day effective dose, a 7-day effective dose and a thyroid dose.

[0063] Based on the same inventive concept, in a second aspect, an embodiment of the present application further provides a device for determining an emergency plan area, the device comprising:

[0064] The first acquisition module is used to acquire n meteorological sequences corresponding to K target types as target meteorological sequences, and obtain K×n target meteorological sequences, wherein the meteorological sequences include values ​​of V meteorological variables at T hours at the meteorological stations corresponding to the target nuclear power plant sites, and n, K, V and T are all positive integers;

[0065] The first determination module is connected to the first acquisition module and is used to determine, for the first type of accident, a first dose result set of C doses at D target distances and P target directions corresponding to the release of the source term of the first type of accident under each target meteorological sequence, so as to obtain K×n groups of first dose result sets, the first type of accident is any one of the accident set, and D, P, and C are all positive integers;

[0066] A second determination module, connected to the first determination module, is used to determine a second dose result set of C doses at D target distances and P target orientations from the K×n groups of first dose result sets;

[0067] a third determination module connected to the second determination module, configured to determine, for each of the D target distances in the second dose result set, a third dose result set of C doses corresponding to the target distance from the second dose result set according to the dose result sets of C doses at P target orientations corresponding to the target distance in the second dose result set, so as to obtain a C×D group of third dose result sets;

[0068] The fourth determination module is connected to the third determination module and is used to determine the emergency plan area applicable to the first type of accident according to the third dose result set of the C×D group and the preset target dose result set.

[0069] In some embodiments of the second aspect, the device further comprises:

[0070] The second acquisition module is used to acquire Q meteorological sequences of the target nuclear power plant site, wherein the Q meteorological sequences are meteorological sequences for each hour of the whole year; Q is a positive integer;

[0071] The first classification module is connected to the second acquisition module and is used to classify Q meteorological sequences to obtain K target types, each target type includes multiple meteorological sequences, and K is less than Q.

[0072] In some embodiments of the second aspect, the first classification module comprises:

[0073] A selection submodule, connected to the second acquisition module, is used to select K meteorological sequence features from the Q meteorological sequence features corresponding to the Q meteorological sequences as K initial classification centers;

[0074] A first determination submodule is connected to the selection submodule and is used to determine the first similarity between the meteorological sequence features and each initial classification center for the remaining QK meteorological sequence features among the Q meteorological sequence features;

[0075] The first classification submodule is connected to the selection submodule and is used to classify the QK meteorological sequences into K initial classification centers according to the first similarity to obtain K target types.

[0076] According to the method and device for determining the emergency plan area provided in the embodiment of the present application, by obtaining n meteorological sequences corresponding to K target types as target meteorological sequences, K×n target meteorological sequences are obtained. On the one hand, since n meteorological sequences are extracted from each target type as target meteorological sequences, the obtained K×n target meteorological sequences can cover various meteorological conditions, thereby ensuring the accuracy of the emergency plan area determination; on the other hand, compared with the prior art of taking all meteorological sequences as target meteorological sequences, the number of target meteorological sequences can be reduced, thereby reducing the amount of calculation, time spent and storage space in the entire emergency plan area determination process.

[0077] In summary, the embodiments of the present application can reduce the amount of calculation, the time spent, and the storage space while ensuring the accuracy of the emergency plan area determination. BRIEF DESCRIPTION OF THE DRAWINGS

[0078] Figure 1 A schematic diagram showing a flow chart of a method for determining an emergency planning area provided in an embodiment of the present application;

[0079] Figure 2A schematic diagram showing the 48-hour changes of various meteorological variables after superposition of all meteorological sequences of the first category provided in an embodiment of the present application;

[0080] Figure 3 A schematic diagram showing the 48-hour changes of various meteorological variables after superposition of all meteorological sequences of the 11th category provided in an embodiment of the present application;

[0081] Figure 4 A schematic diagram showing the 48-hour changes of various meteorological variables after superposition of all meteorological sequences of the 21st category provided in an embodiment of the present application;

[0082] Figure 5 A schematic diagram showing the 48-hour changes of various meteorological variables after superposition of all meteorological sequences of the 31st category provided in an embodiment of the present application;

[0083] Figure 6 Another schematic flow chart of a method for determining an emergency planning area provided in an embodiment of the present application is shown;

[0084] Figure 7 A structural schematic diagram of a device for determining an emergency plan area provided in an embodiment of the present application is shown. DETAILED DESCRIPTION

[0085] In order to enable those skilled in the art to better understand the technical solution of the present application, the present application is further described in detail below with reference to the accompanying drawings and embodiments.

[0086] The features and exemplary embodiments of various aspects of the present application will be described in detail below. In order to make the purpose, technical solutions and advantages of the present application clearer, the present application will be further described in detail below in conjunction with the accompanying drawings and specific embodiments. It should be understood that the specific embodiments described herein are only configured to explain the present application and are not configured to limit the present application. For those skilled in the art, the present application can be implemented without the need for some of these specific details. The following description of the embodiments is only to provide a better understanding of the present application by illustrating the examples of the present application.

[0087] It should be noted that, in this article, relational terms such as first and second, etc. are only used to distinguish one entity or operation from another entity or operation, and do not necessarily require or imply any such actual relationship or order between these entities or operations. Moreover, the terms "include", "comprise" or any other variants thereof are intended to cover non-exclusive inclusion, so that a process, method, article or device including a series of elements includes not only those elements, but also other elements not explicitly listed, or also includes elements inherent to such process, method, article or device. In the absence of further restrictions, the elements defined by the statement "include..." do not exclude the presence of other identical elements in the process, method, article or device including the elements.

[0088] It should be understood that the term "and / or" used in this article is only a description of the association relationship of associated objects, indicating that there can be three relationships. For example, A and / or B can represent: A exists alone, A and B exist at the same time, and B exists alone. In addition, the character " / " in this article generally indicates that the associated objects before and after are in an "or" relationship.

[0089] Example 1

[0090] First, a method for determining an emergency planning area provided in an embodiment of the present application is described with reference to the accompanying drawings.

[0091] Figure 1 A schematic flow chart of a method for determining an emergency plan area provided in an embodiment of the present application is shown.

[0092] The execution subject of the method for determining the emergency plan area provided in the embodiment of the present application may be an electronic device and a device for determining the emergency plan area, etc., that is, the method for determining the emergency plan area provided in the embodiment of the present application may be executed by an electronic device and a device for determining the emergency plan area, etc. The following description is made by taking the execution subject as an electronic device as an example.

[0093] like Figure 1 As shown, the method for determining the emergency plan area provided in the embodiment of the present application may include steps S11 to S15.

[0094] S11. Obtain n meteorological sequences corresponding to K target types as target meteorological sequences, and obtain K×n target meteorological sequences. The meteorological sequences include the values ​​of V meteorological variables at T hours at the meteorological stations corresponding to the target nuclear power plant sites, and n, K, V and T are all positive integers.

[0095] S12. For the first type of accident, determine the first dose result set of C doses at D target distances and P target directions corresponding to the release of the source term of the first type of accident under each target meteorological sequence, so as to obtain K×n groups of first dose result sets. The first type of accident is any one in the accident set, and D, P, and C are all positive integers.

[0096] S13. Determine a second dose result set of C doses at D target distances and P target positions from the K×n groups of first dose result sets.

[0097] S14. For each of the D target distances in the second dose result set, determine a third dose result set of C doses corresponding to the target distance from the second dose result set according to the dose result set of C doses at P target orientations corresponding to the target distance in the second dose result set, so as to obtain a C×D group of third dose result sets.

[0098] S15. Determine an emergency plan area applicable to the first type of accident based on the third dose result set of group C×D and the preset target dose result set.

[0099] According to the method for determining the emergency planning area provided in an embodiment of the present application, by obtaining n meteorological sequences corresponding to K target types as target meteorological sequences, K×n target meteorological sequences are obtained. On the one hand, since n meteorological sequences are extracted from each target type as target meteorological sequences, the obtained K×n target meteorological sequences can cover various meteorological conditions, thereby ensuring the accuracy of the emergency planning area determination; on the other hand, compared with the prior art of taking all meteorological sequences as target meteorological sequences, the number of target meteorological sequences can be reduced, thereby reducing the amount of calculation, time spent and storage space in the entire emergency planning area determination process.

[0100] The specific implementation methods of the above steps are introduced below.

[0101] In S11, each of the K target types may include at least n meteorological sequences.

[0102] Each meteorological sequence corresponds to a meteorological sequence feature. The meteorological sequence may be a sequence of meteorological element values ​​recorded at certain time intervals. The meteorological sequence feature may be characteristic information obtained after statistical analysis and description of the data in the meteorological sequence.

[0103] Exemplarily, the meteorological sequence features corresponding to the meteorological sequence can be used to characterize the changes in V meteorological variables at the meteorological station corresponding to the target nuclear power plant site in T hours.

[0104] Exemplarily, the meteorological station corresponding to the target nuclear power plant site may be the meteorological station where the target nuclear power plant site is located. The target nuclear power plant site may be any nuclear power plant site.

[0105] Exemplarily, the electronic device stores multiple meteorological sequences corresponding to K target types. For each target type, the electronic device can randomly obtain n meteorological sequences corresponding to it from the multiple meteorological sequences corresponding to it as the target meteorological sequence. Alternatively, for each target type, the electronic device can obtain n meteorological sequences with odd sequence numbers from the multiple meteorological sequences corresponding to it as the target meteorological sequence.

[0106] In some embodiments, the V meteorological variables include at least one of precipitation, stability, wind speed in the U direction, and wind speed in the V direction.

[0107] In this embodiment, by considering at least one meteorological variable among precipitation, stability, wind speed in U direction and wind speed in V direction, the emergency plan area determined according to the target meteorological sequence can effectively respond to various meteorological disasters and emergencies.

[0108] Precipitation is the phenomenon that water in the atmosphere falls to the ground from clouds or the atmosphere in liquid or solid form. Precipitation can be in the form of rain, snow, hail, sleet, and sleet.

[0109] Stability can be a measure of whether the disturbance superimposed on the atmospheric background field can be strengthened over time. It also refers to the degree to which a certain air mass in the air generates buoyancy due to the intensity difference in density, temperature and velocity with the surrounding air, causing it to rise or fall with acceleration. The ability of the atmosphere to suppress the vertical movement of air is called atmospheric stability.

[0110] The wind speed in the U direction usually represents the east-west wind speed component.

[0111] The wind speed in the V direction usually represents the north-south wind speed component.

[0112] As an example, V meteorological variables may include precipitation.

[0113] As another example, V meteorological variables include precipitation and stability.

[0114] As yet another example, the V meteorological variables may include stability, wind speed in the U direction, and wind speed in the V direction.

[0115] As yet another example, the V meteorological variables may include precipitation, stability, wind speed in the U direction, and wind speed in the V direction.

[0116] In some embodiments, the C dose includes at least one of a 2-day effective dose, a 7-day effective dose, and a thyroid dose.

[0117] In this embodiment, by considering at least one of the 2-day effective dose, the 7-day effective dose and the thyroid dose, the emergency plan area determined according to the target meteorological sequence can effectively cope with various radiations.

[0118] The 2-day effective dose refers to the effective dose received by the human body in a short period of time (usually within two days) in the event of a radiation accident or radiation exposure. It can be used to evaluate the impact of short-term radiation exposure on the human body and is usually used to quickly evaluate the radiation exposure of people after a radiation accident.

[0119] The 7-day effective dose refers to the effective dose received by the human body over a longer period of time (usually within seven days) under radiation exposure. It can be used to evaluate the impact of radiation exposure over a longer period of time on the human body and is usually used for long-term monitoring and evaluation after radiation accidents.

[0120] Thyroid dose is the amount of radiation received by the thyroid gland.

[0121] In some examples, dose C may include 2 days of effective doses.

[0122] In other examples, dose C may include a 7-day effective dose and a thyroid dose.

[0123] In yet other examples, type C doses may include a 2-day effective dose, a 7-day effective dose, and a thyroid dose.

[0124] It can be understood that the number of data in a meteorological series is V×T.

[0125] It should be noted that the values ​​of n, K, V and T can be set according to actual conditions and are not limited here. For example, n can be 50, K can be 50, V can be 4, T can be 48, etc.

[0126] In some embodiments, before obtaining n meteorological sequences corresponding to K target types respectively as target meteorological sequences and obtaining K×n target meteorological sequences, the method may further include:

[0127] Obtain Q meteorological sequences of the target nuclear power plant site, where Q meteorological sequences are meteorological sequences for each hour of the whole year; Q is a positive integer;

[0128] Q meteorological sequences are classified to obtain K target types, each target type includes multiple meteorological sequences, and K is less than Q.

[0129] In this embodiment, by obtaining Q meteorological sequences of the target nuclear power plant site and classifying the Q meteorological sequences, K target types are obtained, which can provide a basis for the subsequent acquisition of n meteorological sequences corresponding to the K target types as target meteorological sequences. In addition, since the Q meteorological sequences are meteorological sequences for each hour of the whole year, that is, there are 365*24=8760 meteorological sequences in total, the K target types obtained by classifying the Q meteorological sequences can cover various meteorological conditions, and then n meteorological sequences are extracted from each target type as target meteorological sequences, and the K×n target meteorological sequences obtained can cover various meteorological conditions, thereby ensuring the accuracy of the emergency plan area determination.

[0130] For example, the electronic device may store Q meteorological sequences of the target nuclear power plant site, so the Q meteorological sequences of the target nuclear power plant site may be directly called from the electronic device. Alternatively, the Q meteorological sequences of the target nuclear power plant site may be obtained in response to the Q meteorological sequences of the target nuclear power plant site input by the user based on the electronic device.

[0131] In some examples, Q meteorological sequences are classified to obtain K target types, including:

[0132] Select K meteorological sequence features from the Q meteorological sequence features that correspond one to one to the Q meteorological sequences as K initial classification centers;

[0133] For the remaining QK meteorological sequence features among the Q meteorological sequence features, respectively determine the first similarity between the meteorological sequence features and each initial classification center;

[0134] According to the first similarity, the QK meteorological sequences are classified into K initial classification centers respectively to obtain K target types.

[0135] In this embodiment, by respectively determining the first similarity between the meteorological sequence characteristics and each initial classification center, and classifying QK meteorological sequences into K initial classification centers according to the first similarity, K target types are obtained, thereby improving the rationality and accuracy of meteorological sequence classification.

[0136] Exemplarily, the electronic device may randomly select K meteorological sequence features from Q meteorological sequence features as K initial classification centers, or may select the first K meteorological sequence features with even numbers as K initial classification centers, which is not limited here.

[0137] In some examples, for the remaining QK meteorological sequence features among the Q meteorological sequence features, respectively determining the first similarity between the meteorological sequence features and each initial classification center may include:

[0138] For the remaining QK meteorological sequence features among the Q meteorological sequence features, the first similarity between the meteorological sequence features and each initial classification center is determined respectively according to the target similarity algorithm;

[0139] The target similarity algorithm includes any one of Euclidean distance, cosine similarity, edit distance and Manhattan distance.

[0140] In this embodiment, according to any similarity algorithm among Euclidean distance, cosine similarity, edit distance and Manhattan distance, the first similarity between the meteorological sequence characteristics and each initial classification center can be quickly and accurately determined, providing a basis for subsequently classifying QK meteorological sequences into K initial classification centers according to the first similarity to obtain K target types.

[0141] As an example, the target similarity algorithm may include Euclidean distance. The calculation formula of Euclidean distance may include formula (2).

[0142]

[0143] Among them, d can represent the first similarity, x vt It can represent the value of the vth meteorological variable in the meteorological sequence at the tth hour; y vt It can represent the value of the vth meteorological variable at the initial classification center at the tth hour; x VY It can represent the value of the Vth meteorological variable in the meteorological sequence at the Tth hour; VT It can represent the value of the Vth meteorological variable at the initial classification center at the Tth hour.

[0144] As another example, the target similarity algorithm includes cosine similarity.

[0145] As yet another example, the target similarity algorithm may include edit distance.

[0146] As yet another example, the target similarity algorithm may include Manhattan distance.

[0147] It is worth noting that for each meteorological sequence, it corresponds to K first similarities, that is, the first similarity between it and each initial classification center.

[0148] It should be noted that the similarity algorithm in the embodiment of the present application may also include other algorithms that can calculate the first similarity between the meteorological sequence and the initial classification center, which is not limited here.

[0149] In some examples, according to the first similarity, QK meteorological sequences are respectively classified into K initial classification centers to obtain K target types, including:

[0150] For QK meteorological sequences, the meteorological sequences are classified into the target type where the initial classification center corresponding to the maximum first similarity is located, and K target types are obtained.

[0151] In this embodiment, by classifying the meteorological sequence into the target type where the initial classification center corresponding to the maximum first similarity is located, K target types are obtained, which can improve the rationality and accuracy of the meteorological sequence classification.

[0152] Exemplarily, the K target types include the first target type where the first initial classification center is located, the second target type where the second initial classification center is located, and the third target type where the third initial classification center is located. The similarities between a certain meteorological sequence and the first initial classification center, the second initial classification center, and the third initial classification center are 0.98, 0.87, and 0.94, respectively. The meteorological sequence is classified into the first target type where the first initial classification center of the island is located.

[0153] It should be noted that the larger the Euclidean distance, the smaller the similarity; the smaller the Euclidean distance, the greater the similarity. The meteorological sequence is classified into the target type where the initial classification center corresponding to the maximum first similarity is located, and K target types are obtained. The meteorological sequence can be classified into the target type where the initial classification center corresponding to the minimum Euclidean distance is located, and K target types are obtained.

[0154] In some other examples, according to the first similarity, QK meteorological sequences are respectively classified into K initial classification centers to obtain K target types, including:

[0155] S21. According to the first similarity, the QK meteorological sequences are respectively classified into K initial classification centers to obtain K initial types;

[0156] S22, for each of the K initial types, taking the average value of multiple meteorological sequence features corresponding to the initial type as the initial classification center, to obtain K initial classification centers;

[0157] S23. Repeat step S21 and step S22 on an iterative basis, and count the number of executions until the number of executions is greater than or equal to a preset number of executions, and determine the K initial types finally obtained after the iteration as K target types.

[0158] In this embodiment, by repeatedly executing step S21 and step S22 on an iteration basis and counting the number of executions until the number of executions is greater than or equal to the preset number of executions, the K initial types finally obtained after the iteration are determined as K target types, which can further improve the accuracy of target type determination.

[0159] In S21, for the QK meteorological sequences, the meteorological sequences are classified into the initial types where the initial classification centers corresponding to the maximum first similarities are located, and K initial types are obtained.

[0160] In S22, illustratively, the numerical values ​​of the multiple meteorological sequences corresponding to the initial type may include a1, a2, a3, a4, a5, a6 ... a50; and (a1+a2+a3+ ... a50) / 50 may be used as the initial classification center.

[0161] In S23, the preset execution times can be set to any number according to actual conditions, and is not limited here. For example, the preset execution times can be 3 times, 5 times, 10 times, etc.

[0162] In S12, after obtaining n meteorological sequences corresponding to K target types as target meteorological sequences and obtaining K×n target meteorological sequences, the electronic device can also determine, for the first type of accident, under each target meteorological sequence, a first dose result set of C types of doses at D target distances and P target orientations corresponding to the release of the source item of the first type of accident, so as to obtain K×n groups of first dose result sets.

[0163] Category 1 accidents can be any type of accident that may occur in a nuclear power plant.

[0164] The release source term can be the amount or rate of radioactive material released into the environment in a nuclear accident or incident.

[0165] The accident set may include the set of all types of accidents that may occur in a nuclear power plant.

[0166] For example, any target distance among the D target distances can be expressed as d m ∈{d min , …, d max}, d min is the minimum target distance to be considered, d max is the maximum target distance to be considered.

[0167] It should be noted that the first dose result set may include dose results of C doses at D target distances and P target orientations. Each target distance in the D target distances corresponds to P target orientations, and each target orientation in the P target orientations corresponds to C doses. The values ​​of D, P, and C can be set according to actual conditions and are not limited here. For example, D is 20, P is 16, and C is 3.

[0168] In some embodiments, determining a first dose result set of C doses at D target distances and P target positions corresponding to a first type of accident release source term under each target meteorological sequence includes:

[0169] Determine the initial dose result set of C doses at Y distances and W azimuths corresponding to the first type of accident release source term under each target meteorological sequence, where 0 < D ≤ Y, 0 < P ≤ W, and both Y and W are integers;

[0170] Determine D target distances related to the emergency planning zone from the Y distances, and determine P target azimuths related to the emergency planning zone from the W azimuths;

[0171] Determine the first dose result set of C doses at D target distances and P target azimuths from the initial dose result set.

[0172] In this embodiment, by determining the first dose result set of C doses at D target distances and P target azimuths from the initial dose result set, since both the D target distances and the P target azimuths are related to the emergency planning zone, it can provide an accurate data basis for subsequent determination of the emergency planning zone, and thus improve the accuracy of determining the emergency planning zone.

[0173] It can be understood that the initial dose result set may include the dose results of C doses at Y distances and W azimuths corresponding to the first type of accident release source term. Each set of first dose result sets may include the dose results of C doses at D target distances and P target azimuths corresponding to the first type of accident release source term.

[0174] It should be noted that both Y and W can be set according to the actual situation and are not limited here. For example, Y can be 50 and W can be 20.

[0175] In S13, after the electronic device determines the first dose result set of C doses at D target distances and P target azimuths corresponding to the first type of accident release source term under each target meteorological sequence to obtain K × n sets of first dose result sets, it can also determine the second dose result set of C doses at D target distances and P target azimuths from the K × n sets of first dose result sets.

[0176] The second dose result set may include the dose results of C doses at D target distances and P target azimuths.

[0177] In some embodiments, determining the second dose result set of C doses at D target distances and P target azimuths from the K × n sets of first dose result sets may include:

[0178] For each target meteorological sequence, obtain the first weight value corresponding to the target meteorological sequence;

[0179] For each target meteorological sequence, take the product of the first weight value and the first dose result set corresponding to the target meteorological sequence as the fourth dose result set to obtain K × n sets of fourth dose result sets;

[0180] Sorting the fourth dose result set of K×n groups to obtain a first sorting;

[0181] The fourth dose result set located at a preset percentile in the first sorting is used as the second dose result set.

[0182] In this embodiment, for each target meteorological sequence, the product of the first weight value and the first dose result set corresponding to the target meteorological sequence is used as the fourth dose result set, and then the K×n groups of fourth dose result sets are sorted to obtain the first sorting, and then the fourth dose result set located at a preset percentile in the first sorting is used as the second dose result set. Compared with using all first dose result sets as the second dose result set, the amount of calculation, time spent and storage space in the process of determining the entire emergency plan area can be reduced.

[0183] In some examples, the first weight values ​​corresponding to the target meteorological sequences under the same target type are the same;

[0184] Obtaining a first weight value corresponding to the target meteorological sequence may include:

[0185] Acquire a first quantity and a total quantity of meteorological sequences of the target nuclear power plant site, wherein the first quantity is the quantity of meteorological sequences included in the target type corresponding to the target meteorological sequence;

[0186] Substitute the first quantity and the total quantity into formula (1) to calculate and obtain the first weight value corresponding to the target meteorological sequence;

[0187] Formula (1) includes:

[0188]

[0189] Wherein, M represents the first weight value; N k represents a first quantity; N k Indicates the total quantity.

[0190] In this embodiment, by substituting the first number and the total number into formula (1), the first weight value corresponding to the target meteorological sequence can be determined quickly and accurately.

[0191] For example, the electronic device may store the first quantity and the total quantity of the meteorological sequence of the target nuclear power plant site, and then directly call the first quantity and the total quantity. Alternatively, the first quantity and the total quantity may be obtained in response to the first quantity and the total quantity input by the user based on the electronic device.

[0192] In the embodiment of the present application, for each target meteorological sequence, the first dose result sets of C types of doses at D distances and P directions of concern for the emergency plan area division are weighted and sorted respectively, and sorted D×P×C times in total.

[0193] Exemplarily, the K×n groups of fourth dose result sets may be sorted in order from small to large, or from large to small, to obtain a first sorting.

[0194] The preset percentile can be set according to actual conditions and is not limited here. For example, the preset percentile can be 70%.

[0195] In S14, after the electronic device determines a second dose result set of C doses at D target distances and P target orientations from the K×n group of first dose result sets, it can also determine, for each of the D target distances in the second dose result set, a third dose result set of C doses corresponding to the target distance from the second dose result set based on the dose result set of C doses at P target orientations corresponding to the target distance in the second dose result set, so as to obtain a C×D group of third dose result sets.

[0196] Each third dose result set may include dose results of C doses corresponding to the target distance.

[0197] In some embodiments, according to the dose result set of C doses at P target positions corresponding to the target distances in the second dose result set, determining a third dose result set of C doses corresponding to the target distances from the second dose result may include:

[0198] Determine from the second dose result set a fifth dose result set of C doses corresponding to each target orientation at the target distance;

[0199] The fifth dose result set with the largest value is determined as the third dose result set of C doses corresponding to the target distance.

[0200] In this embodiment, by determining the fifth dose result set with the largest value as the third dose result set of C doses corresponding to the target distance, the emergency planning area determined subsequently can meet the fifth dose result set with the largest value, and naturally can meet other fifth dose result sets, thereby improving the applicability of the emergency planning area.

[0201] That is to say, a certain target distance corresponds to P target positions, each target position corresponds to the fifth dose result set of C doses, that is, the P target positions correspond one-to-one to P fifth dose result sets, and the fifth dose result set with the largest value among the P fifth dose result sets is taken as the third dose result set of C doses corresponding to the target distance.

[0202] In S15, the electronic device determines a third dose result set of C types of doses corresponding to the target distance from the second dose result set according to the dose result set of C types of doses at P target orientations corresponding to the target distance in the second dose result set for each target distance among the D target distances in the second dose result set, so as to obtain a C×D group third dose result set. Then, the electronic device can also determine an emergency plan area suitable for the first type of accident according to the C×D group third dose result set and the preset target dose result set.

[0203] The emergency planning area can be an area established in advance around a nuclear facility with emergency plans and emergency preparations so that timely and effective protective actions can be taken to protect the public when an accident occurs at a nuclear facility.

[0204] In some embodiments, determining the emergency plan area applicable to the first type of accident according to C×D third dose result sets and a preset target dose result set may include:

[0205] For each target distance, according to the third dose result set of C doses corresponding to the target distance and the preset target dose result set, C initial radius ranges corresponding to the C doses are determined to obtain C×D initial radius ranges;

[0206] For each target distance, determine a first radius range from C initial radius ranges corresponding to the target distance to obtain D first radius ranges;

[0207] For each target distance, the intersection of D first radius ranges is used as the target radius range of the emergency planning area for the first type of accident;

[0208] The area formed by taking the center of the reactor in the target nuclear power plant as the center of the circle and the maximum target distance within the target radius as the radius is determined as the emergency planning area suitable for the first type of accidents.

[0209] In this embodiment, C×D initial radius ranges are first determined, and then D first radius ranges are determined from the C×D initial radius ranges. The intersection of the D first radius ranges is then used as the target radius range of the emergency planning area for the first type of accident. Subsequently, an area formed with the center of the reactor in the target nuclear power plant as the center of the circle and the maximum target distance within the target radius range as the radius is determined as the emergency planning area suitable for the first type of accident. This allows the emergency planning area suitable for the first type of accident to be accurately determined.

[0210] The target dose result set may include a common optimized intervention level or action level for C doses at D target distances.

[0211] A preset target dose result set may be stored in the electronic device so that the target dose result set can be directly called later.

[0212] In some examples, according to the third dose result set of C doses corresponding to the target distance and the preset target dose result set, C initial radius ranges corresponding to the C doses are determined, including:

[0213] For the first dose result corresponding to each dose in the third dose result set, when the first dose result is less than or equal to the target dose result corresponding to the dose in the preset target dose result set, the initial radius range is determined to be greater than 0 and less than or equal to the target distance corresponding to the first dose result, or, when the first dose result is greater than the target dose result, the initial radius range is determined to be greater than the target distance corresponding to the first dose result, so as to determine C initial radius ranges corresponding to the C doses one by one.

[0214] In this embodiment, by comparing the first dose result with the target dose result, the initial radius range can be determined quickly and accurately.

[0215] In this embodiment, when the first dose result is less than or equal to the target dose result, it means that the public at the location does not need to take shelter, evacuate, or take iodine protection actions, so the initial radius range is determined to be greater than 0 and less than or equal to the target distance corresponding to the first dose result. When the first dose result is greater than the target dose result, it means that the public at the location should take shelter, evacuate, or take iodine protection actions, so the initial radius range is determined to be greater than the target distance corresponding to the first dose result.

[0216] It should be noted that the third dose result set may include the first dose result corresponding to each of the C doses, and the target dose result set may include the target dose result corresponding to each of the C doses. The target dose result compared with the first dose result corresponds to the same dose as the first dose result, for example, the first dose result corresponding to the c1 dose is less than or equal to the target dose result corresponding to the c1 dose.

[0217] Exemplarily, when the target distance is 2km, the third dose result set includes a 2-day effective dose of 15mSv, a 7-day effective dose of 40mSV, and a thyroid dose of 150mGy; the target dose result set may include a 2-day effective dose of 10mSV, a 7-day effective dose of 50mSV, and a thyroid dose of 100mGy. Among them, a 2-day effective dose of 15mSv, a 7-day effective dose of 40mSV, and a thyroid dose of 150mGy may be three first dose results; a 2-day effective dose of 10mSV, a 7-day effective dose of 50mSV, and a thyroid dose of 100mGy may be three target dose results. The initial radius range corresponding to the 2-day effective dose is greater than 2km, the initial radius range corresponding to the 7-day effective dose is less than or equal to 2km, and the initial radius range corresponding to the thyroid dose is greater than 2km.

[0218] In some examples, determining the first radius range from C initial radius ranges corresponding to the target distance includes:

[0219] When the C initial radius ranges corresponding to the target distances are all the same, the initial radius range is determined as the first radius range;

[0220] Alternatively, when there are at least two different initial radius ranges among the C initial radius ranges corresponding to the target distance, the second radius range among the C initial radius ranges is determined as the first radius range;

[0221] Among them, any value within the second radius range is greater than or equal to all values ​​within the other C-1 initial radius ranges except the second radius range.

[0222] In this embodiment, the first radius range can be quickly determined in the above manner, providing a basis for subsequent determination of the emergency plan area.

[0223] As an example, if the C initial radius ranges are all greater than 2 km, then the radius range greater than 2 km is determined as the first radius range.

[0224] As another example, if c of the C initial radius ranges are all greater than 2 km, and Cc initial radius ranges are all less than or equal to 2 km, then greater than 2 km is determined as the first radius range, where c is an integer greater than 1 and can be set according to actual conditions.

[0225] In some embodiments, after determining the emergency plan area applicable to the first type of accident according to C×D third dose result sets and the preset target dose result set, the method further includes:

[0226] S16. All accidents in the accident set except the first category accidents are regarded as the first category accidents;

[0227] S17, repeating steps S12 to S16 until an emergency plan area corresponding to each accident in the accident set is obtained;

[0228] S18. Determine the largest emergency planning area among all emergency planning areas as the emergency planning area applicable to the accident collection.

[0229] In this embodiment, the emergency planning area with the largest area among all emergency planning areas is used as the emergency planning area applicable to the accident set. Since the emergency planning area with the largest area can cover the emergency planning areas corresponding to other accidents, it is ensured that the safety of personnel and the environment can be effectively protected in all accidents in the accident set.

[0230] In order to better understand the method for determining the emergency plan area provided in the embodiment of the present application, it is described below in conjunction with a specific embodiment.

[0231] Taking a certain nuclear power plant site (i.e., the target nuclear power plant site) as an example, the division of the smoke plume emergency planning area (i.e., the emergency planning area) includes the following steps:

[0232] (1) A serious accident (i.e., Category I accident) that has been demonstrated to represent most serious accident conditions is selected. Considering 8760 meteorological series at a nuclear power plant site in 2018, the 48-hour (i.e., T-hour) changes in four meteorological variables, namely, precipitation, stability, wind speed in the U direction, and wind speed in the V direction (i.e., V meteorological variables) at the meteorological station are taken as meteorological series features (i.e., meteorological series). The number of data for a meteorological series feature is 4×48;

[0233] (2) Randomly select 50 (i.e., K) meteorological sequence features from the annual meteorological sequence features (i.e., Q meteorological sequence features) as the initial classification centers;

[0234] (3) Calculate the Euclidean distance between each meteorological sequence feature and 50 classification centers (i.e., initial classification centers) in turn, and classify each meteorological sequence into the class (i.e., target type) where the classification center with the smallest Euclidean distance is located. Take the mean of all meteorological sequence features of each class in turn (i.e., the average of multiple meteorological sequence features corresponding to the initial type) to re-obtain the classification center (initial classification center) of the class, and repeat this step 50 times (i.e., the preset number of executions).

[0235] Figure 2-Figure 5 The 48-hour changes of each meteorological variable after all meteorological sequences of categories 1, 11, 21, and 31 are superimposed after 50 repetitions are given. The four sub-graphs of each graph are the 48-hour changes of precipitation (mm / h), atmospheric stability (i.e. stability), wind speed in the U direction (m / s), and wind speed in the V direction (m / s). The horizontal axis of each sub-graph corresponds to the time (h) of the meteorological sequence (i.e. meteorological sequence) from the starting time, and the vertical axis is the corresponding meteorological variable value, that is, Figures 2 to 5 The vertical axis of the upper left sub-image is the precipitation value. Figures 2 to 5 The vertical axis of the upper right sub-graph is the atmospheric stability value. Figures 2 to 5 The vertical coordinate of the lower left sub-graph is the wind speed value in the U direction. Figures 2 to 5 The vertical coordinate of the lower right sub-graph is the wind speed value in the V direction. Figure 2-5 As shown in the figure, various meteorological sequences (i.e., each meteorological variable of meteorological sequence 0 has a darker overlapping part, showing an obvious classification center (initial classification center)), and the classification centers of different categories are different. Therefore, it can be considered that meteorological sequences classified into the same category have similar meteorological characteristics, and meteorological sequences of different categories have different characteristics.

[0236] (4) 10 (i.e., n) meteorological sequences are selected from each class (i.e., target type) as samples (i.e., target meteorological sequences), and the total number of meteorological sequence samples selected is 50 × 10 (i.e., K × n target meteorological sequences);

[0237] (5) According to the meteorological sequence samples drawn, the 2-day effective dose, 7-day effective dose and thyroid dose (i.e., C-type dose) after the release of the source term of the severe accident under different meteorological conditions are calculated in turn, with a total of 50×10 groups of calculation results (i.e., K×n groups of first dose result sets). The calculation results (i.e., K×n groups of first dose result sets) at 0.5 km, 1 km, 1.5 km, 2 km, …, 10 km (i.e., D target distances) of the nuclear power plant site that need to be concerned about the emergency planning area are weighted and sorted respectively, and a total of 20×16×3 sortings are performed. The weight (i.e., the first weight value) of the calculation result obtained according to the k-th meteorological sequence is Where 10 is the number of meteorological sequences extracted from each class (i.e., n), N k is the total number of meteorological sequences classified into the kth category (i.e., the first number);

[0238] (6) Considering the 2-day effective dose, 7-day effective dose and thyroid dose at the 70% percentile (i.e., the preset percentile) of each direction, they are compared with the universal optimized intervention levels of shelter, evacuation and iodine protection (2-day effective dose 10mSv, 7-day effective dose 50mSv and thyroid dose 100mGy) (i.e., the target dose result set). Taking the dose level at 2km (i.e., the target distance) as an example, the calculated results are 15mSv, 80mSv and 150mGy (i.e., the third dose result set), which are greater than the corresponding universal optimized intervention levels, indicating that shelter, evacuation and iodine protection actions should be taken for the public at this location. Therefore, the range of the smoke plume emergency plan area (i.e., the initial radius range) should be greater than 2km. The same judgment is made for each distance considered, and the radius of the smoke plume emergency plan area under this severe accident condition (i.e., the target radius range) is 3km;

[0239] (7) For each type of design basis accident, most severe accidents and the most severe accident with the most serious consequences (i.e., accidents other than the first type of accident in the accident set) that needs to be considered in demarcating the emergency planning area, steps (5) and (6) are repeated, and the radius of the smoke plume emergency planning area for the nuclear power plant site is finally determined to be 4 km (i.e., the maximum target distance within the target radius corresponding to the largest emergency planning area).

[0240] like Figure 6 As shown, the method for determining the emergency plan area provided in the embodiment of the present application may include steps S101 to S114.

[0241] S101. For type A accidents (i.e., the accident set), clarify the B intervention levels or action levels (i.e., the target dose result set) that need to be considered in the emergency plan area;

[0242] S102. Extract the characteristics of the annual meteorological sequence;

[0243] S103. Randomly select K meteorological sequence characteristics as the initial stratification centers (i.e., the initial classification centers);

[0244] S104. Assign each meteorological sequence to the class where the stratification center with the smallest Euclidean distance is located (i.e., classify the meteorological sequence into the initial type where the initial classification center with the maximum first similarity is located);

[0245] S105. Take the mean value to re-obtain the class center (i.e., take the average value of the multiple meteorological sequence characteristics corresponding to the initial type as the initial classification center);

[0246] S106. Extract n meteorological sequences from each class as samples (i.e., the target meteorological sequences);

[0247] Among them, steps S105 and S106 are looped j times (i.e., the preset execution times);

[0248] S107. For the release source term of a certain accident a (i.e., the first type of accident), calculate the results under the selected meteorological conditions; calculate a total of K×n times;

[0249] S108. Sort the calculation results for each distance (i.e., D target distances) and each azimuth (i.e., P target azimuths) of the b dose levels (i.e., C doses) that need to be considered after the release source term of accident a, and sort a total of b×D×P times (i.e., C×D×P times);

[0250] S109. Take the maximum dose level D dm (x) at the specified percentile x% (i.e., the preset percentile) of the calculation results as the result at distance d m (i.e., determine the fifth dose result set with the largest value as the third dose result set of the C doses corresponding to the target distance), considering a certain action level or intervention level b (i.e., the target dose result), where d m represents the mth target distance;

[0251] S110. Judge whether D dm (x) < b holds;

[0252] If it does not hold and m = m + 1 & m ≤ max, repeat step S110;

[0253] If it holds and m = m - 1 & m ≥ min, execute step S111;

[0254] S111, judgement D dm (x) ≥ b whether it holds;

[0255] If not, and m=m-1&m≥min, repeat step S111;

[0256] If true, m=m+1, and execute step S112;

[0257] S112. Determine the scope of the emergency planning area r of the nuclear power plant site after the release of the source term of accident a and considering the action level or intervention level b. ab (i.e. determine the initial radius range corresponding to a dose);

[0258] Wherein, step S108 to step S112 are calculated B times in total;

[0259] S113. Consider the emergency planning area for all dose levels B related to accident a, and take the maximum value as r a (i.e., determine C initial radius ranges corresponding to C doses);

[0260] S114. Consider all accident types A (i.e., accident set), and take the maximum range to determine the emergency planning area range r (i.e., target radius range) of the nuclear power plant site.

[0261] The embodiment of the present application proposes a method for determining an emergency plan area, which can greatly save calculation costs. The method starts from the annual meteorological sequence of the plant site, extracts the meteorological sequence data features (i.e., meteorological sequence features) required for calculating the dose results (i.e., the dose result set) in sequence, and uses the Euclidean distance between the meteorological sequence features as the basis for determining the similarity of the meteorological sequences. The annual meteorological sequence is classified in this way, and meteorological sequence samples (i.e., target meteorological sequences) are sequentially extracted from each category to fully consider various meteorological conditions throughout the year. The extracted meteorological sequence samples represent the meteorological sequence throughout the year, which is used as meteorological input to obtain all dose calculation results (i.e., the fourth dose result set) under the sampled meteorological sequence (i.e., the target meteorological sequence) and sort them. The dose levels of various types at different distances from the source release point are compared with the general optimized intervention level or action level to determine the scope of the emergency plan area.

[0262] The division of emergency planning areas needs to consider the dose impact on the public under various meteorological conditions after the accident release source items. The embodiment of the present application proposes the idea of ​​meteorological sequence classification sampling, which classifies the meteorological sequence features throughout the year according to the Euclidean distance between them, and extracts meteorological sequences from each category as meteorological inputs for calculation. While reducing the calculation time, the extracted meteorological sequences can also cover various meteorological conditions throughout the year. The embodiment of the present application can provide ideas for optimizing the division of emergency planning areas in nuclear power plants.

[0263] Example 2

[0264] Next, a device for determining an emergency plan area provided in an embodiment of the present application is introduced.

[0265] like Figure 7 As shown, the device for determining the emergency plan area provided in the embodiment of the present application may include a first acquisition module 710 , a first determination module 720 , a second determination module 730 , a third determination module 740 and a fourth determination module 750 .

[0266] The first acquisition module 710 is used to acquire n meteorological sequences corresponding to K target types as target meteorological sequences, and obtain K×n target meteorological sequences, wherein the meteorological sequences include values ​​of V meteorological variables at T hours at the meteorological stations corresponding to the target nuclear power plant sites, and n, K, V and T are all positive integers;

[0267] The first determination module 720 is connected to the first acquisition module 710 and is used to determine, for the first type of accident, a first dose result set of C doses at D target distances and P target directions corresponding to the release of the source term of the first type of accident under each target meteorological sequence, so as to obtain K×n groups of first dose result sets, the first type of accident is any one of the accident set, and D, P, and C are all positive integers;

[0268] The second determination module 730 is connected to the first determination module 720 and is used to determine a second dose result set of C doses at D target distances and P target orientations from the K×n groups of first dose result sets;

[0269] The third determination module 740 is connected to the second determination module 730 and is configured to determine, for each of the D target distances in the second dose result set, a third dose result set of C doses corresponding to the target distance from the second dose result set according to the dose result sets of C doses at P target orientations corresponding to the target distance in the second dose result set, so as to obtain a C×D group of third dose result sets;

[0270] The fourth determination module 750 is connected to the third determination module 740 and is used to determine the emergency plan area applicable to the first type of accident according to the third dose result set of the C×D group and the preset target dose result set.

[0271] According to the device for determining the emergency plan area provided in the embodiment of the present application, by acquiring n meteorological sequences corresponding to K target types respectively as target meteorological sequences, K×n target meteorological sequences are obtained. On the one hand, since n meteorological sequences are extracted from each target type as target meteorological sequences, the obtained K×n target meteorological sequences can cover various meteorological conditions, thereby ensuring the accuracy of the emergency plan area determination; on the other hand, compared with the prior art of taking all meteorological sequences as target meteorological sequences, the number of target meteorological sequences can be reduced, thereby reducing the amount of calculation, time spent and storage space in the entire emergency plan area determination process.

[0272] In some embodiments, the device further comprises:

[0273] The second acquisition module is used to acquire Q meteorological sequences of the target nuclear power plant site, wherein the Q meteorological sequences are meteorological sequences for each hour of the whole year; Q is a positive integer;

[0274] The first classification module is connected to the second acquisition module and is used to classify Q meteorological sequences to obtain K target types, each target type includes multiple meteorological sequences, and K is less than Q.

[0275] In some embodiments, the first classification module comprises:

[0276] A selection submodule, connected to the second acquisition module, is used to select K meteorological sequence features from the Q meteorological sequence features corresponding to the Q meteorological sequences as K initial classification centers;

[0277] A first determination submodule is connected to the selection submodule and is used to determine the first similarity between the meteorological sequence features and each initial classification center for the remaining QK meteorological sequence features among the Q meteorological sequence features;

[0278] The first classification submodule is connected to the selection submodule and is used to classify the QK meteorological sequences into K initial classification centers according to the first similarity to obtain K target types.

[0279] In some implementations, the first determination submodule may be specifically used to:

[0280] For the remaining QK meteorological sequence features among the Q meteorological sequence features, the first similarity between the meteorological sequence features and each initial classification center is determined respectively according to the target similarity algorithm;

[0281] The target similarity algorithm includes any one of Euclidean distance, cosine similarity, edit distance and Manhattan distance.

[0282] In some implementations, the first classification submodule may be specifically used to:

[0283] For Q - K meteorological sequences, classify the meteorological sequences into the target types where the initial classification centers corresponding to their maximum first similarity are located, to obtain K target types;

[0284] Alternatively, according to the first similarity, classify the Q - K meteorological sequences into K initial classification centers respectively, to obtain K target types, including:

[0285] S21. According to the first similarity, classify the Q - K meteorological sequences into K initial classification centers respectively, to obtain K initial types;

[0286] S22. For each initial type among the K initial types, take the average value of the multiple meteorological sequence features corresponding to the initial type as the initial classification center, to obtain K initial classification centers;

[0287] S23. On the basis of iteration, repeat steps S21 and S22, and count the number of executions. Until the number of executions is greater than or equal to the preset number of executions, determine the K initial types finally obtained after iteration as the K target types.

[0288] In some embodiments, the first determination module 720 may be specifically configured to:

[0289] Determine, for each target meteorological sequence, an initial dose result set of C doses at Y distances and W azimuths corresponding to the first - type accident release source term, where 0 < D ≤ Y, 0 < P ≤ W, and both Y and W are integers;

[0290] Determine D target distances related to the emergency planning zone from the Y distances, and determine P target azimuths related to the emergency planning zone from the W azimuths;

[0291] Determine a first dose result set of C doses at D target distances and P target azimuths from the initial dose result set.

[0292] In some embodiments, the second determination module 730 may be specifically configured to:

[0293] An acquisition sub - module, configured to, for each target meteorological sequence, acquire a first weight value corresponding to the target meteorological sequence;

[0294] A first obtaining sub - module, connected to the acquisition sub - module, configured to, for each target meteorological sequence, take the product of the first weight value and the first dose result set corresponding to the target meteorological sequence as the fourth dose result set, so as to obtain K×n groups of fourth dose result sets;

[0295] A second obtaining sub - module, connected to the first obtaining sub - module, configured to sort the K×n groups of fourth dose result sets to obtain a first sorting;

[0296] The first one is used as a submodule, connected to the second obtaining submodule, and is used to take the fourth dose result set located at a preset percentile in the first sorting as the second dose result set.

[0297] In some implementations, the first weight values ​​corresponding to the target meteorological sequences under the same target type are the same;

[0298] Get submodules, which can be used specifically for:

[0299] Acquire a first quantity and a total quantity of meteorological sequences of the target nuclear power plant site, wherein the first quantity is the quantity of meteorological sequences included in the target type corresponding to the target meteorological sequence;

[0300] Substitute the first quantity and the total quantity into formula (1) to calculate and obtain the first weight value corresponding to the target meteorological sequence;

[0301] Formula (1) includes:

[0302]

[0303] Wherein, M represents the first weight value; N k represents a first quantity; N k Indicates the total quantity.

[0304] In some implementations, the third determination module 740 may be specifically configured to:

[0305] Determine from the second dose result set a fifth dose result set of C doses corresponding to each target orientation at the target distance;

[0306] The fifth dose result set with the largest value is determined as the third dose result set of C doses corresponding to the target distance.

[0307] In some implementations, the fourth determination module 750 may include:

[0308] The second determination submodule is used to determine, for each target distance, C initial radius ranges corresponding to the C types of doses according to the third dose result set of the C types of doses corresponding to the target distance and the preset target dose result set, so as to obtain C×D initial radius ranges;

[0309] A third determination submodule is used to determine, for each target distance, a first radius range from C initial radius ranges corresponding to the target distance, so as to obtain D first radius ranges;

[0310] The second is a submodule for taking the intersection of D first radius ranges as the target radius range of the emergency planning area for the first type of accident for each target distance;

[0311] The fourth determination submodule is used to determine an area formed by taking the center of the reactor in the target nuclear power plant as the center of the circle and the maximum target distance within the target radius as the radius as the emergency plan area suitable for the first type of accident.

[0312] In some implementations, the second determining submodule may be specifically configured to:

[0313] For the first dose result corresponding to each dose in the third dose result set, when the first dose result is less than or equal to the target dose result corresponding to the dose in the preset target dose result set, the initial radius range is determined to be greater than 0 and less than or equal to the target distance corresponding to the first dose result, or, when the first dose result is greater than the target dose result, the initial radius range is determined to be greater than the target distance corresponding to the first dose result, so as to determine C initial radius ranges corresponding to the C doses one by one.

[0314] In some implementations, the third determining submodule may be specifically used to:

[0315] When the C initial radius ranges corresponding to the target distances are all the same, the initial radius range is determined as the first radius range;

[0316] Alternatively, when there are at least two different initial radius ranges among the C initial radius ranges corresponding to the target distance, the second radius range among the C initial radius ranges is determined as the first radius range;

[0317] Among them, any value within the second radius range is greater than or equal to all values ​​within the other C-1 initial radius ranges except the second radius range.

[0318] In some embodiments, the device may further include:

[0319] As a module, it is used to treat all accidents in the accident set except the first-category accidents as the first-category accidents;

[0320] A repeat execution module, connected to the as module, for repeatedly executing steps S12 to S16 until an emergency plan area corresponding to each accident in the accident set is obtained;

[0321] The fifth determination module is connected to the repeated execution module and is used to determine the emergency plan area with the largest area among all the emergency plan areas as the emergency plan area suitable for the accident set.

[0322] In some embodiments of the first aspect, the V meteorological variables include at least one of precipitation, stability, wind speed in the U direction, and wind speed in the V direction;

[0323] And / or, type C dose includes at least one of a 2-day effective dose, a 7-day effective dose and a thyroid dose.

[0324] The device for determining the emergency plan area provided in the embodiment of the present application has the beneficial effects and implementation methods of the method for determining the emergency plan area provided in Example 1 of the present application. For details, please refer to the specific description of the method for determining the emergency plan area in the above Example 1, which will not be repeated in this embodiment.

[0325] It is to be understood that the above embodiments are merely exemplary embodiments used to illustrate the principles of the present application, but the present application is not limited thereto. For those skilled in the art, various modifications and improvements can be made without departing from the spirit and substance of the present application, and these modifications and improvements are also considered to be within the scope of protection of the present application.

Claims

1. A method for determining an emergency planning area, characterized in that: include: S11, obtaining n meteorological sequences corresponding to K target types respectively as target meteorological sequences, and obtaining K×n target meteorological sequences, wherein the meteorological sequences include values ​​of V meteorological variables at T hours at the meteorological stations corresponding to the target nuclear power plant sites, and n, K, V and T are all positive integers; S12. For the first type of accident, determine the first dose result set of C doses at D target distances and P target directions corresponding to the release of the source term of the first type of accident under each target meteorological sequence, so as to obtain K×n groups of the first dose result sets, where the first type of accident is any one of the accident set, and D, P, and C are all positive integers; S13, determining a second dose result set of C doses at D target distances and P target positions from the K×n groups of the first dose result sets; S14. For each target distance among the D target distances in the second dose result set, according to the dose result sets of C doses at P target orientations corresponding to the target distance in the second dose result set, determine a third dose result set of C doses corresponding to the target distance from the second dose result set, so as to obtain a C×D group of the third dose result sets; S15. Determine an emergency plan area applicable to the first type of accident based on the third dose result set of group C×D and a preset target dose result set.

2. The method according to claim 1, characterized in that Before acquiring n meteorological sequences corresponding to K target types respectively as target meteorological sequences to obtain K×n target meteorological sequences, the method further includes: Obtain Q meteorological sequences of the target nuclear power plant site, where the Q meteorological sequences are meteorological sequences for each hour throughout the year, and Q is a positive integer; Q meteorological sequences are classified to obtain K target types, each target type includes multiple meteorological sequences, and K is less than Q.

3. The method according to claim 2, characterized in that The Q meteorological sequences are classified to obtain K target types, including: Select K meteorological sequence features from the Q meteorological sequence features that correspond one to one to the Q meteorological sequences as K initial classification centers; For the remaining QK meteorological sequence features among the Q meteorological sequence features, respectively determine the first similarity between the meteorological sequence features and each initial classification center; According to the first similarity, QK meteorological sequences are respectively classified into K initial classification centers to obtain K target types.

4. The method according to claim 3, characterized in that: The step of determining the first similarity between the remaining QK meteorological sequence features and each initial classification center respectively for the meteorological sequence features among the Q meteorological sequence features includes: For the remaining QK meteorological sequence features among the Q meteorological sequence features, respectively determine the first similarity between the meteorological sequence features and each initial classification center according to the target similarity algorithm; The target similarity algorithm includes any one of Euclidean distance, cosine similarity, edit distance and Manhattan distance.

5. The method according to claim 3 or 4, characterized in that: According to the first similarity, the QK meteorological sequences are respectively classified into K initial classification centers to obtain K target types, including: For Q - K meteorological sequences, classify the meteorological sequences into the target type where the initial classification center corresponding to their maximum first similarity is located, to obtain K target types; Alternatively, the step of classifying the Q - K meteorological sequences into K initial classification centers respectively according to the first similarity to obtain K target types includes: S21. Classify the Q - K meteorological sequences into K initial classification centers respectively according to the first similarity to obtain K initial types; S22. For each initial type among the K initial types, use the average value of the multiple meteorological sequence features corresponding to the initial type as the initial classification center to obtain K initial classification centers; S23. Repeat steps S21 and S22 iteratively and count the number of executions until the number of executions is greater than or equal to the preset number of executions, and determine the K initial types finally obtained after iteration as K target types.

6. The method according to any one of claims 1 to 4, characterized in that The step of determining the first dose result set of C doses at D target distances and P target azimuths after the first - type accident release source term under each target meteorological sequence includes: Determine the initial dose result set of C doses at Y distances and W azimuths after the first - type accident release source term under each target meteorological sequence, where 0 < D ≤ Y, 0 < P ≤ W, and both Y and W are integers; Determine D target distances related to the emergency planning zone from the Y distances, and determine P target azimuths related to the emergency planning zone from the W azimuths; Determine the first dose result set of C doses at D target distances and P target azimuths from the initial dose result set.

7. The method according to any one of claims 1 to 4, characterized in that The step of determining the second dose result set of C doses at D target distances and P target azimuths from the K×n groups of the first dose result sets includes: For each target meteorological sequence, obtain the first weight value corresponding to the target meteorological sequence; For each target meteorological sequence, use the product of the first weight value and the first dose result set corresponding to the target meteorological sequence as the fourth dose result set to obtain K×n groups of the fourth dose result sets; Sort the K×n groups of the fourth dose result sets to obtain the first sorting; Use the fourth dose result set at the preset quantile in the first sorting as the second dose result set.

8. The method according to claim 7, characterized in that The first weight values corresponding to the target meteorological sequences under the same target type are the same; The step of obtaining the first weight value corresponding to the target meteorological sequence includes: obtaining the first quantity and the total number of meteorological sequences at the target nuclear power plant site, where the first quantity is the number of meteorological sequences included in the target type corresponding to the target meteorological sequence; Substitute the first quantity and the total number into formula (1) for calculation to obtain the first weight value corresponding to the target meteorological sequence; The formula (1) includes: Wherein, M represents the first weight value; N k represents a first quantity; Indicates the total quantity.

9. The method according to claim 1, characterized in that: The step of determining a third dose result set of C doses corresponding to the target distance from the second dose result according to the dose result set of C doses at P target orientations corresponding to the target distance in the second dose result set includes: determining a fifth dose result set of C doses corresponding to each target orientation at the target distance from the second dose result set; The fifth dose result set with the largest value is determined as the third dose result set of C doses corresponding to the target distance.

10. The method according to claim 1, characterized in that Determining the emergency plan area applicable to the first type of accident according to the C×D third dose result sets and the preset target dose result set includes: For each target distance, according to the third dose result set of C doses corresponding to the target distance and the preset target dose result set, C initial radius ranges corresponding to the C doses are determined to obtain C×D initial radius ranges; For each target distance, determining a first radius range from C initial radius ranges corresponding to the target distance to obtain D first radius ranges; For each target distance, the intersection of D first radius ranges is used as the target radius range of the emergency planning area for the first type of accident; An area formed by taking the center of the reactor in the target nuclear power plant as the center of the circle and the maximum target distance within the target radius as the radius is determined as the emergency planning area suitable for the first type of accident.

11. The method according to claim 10, characterized in that The determining of C initial radius ranges corresponding to the C doses in one-to-one correspondence according to the third dose result set of the C doses corresponding to the target distance and the preset target dose result set includes: For the first dose result corresponding to each dose in the third dose result set, when the first dose result is less than or equal to the target dose result corresponding to the dose in the preset target dose result set, the initial radius range is determined to be greater than 0 and less than or equal to the target distance corresponding to the first dose result, or, when the first dose result is greater than the target dose result, the initial radius range is determined to be greater than the target distance corresponding to the first dose result, so as to determine C initial radius ranges corresponding to C doses one by one.

12. The method according to claim 10, characterized in that Determining a first radius range from the C initial radius ranges corresponding to the target distance includes: When the C initial radius ranges corresponding to the target distances are all the same, determining the initial radius range as the first radius range; Alternatively, when at least two of the C initial radius ranges corresponding to the target distance are different, the second radius range among the C initial radius ranges is determined as the first radius range; Among them, any value within the second radius range is greater than or equal to all values ​​within the other C-1 initial radius ranges except the second radius range.

13. The method according to claim 1, characterized in that After determining the emergency plan area applicable to the first type of accident according to the C×D third dose result sets and the preset target dose result set, the method further includes: S16. taking all accidents in the accident set except the first category of accidents as the first category of accidents; S17, repeating steps S12 to S16 until an emergency plan area corresponding to each accident in the accident set is obtained; S18. Determine the largest emergency planning area among all emergency planning areas as the emergency planning area applicable to the accident set.

14. The method according to any one of claims 1 to 4 and 8 to 14, characterized in that The meteorological variables of type V include at least one of precipitation, stability, wind speed in the U direction and wind speed in the V direction; And / or, type C dose includes at least one of a 2-day effective dose, a 7-day effective dose and a thyroid dose.

15. A device for determining an emergency planning area, characterized in that: include: The first acquisition module is used to acquire n meteorological sequences corresponding to K target types respectively as target meteorological sequences, and obtain K×n target meteorological sequences, wherein the meteorological sequences include values ​​of V meteorological variables at T hours at the meteorological stations corresponding to the target nuclear power plant sites, and n, K, V and T are all positive integers; A first determination module is connected to the first acquisition module, and is used to determine, for a first type of accident, a first dose result set of C doses at D target distances and P target directions corresponding to the release of the source term of the first type of accident under each target meteorological sequence, so as to obtain K×n groups of the first dose result sets, wherein the first type of accident is any one of the accident set, and D, P, and C are all positive integers; A second determination module, connected to the first determination module, is used to determine a second dose result set of C doses at D target distances and P target orientations from K×n groups of the first dose result sets; a third determining module, connected to the second determining module, configured to determine, for each target distance among the D target distances in the second dose result set, a third dose result set of C doses corresponding to the target distance from the second dose result set according to dose result sets of C doses at P target orientations corresponding to the target distance in the second dose result set, so as to obtain a C×D group of the third dose result sets; The fourth determination module is connected to the third determination module and is used to determine the emergency plan area applicable to the first type of accident according to the third dose result set of the C×D group and the preset target dose result set.

16. The device according to claim 15, characterized in that The device also includes: The second acquisition module is used to acquire Q meteorological sequences of the target nuclear power plant site, where Q is a positive integer; The first classification module is connected to the second acquisition module and is used to classify Q meteorological sequences to obtain K target types, each target type includes multiple meteorological sequences, wherein the Q meteorological sequences are meteorological sequences for each hour throughout the year; K is less than Q.

17. The device according to claim 15, characterized in that The first classification module comprises: A selection submodule, connected to the second acquisition module, for selecting K meteorological sequence features from the Q meteorological sequence features corresponding to the Q meteorological sequences as K initial classification centers; A first determination submodule, connected to the selection submodule, is used to determine the first similarity between the remaining QK meteorological sequence features and each initial classification center respectively for the remaining QK meteorological sequence features in the Q meteorological sequence features; The first classification submodule is connected to the selection submodule and is used to classify the QK meteorological sequences into K initial classification centers according to the first similarity to obtain K target types.

Images