Earthquake and electromagnetic inversion data processing method based on composite trigonometric function

Abstract

The invention discloses an earthquake and electromagnetic inversion data processing method based on a composite trigonometric function. The method comprises the following steps: acquiring first seismic inversion data and first electromagnetic inversion data of a work area; performing data matching on the first seismic inversion data and the first electromagnetic inversion data based on a preset rule to obtain second seismic inversion data and second electromagnetic inversion data; performing data transformation on the second seismic inversion data and the second electromagnetic inversion data based on a composite trigonometric function to obtain third seismic inversion data and third electromagnetic inversion data; and carrying out data fusion on the third seismic inversion data and the third electromagnetic inversion data so as to explain the geological structure of the work area based on the inversion data obtained after fusion. According to the technical scheme provided by the embodiment of the invention, two different types of inversion data can be quickly fused while original structural feature information in the inversion data is reserved, so that the accuracy and reliability of geological structure interpretation are improved.

Description

Technical Field

[0001] The present invention relates to the field of geophysical inversion technology, and in particular, to a method, device, electronic device and storage medium for processing seismic and electromagnetic inversion data based on composite trigonometric functions. Background Art

[0002] The Earth is an extremely large and very complex integrated system, and it is necessary to obtain and utilize Earth information to study it. Geophysical exploration, abbreviated as geophysical prospecting, refers to detecting the lithology of strata, geological structures and the distribution of underground ore bodies by studying and observing the changes of various geophysical fields. Generally speaking, different geophysical exploration methods are based on different physical parameters, and the information content of the single geophysical inversion data obtained is limited, and there are obvious defects. For example, the seismic exploration method has low resolution in special structures such as igneous rocks and salt domes, and the electromagnetic exploration method has poor lateral resolution, resulting in one-sidedness and limitation in the interpretation results of the inversion data. Summary of the Invention

[0003] The present invention provides a method, device, electronic device and storage medium for processing seismic and electromagnetic inversion data based on composite trigonometric functions, which can quickly fuse two different types of inversion data while retaining the original structural feature information in the inversion data, thereby improving the accuracy and reliability of geological structure interpretation.

[0004] According to one aspect of the present invention, there is provided a method for processing seismic and electromagnetic inversion data based on composite trigonometric functions, the method comprising:

[0005] Obtaining first seismic inversion data and first electromagnetic inversion data of a work area; wherein, the first seismic inversion data includes first spatial coordinates and the velocity corresponding to the first spatial coordinates, and the first electromagnetic inversion data includes second spatial coordinates and the resistivity corresponding to the second spatial coordinates;

[0006] Performing data matching on the first seismic inversion data and the first electromagnetic inversion data based on a preset rule to obtain second seismic inversion data and second electromagnetic inversion data;

[0007] Performing data transformation on the second seismic inversion data and the second electromagnetic inversion data respectively based on composite trigonometric functions to obtain third seismic inversion data and third electromagnetic inversion data;

[0008] Performing data fusion on the third seismic inversion data and the third electromagnetic inversion data to interpret the geological structure of the work area based on the fused inversion data.

[0009] According to another aspect of the present invention, there is provided a device for processing seismic and electromagnetic inversion data based on composite trigonometric functions, the device comprising:

[0010] An inversion data acquisition module, configured to acquire first seismic inversion data and first electromagnetic inversion data of a work area; wherein, the first seismic inversion data includes first spatial coordinates and the velocity corresponding to the first spatial coordinates, and the first electromagnetic inversion data includes second spatial coordinates and the resistivity corresponding to the second spatial coordinates;

[0011] An inversion data matching module, configured to perform data matching on the first seismic inversion data and the first electromagnetic inversion data based on a preset rule to obtain second seismic inversion data and second electromagnetic inversion data;

[0012] An inversion data transformation module, configured to perform data transformation on the second seismic inversion data and the second electromagnetic inversion data respectively based on a composite trigonometric function to obtain third seismic inversion data and third electromagnetic inversion data;

[0013] An inversion data fusion module, configured to perform data fusion on the third seismic inversion data and the third electromagnetic inversion data to interpret the geological structure of the work area based on the fused inversion data.

[0014] According to another aspect of the present invention, there is provided an electronic device, which includes:

[0015] At least one processor; and

[0016] A memory communicatively connected to the at least one processor; wherein,

[0017] The memory stores a computer program executable by the at least one processor, and when the computer program is executed by the at least one processor, the at least one processor can execute the method for processing seismic and electromagnetic inversion data based on a composite trigonometric function according to any embodiment of the present invention.

[0018] According to another aspect of the present invention, there is provided a computer-readable storage medium, which stores computer instructions for causing a processor to implement the method for processing seismic and electromagnetic inversion data based on a composite trigonometric function according to any embodiment of the present invention when executed.

[0019] In the technical solution of the embodiment of the present invention, the first seismic inversion data and the first electromagnetic inversion data of the work area are obtained; wherein, the first seismic inversion data includes the first spatial coordinates and the velocity corresponding to the first spatial coordinates, and the first electromagnetic inversion data includes the second spatial coordinates and the resistivity corresponding to the second spatial coordinates; based on a preset rule, the first seismic inversion data and the first electromagnetic inversion data are subjected to data matching to obtain the second seismic inversion data and the second electromagnetic inversion data; based on the composite trigonometric function, the second seismic inversion data and the second electromagnetic inversion data are respectively subjected to data transformation to obtain the third seismic inversion data and the third electromagnetic inversion data; the third seismic inversion data and the third electromagnetic inversion data are subjected to data fusion to interpret the geological structure of the work area based on the inversion data obtained after fusion. In the technical solution of the embodiment of the present invention, based on the composite trigonometric function, data transformation is performed on the second seismic inversion data and the second electromagnetic inversion data obtained by data matching, which can retain the original structural feature information in the inversion data while enabling two different types of inversion data to be quickly fused, thereby improving the accuracy and reliability of geological structure interpretation.

