A method for integrating and sharing multi-source seismic data

Abstract

The invention discloses a method for integrating and sharing multi-source seismic data, which belongs to the field of data processing and analysis, including collecting seismic data acquired from the same seismic data source into the same data group according to all data sources of monitored seismic data, obtaining a seismic source data set; screening and classifying seismic data acquired from different seismic data sources according to data encoding formats, obtaining data groups of different target seismic data types; extracting data features of seismic data in the data group, calculating seismic response coefficients by seismic response degrees characterized by the data features, and integrating the seismic response coefficients to obtain seismic response data packets; utilizing data transmission times of different transmission paths and the number of sub-data packets that can be transmitted, sequentially sending sub-data packets in a sub-data packet sharing sequence to a shared user end. The invention ensures the wide dissemination and effective use of seismic data, and has significant advantages and application prospects in promoting the progress of seismic science research.

Description

Technical Field

[0001] The invention relates to the field of data processing and analysis, and in particular to a method for integrating and sharing multi-source seismic data. Background Art

[0002] Earthquakes, as one of the natural disasters, pose a huge threat to human life and property. Improving the accuracy and timeliness of earthquake prediction is an important means to reduce the losses caused by earthquake disasters. Traditional earthquake prediction methods mainly rely on earthquake monitoring data from a single source, and their prediction accuracy and reliability are limited. In recent years, with the continuous advancement of earthquake monitoring technology and the continuous enrichment of data collection methods, the sources and types of earthquake data have increased, including earthquake monitoring data, geophysical measurement data, geographic information data, etc. However, these multi-source earthquake data have large differences in format, resolution, spatiotemporal distribution, etc. How to effectively integrate and share these multi-source earthquake data and comprehensively utilize earthquake data from multiple sources has become a difficulty in current earthquake research. Therefore, it is urgent to propose a method for integrating and sharing multi-source earthquake data. Summary of the invention

[0003] In view of the above-mentioned deficiencies in the prior art, the present invention provides a method for integrating and sharing multi-source seismic data, which comprehensively utilizes seismic data from multiple sources and realizes the integration and sharing of multi-source data.

[0004] In order to achieve the above-mentioned object of the invention, the technical solution adopted by the present invention is:

[0005] A method for integrating and sharing multi-source seismic data is provided, comprising:

[0006] Step S1: According to all data sources of monitored seismic data, seismic data acquired from the same seismic data source are collected into the same data group to obtain a seismic source data set. ;in, For the n The first m earthquake data, n is the number of the earthquake data source, m is the number of the earthquake data;

[0007] S2: define the data encoding format of the target seismic data type to be acquired, screen the seismic data acquired from different seismic data sources according to the data encoding format, classify the seismic data according to the size of the storage space occupied by the seismic data, classify the seismic data acquired from different seismic data sources into the corresponding target seismic data type, and obtain data groups of different target seismic data types;

[0008] S3: extracting data features of seismic data in the data group according to the data group of different target seismic data types, calculating seismic response coefficients according to the seismic response degree represented by the data features, and integrating the seismic response coefficients to obtain seismic response data packets;

[0009] S4: Pack each seismic response coefficient in the seismic response data packet into several sub-packets, sort them in sequence according to the size of the seismic response coefficient value to form a sub-packet sharing sequence, and use the data sending time corresponding to different transmission paths linked to the shared user end and the number of sub-packets that can be sent to send the sub-packets in the sub-packet sharing sequence to the shared user end in sequence.

[0010] Further, step S2 includes:

[0011] Step S21: define the data encoding format of the target earthquake data type to be acquired, according to the earthquake data Data encoding format, calculation of earthquake data The encoding difference coefficient between the data encoding format of and the data encoding format of each target seismic data type ;

[0012] ;

[0013] in, For earthquake data The number of data encoding symbols that are different from the data encoding symbols of the target seismic data type, For earthquake data Data encoding symbol, For the i Data encoding symbol of target seismic data type, i is the target earthquake data type, The number of different thresholds for data encoding symbols, l a For earthquake data The data encoding length, l i For the i The data encoding length of each target seismic data type, is the difference threshold of data encoding length, is the natural exponential function;

[0014] Step S22: Obtain earthquake data A dataset of coded difference coefficients for all target seismic data types , set the threshold of the encoding difference coefficient , if exists , then the earthquake data is determined is the target seismic data type required, otherwise, determine the seismic data The target earthquake data type is not needed and is deleted;

[0015] Step S23: Collect the earthquake source data After deleting the seismic data of the target seismic data type that is not needed, the remaining seismic data set is obtained. , u is the number of earthquake data deleted, The remaining earthquake data set m - u earthquake data;

