Method and apparatus for constructing three-dimensional model of seismogenic fault based on aftershock relocation data

Abstract

The present invention relates to a method and apparatus for constructing a three-dimensional model of a seismogenic fault based on aftershock relocation data, and relates to the technical field of model construction. The method includes: obtaining real-time aftershock relocation data; retrieving the previous aftershock relocation data based on the real-time aftershock relocation data; determining the influence value of aftershock data according to the previous aftershock relocation data; determining the adjusted aftershock relocation data based on the influence value of aftershock data and the real-time aftershock relocation data; determining the final three-dimensional model of the seismogenic fault according to the adjusted aftershock relocation data, and outputting the final three-dimensional model of the seismogenic fault. The present invention has the effect of improving the accuracy of the constructed three-dimensional model.

Description

Technical Field

[0001] The present invention relates to the technical field of model construction, and in particular to a method and device for constructing a three-dimensional model of a seismogenic fault based on aftershock relocation data. Background Art

[0002] A seismogenic fault refers to a fault formed when, under the drive of factors such as plate movement, the crustal rocks are subjected to strong pressure, tension or shear force exceeding the strength limit of the rocks, and the stress accumulated over a long time is suddenly released at the fault, triggering an earthquake and causing the rocks to rupture.

[0003] In order to accurately evaluate the risk of earthquake disasters and better carry out prevention, it is generally necessary to construct a three-dimensional model of the seismogenic fault, so as to understand the spatial distribution, geometric shape and shallow and deep structures of the seismogenic fault, which is convenient for scientific research and prevention guidance by scientific researchers. At present, it is generally adopted to obtain aftershock relocation data, determine the two-dimensional interpreted fault line through three-dimensional slicing and perform spatial interpolation to generate an initial three-dimensional model, and adjust the seismogenic fault plane of the initial three-dimensional model to be close to the middle of the specified large aftershocks, so that the distances from multiple specified large aftershocks to the seismogenic fault plane are within a preset distance, thereby obtaining the final three-dimensional model.

[0004] Since when constructing a three-dimensional model of a seismogenic fault, it is necessary to obtain aftershock relocation data, and after an earthquake, the geological structure has changed, resulting in certain errors in the collected aftershock relocation data, and further reducing the accuracy of the finally constructed three-dimensional model. Summary of the Invention

[0005] In order to improve the accuracy of the constructed three-dimensional model, the present invention provides a method and device for constructing a three-dimensional model of a seismogenic fault based on aftershock relocation data.

[0006] In a first aspect, the present invention provides a method for constructing a three-dimensional model of a seismogenic fault based on aftershock relocation data, adopting the following technical solution:

[0007] A method for constructing a three-dimensional model of a seismogenic fault based on aftershock relocation data includes:

[0008] Obtaining real-time aftershock relocation data;

[0009] Retrieving the previous aftershock relocation data based on the real-time aftershock relocation data;

[0010] Determining the aftershock data influence value according to the previous aftershock relocation data;

[0011] Determining the aftershock relocation adjustment data according to the aftershock data influence value and the real-time aftershock relocation data;

[0012] Determine the final three-dimensional model of the seismogenic fault according to the adjusted data of aftershock relocation, and output the final three-dimensional model of the seismogenic fault.

[0013] Optionally, the method for determining the influence value of aftershock data includes:

[0014] Retrieve the previous aftershock position points based on the previous aftershock relocation data;

[0015] Calculate the distance value between the previous aftershock position points and the preset data acquisition position points and use it as the previous position distance value;

[0016] Determine the influence value of the previous position distance corresponding to the previous position distance value according to the corresponding relationship between the previous position distance value and the preset influence value of the previous position distance;

[0017] Retrieve the surrounding environment information around the previous position based on the previous aftershock position points;

[0018] Determine the environmental change information according to the surrounding environment information of the previous position;

[0019] Determine the influence value of environmental change according to the environmental change information;

[0020] Calculate the sum value between the influence value of the previous position distance and the influence value of environmental change and use it as the influence value of aftershock data.

[0021] Optionally, the method for determining the influence value of environmental change includes:

[0022] Retrieve the environmental type information based on the surrounding environment information of the previous position;

[0023] Determine the corresponding type area reference value for the environmental type information according to the corresponding relationship between the environmental type information and the preset type area reference value;

[0024] Retrieve the changed area value based on the environmental change information;

[0025] Determine whether the changed area value is greater than the type area reference value;

[0026] If it is, calculate the difference between the changed area value and the type area reference value and use it as the area deviation value;

[0027] Determine the area deviation influence value according to the area deviation value and the environmental type information, and use the area deviation influence value as the influence value of environmental change;

[0028] If it is not, determine the changed area influence value according to the previous aftershock relocation data and the changed area value, and use the changed area influence value as the influence value of environmental change.

[0029] Optionally, the method for determining the area deviation influence value includes:

[0030] Obtain the current time point;

[0031] Relocate the last data based on aftershocks and retrieve the last data time point;

[0032] Calculate the time period between the last data time point and the current time point and use it as the elapsed time value of the last data;

[0033] According to the correspondence between the environmental type information and the preset type time change reference value, determine the type time change reference value corresponding to the environmental type information;

[0034] Calculate the product value between the type time change reference value and the elapsed time value of the last data and use it as the type time change estimated value;

[0035] Determine the time change influence value based on the area deviation value and the type time change estimated value, and use the time change influence value as the area deviation influence value.

[0036] Optionally, the method for determining the time change influence value includes:

[0037] According to the correspondence between the environmental type information and the preset type deviation reference influence value, determine the type deviation reference influence value corresponding to the environmental type information;

[0038] Determine whether the area deviation value is less than the type time change estimated value;

[0039] If yes, calculate the product value between the area deviation value and the type deviation reference influence value and use it as the time change influence value;

