Sea area seismic activity parameter calculation method and system

Abstract

The invention provides a sea area seismic activity parameter calculation method and system, and is suitable for the technical field of data processing, and the method comprises the steps: obtaining seismic event information and geological structure information; generating seismic statistical region information according to the seismic event information, the geological structure information, preset neighborhood radius information and preset region position point quantity information; and calculating seismic activity parameters according to the seismic statistical region information, a preset weight coefficient matrix, a preset offset matrix, a preset space mapping function and a preset threshold value of mapping iterations. According to the method, the seismic event information and the geological structure information are analyzed and calculated, the seismic statistical regions are divided, the data of the plurality of seismic statistical regions are mapped to the high-dimensional space for depth feature extraction and prediction calculation, and the seismic activity parameters are further solved through the prediction calculation result, so that the seismic activity of the seismic data is improved. Therefore, the accuracy and the effectiveness of calculating the seismic activity parameters are ensured.

Description

Technical Field

[0001] The present application belongs to the field of data processing technology, and in particular, relates to a method and system for calculating sea area seismic activity parameters. Background Art

[0002] my country has a vast sea area, with both intraplate earthquakes and plate margin earthquakes, and frequent earthquake activities. With the development of national economic construction, the types and number of construction projects in sea and coastal areas have increased. Once an earthquake occurs, it will cause huge earthquake disaster losses to coastal areas and sea and coastal construction project areas. Therefore, it is necessary to conduct marine earthquake analysis and prepare and protect before marine earthquake disasters occur. In the marine earthquake analysis work, the calculation of seismic activity parameters is crucial. The calculated activity parameters can predict indicators such as the danger of seismic activity, thereby providing reliable data support for prevention and post-disaster reconstruction.

[0003] In the existing technology, it is usually necessary to rely on experts to manually judge the earthquake statistical area, and then use the hierarchical analysis method or the least squares method to fit the data of the land earthquake statistical area, so that the fitting results are used as seismic activity parameters for subsequent seismic activity analysis.

[0004] However, in the existing technology, large errors are easily generated when processing high-dimensional and complex seismic data, and the accuracy and reliability of seismic activity analysis and calculation cannot be ensured. Furthermore, the existing seismic zoning map only involves land areas, lacks seismic zoning for sea areas and models suitable for sea area seismic activity analysis, and cannot calculate sea area seismic activity parameters through existing models, so it is impossible to analyze sea area seismic activity. Therefore, when an earthquake disaster occurs in the sea area, the timeliness of pre-disaster prevention and post-disaster reconstruction will not be guaranteed, causing huge economic losses to coastal areas. Summary of the invention

[0005] In view of this, an embodiment of the present application provides a method and system for calculating sea area seismic activity parameters, aiming to solve the problem in the prior art that the calculation errors of sea area seismic activity parameters are large, the calculation results have low accuracy and reliability, and thus cannot be used for effective analysis of sea area seismic activity.

[0006] A first aspect of an embodiment of the present application provides a method for calculating sea area seismicity parameters, comprising: Obtain information on earthquake events and geological structures; Generate earthquake statistical regional information based on the earthquake event information, geological structure information, preset neighborhood radius information and preset regional location point quantity information; The seismic activity parameters are calculated according to the earthquake statistical area information, a preset weight coefficient matrix, a preset bias matrix, a preset spatial mapping function and a preset mapping iteration number threshold.

[0007] A second aspect of an embodiment of the present application provides a system for calculating sea area seismicity parameters, comprising: An information acquisition module is used to obtain earthquake event information and geological structure information; An earthquake statistical region information generating module, configured to generate earthquake statistical region information according to the earthquake event information, geological structure information, preset neighborhood radius information, and preset region location point quantity information; and The sea area seismic activity parameter calculation module is used to calculate the seismic activity parameters according to the seismic statistical area information, a preset weight coefficient matrix, a preset bias matrix, a preset spatial mapping function and a preset mapping iteration number threshold.

[0008] The third aspect of an embodiment of the present application provides a terminal device, which includes a memory and a processor, wherein the memory stores a computer program that can be run on the processor, and when the processor executes the computer program, it implements the steps of the method for calculating sea area seismic activity parameters as described in the first aspect above.

[0009] The fourth aspect of the embodiments of the present application provides a computer-readable storage medium, comprising: a computer program stored therein, which, when executed by a processor, implements the steps of the method for calculating sea area seismic activity parameters as described in the first aspect above.

[0010] Compared with the prior art, the embodiments of the present application have the following beneficial effects: by analyzing and calculating complex and large earthquake event information and geological structure information, dividing earthquake statistical areas by data feature analysis of multiple information, reducing data complexity, and then accurately capturing the complex relationship between earthquake influencing factors between regions by mapping the data of multiple earthquake statistical areas into a high-dimensional space for deep feature extraction and predictive calculation, and further solving the earthquake activity parameters through the predictive calculation results, so that the calculated earthquake activity parameters can adapt to different earthquake activity modes, so as to ensure the accuracy and effectiveness of the calculated earthquake activity parameters, thereby improving the accuracy and reliability of the analysis and calculation of the earthquake activity in the sea area, and providing scientific and effective data support for the planning and construction of marine facilities and coastal cities. BRIEF DESCRIPTION OF THE DRAWINGS

[0011] In order to more clearly illustrate the technical solutions in the embodiments of the present application, the drawings required for use in the embodiments or the description of the prior art will be briefly introduced below. Obviously, the drawings described below are only some embodiments of the present application. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying creative labor.