[0016] Step S24: Extract the remaining earthquake data set Each earthquake data Coding difference coefficient between the data coding format of the target seismic data type and the target seismic data type ;

[0017] Step S25: Calculate earthquake data The difference between the size of the storage space occupied and the size of the storage space occupied by each target seismic data type , get the storage space difference data set , For earthquake data The amount of storage space occupied, For the i The size of the storage space occupied by each target earthquake data type;

[0018] Step S26: Calculate earthquake data based on the difference in storage space size and the encoding difference coefficient Data difference coefficient between the target seismic data type:

[0019] ;in, is the weight of the impact of the difference in storage space occupied by data on data differences, Weights of impact of differences in data coding on data differences;

[0020] Step S27: Obtaining a data set of data difference coefficients , and filter out the minimum value , minimum The corresponding target earthquake data type is earthquake data The type of earthquake data;

[0021] Step S28: Repeat steps S24-S27 to collect the seismic data Each earthquake data All are classified to form seismic data groups of different target seismic data types;

[0022] Step S29: Repeat steps S21-S28 to classify the seismic data from different seismic data sources according to different target seismic data types to form data sets of different target seismic data types. ,in, For the i The target seismic data type v earthquake data, v After earthquake data classification i The earthquake data number in the target earthquake data type.

[0023] Further, step S3 includes:

[0024] Step S31: According to the data set According to the target earthquake data type, a data feature extraction model corresponding to the target earthquake data type is constructed to obtain the data features of each earthquake data;

[0025] ;

[0026] in, is the data feature extraction function, For earthquake data About time t and space s No. k data features, k is the number of the data feature;

[0027] Step S32: Setting the normal value of each data feature in the target earthquake data type , calculate the characteristic coefficients of the data features ; Characteristic coefficient Characterizing Data The degree of earthquake response reflected;

[0028] Step S33: Obtain earthquake data The characteristic coefficient , and the characteristic coefficient The maximum value in and minimum As for earthquake data In time t and space s The seismic response coefficient on ;

[0029] Step S34: Calculate the data set The seismic response coefficient of each earthquake data in time t ,space s and magnitude-repeated earthquake response coefficients, and obtains Seismic response coefficient data set;

[0030] Step S35: Repeat steps S31-S34 to obtain the seismic response coefficient data set corresponding to each data set of the target seismic data type, and compare the time intervals between the seismic response coefficient data sets. t ,space s The earthquake response coefficients that are repeated in terms of size are deleted, and then all the retained earthquake response coefficients are integrated into one earthquake response data package.

[0031] Further, step S4 includes:

[0032] Step S41: Pack each seismic response coefficient in the seismic response data packet into a plurality of sub-data packets, and sort them in order according to the magnitude of the seismic response coefficient value to form a sub-data packet sharing sequence, wherein the sub-data packet with the largest seismic response coefficient value in the sub-data packet sharing sequence is preferentially shared with the user end, and the sub-data packet with the smallest seismic response coefficient value is shared with the user end last;

[0033] Step S42: Filter out the transmission path with the shortest data transmission time according to the data transmission times corresponding to the different transmission paths linked to the shared user terminal. W , and the number of sub-packets that can be sent M , the sub-packet shared sequence M Sub-packets use the transmission path W Transmit to the shared user end;

[0034] Step S43: Afterwards, according to the data transmission time of the remaining transmission paths and the number of sub-data packets that can be transmitted, the remaining sub-data packets in the sub-data packet sharing sequence are sequentially transmitted to the shared user end.

[0035] The beneficial effects of the present invention are as follows: the present invention can realize efficient processing and accurate analysis of seismic data from multiple sources, and integrate seismic data from different sources, thereby reducing the complexity of seismic data types from different sources, while making seismic data more integrated, facilitating data transmission and sharing, increasing the transmission efficiency of data related to the degree of earthquake response, greatly improving the accuracy of earthquake prediction, reducing the difficulty of subsequent user-side seismic data analysis, ensuring the wide dissemination and effective use of seismic data, and having significant advantages and application prospects in promoting the progress of earthquake science research. BRIEF DESCRIPTION OF THE DRAWINGS

[0036] Figure 1 Flowchart of a method for sharing integrated multi-source seismic data. DETAILED DESCRIPTION

[0037] The specific implementation modes of the present invention are described below so that those skilled in the art can understand the present invention. However, it should be clear that the present invention is not limited to the scope of the specific implementation modes. For those of ordinary skill in the art, as long as various changes are within the spirit and scope of the present invention as defined and determined by the attached claims, these changes are obvious, and all inventions and creations utilizing the concept of the present invention are protected.