[0040] If no, calculate the difference between the area deviation value and the type time change estimated value and use it as the deviation abnormal value;

[0041] According to the correspondence between the deviation abnormal value and the preset abnormal deviation influence value, determine the abnormal deviation influence value corresponding to the deviation abnormal value;

[0042] According to the correspondence between the type time change estimated value and the preset change estimation reference influence value, determine the change estimation reference influence value corresponding to the type time change estimated value;

[0043] Calculate the sum value between the change estimation reference influence value and the abnormal deviation influence value and use it as the time change influence value.

[0044] Optionally, the method for determining the changed area influence value includes:

[0045] Relocate the last data based on aftershocks and retrieve the last data time point;

[0046] Determine the last location weather information according to the last data time point and the aftershock last location point;

[0047] Determine the changing weather reference area value and the changing weather impact value based on the weather information at the last location;

[0048] Based on the corresponding relationship between the changing area value and the preset initial impact value of the changing area, determine the initial impact value of the changing area corresponding to the changing area value;

[0049] Determine whether the changing area value is less than the changing weather reference area value;

[0050] If it is, use the initial impact value of the changing area as the impact value of the changing area;

[0051] If it is not, calculate the difference between the changing area value and the changing weather reference area value and use it as the area anomaly value;

[0052] Calculate the product value between the area anomaly value and the changing weather impact value and use it as the area anomaly impact value;

[0053] Calculate the sum value between the area anomaly impact value and the initial impact value of the changing area and use it as the adjusted impact value of the changing area, and use the adjusted impact value of the changing area as the impact value of the changing area.

[0054] Optionally, the method for determining the changing weather reference area value and the changing weather impact value includes:

[0055] Retrieve the weather type information and the weather detection value based on the weather information at the last location;

[0056] Based on the corresponding relationship between the weather type information and the preset weather type reference area value, determine the weather type reference area value corresponding to the weather type information;

[0057] Calculate the product value between the weather detection value and the weather type reference area value and use it as the changing weather reference area value;

[0058] Based on the corresponding relationship between the weather type information and the preset weather type reference interval, determine the weather type reference interval corresponding to the weather type information;

[0059] Determine whether the weather detection value is within the weather type reference interval;

[0060] If it is, based on the corresponding relationship between the weather detection value and the preset weather detection impact value, determine the weather detection impact value corresponding to the weather detection value, and use the weather detection impact value as the changing weather impact value;

[0061] If it is not, determine the weather detection anomaly impact value based on the weather detection value and the weather type reference interval, and use the weather detection anomaly impact value as the changing weather impact value.

[0062] Optionally, the method for determining the weather detection anomaly influence value includes:

[0063] Calculate the difference between the weather detection value and the weather type reference interval and use it as the interval end deviation value;

[0064] Perform curve fitting based on the interval end deviation value to form an end deviation curve;

[0065] Analyze the deviation between the end deviation curve and the preset end deviation reference curve to form curve deviation information;

[0066] Retrieve the curve deviation position points and curve deviation values based on the curve deviation information;

[0067] Calculate the distance value between two adjacent curve deviation position points and use it as the curve adjacent deviation distance value;

[0068] According to the correspondence between the curve adjacent deviation distance value and the preset deviation distance unit influence value, determine the deviation distance unit influence value corresponding to the curve adjacent deviation distance value;

[0069] Calculate the product value between the deviation distance unit influence value and the curve deviation value and use it as the weather detection anomaly influence value.

[0070] In a second aspect, the present invention provides a device for constructing a three-dimensional model of a seismogenic fault based on aftershock relocation data, adopting the following technical solution:

[0071] A device for constructing a three-dimensional model of a seismogenic fault based on aftershock relocation data, comprising:

[0072] An acquisition module for acquiring real-time aftershock relocation data and the current time point;

[0073] A memory for storing the method for constructing a three-dimensional model of a seismogenic fault based on aftershock relocation data as described in any one of the first aspect;

[0074] A processor for loading and executing the program in the memory.

[0075] In summary, the present invention includes at least one of the following beneficial technical effects:

[0076] 1. By acquiring real-time aftershock relocation data and retrieving the previous aftershock relocation data to determine the aftershock data influence value, determining the aftershock relocation adjustment data through the aftershock data influence value and the real-time aftershock relocation data, and determining and outputting the final three-dimensional model of the seismogenic fault through the aftershock relocation adjustment data, the obtained three-dimensional model is affected by the previous aftershocks, thereby improving the accuracy of the constructed three-dimensional model;

[0077] 2. Retrieve the previous aftershock location points from the previous data through aftershock relocation, calculate the distance value from the previous location, then query and determine the influence value of the distance from the previous location through the distance value from the previous location. Retrieve the environmental information around the previous location through the previous aftershock location points to determine the environmental change information, determine the influence value of the environmental change through the environmental change information, and calculate the sum of the influence value of the distance from the previous location and the influence value of the environmental change as the influence value of the aftershock data, thereby improving the accuracy of the obtained influence value of the aftershock data;

[0078] 3. Retrieve the environmental type information through the environmental information around the previous location, query and determine the reference value of the type area, retrieve the changed area value through the environmental change information, and determine whether the changed area value is greater than the reference value of the type area. When it is greater, calculate the area deviation value and determine the influence value of the area deviation based on the environmental type information as the influence value of the environmental change. When it is not greater, determine the influence value of the changed area through the aftershock relocation of the previous data and the changed area value as the influence value of the environmental change, thereby improving the accuracy of the obtained influence value of the environmental change. Brief Description of the Drawings

[0079] Figure 1 is the flowchart of the method for constructing a three-dimensional model of the seismogenic fault based on aftershock relocation data in an embodiment of the present application;

[0080] Figure 2 is the flowchart of the method for determining the influence value of aftershock data in an embodiment of the present application;