[0020] It should be understood that the content described in this part is not intended to identify the key or important features of the embodiments of the present invention, nor is it used to limit the scope of the present invention. Other features of the present invention will become easily understood through the following description. BRIEF DESCRIPTION OF THE DRAWINGS

[0021] In order to more clearly illustrate the technical solutions in the embodiments of the present invention, the following will briefly introduce the drawings required for the description of the embodiments. Obviously, the drawings in the following description are only some embodiments of the present invention, and those of ordinary skill in the art can obtain other drawings without creative efforts based on these drawings.

[0022] Figure 1 is a flowchart of a method for processing seismic and electromagnetic inversion data based on a composite trigonometric function according to Embodiment 1 of the present invention;

[0023] Figure 2 is a flowchart of a method for processing seismic and electromagnetic inversion data based on a composite trigonometric function according to Embodiment 2 of the present invention;

[0024] Figure 3 A schematic diagram of the first seismic inversion data provided by Embodiment 2 of the present invention;

[0025] Figure 4 A schematic diagram of the first electromagnetic inversion data provided by Embodiment 2 of the present invention;

[0026] Figure 5Schematic diagram of a third seismic inversion data provided by the second embodiment of the present invention;

[0027] Figure 6 Schematic diagram of a third electromagnetic inversion data provided by the second embodiment of the present invention;

[0028] Figure 7 Schematic diagram of the fused inversion data provided by the second embodiment of the present invention;

[0029] Figure 8 Schematic structural diagram of a seismic and electromagnetic inversion data processing device based on a composite trigonometric function provided by the third embodiment of the present invention;

[0030] Figure 9 Schematic structural diagram of an electronic device for implementing the seismic and electromagnetic inversion data processing method based on a composite trigonometric function of the embodiment of the present invention. Detailed implementation manners

[0031] In order to enable those skilled in the art to better understand the solution of the present invention, the technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings in the embodiments of the present invention. Obviously, the described embodiments are only a part of the embodiments of the present invention, rather than all the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those of ordinary skill in the art without creative efforts shall fall within the protection scope of the present invention.

[0032] It should be noted that the terms "first", "second", etc. in the specification and claims of the present invention and the above drawings are used to distinguish similar objects, and do not necessarily need to be used to describe a specific order or sequence. It should be understood that such data can be interchanged under appropriate circumstances so that the embodiments of the present invention described herein can be implemented in an order other than those illustrated or described herein. In addition, the terms "comprising" and "having" and any variations thereof are intended to cover non-exclusive inclusion. For example, a process, method, system, product or device comprising a series of steps or units does not necessarily have to be limited to those steps or units clearly listed, but may include other steps or units not clearly listed or inherent to these processes, methods, products or devices.

[0033] Embodiment 1

[0034] Figure 1The following is a flowchart of a method for processing seismic and electromagnetic inversion data based on composite trigonometric functions provided by Embodiment 1 of the present invention. This embodiment is applicable to the situation of processing seismic and electromagnetic inversion data of geophysics. This method can be executed by a device for processing seismic and electromagnetic inversion data based on composite trigonometric functions, which can be implemented in the form of hardware and / or software, and can be configured in an electronic device. As Figure 1 shown, the method includes:

[0035] S110. Obtain the first seismic inversion data and the first electromagnetic inversion data of the work area; wherein, the first seismic inversion data includes the first spatial coordinates and the velocity corresponding to the first spatial coordinates, and the first electromagnetic inversion data includes the second spatial coordinates and the resistivity corresponding to the second spatial coordinates.

[0036] Among them, the seismic inversion data is obtained in the process of imaging or solving the underground rock formation by combining the seismic data observed on the surface with other known information. In the embodiment of the present invention, the first seismic inversion data refers to the velocity data of the known spatial coordinates in the seismic inversion data. Specifically, the first seismic inversion data includes the first spatial coordinates and the velocity corresponding to the first spatial coordinates.

[0037] Among them, the electromagnetic inversion data is obtained in the process of imaging or solving the electromagnetic properties or structure of the underground medium by combining the electromagnetic field data observed on the surface or in the air with the known electromagnetic theory and geological laws. In the embodiment of the present invention, the first electromagnetic inversion data refers to the resistivity data of the known spatial coordinates in the electromagnetic inversion data. Specifically, the first electromagnetic inversion data includes the second spatial coordinates and the resistivity corresponding to the second spatial coordinates.

[0038] In the embodiment of the present invention, the first seismic inversion data can be obtained from the existing seismic inversion data of the work area, and the first electromagnetic inversion data can be obtained from the existing electromagnetic inversion data of the work area, providing reliable basic data for the subsequent interpretation of the geological structure of the work area. By obtaining two types of inversion data, the physical characteristics of the geological structure of the work area can be more comprehensively reflected.