[0038] like Figure 1 As shown, a method for integrating and sharing multi-source seismic data includes:

[0039] Step S1: According to all data sources of monitored seismic data, seismic data acquired from the same seismic data source are collected into the same data group to obtain a seismic source data set. ;in, For the n The first m earthquake data, n is the number of the earthquake data source, m The number of earthquake data.

[0040] In recent years, with the continuous development of earthquake data monitoring methods, the increase in monitoring methods will increase the complexity of earthquake data, such as video data, different sensor data, earthquake monitoring networks, geophysical measurement equipment, geographic information systems, etc. to obtain multi-source earthquake data. There are obvious differences in data type and data size of earthquake data from different sources. Therefore, the above characteristics can be used to screen and classify earthquake data.

[0041] S2: Define the data encoding format of the target seismic data type to be acquired, filter the seismic data acquired from different seismic data sources according to the data encoding format, and classify the seismic data based on the size of the storage space occupied by the seismic data, classify the seismic data acquired from different seismic data sources into corresponding target seismic data types, and obtain data groups of different target seismic data types.

[0042] Different data need to be encoded during transmission or processing. Different types of data have very different encoding methods, including differences in encoding symbols and encoding lengths. The present invention uses data encoding formats to screen and classify data.

[0043] Step S2 specifically includes:

[0044] Step S21: define the data encoding format of the target earthquake data type to be acquired, according to the earthquake data Data encoding format, calculation of earthquake data The encoding difference coefficient between the data encoding format of and the data encoding format of each target seismic data type ;

[0045] ;

[0046] in, For earthquake data The number of data encoding symbols that are different from the data encoding symbols of the target seismic data type, For earthquake data Data encoding symbol, For the i Data encoding symbol of target seismic data type, i is the target earthquake data type, The threshold value of the number of different data encoding symbols indicates the maximum value allowed for the number of different data encoding symbols. If the threshold is exceeded, , the data type difference is too large, l a For earthquake data The data encoding length, l i For the i The data encoding length of each target seismic data type, The difference threshold of the data encoding length indicates the maximum value allowed for the difference of the data encoding length. If it exceeds the threshold , which also indicates that the data type difference is too large. is the natural exponential function;

[0047] Step S22: Obtain earthquake data A dataset of coded difference coefficients for all target seismic data types , set the threshold of the encoding difference coefficient , threshold Indicates the maximum value allowed for the encoding difference coefficient, if any , then the earthquake data is determined is the target seismic data type required, otherwise, determine the seismic data The target earthquake data type is not needed and is deleted;

[0048] Step S23: Collect the earthquake source data After deleting the seismic data of the target seismic data type that is not needed, the remaining seismic data set is obtained. , u is the number of earthquake data deleted, The remaining earthquake data set m - u earthquake data;

[0049] Step S24: Extract the remaining earthquake data set Each earthquake data Coding difference coefficient between the data coding format of the target seismic data type and the target seismic data type ;

[0050] Step S25: Calculate earthquake data The difference between the size of the storage space occupied and the size of the storage space occupied by each target seismic data type , get the storage space difference data set , For earthquake data The amount of storage space occupied, For the i The size of the storage space occupied by each target earthquake data type;

[0051] The difference in storage space is also a typical feature that reflects the difference in data types. Digital data obviously takes up less storage space than video data.

[0052] Step S26: Calculate earthquake data based on the difference in storage space size and the encoding difference coefficient Data difference coefficient between the target seismic data type:

[0053] ;in, is the weight of the impact of the difference in storage space occupied by data on data differences, is the weight of the impact of data encoding differences on data differences. Generally, , Take 0.6, Take 0.4;

[0054] Step S27: Obtaining a data set of data difference coefficients , and filter out the minimum value , minimum The corresponding target earthquake data type is earthquake data The type of earthquake data;

[0055] Step S28: Repeat steps S24-S27 to collect the seismic data Each earthquake data All are classified to form seismic data groups of different target seismic data types;

[0056] Step S29: Repeat steps S21-S28 to classify the seismic data from different seismic data sources according to different target seismic data types to form data sets of different target seismic data types. ,in, For the i The target seismic data type vearthquake data, v After earthquake data classification i The earthquake data number in the target earthquake data type.

[0057] S3: According to the data groups of different target seismic data types, the data features of the seismic data in the data groups are extracted, the seismic response coefficients are calculated by the seismic response degree represented by the data features, and the seismic response coefficients are integrated to obtain the seismic response data packets.