[0081] Figure 3 is the flowchart of the method for determining the influence value of environmental change in an embodiment of the present application;

[0082] Figure 4 is the flowchart of the method for determining the influence value of area deviation in an embodiment of the present application;

[0083] Figure 5 is the flowchart of the method for determining the influence value of time change in an embodiment of the present application;

[0084] Figure 6 is the flowchart of the method for determining the influence value of the changed area in an embodiment of the present application;

[0085] Figure 7 is the flowchart of the method for determining the reference value of the changed weather area and the influence value of the changed weather in an embodiment of the present application;

[0086] Figure 8 is the flowchart of the method for determining the influence value of abnormal weather detection in an embodiment of the present application. Detailed Embodiments

[0087] The present invention will be further described in detail below with reference to the drawings and embodiments.

[0088] A method for constructing a three-dimensional model of a seismogenic fault based on aftershock relocation data, which obtains the real-time data of aftershock relocation and the current time point, thereby determining the influence value of aftershock data generated by environmental and other factors, and then determines the adjusted data of aftershock relocation through the influence value of aftershock data and the real-time data of aftershock relocation, and determines and outputs the final three-dimensional model of the seismogenic fault through the adjusted data of aftershock relocation, so that the obtained three-dimensional model is affected by the previous aftershock, thereby improving the accuracy of the constructed three-dimensional model.

[0089] Refer to Figure 1 , an embodiment of the present invention discloses a method for constructing a three-dimensional model of a seismogenic fault based on aftershock relocation data, which includes:

[0090] Step S100: Obtain the real-time data of aftershock relocation.

[0091] Among them, the real-time data of aftershock relocation refers to the data obtained by real-time monitoring of the location, depth, magnitude and other information of aftershocks after the main shock occurs. The real-time data of aftershock relocation is obtained through seismic detection stations.

[0092] Step S101: Retrieve the previous data of aftershock relocation based on the real-time data of aftershock relocation.

[0093] Among them, the previous data of aftershock relocation refers to the data corresponding to the previous aftershock of the aftershock detected at the current time. The previous data of aftershock relocation is retrieved through the real-time data of aftershock relocation for subsequent use.

[0094] Step S102: Determine the influence value of aftershock data according to the previous data of aftershock relocation.

[0095] Among them, the influence value of aftershock data refers to the influence degree value corresponding to the influence of the previous aftershock on the data of the current aftershock. By analyzing the previous data of aftershock relocation, the influence value of aftershock data is determined for subsequent use. The specific steps for determining the influence value of aftershock data refer to Step S200 to Step S206.

[0096] Step S103: Determine the adjusted data of aftershock relocation according to the influence value of aftershock data and the real-time data of aftershock relocation.

[0097] Among them, the adjusted data of aftershock relocation refers to the data obtained by adjusting the aftershock data detected at the current time. The real-time data of aftershock relocation is updated and adjusted through the influence value of aftershock data to form the adjusted data of aftershock relocation for subsequent use.

[0098] Step S104: Determine the final three-dimensional model of the seismogenic fault according to the adjusted data of aftershock relocation, and output the final three-dimensional model of the seismogenic fault.

[0099] Among them, the final three-dimensional model of the seismogenic fault refers to the final three-dimensional model of the fault corresponding to the earthquake. By using the aftershock relocation adjustment data to form a two-dimensional fault line through three-dimensional slicing technology, and then performing spatial interpolation to generate an initial three-dimensional model, and then adjusting according to the distance from the aftershocks of the specified magnitude to the seismogenic fault plane of the initial three-dimensional model to obtain the final three-dimensional model and use it as the final three-dimensional model of the seismogenic fault, and then output the final three-dimensional model of the seismogenic fault, so that the obtained three-dimensional model is affected by the previous aftershocks, thereby improving the accuracy of the constructed three-dimensional model.

[0100] In Figure 1 In step S102 shown, in order to further ensure the rationality of the influence value of the aftershock data, it is necessary to perform a further separate analysis and calculation on the influence value of the aftershock data. Specifically, it is described in detail through Figure 2 the steps shown.

[0101] Referring to Figure 2 , the method for determining the influence value of the aftershock data includes the following steps:

[0102] Step S200: Retrieve the previous position point of the aftershock based on the previous aftershock relocation data.

[0103] Among them, the previous position point of the aftershock refers to the position point corresponding to the previous aftershock. Retrieving the previous position point of the aftershock through the previous aftershock relocation data facilitates subsequent use.

[0104] Step S201: Calculate the distance value between the previous position point of the aftershock and the preset data acquisition position point and use it as the previous position distance value.

[0105] Among them, the data acquisition position point refers to the position point where the seismic detection station for collecting aftershock data is located, and the data acquisition position point is obtained through pre-input. The previous position distance value refers to the distance value between the position where the previous aftershock occurred and the detection position. Calculating the distance value between the previous position point of the aftershock and the preset data acquisition position point and using it as the previous position distance value facilitates subsequent use.

[0106] Step S202: Determine the previous position distance influence value corresponding to the previous position distance value according to the corresponding relationship between the previous position distance value and the preset previous position distance influence value.

[0107] Among them, the influence value of the last position distance refers to the degree of influence of the distance between the position of the last aftershock and the detection position on the detection data. Different last position distance values correspond to different influence values of the last position distance. The influence value of the last position distance is obtained by querying from a database that stores different last position distance values and their corresponding influence values of the last position distance. This database is obtained through pre-input. Querying and determining the influence value of the last position distance through the last position distance value facilitates subsequent use.

[0108] Step S203: Retrieve the environmental information around the last position based on the last position point of the aftershock.

[0109] Among them, the environmental information around the last position refers to the environmental information located around the position of the last aftershock. The environmental information around the last position is obtained by querying and retrieving from a database that pre-stores the real-time environment around this position through the last position point of the aftershock, facilitating subsequent use.