[0012] Figure 1 It is a schematic diagram of the implementation process of the method for calculating the sea area seismic activity parameters provided in the first embodiment of the present application; Figure 2 This is a schematic diagram of the implementation process of the method for calculating sea area seismic activity parameters provided in Example 2 of the present application; Figure 3 It is a schematic diagram of the implementation flow of the method for calculating sea area seismic activity parameters provided in Example 3 of the present application; Figure 4 This is a schematic diagram of the implementation process of the method for calculating sea area seismic activity parameters provided in Example 4 of the present application; Figure 5 This is a schematic diagram of the implementation process of the method for calculating sea area seismic activity parameters provided in Example 5 of the present application; Figure 6 This is a schematic diagram of the implementation process of the method for calculating sea area seismic activity parameters provided in Example 6 of the present application; Figure 7 It is a structural schematic diagram of a system for calculating sea area seismic activity parameters provided in an embodiment of the present application; Figure 8 It is a schematic diagram of a terminal device provided in an embodiment of the present application. DETAILED DESCRIPTION

[0013] In the following description, specific details such as specific system structures, technologies, etc. are provided for the purpose of illustration rather than limitation, so as to provide a thorough understanding of the embodiments of the present application. However, it should be clear to those skilled in the art that the present application may also be implemented in other embodiments without these specific details. In other cases, detailed descriptions of well-known systems, systems, circuits, and methods are omitted to prevent unnecessary details from obstructing the description of the present application.

[0014] In order to illustrate the technical solution described in this application, a specific embodiment is provided below for illustration.

[0015] Figure 1 The following is a flowchart of the method for calculating the sea area seismic activity parameters provided in the first embodiment of the present application, which is described in detail as follows: Step S101, obtaining earthquake event information and geological structure information.

[0016] In this embodiment, earthquake event information may be information such as the epicenter location, magnitude, and time of earthquake occurrence, which can reflect the spatial distribution, intensity, and time regularity of earthquake occurrence. It can be obtained through earthquake data institutions such as the Global Seismic Network Center and the China Earthquake Network Center; geological structure information may include information such as the location, length, strike, activity nature, and activity rate of faults, which can be obtained through field survey data or relevant research literature of geological survey departments, or can be determined by geological staff using methods such as address mapping or geophysical exploration and then input into a computer to complete the acquisition.

[0017] Step S102, generating earthquake statistical regional information according to the earthquake event information, geological structure information, preset neighborhood radius information and preset regional location point quantity information.

[0018] In this embodiment, the preset neighborhood radius information and the preset regional location point quantity information may be manually set, and may be calculation parameters when performing the division processing of earthquake statistical regions according to earthquake event information and geological structure information. The scope of each earthquake statistical region may be determined by the neighborhood radius information and the regional location point quantity information, and each earthquake statistical region may be obtained by clustering calculation. The earthquake statistical region may be used to further calculate and divide the potential source region.

[0019] Step S103, calculating seismicity parameters according to the earthquake statistical area information, a preset weight coefficient matrix, a preset bias matrix, a preset spatial mapping function and a preset mapping iteration number threshold.

[0020] In this embodiment, the preset weight coefficient matrix, the preset bias matrix, the preset space mapping function and the preset mapping iteration number threshold can be set manually, wherein the preset space mapping function can be designed based on an exponential function or a hyperbolic tangent function. It can be that each earthquake statistical area is weighted and summed by the weight coefficient matrix first, and then the weighted summation result is summed with the bias matrix, and then the summation result is used as the independent variable of the space mapping function, and the calculated function value is used as the intermediate variable, and then the intermediate variable is recalculated by the weight coefficient matrix and the bias matrix, and the calculated result is transformed and calculated as the independent variable of the space mapping function, so as to realize iterative calculation, and the number of iterations can be determined by the mapping iteration number threshold. The seismic activity parameter can be a magnitude frequency characterization parameter and an earthquake occurrence rate characterization parameter, wherein the magnitude frequency characterization parameter of the earthquake event can be a magnitude frequency characterization parameter in the Gutenberg-Richter law, which is used to describe the logarithmic linear relationship between the number of occurrences of earthquakes of a specific magnitude and the magnitude, and the earthquake occurrence rate characterization parameter can be used to represent the annual occurrence rate of earthquakes greater than or equal to a specific earthquake magnitude.

[0021] The method for calculating the parameters of marine seismic activity provided in the embodiment of the present application analyzes and calculates the complex and huge information on earthquake events and geological structures, divides the seismic statistical areas by performing data feature analysis on multiple information, reduces the complexity of data, and then accurately captures the complex relationship between earthquake influencing factors in various areas by mapping the data of multiple seismic statistical areas into a high-dimensional space for deep feature extraction and predictive calculation. The seismic activity parameters are further solved through the predictive calculation results, so that the calculated seismic activity parameters can adapt to different seismic activity modes to ensure the accuracy and effectiveness of the calculated seismic activity parameters, thereby improving the accuracy and reliability of the analysis and calculation of marine seismic activity, and providing scientific and effective data support for the planning and construction of marine facilities and coastal cities.

[0022] Figure 2 The flowchart of the method for calculating the sea area seismic activity parameters provided in the second embodiment of the present application is shown. The difference between the method and the first embodiment is as follows: The earthquake event information includes magnitude information and epicenter location information; The step S102 specifically includes: Step S201, normalizing the magnitude information to obtain magnitude characterization information.

[0023] In this embodiment, it can be understood that in earthquake research work, the earthquake magnitude information collected usually has different magnitude measurement standards and magnitude differences, and needs to be normalized to unify the dimensions of various magnitude information so that magnitude data from different sources and different types are comparable, ensuring the effectiveness of subsequent calculations. The normalization method can be linear normalization or Z-score normalization. The result after normalization is the magnitude characterization information.

[0024] Step S202, calculating the logical spatial distance of each of the magnitude characterization information to obtain the magnitude logical distance.

[0025] In this embodiment, the logical spatial distance may be the Euclidean distance. The magnitude logical distance is obtained by calculating the Euclidean distance between each magnitude characterization information. Step S203, calculating the logical spatial distance of each epicenter position information to obtain the position logical distance.