[0039] S120. Perform data matching on the first seismic inversion data and the first electromagnetic inversion data based on a preset rule to obtain the second seismic inversion data and the second electromagnetic inversion data.

[0040] Among them, the preset rule is used to perform data matching on the first seismic inversion data and the first electromagnetic inversion data, so that the number and interval of data points in the first seismic inversion data and the first electromagnetic inversion data after the matching process are exactly the same.

[0041] In an embodiment of the present invention, the first seismic inversion data and the first electromagnetic inversion data can be subjected to data matching based on a preset rule to obtain the second seismic inversion data and the second electromagnetic inversion data. Specifically, based on the preset rule, according to the first spatial coordinate and the second spatial coordinate, the first seismic inversion data and the first electromagnetic inversion data can be subjected to matching processing, and the mismatched or abnormal data points are excluded during the processing, and the first seismic inversion data obtained after the processing is used as the second seismic inversion data, and the first electromagnetic inversion data obtained after the processing is used as the second electromagnetic inversion data. By performing data matching on the first seismic inversion data and the first electromagnetic inversion data based on the preset rule, the consistency of the obtained second seismic inversion data and the second electromagnetic inversion data in the spatial coordinate can be ensured, the overall quality of the inversion data can be improved, and a data basis can be provided for subsequent data fusion and geological interpretation.

[0042] S130. Perform data transformation on the second seismic inversion data and the second electromagnetic inversion data respectively based on a composite trigonometric function to obtain the third seismic inversion data and the third electromagnetic inversion data.

[0043] Among them, the composite trigonometric function is a function obtained by operating (such as addition, subtraction, multiplication, division, composition, etc.) on two or more basic trigonometric functions (such as sine function, cosine function, tangent function, etc.).

[0044] In an embodiment of the present invention, the second seismic inversion data and the second electromagnetic inversion data can be subjected to data transformation respectively based on a composite trigonometric function to obtain the third seismic inversion data and the third electromagnetic inversion data. Specifically, the second seismic inversion data can be used as an independent variable to input into the composite trigonometric function, and the final output of the composite trigonometric function is calculated as the third seismic inversion data. The second electromagnetic inversion data is used as an independent variable to input into the composite trigonometric function, and the final output of the composite trigonometric function is calculated as the third electromagnetic inversion data. By performing data transformation on the second seismic inversion data and the second electromagnetic inversion data respectively based on the composite trigonometric function to obtain the third seismic inversion data and the third electromagnetic inversion data, the characteristics of the periodic description ability, smoothing and filtering effect, and nonlinear transformation ability of the composite trigonometric function can be fully utilized to meet the requirements of retaining structural features, enhancing data compatibility, and improving the fusion efficiency during the data fusion process, so as to facilitate the rapid fusion of two different types of inversion data while retaining and displaying the original structural feature information in the second seismic inversion data and the second electromagnetic inversion data, and the implementation method is simple and efficient.

[0045] S140. Perform data fusion on the third seismic inversion data and the third electromagnetic inversion data to interpret the geological structure of the work area based on the inversion data obtained after fusion.

[0046] In an embodiment of the present invention, the third seismic inversion data and the third electromagnetic inversion data can be fused to interpret the geological structure of the work area based on the fused inversion data. Specifically, after obtaining the third seismic inversion data and the third electromagnetic inversion data, a suitable data fusion method, such as weighted average, principal component analysis, wavelet transform, etc., can be selected to fuse the third seismic inversion data and the third electromagnetic inversion data, and the geological structure of the work area can be interpreted according to the fused inversion data, providing a scientific basis for resource exploration and development. By fusing the third seismic inversion data and the third electromagnetic inversion data, the complementarity of the third seismic inversion data and the third electromagnetic inversion data in reflecting the geophysical characteristics of the work area can be utilized to fully utilize the information of these two types of data and improve the accuracy and reliability of geological interpretation.

[0047] The technical solution of the embodiment of the present invention is to obtain the first seismic inversion data and the first electromagnetic inversion data of the work area; wherein, the first seismic inversion data includes the first spatial coordinates and the velocity corresponding to the first spatial coordinates, and the first electromagnetic inversion data includes the second spatial coordinates and the resistivity corresponding to the second spatial coordinates; based on a preset rule, the first seismic inversion data and the first electromagnetic inversion data are matched to obtain the second seismic inversion data and the second electromagnetic inversion data; based on the composite trigonometric function, the second seismic inversion data and the second electromagnetic inversion data are respectively transformed to obtain the third seismic inversion data and the third electromagnetic inversion data; the third seismic inversion data and the third electromagnetic inversion data are fused to interpret the geological structure of the work area based on the fused inversion data. The technical solution of the embodiment of the present invention can, while retaining the original structural feature information in the inversion data, enable two different types of inversion data to be quickly fused by performing data transformation on the second seismic inversion data and the second electromagnetic inversion data obtained by data matching based on the composite trigonometric function, thereby improving the accuracy and reliability of geological structure interpretation.

[0048] Embodiment 2

[0049] Figure 2 It is a flowchart of a method for processing seismic and electromagnetic inversion data based on a composite trigonometric function provided by Embodiment 2 of the present invention. The embodiment of the present invention is optimized based on the above embodiment, and the solutions not described in detail in the embodiment of the present invention can be seen in the above embodiment. As Figure 2 shown, the method includes:

[0050] S210. Obtain the first seismic inversion data and the first electromagnetic inversion data of the work area.