[0110] Step S204: Determine the environmental change information based on the environmental information around the last position.

[0111] Among them, the environmental change information refers to parameter information such as the change area corresponding to when the environment changes. By retrieving the environmental information around the last position from the environment from the time of the last aftershock to the current time, and comparing the environments at the two times, the environmental change information is obtained, facilitating subsequent use.

[0112] Step S205: Determine the environmental change influence value based on the environmental change information.

[0113] Among them, the environmental change influence value refers to the degree of influence generated when the environment changes. By analyzing the environmental change information, the environmental change influence value is determined, facilitating subsequent use. The specific steps for determining the environmental change influence value refer to Step S300 to Step S306.

[0114] Step S206: Calculate the sum of the influence value of the last position distance and the environmental change influence value and use it as the aftershock data influence value.

[0115] Among them, by calculating the sum of the influence value of the last position distance and the environmental change influence value and using it as the aftershock data influence value, the accuracy of the obtained aftershock data influence value is improved.

[0116] In Figure 2 In Step S205 shown, in order to further ensure the rationality of the environmental change influence value, it is necessary to perform a further separate analysis and calculation on the environmental change influence value. Specifically, it is described in detail through the Figure 3 steps shown.

[0117] Refer toFigure 3 , the method for determining the environmental change influence value includes the following steps:

[0118] Step S300: Retrieve the environmental type information based on the environmental information around the last position.

[0119] Among them, the environmental type information refers to the type information of the environment surrounding the position of the last aftershock. Retrieving the environmental type information through the environmental information around the last position facilitates subsequent use.

[0120] Step S301: Determine the type area reference value corresponding to the environmental type information according to the correspondence between the environmental type information and the preset type area reference value.

[0121] Among them, the type area reference value refers to the reference conversion value of the area change corresponding to the environmental type. Different environmental type information corresponds to different type area reference values. The type area reference value is obtained by querying from a database storing environmental type information and the corresponding type area reference values, and this database is obtained through pre-input. Querying and determining the type area reference value through the environmental type information facilitates subsequent use.

[0122] Step S302: Retrieve the changed area value based on the environmental change information.

[0123] Among them, the changed area value refers to the actual area value that changes when the environment changes. Retrieving the changed area value through the environmental change information facilitates subsequent use.

[0124] Step S303: Determine whether the changed area value is greater than the type area reference value. If it is, execute Step S304; if not, execute Step S306.

[0125] Among them, by judging whether the changed area value is greater than the type area reference value, it is thus judged whether there is a large change in the environment around the last aftershock.

[0126] Step S304: Calculate the difference between the changed area value and the type area reference value and use it as the area deviation value.

[0127] Among them, the area deviation value refers to the deviation value corresponding when the actual area of the change in the surrounding environment exceeds the reference change area. When the changed area value is greater than the type area reference value, it indicates that there is a large change in the environment around the last aftershock at this time. Therefore, calculate the difference between the changed area value and the type area reference value and use it as the area deviation value for subsequent use.

[0128] Step S305: Determine the area deviation influence value according to the area deviation value and the environmental type information, and use the area deviation influence value as the environmental change influence value.

[0129] Among them, the area deviation influence value refers to the influence degree value generated when there is a deviation in the changed area. By analyzing the area deviation value and the environmental type information, the area deviation influence value is determined, and the area deviation influence value is used as the environmental change influence value to improve the accuracy of the obtained environmental change influence value. The specific determination steps of the area deviation influence value refer to steps S400 to S405.

[0130] Step S306: Determine the changed area influence value according to the aftershock relocation previous data and the changed area value, and use the changed area influence value as the environmental change influence value.

[0131] Among them, the changed area influence value refers to the influence degree value generated when there is no deviation in the changed area. When the changed area value is not greater than the type area reference value, it indicates that there has been no significant change in the environment around the previous aftershock at this time. Therefore, by analyzing the aftershock relocation previous data and the changed area value, the changed area influence value is determined, and the changed area influence value is used as the environmental change influence value to improve the accuracy of the obtained environmental change influence value. The specific determination steps of the changed area influence value refer to steps S600 to S608.

[0132] In Figure 3 In step S305 shown, in order to further ensure the rationality of the area deviation influence value, it is necessary to conduct a further separate analysis and calculation of the area deviation influence value, specifically through Figure 4 the steps shown for detailed description.

[0133] Referring to Figure 4 , the determination method of the area deviation influence value includes the following steps:

[0134] Step S400: Obtain the current time point.

[0135] Among them, the current time point refers to the time point corresponding to the current time, and the current time point is obtained by querying the database for timing through the current time.

[0136] Step S401: Based on the aftershock relocation previous data, retrieve the previous data time point.

[0137] Among them, the previous data time point refers to the time point corresponding to when the aftershock relocation previous data is detected. By querying the time point corresponding to when the aftershock relocation previous data is detected and retrieving it as the previous data time point, it is convenient for subsequent use.

[0138] Step S402: Calculate the time period between the previous data time point and the current time point and use it as the previous data elapsed time value.

[0139] Among them, the time elapsed since the last data refers to the time elapsed from the last data to the current time. By calculating the time period between the time point of the last data and the current time point and using it as the time elapsed since the last data, it is convenient for subsequent use.

[0140] Step S403: According to the correspondence between the environmental type information and the preset type time change reference value, determine the type time change reference value corresponding to the environmental type information.

[0141] Among them, the type time change reference value refers to the reference value of the area change of the environmental type after passing through a unit time. Different environmental type information corresponds to different type time change reference values. The type time change reference value is obtained by querying from a database that stores different environmental type information and the corresponding type time change reference values, and this database is obtained through pre-input. Querying and obtaining the type time change reference value through the environmental type information is convenient for subsequent use.

[0142] Step S404: Calculate the product value between the type time change reference value and the time elapsed since the last data and use it as the predicted type time change value.