[0026] In this embodiment, the logical spatial distance may be the Euclidean distance, and the position logical distance is obtained by calculating the Euclidean distance between each epicenter position information.

[0027] Step S204, calculating the logical spatial distance of each piece of geological structure information to obtain the geological structure logical distance.

[0028] In this embodiment, the logical spatial distance may be the Euclidean distance. The geological structure logical distance is obtained by calculating the Euclidean distance between each geological structure information.

[0029] Step S205, calculating the comprehensive distance of earthquake information according to the magnitude logical distance, location logical distance, geological structure logical distance and preset weight coefficient.

[0030] In this embodiment, the preset weight coefficient may be manually set, and the magnitude logical distance, location logical distance, and geological structure logical distance may be weighted and summed by the preset weight coefficient, and the calculation result is the comprehensive distance of earthquake information.

[0031] Step S206, generating a plurality of earthquake characterization information points according to the magnitude characterization information, epicenter location information, geological structure information and comprehensive distance of earthquake information.

[0032] In this embodiment, the magnitude characterization information, epicenter location information, geological structure information and the comprehensive distance of earthquake information can be used as matrix elements, and multiple earthquake characterization information matrices are generated through multiple groups of magnitude characterization information, epicenter location information, geological structure information and the comprehensive distance of earthquake information. The multiple matrices are stored in a single logical storage space as multiple earthquake characterization information points. An earthquake characterization information point represents a matrix generated by the magnitude characterization information, epicenter location information, geological structure information and the comprehensive distance of earthquake information.

[0033] Step S207, determine whether the comprehensive distance of the earthquake information is less than or equal to the preset neighborhood radius information; if so, generate neighborhood set information of each earthquake characterization information point according to the earthquake characterization information point corresponding to the comprehensive distance of the earthquake information; if not, do not process the earthquake characterization information point corresponding to the comprehensive distance of the earthquake information, and do not generate neighborhood set information.

[0034] In this embodiment, one of the earthquake characterization information points may be randomly selected as a point in the main domain. When the earthquake information comprehensive distance is less than or equal to the preset neighborhood radius information, it means that the two earthquake characterization information points corresponding to the earthquake information comprehensive distance are close to each other. If one of the earthquake characterization information points is used as a point in the main domain, the other earthquake characterization information point is used as a point in the neighborhood, so that all points whose earthquake information comprehensive distances with the points in the main domain are less than or equal to the preset neighborhood radius information are counted to obtain all points in the neighborhood, and then the neighborhood set information of each earthquake characterization information point is generated based on all points in the neighborhood. When the earthquake information comprehensive distance is greater than the preset neighborhood radius information, it means that the two earthquake characterization information points corresponding to the earthquake information comprehensive distance are far from each other, and the earthquake characterization information point is not used as a point in the neighborhood, so the earthquake characterization information point corresponding to the earthquake information comprehensive distance is not processed, and the neighborhood set information is not generated.

[0035] Step S208, generating earthquake statistical area information according to the neighborhood set information and preset area location point quantity information.

[0036] In this embodiment, the preset regional location point quantity information can be manually set to constrain the scope of each earthquake statistical regional information. It can be divided into each earthquake statistical regional information according to the regional location point quantity information by counting the earthquake characterization information points in each neighborhood set information.

[0037] The method for calculating marine seismic activity parameters provided in the embodiment of the present application calculates the comprehensive logical distance by integrating earthquake magnitude information and geological structure information, so as to more comprehensively reflect the similarities and correlations between various earthquake points, and help to effectively extract the features between various earthquake data in the subsequent calculation process. By performing statistical analysis on the number of neighborhood points, areas with frequent seismic activities can be quickly located, thereby increasing the accuracy and effectiveness of the division of earthquake statistical areas, so as to improve the accuracy and robustness of the calculation of marine seismic activity parameters and provide more comprehensive information for the analysis of marine seismic activity.

[0038] Figure 3 The flowchart of the method for calculating the sea area seismic activity parameters provided in the third embodiment of the present application is shown. The difference between the method and the second embodiment is that the step S208 specifically includes: Step S301, counting the number of information points in each of the neighborhood set information to obtain neighborhood set point quantity information.

[0039] In this embodiment, the number of earthquake characterization information points in each neighborhood set information is counted, and this number is the neighborhood set point quantity information.

[0040] Step S302, determine whether the number information of the neighborhood collection points is greater than or equal to the preset number information of regional location points; if so, proceed to step S303; if not, proceed to step S304.

[0041] In this embodiment, the preset regional location point quantity information can be manually set. When the neighborhood set point quantity information is greater than or equal to the preset regional location point quantity information, it means that the earthquake characterization information points in the neighborhood set corresponding to the points in the main domain are strongly correlated with each other, and the points in the main domain are frequently seismically active, and the earthquake characterization information points in the neighborhood set need to be used as core points, and multiple core points constitute the core point information. When the neighborhood set point quantity information is less than the preset regional location point quantity information, it means that the earthquake points corresponding to the neighborhood set are weakly correlated with each other. When the neighborhood set point quantity information is equal to zero, it means that the points in the main domain corresponding to the neighborhood set are of little reference significance for seismic activity analysis, and the points in the main domain need to be discarded; when the neighborhood set point quantity information is not equal to zero, the earthquake characterization information points in the neighborhood set are used as boundary points, and multiple boundary points constitute the boundary point information.

[0042] Step S303: taking the earthquake characterization information points corresponding to the neighborhood set point quantity information as core point information.

[0043] In this embodiment, when the number of neighborhood set point information is greater than or equal to the preset number of regional location point information, it means that the earthquake characterization information points in the neighborhood set corresponding to the points in the main domain are strongly correlated with each other, and the seismic activity of the points in the main domain is frequent. It is necessary to take the earthquake characterization information points in the neighborhood set as core points, and multiple core points constitute the core point information.