[0058] Step S3 specifically includes:

[0059] Step S31: According to the data set According to the target earthquake data type, a data feature extraction model corresponding to the target earthquake data type is constructed to obtain the data features of each earthquake data;

[0060] ;

[0061] in, is the data feature extraction function, For earthquake data About time t and space s No. k data features, k is the number of the data feature;

[0062] Step S32: Setting the normal value of each data feature in the target earthquake data type , calculate the characteristic coefficients of the data features ; Characteristic coefficient Characterizing Data The degree of earthquake response reflected;

[0063] Step S33: Obtain earthquake data The characteristic coefficient , and the characteristic coefficient The maximum value in and minimum As for earthquake data In time t and space s The seismic response coefficient on ;

[0064] Step S34: Calculate the data set The seismic response coefficient of each earthquake data in time t ,space s and magnitude-repeated earthquake response coefficients, and obtains Seismic response coefficient data set;

[0065] Step S35: Repeat steps S31-S34 to obtain the seismic response coefficient data set corresponding to each data set of the target seismic data type, and compare the time intervals between the seismic response coefficient data sets. t ,space s The earthquake response coefficients that are repeated in terms of size are deleted, and then all the retained earthquake response coefficients are integrated into one earthquake response data package.

[0066] S4: Pack each seismic response coefficient in the seismic response data packet into several sub-packets, sort them in sequence according to the size of the seismic response coefficient value to form a sub-packet sharing sequence, and use the data sending time corresponding to different transmission paths linked to the shared user end and the number of sub-packets that can be sent to send the sub-packets in the sub-packet sharing sequence to the shared user end in sequence.

[0067] Step S4 specifically includes:

[0068] Step S41: Pack each seismic response coefficient in the seismic response data packet into a plurality of sub-data packets, and sort them in order according to the magnitude of the seismic response coefficient value to form a sub-data packet sharing sequence, wherein the sub-data packet with the largest seismic response coefficient value in the sub-data packet sharing sequence is preferentially shared with the user end, and the sub-data packet with the smallest seismic response coefficient value is shared with the user end last;

[0069] Step S42: Filter out the transmission path with the shortest data transmission time according to the data transmission times corresponding to the different transmission paths linked to the shared user terminal. W , and the number of sub-packets that can be sent M , the sub-packet shared sequence M Sub-packets use the transmission path W Transmit to the shared user end; ensure that the maximum earthquake response can reach the shared user end as quickly as possible, according to the time corresponding to the earthquake response coefficient t and space s Execute earthquake warning command; of course, the transmission path W The more sub-data packets that can be sent, the better, which is more conducive to the shared user end to judge the degree of earthquake damage or the corresponding earthquake level according to the earthquake response coefficient.

[0070] Step S43: Afterwards, according to the data transmission time of the remaining transmission paths and the number of sub-data packets that can be transmitted, the remaining sub-data packets in the sub-data packet sharing sequence are sequentially transmitted to the shared user end.

[0071] The present invention can realize efficient processing and accurate analysis of seismic data from multiple sources, integrate seismic data from different sources, reduce the complexity of seismic data types from different sources, and at the same time make seismic data more integrated, facilitate data transmission and sharing, increase the transmission efficiency of data related to earthquake response degree, greatly improve the accuracy of earthquake prediction, reduce the difficulty of subsequent user-side seismic data analysis, ensure the wide dissemination and effective use of seismic data, and have significant advantages and application prospects in promoting the progress of earthquake science research.

Claims