[0143] Among them, the predicted type time change value refers to the predicted value of the area change of the environmental type based on the time elapsed since the last data collection. By calculating the product value between the type time change reference value and the time elapsed since the last data and using it as the predicted type time change value, it is convenient for subsequent use.

[0144] Step S405: Determine the time change influence value based on the area deviation value and the predicted type time change value, and use the time change influence value as the area deviation influence value.

[0145] Among them, the time change influence value refers to the influence degree value generated by the area deviation value in the time elapsed since the last data collection. By analyzing the area deviation value and the predicted type time change value, the time change influence value is determined, and the time change influence value is used as the area deviation influence value to improve the accuracy of the obtained area deviation influence value. The specific determination steps of the time change influence value refer to Step S500 to Step S506.

[0146] In Figure 4 the shown Step S405, in order to further ensure the rationality of the time change influence value, it is necessary to perform a further separate analysis and calculation on the time change influence value. Specifically, it is described in detail through Figure 5 the shown steps.

[0147] Referring to Figure 5 , the determination method of the time change influence value includes the following steps:

[0148] Step S500: Determine the type deviation reference impact value corresponding to the environmental type information according to the correspondence between the environmental type information and the preset type deviation reference impact value.

[0149] Among them, the type deviation reference impact value refers to the reference impact degree value generated per unit area of the environmental type change over time. Different environmental type information corresponds to different type deviation reference impact values. The type deviation reference impact value is obtained by querying a database that stores different environmental type information and the corresponding type deviation reference impact values. This database is obtained through pre-input. Querying and determining the type deviation reference impact value through the environmental type information facilitates subsequent use.

[0150] Step S501: Determine whether the area deviation value is less than the estimated type time change value. If yes, execute Step S502; if no, execute Step S503.

[0151] Among them, by judging whether the area deviation value is less than the estimated type time change value, it is thus judged whether the deviation of the changing area is caused by the elapsed time of data collection.

[0152] Step S502: Calculate the product value between the area deviation value and the type deviation reference impact value and use it as the time change impact value.

[0153] Among them, the time change impact value refers to the impact degree value generated by the area deviation value over time. When the area deviation value is less than the estimated type time change value, it indicates that the deviation of the changing area at this time is caused by the elapsed time of data collection. Therefore, calculate the product value between the area deviation value and the type deviation reference impact value and use it as the time change impact value, thereby improving the accuracy of the obtained time change impact value.

[0154] Step S503: Calculate the difference between the area deviation value and the estimated type time change value and use it as the deviation anomaly value.

[0155] Among them, the deviation anomaly value refers to the deviation value between the area deviation value and the estimated changing area over time. When the area deviation value is not less than the estimated type time change value, it indicates that the deviation of the changing area at this time is not caused by the elapsed time of data collection. Therefore, calculate the difference between the area deviation value and the estimated type time change value and use it as the deviation anomaly value, which facilitates subsequent use.

[0156] Step S504: Determine the abnormal deviation impact value corresponding to the deviation anomaly value according to the correspondence between the deviation anomaly value and the preset abnormal deviation impact value.

[0157] Among them, the abnormal deviation influence value refers to the influence degree value generated by the abnormal deviation. Different abnormal deviation values correspond to different abnormal deviation influence values. The abnormal deviation influence value is obtained by querying from a database that stores different abnormal deviation values and their corresponding abnormal deviation influence values. This database is obtained through pre-input. Querying and determining the abnormal deviation influence value through the abnormal deviation value is convenient for subsequent use.

[0158] Step S505: According to the corresponding relationship between the predicted value of the change over time of the category and the preset baseline influence value of the change prediction, determine the baseline influence value of the change prediction corresponding to the predicted value of the change over time of the category.

[0159] Among them, the baseline influence value of the change prediction refers to the baseline influence degree value generated by the predicted change area estimated during the time elapsed since the last data collection according to the environmental category. Different predicted values of the change over time of the category correspond to different baseline influence values of the change prediction. The baseline influence value of the change prediction is obtained by querying from a database that stores different predicted values of the change over time of the category and their corresponding baseline influence values of the change prediction. This database is obtained through pre-input. Querying and determining the baseline influence value of the change prediction through the predicted value of the change over time of the category is convenient for subsequent use.

[0160] Step S506: Calculate the sum of the baseline influence value of the change prediction and the abnormal deviation influence value and use it as the influence value of the change over time.

[0161] Among them, by calculating the sum of the baseline influence value of the change prediction and the abnormal deviation influence value and using it as the influence value of the change over time, the accuracy of the obtained influence value of the change over time is improved.

[0162] In Figure 3 In step S306 shown, in order to further ensure the rationality of the influence value of the change area, it is necessary to conduct a further separate analysis and calculation of the influence value of the change area. Specifically, it is described in detail through the Figure 6 steps shown.

[0163] Referring to Figure 6 , the method for determining the influence value of the change area includes the following steps:

[0164] Step S600: Based on the aftershock relocation, retrieve the last data time point of the last data.

[0165] Step S601: Determine the weather information at the last location according to the last data time point and the last location point of the aftershock.

[0166] Among them, the weather information of the previous position refers to the weather information corresponding to the position of the previous aftershock at the previous data time point. By retrieving the real-time weather conditions corresponding to the previous position point of the aftershock and selecting the weather conditions consistent with the previous data time point as the weather information of the previous position, it is convenient for subsequent use.

[0167] Step S602: Determine the reference area value of the changing weather and the influence value of the changing weather according to the weather information of the previous position.

[0168] Among them, the reference area value of the changing weather refers to the minimum area value corresponding to the change in area caused by the weather, and the influence value of the changing weather refers to the influence degree value generated by the unit changing area under the weather corresponding to the weather information of the previous position. By analyzing the weather information of the previous position, the reference area value of the changing weather and the influence value of the changing weather are determined, which is convenient for subsequent use. The specific determination steps of the reference area value of the changing weather and the influence value of the changing weather refer to Step S700 to Step S706.