[0051] Exemplarily, Figure 3 shows a schematic diagram of a first seismic inversion data, Figure 4 shows a schematic diagram of a first electromagnetic inversion data.

[0052] S220. Align the first seismic inversion data and the first electromagnetic inversion data in the spatial domain based on the first spatial coordinate and the second spatial coordinate.

[0053] In the embodiment of the present invention, the first seismic inversion data and the first electromagnetic inversion data can be aligned in the spatial domain based on the first spatial coordinate and the second spatial coordinate. Specifically, due to factors such as data acquisition methods, instrument settings, or geographical locations, there may be differences in the coordinate systems of the first spatial coordinate and the second spatial coordinate. Therefore, it is necessary to first determine whether the coordinate systems of the two are the same. If they are different, it is necessary to determine the conversion relationship between the two coordinate systems and unify them into one coordinate system through the conversion relationship. After that, the first seismic inversion data and the first electromagnetic inversion data are aligned in the spatial domain in the same coordinate system. By aligning the first seismic inversion data and the first electromagnetic inversion data in the spatial domain based on the first spatial coordinate and the second spatial coordinate, it can be ensured that the aligned first seismic inversion data and the first electromagnetic inversion data are superimposed or compared at the same spatial position, facilitating comprehensive analysis.

[0054] S230. Perform interpolation processing on the first seismic inversion data and the first electromagnetic inversion data after spatial domain alignment respectively to obtain the second seismic inversion data and the second electromagnetic inversion data.

[0055] In the embodiment of the present invention, since the acquisition densities and resolutions of the seismic inversion data and the electromagnetic inversion data may be different, after the first seismic inversion data and the first electromagnetic inversion data are aligned in the spatial domain, it is also necessary to perform interpolation processing on the first seismic inversion data and the first electromagnetic inversion data after spatial domain alignment respectively to obtain the second seismic inversion data and the second electromagnetic inversion data to ensure that the two types of inversion data are on the same spatial grid.

[0056] Optionally, performing interpolation processing on the first seismic inversion data and the first electromagnetic inversion data after spatial domain alignment respectively to obtain the second seismic inversion data and the second electromagnetic inversion data includes: determining the range and interval of the interpolation processing based on the first spatial coordinate and the second spatial coordinate, and generating an interpolation grid based on the range and the interval; mapping the first seismic inversion data after spatial domain alignment to the interpolation grid through an inverse distance interpolation algorithm to obtain the second seismic inversion data; mapping the first electromagnetic inversion data after spatial domain alignment to the interpolation grid through an inverse distance interpolation algorithm to obtain the second electromagnetic inversion data.

[0057] Among them, the inverse distance interpolation algorithm, that is, the inverse distance weighted interpolation algorithm, is an interpolation method based on the distance between discrete points. Its basic idea is that the value of an unknown point is obtained by weighting the values of the known points closest to it, and the weight is proportional to the reciprocal of the distance.

[0058] In an embodiment of the present invention, the range and interval of interpolation processing can be determined based on the first spatial coordinate and the second spatial coordinate. The range of interpolation processing can be determined according to the maximum and minimum values in the first spatial coordinate and the second spatial coordinate. The interval of interpolation processing needs to select an appropriate interval to include all data points including the first spatial coordinate and the second spatial coordinate. After determining the range and interval of interpolation processing, an interpolation network can be generated according to the range and interval of interpolation processing. After that, the first seismic inversion data aligned in the spatial domain can be mapped to the interpolation grid through the inverse distance interpolation algorithm to complete the interpolation processing of the first seismic inversion data aligned in the spatial domain, and the second seismic inversion data can be obtained. The second electromagnetic inversion data can be obtained in the same way. Thus, it can be ensured that the number and interval of corresponding data points of the velocity data in the second seismic inversion data and the resistivity data in the second electromagnetic inversion data are exactly the same, and the matching processing of the velocity data and the resistivity data can be realized.

[0059] Optionally, after respectively performing interpolation processing on the first seismic inversion data and the first electromagnetic inversion data aligned in the spatial domain to obtain the second seismic inversion data and the second electromagnetic inversion data, it further includes: unifying the formats of the second seismic inversion data and the second electromagnetic inversion data into a file format for storing seismic record data.

[0060] In an embodiment of the present invention, after obtaining the second seismic inversion data and the second electromagnetic inversion data, the formats of the second seismic inversion data and the second electromagnetic inversion data can be unified into a file format for storing seismic record data, such as the SGY format, so as to facilitate subsequent data transformation and data fusion processing.

[0061] S240. Perform preprocessing on the second seismic inversion data and the second electromagnetic inversion data; wherein, the preprocessing includes at least one of denoising, filtering, and normalization.

[0062] In an embodiment of the present invention, before performing data transformation on the second seismic inversion data and the second electromagnetic inversion data, preprocessing operations including at least one of denoising, filtering, and normalization can be performed on the second seismic inversion data and the second electromagnetic inversion data to improve the data quality of the second seismic inversion data and the second electromagnetic inversion data and ensure the accuracy and reliability of the data.

[0063] S250. Input the preprocessed second seismic inversion data into a composite trigonometric function to obtain the third seismic inversion data, and input the preprocessed second electromagnetic inversion data into a composite trigonometric function to obtain the third electromagnetic inversion data.