[0044] Step S304, determining whether the neighborhood set point quantity information is not zero; if so, proceeding to step S305; if not, discarding the earthquake characterization information point corresponding to the neighborhood set quantity information.

[0045] In this embodiment, when the number information of neighborhood set points is less than the preset number information of regional location points, it means that the earthquake points corresponding to the neighborhood set are weakly correlated with each other. When the number information of neighborhood set points is equal to zero, it means that the points in the main domain corresponding to the neighborhood set have little reference significance for earthquake activity analysis, and the points in the main domain need to be discarded.

[0046] Step S305: taking the earthquake characterization information points corresponding to the neighborhood set point quantity information as boundary point information.

[0047] In this embodiment, when the number of neighborhood set points is not equal to zero, the earthquake characterization information points in the neighborhood set are used as boundary points, and multiple boundary points constitute boundary point information.

[0048] Step S306: Generate a seismic statistical regional information representation set based on the core point information, boundary point information and neighborhood set information.

[0049] In this embodiment, starting from any core point in the core point information, the neighborhood set of the core point can be determined, and then the core point can be diffused through the neighborhood set in a recursive manner until all core points are assigned to the neighborhood set and can no longer be diffused. The neighborhood set information finally generated is used as the information representation set of each earthquake statistical area.

[0050] Step S307: Generate earthquake statistical area information according to the earthquake statistical area information representation set.

[0051] In this embodiment, it can be understood that the earthquake statistical area information representation set includes multiple matrix elements, the values ​​of the matrix elements are extracted, and multiple matrices are regenerated according to the extracted values. The multiple matrices are the earthquake statistical area information, and the earthquake statistical area information includes earthquake event information and geological structure information of each earthquake statistical area.

[0052] The method for calculating the sea area seismic activity parameters provided in the embodiment of the present application performs set division based on the density of seismic data points, can adaptively identify seismic statistical areas of different shapes and sizes, and effectively filter out redundant data such as accidental small earthquakes or data measurement errors by determining core point information and boundary point information, thereby increasing the accuracy and reliability of seismic statistical area division and providing effective data support for the calculation of seismic activity parameters.

[0053] Figure 4 The flowchart of the method for calculating the sea area seismic activity parameters provided by the fourth embodiment of the present application is shown. The difference between the method and the first embodiment is that the step S103 specifically includes: Step S401, generating an earthquake statistical parameter matrix according to the earthquake statistical area information.

[0054] In this embodiment, each value in multiple earthquake statistical area information may be extracted to generate one-dimensional data, the dimension of the matrix is ​​set, and the one-dimensional array is converted into multiple earthquake statistical parameter matrices according to the dimension of the matrix.

[0055] Step S402, calculating a seismicity characterization variable matrix according to the earthquake statistical parameter matrix, a preset weight coefficient matrix and a preset bias matrix.

[0056] In this embodiment, the earthquake statistical parameter matrix may be first weighted and summed using a weight coefficient matrix, and then the weighted summation result is added to the bias matrix, and the addition result is used as the earthquake activity characterization variable matrix.

[0057] Step S403: Using the seismicity characterization variable matrix as an independent variable of a preset spatial mapping function to generate a seismicity characterization spatial mapping variable matrix.

[0058] In this embodiment, the seismic activity characterization variable matrix is ​​used as the independent variable of the preset space mapping function, and the function value is calculated as the seismic activity characterization space mapping variable matrix.

[0059] Step S404, calculating seismicity parameters according to the seismicity characterization space mapping variable matrix and a preset mapping iteration number threshold.

[0060] In this embodiment, the seismic activity characterization space mapping variable matrix can be calculated through the weight coefficient matrix and the bias matrix, and then the calculation result is used as the independent variable of the space mapping function to calculate the function value, thereby realizing one iteration. The number of iterations can be determined by the mapping iteration number threshold, and the results obtained after multiple iterative calculations are output as seismic activity parameters.

[0061] The method for calculating marine seismic activity parameters provided in the embodiment of the present application performs nonlinear spatial mapping transformation on seismic statistical area data through a preset spatial mapping function, performs offset transformation on seismic statistical area data in nonlinear space through a preset bias function, fully extracts and transforms nonlinear features in the seismic statistical area, increases the spatial distance between each nonlinear feature, and enables a variety of complex features in the seismic statistical area data to be extracted and analyzed, thereby avoiding the loss of effective feature information and improving the accuracy of calculating marine seismic activity parameters.

[0062] Figure 5 The flowchart of the method for calculating the sea area seismic activity parameters provided in the fifth embodiment of the present application is shown. The difference between the method and the fourth embodiment is that the step S404 specifically includes: Step S501, counting the number of times the seismicity characterization space mapping variable matrix is ​​generated to obtain the number of mapping iterations.

[0063] In this embodiment, each time the spatial mapping variable matrix for characterizing seismic activity is generated, it indicates that the spatial mapping transformation has been performed once. Therefore, the number of times the spatial mapping variable matrix for characterizing seismic activity is generated can be counted to calculate the number of mapping iterations to determine whether the spatial mapping transformation calculation needs to be terminated.

[0064] Step S502, determining whether the mapping iteration number is less than a preset mapping iteration number threshold; if so, proceeding to step S503; if not, proceeding to step S504.

[0065] In this embodiment, the preset mapping iteration number threshold can be manually set. When the mapping iteration number is less than the preset mapping iteration number threshold, it means that the feature extraction and analysis of the seismic statistical area data is still insufficient, and it is necessary to continue iterative calculation to ensure that the seismic statistical area data is deeply feature extracted and analyzed; when the mapping iteration number is greater than or equal to the preset mapping iteration number threshold, it means that the feature extraction and analysis of the seismic statistical area data is sufficient, and it is not necessary to continue iterative calculation. The seismic statistical area data has been deeply feature extracted and analyzed and can be used to output the calculation results.