[0169] Step S603: Determine the initial influence value of the changing area corresponding to the changing area value according to the corresponding relationship between the changing area value and the preset initial influence value of the changing area.

[0170] Among them, the initial influence value of the changing area refers to the initial influence degree value corresponding to the influence generated by the changing area value. Different changing area values correspond to different initial influence values of the changing area. The initial influence value of the changing area is obtained by querying from a database storing different changing area values and the corresponding initial influence values of the changing area, and this database is obtained through pre-input. Determining the initial influence value of the changing area by querying the changing area value is convenient for subsequent use.

[0171] Step S604: Determine whether the changing area value is less than the reference area value of the changing weather. If it is, execute Step S605; if not, execute Step S606.

[0172] Among them, by judging whether the changing area value is less than the reference area value of the changing weather, it is judged whether the changing area value is caused by the weather.

[0173] Step S605: Take the initial influence value of the changing area as the influence value of the changing area.

[0174] Among them, when the changing area value is less than the reference area value of the changing weather, it indicates that the changing area value is not caused by the weather at this time. Therefore, the initial influence value of the changing area is taken as the influence value of the changing area, thereby improving the accuracy of the obtained influence value of the changing area.

[0175] Step S606: Calculate the difference between the changing area value and the reference area value of the changing weather and take it as the area anomaly value.

[0176] Among them, the area outlier refers to the abnormal deviation value between the changed area value and the changed area caused by the weather. When the changed area value is not less than the changed weather reference area value, it indicates that the changed area value is caused by the weather at this time. Therefore, the difference between the changed area value and the changed weather reference area value is calculated and used as the area outlier for subsequent use.

[0177] Step S607: Calculate the product value of the area outlier and the changed weather impact value and use it as the area outlier impact value.

[0178] Among them, the area outlier impact value refers to the impact degree value generated by the area outlier under the influence of the weather. By calculating the product value of the area outlier and the changed weather impact value and using it as the area outlier impact value, it is convenient for subsequent use.

[0179] Step S608: Calculate the sum value of the area outlier impact value and the initial changed area impact value and use it as the adjusted changed area impact value, and use the adjusted changed area impact value as the changed area impact value.

[0180] Among them, the adjusted changed area impact value refers to the impact degree value after adjusting the impact on the changed area value. By calculating the sum value of the area outlier impact value and the initial changed area impact value and using it as the adjusted changed area impact value, and using the adjusted changed area impact value as the changed area impact value, the accuracy of the obtained changed area impact value can be improved.

[0181] In Figure 6 In step S602 shown, in order to further ensure the rationality of the changed weather reference area value and the changed weather impact value, it is necessary to conduct a further separate analysis and calculation on the changed weather reference area value and the changed weather impact value. Specifically, it is described in detail through the Figure 7 steps shown.

[0182] Referring to Figure 7 , the determination method of the changed weather reference area value and the changed weather impact value includes the following steps:

[0183] Step S700: Retrieve the weather type information and the weather detection value based on the previous location weather information.

[0184] Among them, the weather type information refers to the type information to which the weather belongs. The weather type information can be types such as sunny, rainy, snowy, etc. The weather detection value refers to the intensity value corresponding to the weather under this type. Retrieving the weather type information and the weather detection value through the previous location weather information is convenient for subsequent use.

[0185] Step S701: Determine the weather type reference area value corresponding to the weather type information according to the correspondence between the weather type information and the preset weather type reference area value.

[0186] Among them, the weather type reference area value refers to the reference value of the area change generated by the weather of the corresponding type under unit intensity. Different weather type information corresponds to different weather type reference area values. The weather type reference area value is obtained by querying from a database that stores different weather type information and the corresponding weather type reference area values, and this database is obtained through pre-input. Querying and determining the weather type reference area value through the weather type information is convenient for subsequent use.

[0187] Step S702: Calculate the product value between the weather detection value and the weather type reference area value and use it as the variable weather reference area value.

[0188] Among them, by calculating the product value between the weather detection value and the weather type reference area value and using it as the variable weather reference area value, the accuracy of the obtained variable weather reference area value is improved.

[0189] Step S703: Determine the weather type reference interval corresponding to the weather type information according to the correspondence between the weather type information and the preset weather type reference interval.

[0190] Among them, the weather type reference interval refers to the reference interval corresponding to the area of normal change of the weather of the corresponding type. Different weather type information corresponds to different weather type reference intervals. The weather type reference interval is obtained by querying from a database that stores different weather type information and the corresponding weather type reference intervals, and this database is obtained through pre-input. Querying and determining the weather type reference interval through the weather type information is convenient for subsequent use.

[0191] Step S704: Determine whether the weather detection value is within the weather type reference interval. If it is, execute Step S705; if not, execute Step S706.

[0192] Among them, by judging whether the weather detection value is within the weather type reference interval, it is thus judged whether the weather detection value is abnormal.

[0193] Step S705: Determine the weather detection impact value corresponding to the weather detection value according to the correspondence between the weather detection value and the preset weather detection impact value, and use the weather detection impact value as the variable weather impact value.

[0194] Among them, the weather detection impact value refers to the impact degree value corresponding to when the weather detection value has an impact. Different weather detection values correspond to different weather detection impact values. The weather detection impact value is obtained by querying from a database that stores different weather detection values and their corresponding weather detection impact values. This database is obtained through pre-input. When the weather detection value is within the weather type reference interval, it indicates that the weather detection value is normal at this time. Therefore, the weather detection impact value is determined by querying the weather detection value, and the weather detection impact value is used as the variable weather impact value, thereby improving the accuracy of the obtained variable weather impact value.