[0064] Optionally, input the preprocessed second seismic inversion data into a composite trigonometric function to obtain third seismic inversion data, and input the preprocessed second electromagnetic inversion data into the composite trigonometric function to obtain third electromagnetic inversion data, including: determining an intermediate output of the composite trigonometric function based on the input inversion data, the sine value of the input inversion data, and the cosine value of the input inversion data; determining the target inversion data based on the calculated intermediate output, the maximum value of the calculated intermediate output, and the minimum value of the calculated intermediate output; wherein, when the input inversion data is the preprocessed second seismic inversion data, the target inversion data is the third seismic inversion data, and when the input inversion data is the preprocessed second electromagnetic inversion data, the target inversion data is the third electromagnetic inversion data.

[0065] In an embodiment of the present invention, the input inversion data, that is, the preprocessed second seismic inversion data or the preprocessed second electromagnetic inversion data, can be input into the inversion data and used as an independent variable to input into the composite trigonometric function. First, determine the intermediate output of the composite trigonometric function based on the input inversion data, the sine value of the input inversion data, and the cosine value of the input inversion data, and then determine the target inversion data, that is, the third seismic inversion data or the third electromagnetic inversion data, based on the calculated intermediate output, the maximum value of the calculated intermediate output, and the minimum value of the calculated intermediate output. It should be noted that when the input inversion data is the preprocessed second seismic inversion data, the target inversion data is the third seismic inversion data, and when the input inversion data is the preprocessed second electromagnetic inversion data, the target inversion data is the third electromagnetic inversion data. Exemplarily, Figure 5 shows a schematic diagram of a kind of third seismic inversion data, Figure 6 shows a schematic diagram of a kind of third electromagnetic inversion data.

[0066] Optionally, the composite trigonometric function is:

[0067]

[0068] wherein, x is the input inversion data, y is the calculated intermediate output, y min is the minimum value of the calculated intermediate output, y max is the maximum value of the calculated intermediate output, and d is the target inversion data.

[0069] S260. Take the square root of the product of the third seismic inversion data and the third electromagnetic inversion data to obtain the fused inversion data.

[0070] In an embodiment of the present invention, the product of the third seismic inversion data and the third electromagnetic inversion data can be calculated first, and then the square root of the product of the third seismic inversion data and the third electromagnetic inversion data can be taken to complete the data fusion of the third seismic inversion data and the third electromagnetic inversion data, and obtain the fused inversion data. Exemplarily,Figure 7 A schematic diagram of the fused inversion data is shown. Optionally, taking the square root of the product of the third seismic inversion data and the third electromagnetic inversion data gives the fused inversion data, and its calculation formula is as follows:

[0071]

[0072] where d 3 is the fused inversion data, d 1 is the third seismic inversion data, and d 2 is the third electromagnetic inversion data.

[0073] S270. Interpret at least one of the structural characteristics, lithology distribution, and reservoir characteristics of the work area based on the fused inversion data.

[0074] In the embodiments of the present invention, after obtaining the fused inversion data, based on the information of the comprehensive seismic inversion data and electromagnetic inversion data in the fused inversion data, and using the complementarity of the seismic inversion data and electromagnetic inversion data in reflecting the geophysical characteristics of the work area, the accuracy and reliability of the interpretation of at least one of the structural characteristics, lithology distribution, and reservoir characteristics of the work area can be improved, providing a scientific basis for resource exploration and development.

[0075] The technical solution of the embodiments of the present invention is to obtain the first seismic inversion data and the first electromagnetic inversion data of the work area; perform spatial domain alignment on the first seismic inversion data and the first electromagnetic inversion data based on the first spatial coordinate and the second spatial coordinate; perform interpolation processing on the spatially domain-aligned first seismic inversion data and first electromagnetic inversion data respectively to obtain the second seismic inversion data and the second electromagnetic inversion data; perform preprocessing on the second seismic inversion data and the second electromagnetic inversion data; where the preprocessing includes at least one of denoising, filtering, and normalization; input the preprocessed second seismic inversion data into a composite trigonometric function to obtain the third seismic inversion data, and input the preprocessed second electromagnetic inversion data into a composite trigonometric function to obtain the third electromagnetic inversion data; take the square root of the product of the third seismic inversion data and the third electromagnetic inversion data to obtain the fused inversion data; interpret at least one of the structural characteristics, lithology distribution, and reservoir characteristics of the work area based on the fused inversion data. The technical solution of the embodiments of the present invention can perform data transformation on the second seismic inversion data and the second electromagnetic inversion data obtained by data matching based on the composite trigonometric function, and while retaining the original structural characteristic information in the inversion data, enable the two different types of inversion data to be quickly fused, thereby improving the accuracy and reliability of geological structure interpretation and providing a scientific basis for resource exploration and development.

[0076] Embodiment III

[0077] Figure 8The structural schematic diagram of a seismic and electromagnetic inversion data processing device provided in Embodiment 3 of the present invention is as follows. As Figure 8 shown, the device includes:

[0078] An inversion data acquisition module 310, configured to acquire first seismic inversion data and first electromagnetic inversion data of a work area; wherein, the first seismic inversion data includes first spatial coordinates and the velocity corresponding to the first spatial coordinates, and the first electromagnetic inversion data includes second spatial coordinates and the resistivity corresponding to the second spatial coordinates;

[0079] An inversion data matching module 320, configured to perform data matching on the first seismic inversion data and the first electromagnetic inversion data based on a preset rule to obtain second seismic inversion data and second electromagnetic inversion data;