[0066] Step S503, using the seismic activity characterization space mapping variable matrix as the seismic statistical parameter matrix, and returning to step S402.

[0067] In this embodiment, when the number of mapping iterations is less than a preset mapping iteration threshold, it means that the feature extraction and analysis of the earthquake statistical area data is still insufficient, and iterative calculation needs to be continued to ensure that the earthquake statistical area data is deeply feature extracted and analyzed.

[0068] Step S504, extracting the magnitude characterization information and earthquake occurrence frequency information of the seismic activity characterization space mapping variable matrix.

[0069] In this embodiment, when the number of mapping iterations is greater than or equal to the preset mapping iteration threshold, it means that the feature extraction and analysis of the earthquake statistical area data are sufficient, and there is no need to continue iterative calculation. The earthquake statistical area data has been deeply feature extracted and analyzed and can be used to output calculation results. At this time, it is necessary to extract the magnitude characterization information and the number of earthquake occurrence information in the earthquake activity characterization space mapping variable matrix for subsequent calculation of earthquake activity parameters.

[0070] Step S505, calculating seismic activity parameters according to the magnitude characterization information and the earthquake occurrence frequency information.

[0071] In this embodiment, the seismic activity parameter may be a magnitude frequency characterization parameter and an earthquake occurrence rate, which are used to quantitatively describe the seismic activity characteristics. It may be based on the Gutenberg-Rickert law, and the magnitude characterization information may be fitted by the analytic hierarchy process or the least squares method, and the fitting result may be output as a magnitude frequency characterization parameter. It may be that a plurality of magnitudes are divided into different magnitude intervals, and then the number of earthquake occurrences in each magnitude interval is counted, and then the earthquake occurrence rate in each magnitude interval is calculated.

[0072] The method for calculating marine seismic activity parameters provided in the embodiment of the present application maps the seismic statistical area information into a nonlinear space to fully extract and analyze the data features, fully capturing the complex nonlinear relationship between multiple factors such as magnitude, time, and space in the seismic statistical area information, thereby more accurately fitting the relationship between the magnitude frequency and earthquake occurrence rate and various influencing factors, so as to adapt to the differences in geological conditions and seismic activity characteristics in different seismic statistical areas, so that the calculation results are more in line with the actual laws of seismic activity, and ensure that the calculation of marine seismic activity parameters provides accuracy and comprehensiveness.

[0073] Figure 6 The flowchart of the method for calculating the sea area seismic activity parameters provided in the sixth embodiment of the present application is shown, which differs from the first embodiment described above in that: The seismic activity parameters include magnitude frequency characterization parameters and earthquake occurrence rate characterization parameters; After step S103, the method further includes: Step S601, obtaining the distribution information of potential earthquake source marine facilities, the population distribution information of potential earthquake source area and the surrounding environment information of potential earthquake source.

[0074] In this embodiment, the distribution information of potential earthquake source marine facilities can be the geographical location, quantity, type, scale, function and other information of offshore drilling platforms, submarine cables, marine observation stations, port terminals, offshore wind farms, etc. The population distribution information of potential earthquake source regions can be the number, density and spatial distribution characteristics of the population in the earthquake statistical area and its surrounding areas, and can include the population distribution of cities and villages, the distribution of people of different ages and occupations, etc. The surrounding environmental information of potential earthquake sources can be the natural environment and ecosystem characteristics of the surrounding earthquake statistical area involved in the surrounding environmental information of potential earthquake sources. The natural environment can include topography (such as mountains, plains, coastlines, etc.), geological structure, meteorological conditions, etc.; the ecosystem can include marine ecology (such as coral reefs, fishing grounds, etc.), terrestrial ecology (such as forests, wetlands, etc.). The above information can be retrieved on the Internet. Among them, the potential earthquake source area can be obtained by further dividing the earthquake statistical area.

[0075] Step S602, calculating earthquake risk coefficient information based on the magnitude frequency characterization parameters, earthquake occurrence rate characterization parameters, potential earthquake source marine facility distribution information, potential earthquake source area population distribution information and potential earthquake source surrounding environment information.

[0076] In this embodiment, the various types of information collected can be quantified. For example, give corresponding values ​​and weights to marine facilities of different types and importance; convert population distribution information into quantitative indicators such as population density; score the ecological fragility and terrain impact in the surrounding environment information. In order to eliminate the magnitude differences between different data, the quantified data needs to be standardized, such as using the Min-Max normalized Z-score normalization method. After normalization, the comprehensive influencing factors of magnitude frequency and earthquake occurrence rate, marine facility distribution influencing factors, population distribution influencing factors and surrounding environment influencing factors are obtained. The earthquake risk coefficient can be obtained by summing the comprehensive influencing factors of magnitude frequency and earthquake occurrence rate, marine facility distribution influencing factors, population distribution influencing factors and surrounding environment influencing factors.

[0077] Step S603, determining earthquake risk response level information according to the earthquake risk factor information and a preset risk factor calibration interval.

[0078] In this embodiment, the preset risk factor calibration interval can be set manually, and different risk factor calibration intervals can be set with reference to specific earthquake magnitude classification standards, and each interval corresponds to a specific earthquake risk response level, such as low risk, medium risk, high risk and extremely high risk. The calculated earthquake risk coefficient can be compared with the preset calibration intervals one by one. If the earthquake risk coefficient is in the interval corresponding to the low risk, then the earthquake risk response level is determined to be a low risk level; when the risk coefficient falls in the medium risk interval, the response level is medium risk; if it is in the high risk interval, the response level is high risk.

[0079] Step S604, generating a plurality of earthquake emergency rescue information according to the earthquake risk response level information, the preset earthquake risk level information and the preset emergency response information.