[0195] Step S706: Determine the weather detection abnormal impact value according to the weather detection value and the weather type reference interval, and use the weather detection abnormal impact value as the variable weather impact value.

[0196] Among them, the weather detection abnormal impact value refers to the impact degree value generated when the weather detection value is abnormal. When the weather detection value is not within the weather type reference interval, it indicates that the weather detection value is abnormal at this time. Therefore, by analyzing the weather detection value and the weather type reference interval, the weather detection abnormal impact value is determined, and the weather detection abnormal impact value is used as the variable weather impact value to improve the accuracy of the obtained variable weather impact value. The specific steps for determining the weather detection abnormal impact value refer to Step S800 to Step S806.

[0197] In Figure 7 In the shown Step S706, in order to further ensure the rationality of the weather detection abnormal impact value, it is necessary to conduct a further separate analysis and calculation of the weather detection abnormal impact value. Specifically, it is described in detail through the steps shown in Figure 8 shown.

[0198] Referring to Figure 8 , the method for determining the weather detection abnormal impact value includes the following steps:

[0199] Step S800: Calculate the difference between the weather detection value and the weather type reference interval and use it as the interval end deviation value.

[0200] Among them, the interval end deviation value refers to the deviation value when there is a deviation between the weather detection value and the end of the reference interval. By calculating the difference between the weather detection value and the weather type reference interval and using it as the interval end deviation value, it is convenient for subsequent use.

[0201] Step S801: Perform curve fitting based on the interval end deviation value to form an end deviation curve;

[0202] Among them, the weather detection values of the end deviation curve at different times and the curve generated when there is a deviation at the interval end are obtained by curve fitting the interval end deviation values according to time changes to form the end deviation curve, which is convenient for subsequent use.

[0203] Step S802: Analyze the deviation between the end deviation curve and the preset end deviation reference curve to form curve deviation information.

[0204] Among them, the end deviation reference curve refers to the reference curve corresponding to the situation when there is a deviation at the end. The end deviation reference curve is obtained through pre-input. The curve deviation information refers to the deviation information corresponding to the situation when there is a deviation in the curve. By analyzing the deviation between the end deviation curve and the preset end deviation reference curve, curve deviation information is formed, which is convenient for subsequent use.

[0205] Step S803: Retrieve the curve deviation position points and curve deviation values based on the curve deviation information.

[0206] Among them, the curve deviation position point refers to the position point where the curve has a deviation, and the curve deviation value refers to the deviation value corresponding to the situation when the curve has a deviation. By using the curve deviation information, the curve deviation position points and curve deviation values are retrieved, which is convenient for subsequent use.

[0207] Step S804: Calculate the distance value between two adjacent curve deviation position points and use it as the curve adjacent deviation distance value.

[0208] Among them, the curve adjacent deviation distance value refers to the distance value between two adjacent deviation positions. By calculating the distance value between two adjacent curve deviation position points and using it as the curve adjacent deviation distance value, it is convenient for subsequent use.

[0209] Step S805: Determine the deviation distance unit influence value corresponding to the curve adjacent deviation distance value according to the corresponding relationship between the curve adjacent deviation distance value and the preset deviation distance unit influence value.

[0210] Among them, the deviation distance unit influence value refers to the influence degree value generated by the deviation of the deviation distance unit. Different curve adjacent deviation distance values correspond to different deviation distance unit influence values. The deviation distance unit influence value is obtained by querying from a database storing different curve adjacent deviation distance values and their corresponding deviation distance unit influence values. This database is obtained through pre-input. By querying with the curve adjacent deviation distance value, the deviation distance unit influence value is determined, which is convenient for subsequent use.

[0211] Step S806: Calculate the product value of the deviation distance unit influence value and the curve deviation value and use it as the weather detection anomaly influence value.

[0212] Among them, by calculating the product value between the influence value per unit of deviation distance and the curve deviation value and using it as the influence value of abnormal weather detection, the accuracy of the obtained influence value of abnormal weather detection is improved.

[0213] Based on the same inventive concept, an embodiment of the present invention provides a device for constructing a three-dimensional model of a seismogenic fault based on aftershock relocation data, including:

[0214] An acquisition module, configured to acquire real-time aftershock relocation data and the current time point;

[0215] A memory, configured to store the method for constructing a three-dimensional model of a seismogenic fault based on aftershock relocation data as described above;

[0216] A processor, configured to load and execute the program in the memory.

[0217] Those skilled in the art can clearly understand that for the convenience and brevity of description, only the above division of each functional module is used as an example. In actual applications, the above functions can be allocated to different functional modules according to needs, that is, the internal structure of the device is divided into different functional modules to complete all or part of the functions described above. The specific working processes of the above-described system, device, and unit can refer to the corresponding processes in the foregoing method embodiments and will not be elaborated herein.

[0218] The above is only the preferred embodiment of the present invention, and the protection scope of the present invention is not limited to the above embodiments. All technical solutions falling within the idea of the present invention belong to the protection scope of the present invention. It should be pointed out that for those of ordinary skill in the art, several improvements and refinements made without departing from the principle of the present invention should also be regarded as within the protection scope of the present invention.

Claims