[0080] An inversion data transformation module 330, configured to perform data transformation on the second seismic inversion data and the second electromagnetic inversion data respectively based on a composite trigonometric function to obtain third seismic inversion data and third electromagnetic inversion data;

[0081] An inversion data fusion module 340, configured to perform data fusion on the third seismic inversion data and the third electromagnetic inversion data, so as to interpret the geological structure of the work area based on the inversion data obtained after fusion.

[0082] Optionally, the inversion data matching module 320 includes:

[0083] A spatial domain alignment unit, configured to perform spatial domain alignment on the first seismic inversion data and the first electromagnetic inversion data based on the first spatial coordinates and the second spatial coordinates;

[0084] An interpolation processing unit, configured to perform interpolation processing on the first seismic inversion data and the first electromagnetic inversion data after spatial domain alignment respectively to obtain second seismic inversion data and second electromagnetic inversion data.

[0085] Optionally, the interpolation processing unit includes:

[0086] An interpolation grid generation sub-unit, configured to determine the range and interval of interpolation processing based on the first spatial coordinates and the second spatial coordinates, and generate an interpolation grid based on the range and the interval;

[0087] A first interpolation processing sub-unit, configured to map the first seismic inversion data after spatial domain alignment to the interpolation grid through an inverse distance interpolation algorithm to obtain second seismic inversion data;

[0088] The second interpolation processing subunit is configured to map the first electromagnetic inversion data after spatial domain alignment to the interpolation grid through an inverse distance interpolation algorithm to obtain the second electromagnetic inversion data.

[0089] Optionally, the inversion data transformation module 330 includes:

[0090] The inversion data preprocessing unit is configured to preprocess the second seismic inversion data and the second electromagnetic inversion data; wherein, the preprocessing includes at least one of denoising, filtering, and normalization.

[0091] The inversion data transformation unit is configured to input the preprocessed second seismic inversion data into a composite trigonometric function to obtain the third seismic inversion data, and input the preprocessed second electromagnetic inversion data into a composite trigonometric function to obtain the third electromagnetic inversion data.

[0092] Optionally, the inversion data transformation unit includes:

[0093] The intermediate output determination subunit is configured to determine the intermediate output of the composite trigonometric function according to the input inversion data, the sine value of the input inversion data, and the cosine value of the input inversion data.

[0094] The target inversion data determination subunit is configured to determine the target inversion data according to the calculated intermediate output, the maximum value of the calculated intermediate output, and the minimum value of the calculated intermediate output; wherein, when the input inversion data is the preprocessed second seismic inversion data, the target inversion data is the third seismic inversion data, and when the input inversion data is the preprocessed second electromagnetic inversion data, the target inversion data is the third electromagnetic inversion data.

[0095] Optionally, the inversion data fusion module 340 includes:

[0096] The inversion data fusion unit is configured to take the square root of the product of the third seismic inversion data and the third electromagnetic inversion data to obtain the fused inversion data.

[0097] The geological structure interpretation unit is configured to interpret at least one of the structural characteristics, lithology distribution, and reservoir characteristics of the work area based on the fused inversion data.

[0098] Optionally, the inversion data matching module 320 further includes:

[0099] The format conversion unit is configured to unify the formats of the second seismic inversion data and the second electromagnetic inversion data into a file format for storing seismic record data.

[0100] The seismic and electromagnetic inversion data processing device based on composite trigonometric functions provided by the embodiments of the present invention can execute the seismic and electromagnetic inversion data processing method based on composite trigonometric functions provided by any embodiment of the present invention, and has the corresponding functional modules and beneficial effects for executing the method.

[0101] Embodiment 4

[0102] Figure 9 FIG. shows a schematic structural diagram of an electronic device 10 that can be used to implement the embodiments of the present invention. The electronic device is intended to represent various forms of digital computers, such as, laptop computers, desktop computers, workstations, personal digital assistants, servers, blade servers, mainframe computers, and other suitable computers. The electronic device can also represent various forms of mobile devices, such as, personal digital processors, cellular phones, smart phones, wearable devices (such as helmets, glasses, watches, etc.) and other similar computing devices. The components shown herein, their connections and relationships, and their functions are merely illustrative and are not intended to limit the implementation of the present invention described and / or claimed herein.

[0103] As Figure 9 shown, the electronic device 10 includes at least one processor 11, and a memory communicatively connected to the at least one processor 11, such as a read-only memory (ROM) 12, a random access memory (RAM) 13, etc. The memory stores a computer program executable by the at least one processor. The processor 11 can perform various appropriate actions and processes according to the computer program stored in the read-only memory (ROM) 12 or the computer program loaded from the storage unit 18 into the random access memory (RAM) 13. In the RAM 13, various programs and data required for the operation of the electronic device 10 can also be stored. The processor 11, the ROM 12, and the RAM 13 are connected to each other through a bus 14. The input / output (I / O) interface 15 is also connected to the bus 14.

[0104] Multiple components in the electronic device 10 are connected to the I / O interface 15, including: an input unit 16, such as a keyboard, a mouse, etc.; an output unit 17, such as various types of displays, speakers, etc.; a storage unit 18, such as a magnetic disk, an optical disc, etc.; and a communication unit 19, such as a network card, a modem, a wireless communication transceiver, etc. The communication unit 19 allows the electronic device 10 to exchange information / data with other devices through a computer network such as the Internet and / or various telecommunication networks.