[0080] In this embodiment, the earthquake risk response level can be matched with the preset earthquake risk level. If it is matched to a low risk level, the corresponding rescue information can be extracted from the preset emergency response information. The earthquake emergency rescue information can include the deployment of a small number of rescue personnel, the storage of basic materials, etc.; for medium risk levels, more resources can be invested, and the earthquake emergency rescue information can include information such as increasing the rescue team and replenishing materials; when at a high risk or extremely high risk level, the emergency response information can be more comprehensive and more powerful, and the earthquake emergency rescue information can cover large-scale personnel rescue, professional equipment investment, medical insurance, etc.

[0081] The method for calculating sea earthquake activity parameters provided in the embodiment of the present application comprehensively considers the distribution of potential earthquake source marine facilities, regional population distribution and surrounding environmental information to calculate the risk coefficient, comprehensively and accurately quantifies the impact of sea earthquakes on marine facilities, personnel and the surrounding environment, and quickly and accurately generates multiple emergency rescue information that is adapted to the current earthquake risk situation based on the precise risk coefficient, making the rescue strategy more targeted, providing a solid guarantee for the efficient implementation of rescue operations in the sea area after the earthquake disaster, and comprehensively improving the ability to respond to sea earthquake disasters.

[0082] Corresponding to the method of the above embodiment, Figure 7 A structural block diagram of a system for calculating sea area seismic activity parameters provided in an embodiment of the present application is shown. For ease of explanation, only the parts related to the embodiment of the present application are shown. Figure 7 The exemplary marine seismicity parameter calculation system may be an execution body of the marine seismicity parameter calculation method provided in the aforementioned first embodiment.

[0083] Reference Figure 7 , the sea area seismic activity parameter calculation system includes: Information acquisition module 710, used to acquire earthquake event information and geological structure information; An earthquake statistical region information generating module 720 is used to generate earthquake statistical region information according to the earthquake event information, geological structure information, preset neighborhood radius information and preset region location point quantity information; and The marine seismicity parameter calculation module 730 is used to calculate the seismicity parameters according to the seismic statistical area information, a preset weight coefficient matrix, a preset bias matrix, a preset spatial mapping function and a preset mapping iteration number threshold.

[0084] The process of each module realizing its own function in the marine seismic activity parameter calculation system provided in the embodiment of the present application can be specifically referred to the aforementioned Figure 1 The description of the first embodiment is not repeated here.

[0085] It should be understood that the size of the serial numbers of the steps in the above embodiments does not mean the order of execution. The execution order of each process should be determined by its function and internal logic, and should not constitute any limitation on the implementation process of the embodiments of the present application.

[0086] It should be understood that when used in the present specification and the appended claims, the term "comprising" indicates the presence of described features, wholes, steps, operations, elements and / or components, but does not exclude the presence or addition of one or more other features, wholes, steps, operations, elements, components and / or combinations thereof.

[0087] It should also be understood that the term “and / or” used in the specification and appended claims refers to any and all possible combinations of one or more of the associated listed items, and includes these combinations.

[0088] As used in the specification and appended claims of this application, the term "if" can be interpreted as "when" or "uponce" or "in response to determining" or "in response to detecting", depending on the context. Similarly, the phrase "if it is determined" or "if [described condition or event] is detected" can be interpreted as meaning "uponce it is determined" or "in response to determining" or "uponce [described condition or event] is detected" or "in response to detecting [described condition or event]", depending on the context.

[0089] In addition, in the description of the present specification and the appended claims, the terms "first", "second", "third", etc. are only used to distinguish descriptions, and cannot be understood as indicating or suggesting relative importance. It should also be understood that although the terms "first", "second", etc. are used to describe various elements in some embodiments of the present application in the text, these elements should not be limited by these terms. These terms are only used to distinguish one element from another element. For example, the first table can be named as the second table, and similarly, the second table can be named as the first table without departing from the scope of the various described embodiments. The first table and the second table are both tables, but they are not the same table.

[0090] References to "one embodiment" or "some embodiments" etc. described in the specification of this application mean that one or more embodiments of the present application include specific features, structures or characteristics described in conjunction with the embodiment. Therefore, the statements "in one embodiment", "in some embodiments", "in some other embodiments", "in some other embodiments", etc. that appear in different places in this specification do not necessarily refer to the same embodiment, but mean "one or more but not all embodiments", unless otherwise specifically emphasized in other ways. The terms "including", "comprising", "having" and their variations all mean "including but not limited to", unless otherwise specifically emphasized in other ways.

[0091] The method for calculating sea area seismic activity parameters provided in the embodiments of the present application can be applied to terminal devices such as mobile phones, tablet computers, wearable devices, vehicle-mounted devices, augmented reality (AR) / virtual reality (VR) devices, laptop computers, ultra-mobile personal computers (UMPC), netbooks, personal digital assistants (PDA), etc. The embodiments of the present application do not impose any restrictions on the specific types of terminal devices.

[0092] For example, the terminal device can be a station (STAION, ST) in a WLAN, a cellular phone, a cordless phone, a Session Initiation Protocol (SIP) phone, a Wireless Local Loop (WLL) station, a Personal Digital Assistant (PDA) device, a handheld device with wireless communication function, a computing device or other processing device connected to a wireless modem, a vehicle-mounted device, a vehicle networking terminal, a computer, a laptop computer, a handheld communication device, a handheld computing device, a satellite wireless device, a wireless modem card, a TV set top box (STB), a customer premises equipment (CPE) and / or other devices for communicating on a wireless system and a next-generation communication system, such as a mobile terminal in a 5G network or a mobile terminal in a future evolved Public Land Mobile Network (PLMN) network, etc.