1. A method for constructing a three-dimensional model of a seismic fault based on aftershock relocation data, characterized in that: include: Obtain real-time data on aftershock relocation; Retrieve the last aftershock relocation data based on the aftershock relocation real-time data; Determine the impact value of aftershock data based on the last data of aftershock relocation; Determine aftershock relocation adjustment data based on aftershock data impact value and aftershock relocation real-time data; Determine the final three-dimensional model of the seismogenic fault according to the aftershock relocation adjustment data, and output the final three-dimensional model of the seismogenic fault; Methods for determining the impact value of aftershock data include: Retrieve the last location point of the aftershock based on the last aftershock relocation data; Calculate the distance between the last location of the aftershock and the preset data collection location and use it as the last location distance value; According to the corresponding relationship between the last position distance value and the preset last position distance influence value, the last position distance influence value corresponding to the last position distance value is determined; Based on the last location point of the aftershock, retrieve the surrounding environment information of the last location; Determine the environmental change information based on the surrounding environment information of the last location; Determine the impact value of environmental changes based on environmental change information; Calculate the sum of the impact value of the last position distance and the impact value of the environmental change and use it as the impact value of the aftershock data; Methods for determining the impact value of environmental changes include: Retrieving environment type information based on the last location surrounding environment information; According to the correspondence between the environmental category information and the preset category area reference value, the category area reference value corresponding to the environmental category information is determined; Retrieve the changed area value based on the environmental change information; Determine whether the changed area value is greater than the category area benchmark value; If yes, the difference between the changed area value and the category area benchmark value is calculated and used as the area deviation value; Determine the area deviation impact value according to the area deviation value and the environment type information, and use the area deviation impact value as the environment change impact value; If not, the change area impact value is determined based on the last aftershock relocation data and the change area value, and the change area impact value is used as the environmental change impact value.

2. The method for constructing a three-dimensional model of a seismic fault based on aftershock relocation data according to claim 1, characterized in that: Methods for determining the area deviation impact value include: Get the current time point; Retrieve the last data time point based on aftershock relocation; Calculate the time period between the last data time point and the current time point and use it as the last data elapsed time value; According to the correspondence between the environment type information and the preset type time change reference value, the type time change reference value corresponding to the environment type information is determined; Calculate the product of the time change benchmark value of the category and the time value of the last data and use it as the estimated value of the time change of the category; The time change impact value is determined based on the area deviation value and the estimated value of the time change of the type, and the time change impact value is used as the area deviation impact value.

3. The method for constructing a three-dimensional model of a seismic fault based on aftershock relocation data according to claim 2, characterized in that: Methods for determining the time-varying impact value include: According to the correspondence between the environmental category information and the preset category deviation baseline impact value, the category deviation baseline impact value corresponding to the environmental category information is determined; Determine whether the area deviation value is less than the estimated value of species temporal change; If yes, then the product of the area deviation value and the species deviation benchmark impact value is calculated and used as the time variation impact value; If not, the difference between the area deviation value and the estimated value of species temporal change is calculated and used as the deviation outlier; According to the correspondence between the deviation abnormal value and the preset abnormal deviation impact value, the abnormal deviation impact value corresponding to the deviation abnormal value is determined; According to the correspondence between the estimated value of the time change of the category and the preset estimated change benchmark impact value, the estimated change benchmark impact value corresponding to the estimated value of the time change of the category is determined; The sum of the change estimate baseline impact value and the abnormal deviation impact value is calculated and used as the time change impact value.

4. The method for constructing a three-dimensional model of a seismic fault based on aftershock relocation data according to claim 2, characterized in that: Methods for determining the impact value of the change area include: Retrieve the last data time point based on aftershock relocation; Determine the weather information of the last location based on the time point of the last data and the last location of the aftershock; Determine the weather base area value and weather impact value of the change according to the weather information of the last location; According to the correspondence between the change area value and the preset change area initial impact value, the change area initial impact value corresponding to the change area value is determined; Determine whether the change area value is less than the change weather base area value; If yes, the initial impact value of the changed area is used as the impact value of the changed area; If not, the difference between the changed area value and the changed weather base area value is calculated and used as the area anomaly value; Calculate the product of the area anomaly value and the weather change impact value and use it as the area anomaly impact value; The sum of the area abnormality impact value and the change area initial impact value is calculated and used as the change area impact adjustment value, and the change area impact adjustment value is used as the change area impact value.

5. The method for constructing a three-dimensional model of a seismic fault based on aftershock relocation data according to claim 4, characterized in that: The methods for determining the changing weather base area value and the changing weather impact value include: Retrieve weather type information and weather detection value based on the last location weather information; According to the correspondence between the weather type information and the preset weather type reference area value, the weather type reference area value corresponding to the weather type information is determined; Calculate the product value between the weather detection value and the weather type reference area value and use it as the change weather reference area value; According to the correspondence between the weather type information and the preset weather type reference interval, the weather type reference interval corresponding to the weather type information is determined; Determine whether the weather detection value is within the weather type reference range; If yes, then determine the weather detection impact value corresponding to the weather detection value according to the correspondence between the weather detection value and the preset weather detection impact value, and use the weather detection impact value as the change weather impact value; If not, the weather detection anomaly impact value is determined based on the weather detection value and the weather type reference interval, and the weather detection anomaly impact value is used as the changed weather impact value.

6. The method for constructing a three-dimensional model of a seismic fault based on aftershock relocation data according to claim 5, characterized in that: Methods for determining the impact value of weather detection anomaly include: Calculate the difference between the weather detection value and the weather type reference interval and use it as the interval end deviation value; Performing curve processing based on the end deviation value of the interval to form an end deviation curve; Analyzing the deviation between the end deviation curve and a preset end deviation reference curve to form curve deviation information; Retrieve the curve deviation position point and the curve deviation value based on the curve deviation information; Calculate the distance between two adjacent curve deviation position points and use it as the curve adjacent deviation distance value; According to the corresponding relationship between the curve adjacent deviation distance value and the preset deviation distance unit influence value, the deviation distance unit influence value corresponding to the curve adjacent deviation distance value is determined; The product of the deviation distance unit impact value and the curve deviation value is calculated and used as the weather detection anomaly impact value.

7. A device for constructing a three-dimensional model of a seismic fault based on aftershock relocation data, characterized in that: include: The acquisition module is used to obtain the real-time data of aftershock relocation and the current time point; A memory for storing a program for implementing the method for constructing a three-dimensional model of a seismic fault based on aftershock relocation data according to any one of claims 1 to 6; The processor loads and executes the program in the memory.

Images