[0105] The processor 11 can be various general-purpose and / or special-purpose processing components with processing and computing capabilities. Some examples of the processor 11 include, but are not limited to, a central processing unit (CPU), a graphics processing unit (GPU), various dedicated artificial intelligence (AI) computing chips, various processors running machine learning model algorithms, a digital signal processor (DSP), and any suitable processor, controller, microcontroller, etc. The processor 11 executes the various methods and processes described above, such as the seismic and electromagnetic inversion data processing method based on composite trigonometric functions.

[0106] In some embodiments, the seismic and electromagnetic inversion data processing method based on composite trigonometric functions can be implemented as a computer program tangibly embodied in a computer-readable storage medium, such as the storage unit 18. In some embodiments, part or all of the computer program can be loaded and / or installed onto the electronic device 10 via the ROM 12 and / or the communication unit 19. When the computer program is loaded into the RAM 13 and executed by the processor 11, one or more steps of the seismic and electromagnetic inversion data processing method based on composite trigonometric functions described above can be executed. Alternatively, in other embodiments, the processor 11 can be configured to execute the seismic and electromagnetic inversion data processing method based on composite trigonometric functions by any other suitable means (e.g., by means of firmware).

[0107] Various embodiments of the systems and techniques described above in this document can be implemented in digital electronic circuit systems, integrated circuit systems, field-programmable gate arrays (FPGA), application-specific integrated circuits (ASIC), application-specific standard products (ASSP), system-on-chip systems (SOC), complex programmable logic devices (CPLD), computer hardware, firmware, software, and / or combinations thereof. These various embodiments can include: being implemented in one or more computer programs that can be executed and / or interpreted on a programmable system including at least one programmable processor, which can be a special-purpose or general-purpose programmable processor that receives data and instructions from a storage system, at least one input device, and at least one output device, and transmits the data and instructions to the storage system, the at least one input device, and the at least one output device.

[0108] A computer program for implementing the method of the present invention may be written in any combination of one or more programming languages. These computer programs may be provided to a processor of a general purpose computer, a special purpose computer, or other programmable data processing apparatus, such that the computer programs, when executed by the processor, cause the functions / operations specified in the flowchart and / or block diagram to be implemented. The computer program may be executed entirely on the machine, partly on the machine, as a stand-alone software package partly on the machine and partly on a remote machine, or entirely on the remote machine or server.

[0109] In the context of the present invention, a computer-readable storage medium may be a tangible medium that can contain, or store a computer program for use by or in connection with an instruction execution system, apparatus, or device. The computer-readable storage medium may include, but is not limited to, electronic, magnetic, optical, electromagnetic, infrared, or semiconductor systems, apparatus, or devices, or any suitable combination of the foregoing. Alternatively, the computer-readable storage medium may be a machine-readable signal medium. More specific examples of the machine-readable storage medium would include an electrical connection based on one or more wires, a portable computer diskette, a hard disk, a random access memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or Flash memory), an optical fiber, a portable compact disc read-only memory (CD-ROM), an optical storage device, a magnetic storage device, or any suitable combination of the foregoing.

[0110] In order to provide interaction with a user, the systems and techniques described herein may be implemented on an electronic device having: a display device (e.g., a CRT (cathode ray tube) or LCD (liquid crystal display) monitor) for displaying information to the user; and a keyboard and a pointing device (e.g., a mouse or a trackball) by which the user can provide input to the electronic device. Other kinds of devices may also be used to provide interaction with the user; for example, the feedback provided to the user may be any form of sensory feedback (e.g., visual feedback, auditory feedback, or tactile feedback); and input from the user may be received in any form (including acoustic input, voice input, or tactile input).

[0111] The systems and techniques described herein can be implemented in a computing system including backend components (e.g., as a data server), or a computing system including middleware components (e.g., an application server), or a computing system including frontend components (e.g., a user computer having a graphical user interface or a web browser through which a user can interact with an implementation of the systems and techniques described herein), or a computing system including any combination of such backend components, middleware components, or frontend components. The components of the system can be interconnected to each other by digital data communication in any form or medium (e.g., a communication network). Examples of communication networks include: local area network (LAN), wide area network (WAN), blockchain network, and the Internet.

[0112] A computing system can include a client and a server. The client and the server are generally remote from each other and typically interact through a communication network. The client-server relationship is created by computer programs running on respective computers and having a client-server relationship with each other. The server can be a cloud server, also known as a cloud computing server or a cloud host, which is a host product in the cloud computing service system, and solves the defects of difficult management and weak business scalability existing in traditional physical hosts and VPS services.

[0113] It should be understood that the various forms of processes shown above can be used, with steps reordered, added, or deleted. For example, the steps recited in the present invention can be executed in parallel, sequentially, or in a different order, as long as the desired results of the technical solution of the present invention can be achieved, and no limitation is made herein.

[0114] The above specific embodiments do not constitute a limitation on the protection scope of the present invention. Those skilled in the art should understand that various modifications, combinations, sub-combinations, and substitutions can be made according to design requirements and other factors. Any modifications, equivalent substitutions, and improvements made within the spirit and principles of the present invention shall be included within the protection scope of the present invention.

Claims