[0093] As an example but not limitation, when the terminal device is a wearable device, the wearable device can also be a general term for wearable devices that are intelligently designed and developed using wearable technology for daily wear, such as glasses, gloves, watches, clothing and shoes. A wearable device is a portable device that is worn directly on the body or integrated into the user's clothes or accessories. Wearable devices are not just hardware devices, but also achieve powerful functions through software support, data interaction, and cloud interaction. Broadly speaking, wearable smart devices include those that are fully functional, large in size, and can achieve complete or partial functions without relying on smartphones, such as smart watches or smart glasses, as well as those that only focus on a certain type of application function and need to be used in conjunction with other devices such as smartphones, such as various smart bracelets and smart jewelry for vital sign monitoring.

[0094] Figure 8Schematic diagram of the structure of a terminal device provided by an embodiment of the present application. Figure 8 As shown, the terminal device 8 of this embodiment includes: at least one processor 80 ( Figure 8 Only one is shown in the figure), a memory 81, wherein the memory 81 stores a computer program 82 that can be run on the processor 80. When the processor 80 executes the computer program 82, the steps in the above-mentioned embodiments of the method for calculating the parameters of sea area seismicity are implemented, such as Figure 1 Alternatively, when the processor 80 executes the computer program 82, the functions of each module / unit in the above-mentioned system embodiments are realized, for example Figure 7 Functions of modules 710 to 730 are shown.

[0095] The terminal device 8 may be a computing device such as a desktop computer, a notebook, a PDA, or a cloud server. The terminal device may include, but is not limited to, a processor 80 and a memory 81. Those skilled in the art will appreciate that Figure 8 It is only an example of the terminal device 8 and does not constitute a limitation of the terminal device 8. It may include more or fewer components than shown in the figure, or a combination of certain components, or different components. For example, the terminal device may also include an input sending device, a network access device, a bus, etc.

[0096] The processor 80 may be a central processing unit (CPU), or other general-purpose processors, digital signal processors (DSP), application-specific integrated circuits (ASIC), field-programmable gate arrays (FPGA) or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, etc. A general-purpose processor may be a microprocessor or any conventional processor, etc.

[0097] In some embodiments, the memory 81 may be an internal storage unit of the terminal device 8, such as a hard disk or memory of the terminal device 8. The memory 81 may also be an external storage device of the terminal device 8, such as a plug-in hard disk, a smart media card (SMC), a secure digital (SD) card, a flash card, etc. equipped on the terminal device 8. Further, the memory 81 may also include both an internal storage unit and an external storage device of the terminal device 8. The memory 81 is used to store an operating system, an application program, a boot loader (BootLoader), data, and other programs, such as the program code of the computer program, etc. The memory 81 may also be used to temporarily store data that has been sent or is to be sent.

[0098] In addition, each functional unit in each embodiment of the present application may be integrated into one processing unit, or each unit may exist physically separately, or two or more units may be integrated into one unit. The above-mentioned integrated unit may be implemented in the form of hardware or in the form of software functional units.

[0099] An embodiment of the present application also provides a terminal device, which includes at least one memory, at least one processor, and a computer program stored in the at least one memory and executable on the at least one processor, wherein when the processor executes the computer program, the terminal device implements the steps in any of the above-mentioned method embodiments.

[0100] An embodiment of the present application further provides a computer-readable storage medium, wherein the computer-readable storage medium stores a computer program, and when the computer program is executed by a processor, the steps in the above-mentioned method embodiments can be implemented.

[0101] An embodiment of the present application provides a computer program product. When the computer program product runs on a terminal device, the terminal device can implement the steps in the above-mentioned method embodiments when executing the computer program product.

[0102] If the integrated module / unit is implemented in the form of a software functional unit and sold or used as an independent product, it can be stored in a computer-readable storage medium. Based on this understanding, the present application implements all or part of the processes in the above-mentioned embodiment method, and can also be completed by instructing the relevant hardware through a computer program. The computer program can be stored in a computer-readable storage medium, and the computer program can implement the steps of the above-mentioned various method embodiments when executed by the processor. Among them, the computer program includes computer program code, and the computer program code can be in source code form, object code form, executable file or some intermediate form. The computer-readable medium may include: any entity or device that can carry the computer program code, recording medium, U disk, mobile hard disk, disk, optical disk, computer memory, read-only memory (ROM), random access memory (RAM), electrical carrier signal, telecommunication signal and software distribution medium, etc.

[0103] In the above embodiments, the description of each embodiment has its own emphasis. For parts that are not described or recorded in detail in a certain embodiment, reference can be made to the relevant descriptions of other embodiments.

[0104] Those of ordinary skill in the art will appreciate that the units and algorithm steps of each example described in conjunction with the embodiments disclosed herein can be implemented in electronic hardware, or a combination of computer software and electronic hardware. Whether these functions are performed in hardware or software depends on the specific application and design constraints of the technical solution. Professional and technical personnel can use different methods to implement the described functions for each specific application, but such implementation should not be considered to be beyond the scope of this application.

[0105] The units described as separate components may or may not be physically separated, and the components shown as units may or may not be physical units, that is, they may be located in one place or distributed on multiple network units. Some or all of the units may be selected according to actual needs to achieve the purpose of the solution of this embodiment.

[0106] The embodiments described above are only used to illustrate the technical solutions of the present application, rather than to limit them. Although the present application has been described in detail with reference to the aforementioned embodiments, a person skilled in the art should understand that the technical solutions described in the aforementioned embodiments may still be modified, or some of the technical features may be replaced by equivalents. Such modifications or replacements do not deviate the essence of the corresponding technical solutions from the spirit and scope of the technical solutions of the embodiments of the present application, and should all be included in the protection scope of the present application.

Claims