1. A method for integrating and sharing multi-source seismic data, characterized in that: include: Step S1: According to all data sources of monitored seismic data, seismic data acquired from the same seismic data source are collected into the same data group to obtain a seismic source data set. ;in, For the n The first m earthquake data, n is the number of the earthquake data source, m is the number of the earthquake data; S2: define the data encoding format of the target seismic data type to be acquired, screen the seismic data acquired from different seismic data sources according to the data encoding format, classify the seismic data according to the size of the storage space occupied by the seismic data, classify the seismic data acquired from different seismic data sources into the corresponding target seismic data type, and obtain data groups of different target seismic data types; S3: extracting data features of seismic data in the data group according to the data group of different target seismic data types, calculating seismic response coefficients according to the seismic response degree represented by the data features, and integrating the seismic response coefficients to obtain seismic response data packets; S4: Packing each seismic response coefficient in the seismic response data packet into a plurality of sub-data packets, and sequentially sorting them according to the magnitude of the seismic response coefficient values ​​to form a sub-data packet sharing sequence, and sequentially sending the sub-data packets in the sub-data packet sharing sequence to the shared user end by using the data transmission time corresponding to different transmission paths linked to the shared user end and the number of sub-data packets that can be sent; The step S2 comprises: Step S21: define the data encoding format of the target earthquake data type to be acquired, according to the earthquake data Data encoding format, calculation of earthquake data The encoding difference coefficient between the data encoding format of and the data encoding format of each target seismic data type ; ； in, For earthquake data The number of data encoding symbols that are different from the data encoding symbols of the target seismic data type, For earthquake data Data encoding symbol, For the i The data encoding symbol of the target seismic data type, i is the target earthquake data type, The number of different thresholds for data encoding symbols, l a For earthquake data The data encoding length, l i For the i The data encoding length of each target seismic data type, is the difference threshold of data encoding length, is the natural exponential function; Step S22: Obtain earthquake data A dataset of coded difference coefficients for all target seismic data types , set the threshold of the encoding difference coefficient , if exists , then the earthquake data is determined is the target seismic data type required, otherwise, determine the seismic data The target earthquake data type is not needed and is deleted; Step S23: Collect the earthquake source data After deleting the seismic data of the target seismic data type that is not needed, the remaining seismic data set is obtained. , u is the number of earthquake data deleted, The remaining earthquake data set m - u earthquake data; Step S24: Extract the remaining earthquake data set Each earthquake data Coding difference coefficient between the data coding format of the target seismic data type and the target seismic data type ; Step S25: Calculate earthquake data The difference between the size of the storage space occupied and the size of the storage space occupied by each target seismic data type , get the storage space difference data set , For earthquake data The amount of storage space occupied, For the i The size of the storage space occupied by each target earthquake data type; Step S26: Calculate earthquake data based on the difference in storage space size and the encoding difference coefficient Data difference coefficient between the target seismic data type: ;in, is the weight of the impact of the difference in storage space occupied by data on data differences, Weights of impact of differences in data coding on data differences; Step S27: Obtaining a data set of data difference coefficients , and filter out the minimum value , minimum The corresponding target earthquake data type is earthquake data The type of earthquake data; Step S28: Repeat steps S24-S27 to collect the seismic data Each earthquake data All are classified to form seismic data groups of different target seismic data types; Step S29: Repeat steps S21-S28 to classify the seismic data from different seismic data sources according to different target seismic data types to form data sets of different target seismic data types. ,in, For the i Among the target seismic data types v earthquake data, v After earthquake data classification i The earthquake data number in the target earthquake data type; The step S3 comprises: Step S31: According to the data set According to the target earthquake data type, a data feature extraction model corresponding to the target earthquake data type is constructed to obtain the data features of each earthquake data; ； in, is the data feature extraction function, For earthquake data About time t and space s No. k data features, k is the number of the data feature; Step S32: Setting the normal value of each data feature in the target earthquake data type , calculate the characteristic coefficients of the data features ; Characteristic coefficient Characterizing Data The degree of earthquake response reflected; Step S33: Obtain earthquake data The characteristic coefficient , and the characteristic coefficient The maximum value in and minimum As for earthquake data In time t and space s The seismic response coefficient on ; Step S34: Calculate the data set The seismic response coefficient of each earthquake data in time t ,space s and magnitude-repeated earthquake response coefficients, and obtains Seismic response coefficient data set; Step S35: Repeat steps S31-S34 to obtain the seismic response coefficient data set corresponding to each data set of the target seismic data type, and compare the time intervals between the seismic response coefficient data sets. t ,space s The earthquake response coefficients that are repeated in terms of size are deleted, and then all the retained earthquake response coefficients are integrated into one earthquake response data package.

2. The method for integrating and sharing multi-source seismic data according to claim 1, characterized in that: The step S4 comprises: Step S41: Pack each seismic response coefficient in the seismic response data packet into a plurality of sub-data packets, and sort them in order according to the magnitude of the seismic response coefficient value to form a sub-data packet sharing sequence, wherein the sub-data packet with the largest seismic response coefficient value in the sub-data packet sharing sequence is preferentially shared with the user end, and the sub-data packet with the smallest seismic response coefficient value is shared with the user end last; Step S42: Filter out the transmission path with the shortest data transmission time according to the data transmission times corresponding to the different transmission paths linked to the shared user terminal. W , and the number of subpackets that can be sent M , the sub-packet shared sequence M Sub-packets use the transmission path W Transmit to the shared user end; Step S43: Afterwards, according to the data transmission time of the remaining transmission paths and the number of sub-data packets that can be transmitted, the remaining sub-data packets in the sub-data packet sharing sequence are sequentially transmitted to the shared user end.

Images