1. A method for processing seismic and electromagnetic inversion data based on composite trigonometric functions, characterized in that: The method comprises: Acquire first seismic inversion data and first electromagnetic inversion data of the work area; wherein the first seismic inversion data includes a first spatial coordinate and a velocity corresponding to the first spatial coordinate, and the first electromagnetic inversion data includes a second spatial coordinate and a resistivity corresponding to the second spatial coordinate; Performing data matching on the first seismic inversion data and the first electromagnetic inversion data based on a preset rule to obtain second seismic inversion data and second electromagnetic inversion data; Based on the composite trigonometric function, the second seismic inversion data and the second electromagnetic inversion data are respectively transformed to obtain third seismic inversion data and third electromagnetic inversion data; The third seismic inversion data and the third electromagnetic inversion data are fused to interpret the geological structure of the work area based on the fused inversion data.

2. The method according to claim 1, characterized in that The first seismic inversion data and the first electromagnetic inversion data are matched based on a preset rule to obtain second seismic inversion data and second electromagnetic inversion data, including: Aligning the first seismic inversion data and the first electromagnetic inversion data in a spatial domain based on the first spatial coordinate and the second spatial coordinate; Interpolation processing is performed on the first seismic inversion data and the first electromagnetic inversion data after alignment in the spatial domain to obtain second seismic inversion data and second electromagnetic inversion data.

3. The method according to claim 2, characterized in that Interpolation processing is performed on the first seismic inversion data and the first electromagnetic inversion data after spatial domain alignment to obtain second seismic inversion data and second electromagnetic inversion data, including: determining a range and an interval of interpolation processing based on the first spatial coordinates and the second spatial coordinates, and generating an interpolation grid based on the range and the interval; Mapping the first seismic inversion data after spatial domain alignment onto the interpolation grid by an inverse distance interpolation algorithm to obtain second seismic inversion data; The first electromagnetic inversion data after spatial domain alignment is mapped onto the interpolation grid by an inverse distance interpolation algorithm to obtain second electromagnetic inversion data.

4. The method according to claim 1, characterized in that The second seismic inversion data and the second electromagnetic inversion data are respectively transformed based on a composite trigonometric function to obtain third seismic inversion data and third electromagnetic inversion data, including: Preprocessing the second seismic inversion data and the second electromagnetic inversion data; wherein the preprocessing includes at least one of denoising, filtering and normalization; The preprocessed second seismic inversion data is input into a composite trigonometric function to obtain third seismic inversion data, and the preprocessed second electromagnetic inversion data is input into a composite trigonometric function to obtain third electromagnetic inversion data.

5. The method according to claim 4, characterized in that Inputting the preprocessed second seismic inversion data into a composite trigonometric function to obtain third seismic inversion data, and inputting the preprocessed second electromagnetic inversion data into a composite trigonometric function to obtain third electromagnetic inversion data, including: Determining an intermediate output of the composite trigonometric function based on input inversion data, a sine value of the input inversion data, and a cosine value of the input inversion data; The target inversion data is determined according to the calculated intermediate output, the maximum value of the calculated intermediate output, and the minimum value of the calculated intermediate output; wherein, when the input inversion data is the second seismic inversion data after preprocessing, the target inversion data is the third seismic inversion data, and when the input inversion data is the second electromagnetic inversion data after preprocessing, the target inversion data is the third electromagnetic inversion data.

6. The method according to claim 1, characterized in that The third seismic inversion data and the third electromagnetic inversion data are fused to interpret the geological structure of the work area based on the fused inversion data, including: Taking a square root of the product of the third seismic inversion data and the third electromagnetic inversion data to obtain fused inversion data; At least one of the structural characteristics, lithology distribution and reservoir characteristics of the work area is interpreted based on the fused inversion data.

7. The method according to claim 2, characterized in that After interpolating the first seismic inversion data and the first electromagnetic inversion data after spatial domain alignment to obtain the second seismic inversion data and the second electromagnetic inversion data, the method further includes: The formats of the second seismic inversion data and the second electromagnetic inversion data are unified into a file format for storing seismic record data.

8. A seismic and electromagnetic inversion data processing device based on composite trigonometric functions, characterized in that: The device comprises: An inversion data acquisition module, used to acquire first seismic inversion data and first electromagnetic inversion data of a work area; wherein the first seismic inversion data includes a first spatial coordinate and a velocity corresponding to the first spatial coordinate, and the first electromagnetic inversion data includes a second spatial coordinate and a resistivity corresponding to the second spatial coordinate; An inversion data matching module, used for performing data matching on the first seismic inversion data and the first electromagnetic inversion data based on a preset rule to obtain second seismic inversion data and second electromagnetic inversion data; An inversion data transformation module, used for performing data transformation on the second seismic inversion data and the second electromagnetic inversion data based on a composite trigonometric function to obtain third seismic inversion data and third electromagnetic inversion data; The inversion data fusion module is used to fuse the third seismic inversion data and the third electromagnetic inversion data so as to interpret the geological structure of the work area based on the inversion data obtained after fusion.

9. An electronic device, characterized in that: The electronic device comprises: at least one processor; and a memory communicatively connected to the at least one processor; wherein, The memory stores a computer program executable by the at least one processor, and the computer program is executed by the at least one processor to enable the at least one processor to perform the seismic and electromagnetic inversion data processing method based on complex trigonometric functions as described in any one of claims 1 to 7.

10. A computer-readable storage medium, characterized in that: The computer-readable storage medium stores computer instructions, and the computer instructions are used to enable a processor to implement the seismic and electromagnetic inversion data processing method based on composite trigonometric functions as described in any one of claims 1 to 7 when executed.

Images