1. A method for calculating sea area seismicity parameters, characterized in that: include: Obtain information on earthquake events and geological structures; Generate earthquake statistical regional information based on the earthquake event information, geological structure information, preset neighborhood radius information and preset regional location point quantity information; The seismic activity parameters are calculated according to the earthquake statistical area information, a preset weight coefficient matrix, a preset bias matrix, a preset spatial mapping function and a preset mapping iteration number threshold.

2. The method for calculating sea area seismic activity parameters according to claim 1, characterized in that: The earthquake event information includes magnitude information and epicenter location information; The step of generating earthquake statistical regional information according to the earthquake event information, geological structure information, preset neighborhood radius information and preset regional location point quantity information specifically includes: Normalizing the magnitude information to obtain magnitude representation information; Calculate the logical spatial distance of each magnitude characterization information to obtain the magnitude logical distance; Calculate the logical spatial distance of each epicenter location information to obtain the location logical distance; Calculating the logical spatial distance of each of the geological structure information to obtain the geological structure logical distance; Calculate the comprehensive distance of earthquake information according to the magnitude logical distance, location logical distance, geological structure logical distance and preset weight coefficient; Generate multiple earthquake characterization information points according to the magnitude characterization information, epicenter location information, geological structure information and earthquake information comprehensive distance; When the earthquake information comprehensive distance is less than or equal to the preset neighborhood radius information, generating neighborhood set information of each earthquake characterization information point according to the earthquake characterization information point corresponding to the earthquake information comprehensive distance; Earthquake statistical area information is generated based on the neighborhood set information and preset area location point quantity information.

3. The method for calculating sea area seismic activity parameters according to claim 2, characterized in that: The step of generating earthquake statistical regional information according to the neighborhood set information and the preset regional location point quantity information specifically includes: Counting the number of information points in each of the neighborhood set information to obtain neighborhood set point quantity information; Determine whether the number of neighborhood set points is greater than or equal to the preset number of regional location points; If yes, the earthquake characterization information point corresponding to the neighborhood set point quantity information is used as the core point information; If not, when the number of neighborhood set points is not zero, the earthquake characterization information point corresponding to the number of neighborhood set points is used as the boundary point information; Generate a seismic statistical regional information representation set according to the core point information, boundary point information and neighborhood set information; The seismic statistical area information is generated according to the seismic statistical area information representation set.

4. The method for calculating sea area seismic activity parameters according to claim 1, characterized in that: The step of calculating the seismicity parameter according to the earthquake statistical area information, the preset weight coefficient matrix, the preset bias matrix, the preset spatial mapping function and the preset mapping iteration number threshold specifically includes: Generate an earthquake statistical parameter matrix according to the earthquake statistical area information; Calculating a seismicity characterization variable matrix according to the seismic statistical parameter matrix, a preset weight coefficient matrix, and a preset bias matrix; Using the seismicity characterization variable matrix as an independent variable of a preset spatial mapping function to generate a seismicity characterization spatial mapping variable matrix; The seismicity parameter is calculated according to the seismicity characterization space mapping variable matrix and a preset mapping iteration number threshold.

5. The method for calculating sea area seismic activity parameters according to claim 4, characterized in that: The step of calculating the seismicity parameter according to the seismicity characterization space mapping variable matrix and a preset mapping iteration number threshold specifically includes: Counting the number of times the seismic activity characterization space mapping variable matrix is ​​generated to obtain the number of mapping iterations; Determining whether the mapping iteration number is less than a preset mapping iteration number threshold; If yes, the seismic activity characterization space mapping variable matrix is ​​used as the seismic statistical parameter matrix, and the process returns to the step of calculating the seismic activity characterization variable matrix according to the seismic statistical parameter matrix, the preset weight coefficient matrix and the preset bias matrix; If not, extracting the magnitude characterization information and earthquake occurrence number information of the seismic activity characterization space mapping variable matrix; The seismic activity parameters are calculated based on the magnitude characterization information and the earthquake occurrence frequency information.

6. The method for calculating sea area seismic activity parameters according to claim 1, characterized in that: The seismic activity parameters include magnitude frequency characterization parameters and earthquake occurrence rate characterization parameters.

7. The method for calculating sea area seismic activity parameters according to claim 6, characterized in that: After the step of calculating the seismicity parameters according to the earthquake statistical area information, the preset weight coefficient matrix, the preset bias matrix, the preset spatial mapping function and the preset mapping iteration number threshold, the method further includes: Obtain information on the distribution of marine facilities at potential earthquake sources, population distribution in potential earthquake source areas, and the surrounding environment of potential earthquake sources; Calculate earthquake risk coefficient information based on the magnitude frequency characterization parameters, earthquake occurrence rate characterization parameters, potential earthquake source marine facility distribution information, potential earthquake source regional population distribution information, and potential earthquake source surrounding environment information; Determining earthquake risk response level information according to the earthquake risk factor information and a preset risk factor calibration interval; A plurality of earthquake emergency rescue information are generated according to the earthquake risk response level information, the preset earthquake risk level information and the preset emergency response information.

8. A system for calculating parameters of marine seismic activity, characterized in that: include: An information acquisition module is used to obtain earthquake event information and geological structure information; An earthquake statistical area information generation module, used to generate earthquake statistical area information according to the earthquake event information, geological structure information, preset neighborhood radius information and preset area location point quantity information; as well as The sea area seismic activity parameter calculation module is used to calculate the seismic activity parameters according to the seismic statistical area information, a preset weight coefficient matrix, a preset bias matrix, a preset spatial mapping function and a preset mapping iteration number threshold.

9. A terminal device, characterized in that: The terminal device includes a memory and a processor, wherein the memory stores a computer program that can be run on the processor, and when the processor executes the computer program, the steps of the method according to any one of claims 1 to 7 are implemented.

10. A computer-readable storage medium storing a computer program, characterized in that: When the computer program is executed by a processor, the steps of the method according to any one of claims 1 to 7 are implemented.

Images